Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 7446.
Published online 2017 August 7. doi:  10.1038/s41598-017-07950-x
PMCID: PMC5547144

Associations among oxytocin receptor gene (OXTR) DNA methylation in adulthood, exposure to early life adversity, and childhood trajectories of anxiousness


Recent models propose deoxyribonucleic acid methylation of key neuro-regulatory genes as a molecular mechanism underlying the increased risk of mental disorder associated with early life adversity (ELA). The goal of this study was to examine the association of ELA with oxytocin receptor gene (OXTR) methylation among young adults. Drawing from a 21-year longitudinal cohort, we compared adulthood OXTR methylation frequency of 46 adults (23 males and 23 females) selected for high or low ELA exposure based on childhood socioeconomic status and exposure to physical and sexual abuse during childhood and adolescence. Associations between OXTR methylation and teacher-rated childhood trajectories of anxiousness were also assessed. ELA exposure was associated with one significant CpG site in the first intron among females, but not among males. Similarly, childhood trajectories of anxiousness were related to one significant CpG site within the promoter region among females, but not among males. This study suggests that females might be more sensitive to the impact of ELA on OXTR methylation than males.


Exposure to early life adversity (ELA), including physical and sexual abuse and neglect as well as poor socioeconomic conditions during childhood, are associated with increased risk for anxiety and depression later in life1, 2. A growing body of evidence suggests a biological embedding of early experiences leading to lasting consequences on autonomic, neuroendocrine, immune, and neural function3. Recent models propose deoxyribonucleic acid (DNA) methylation of key neuro-regulatory genes as an underlying molecular mechanism responsible for the increased risk for mental disorders associated with ELA exposure4, 5.

Gene methylation is an epigenetic process impacting transcriptional activity of DNA without altering the actual DNA nucleotide sequence. Gene methylation represents a covalent modification of DNA by the addition of a methyl group onto cytosine rings found within 5′-Cytosine-phosphate-Guanine-3′ (CpG) dinucleotide pairs6. Gene methylation typically influences phenotypic variation by suppressing gene expression via the prevention of transcription factors binding to their cis-acting elements in gene promoters and enhancers, and through the recruitment of methylated-DNA binding proteins that lead to the formation of a closed chromatin configuration blocking access to the DNA6. DNA methylation is the closest epigenetic mark to the gene itself and it is assumed to be the most stable epigenetic modification7.

Early life is a period particularly sensitive to the regulatory effects of epigenetic mechanisms8. A seminal study demonstrated that maternal care influenced DNA methylation of the glucocorticoid receptor (GR) in rats9. Offspring of low licking and grooming mothers had greater methylation of GR exon 17 promoter region, compared to offspring of mothers providing more maternal care. Furthermore, these changes in DNA methylation were associated with GR expression in the hippocampus and physiological stress responses in adulthood9. Other animal studies have replicated the effects of ELA exposure on DNA methylation of other candidate genes. Notably, manipulation of ELA impacted methylation of the ESR1, CRH, BDNF, and AVP genes, providing strong evidence for the impact of ELA on gene methylation913.

Early life adversity also leads to genome-wide methylation modifications. Among rhesus macaques, differential rearing conditions (maternal vs. surrogate peer rearing) led to global methylation differences in as many as ~1300 distinct gene promoters in T-cells and in the prefrontal cortex14. In humans, institutionalized children exhibited differential methylation in 800 gene promoters from whole blood, compared to children raised by their biological parents15. Furthermore, genome-wide promoter methylation profiling revealed that ELA was associated with a distinct epigenetic signature in peripheral blood leukocytes and hippocampal tissues in adulthood1619. Notably, the epigenetic profile was more strongly associated with early life adversity than adversity later in life17, 20, 21. This suggests that ELA exposure leads to genome-wide epigenetic changes in peripheral and central tissues.

ELA has been associated with DNA methylation of candidate genes involved in key neuro-regulatory functions among both healthy individuals and individuals suffering from mental disorders. Indeed, childhood adversity has been associated with greater methylation of different candidate genes, such as NR3C1, FKBP5, SLC6A4, and BDNF, in peripheral blood leukocytes of healthy and depressed children and adults2230. Furthermore, post-mortem analysis of suicide victims’ hippocampal tissues revealed that a history of childhood abuse was related to greater methylation of NR3C1 as well as the ribosomal ribonucleic acid (RNA) genes19, 31, indicating that DNA methylation changes were also observed in central tissues. Importantly, DNA methylation mediated the impact of early adversity on later health outcomes3234.

The oxytocinergic system is a key neurobiological system involved in social behaviors. In animal studies, pharmacological oxytocin (OT) manipulation facilitated social recognition, social bonding, maternal care, aggression, and sexual activity35. In humans, intra-nasal oxytocin administration modulated social cognition and behaviors36. The oxytocinergic system is also associated with anxiety and depression. In animal models, oxytocin receptor gene (OXTR) knockout mice did not benefit from the antidepressant effect of mating behavior, compared to wild type mice37, 38. Further, an oxytocin receptor agonist had antidepressant properties39, 40, while an oxytocin receptor antagonist blocked the buffering effect of social relationships on depressive behavior41, 42. In humans, depression and anxiety have been associated with dysregulated plasma OT4347. Further, OXTR polymorphisms have been associated with stress reactivity and psychological distress4852, especially in the context of ELA exposure5356. Importantly, oxytocin appears to have a sexually dimorphic role, with females showing a greater association between oxytocin and prosocial behaviors than males57.

The oxytocinergic system shows considerable plasticity in response to early life events. In rodents and primate studies, ELA exposure was associated with changes in oxytocin receptor binding and expression in limbic and prefrontal areas, as well as altered cerebrospinal fluid OT levels5868. In humans, individuals with a history of ELA had lower cerebrospinal fluid OT levels69 and exhibited altered cortisol, limbic, and caregiving responses to intra-nasal OT administration7072, compared to individuals without an ELA history. Moreover, children raised in an orphanage had lower OT production during mother-infant interactions, compared to children raised by their biological parents73.

To date, only a few studies have investigated methylation of oxytocinergic genes. A single oxytocin receptor gene (OXTR) localized on chromosome 3 at locus 3p25 has been identified74. Early epigenetic studies using luciferase reporter gene assays showed that methylation of a CpG island within the first intron influenced transcriptional activity75. In a seminal paper, five CpG dinucleotides of this CpG island had significantly higher DNA methylation in patients with autism, compared to healthy control subjects76. This difference was found not only in blood DNA, but also in DNA extracted from brain tissue76, providing some evidence for correspondence of brain and blood DNA methylation for this gene. Subsequent studies examining this same CpG island found that differences in methylation frequency from whole blood were related to individual differences in anorexia nervosa symptoms, child conduct problems, unemotional and callous-unemotional trait, acute psychosocial stress, perception of ambiguous social stimuli, and limbic activation in response to fear7782. Furthermore, studies investigating other CpG sites within OXTR have found significant associations with social anxiety and depressive disorders8385, as well as pessimism and interpersonal distrust86.

A few studies have investigated the association of ELA with OXTR methylation in adulthood. In a rodent study, natural variations in maternal care were associated with altered OXTR DNA methylation in central and peripheral tissues among adult rats87. In humans, an epidemiological study reported that low childhood socioeconomic status (SES) was associated with increased OXTR DNA methylation in non-promoter regions of the gene, but was not related to methylation of CpG sites within the promoter region88. Moreover, poorer maternal care in childhood was associated with greater OXTR methylation within exon 3 in peripheral blood cells89. Furthermore, in a sample of African-Americans with low socio-economic status, early child abuse was associated with higher OXTR DNA methylation of two CpG sites in exon 3 in whole blood, but not at the promoter region90. These data provide preliminary evidence that ELA is associated with methylation of different regions within OXTR.

Although past studies have examined the association of ELA with OXTR methylation in adulthood, they have focused on a limited subset of regulatory regions within OXTR. Notably, past studies have focused on a limited number of CpG sites within the promoter region. Promoter regions are interesting targets because they are enriched in both CpG dinucleotides and binding sites for transcriptional activators and repressors. Methylation of such regions can hinder the proper interaction between DNA and the transcription factors. Disruptions of such interactions bring about inhibition of downstream gene expression91. Furthermore, chromatin immune-precipitation (ChIP)-sequencing indicates that there are two enhancer elements within the 3rd intron, which might be subjected to regulation by DNA methylation92. Enhancer elements are short distal regions involved in the co-regulation of gene transcription, usually acting cis to the promoter elements93. Moreover, past studies have used retrospective and limited assessments of ELA.

The primary aim of this study was to examine the predictive association between prospectively assessed ELA and later peripheral OXTR methylation in adulthood. Various regulatory regions of the gene were assessed in order to identify key CpG sites sensitive to ELA. Given that methylation mediated the impact of ELA on later health outcomes in animal and human studies9, 32, 33, a secondary aim of this study was to evaluate the associations between OXTR methylation and childhood trajectories of anxiousness and disruptiveness. Furthermore, given the sexually dimorphic role of oxytocin in social behaviors57, sex-specific analyses were also conducted. We hypothesized that ELA exposure would be associated with greater OXTR methylation, that OXTR methylation would be related to childhood trajectories of anxiousness and disruptiveness, and that females would show a stronger association between ELA and OXTR methylation.



Participants were recruited from l’Étude longitudinale des enfants de maternelle au Québec (ÉLEMQ), a longitudinal study of 3785 children initially recruited while they were attending kindergarten in francophone schools in Québec. From this larger sample, a randomly selected, representative group of 2000 boys and girls was followed longitudinally. The cohort was followed yearly from ages 6–12, and then in mid-adolescence (mean age = 15), in emerging adulthood (mean age = 21) and in adulthood (mean age = 27). At age 27, participants provided a blood sample for epigenetic analysis. A subset of these participants was selected for exposure to high or low levels of early adversity.

In order to create two groups with differential ELA exposure, 46 participants with available DNA were selected based on their scores on both childhood SES and early exposure to abuse (see below). These two aspects were considered given that childhood socioeconomic status has a different methylation signature than childhood abuse17, 94, suggesting that these two factors might independently predict adulthood methylation profiles. Specifically, two extreme groups were established on the basis of children’s z-scores for both indices calculated on the entire cohort. Participants in the high ELA group (n = 24) scored low on the SES index, but high on the Abuse index, whereas participants in the low ELA group (n = 22) scored high on the SES index, but low on the Abuse index. An equal number of males and females were selected in each group. To reduce genetic admixture, only Caucasian individuals of Western European Ancestry were included in the study. This study was approved by the Hospital Ste-Justine Research Ethics Board. The methods were carried out in accordance with the relevant guidelines and regulations. All participants provided informed consent.

Psychosocial Measures

Childhood Socio-Economic Status (SES)

Childhood socioeconomic status was well characterized in this cohort using five relevant socioeconomic indicators collected prospectively across childhood. Maternal and paternal years of schooling, mean maternal and paternal occupational prestige, and mean family income were assessed at 8 occasions when the participants were aged 6–12, and were then averaged. Confirmatory factor analysis aggregating the 5 indicators was used to create a general childhood SES score for each individual.

Childhood Abuse Index

Exposure to physical and sexual abuse during childhood and adolescence was assessed retrospectively at age 2195. Eight items of the Parent–Child Conflict Tactics Scale96 were used to evaluate child physical abuse by a mother or father figure. Participants indicated how often they experienced each of the following items during their childhood: severe physical abuse (4 items; e.g., hit you with an object, brutally threw you against a wall) and very severe physical abuse (4 items; e.g., threatened you with a weapon, beat you up over and over). Five questions regarding childhood sexual abuse were adapted from the Adverse Childhood Experiences Questionnaire97 and from the Sexually Victimized Children Questionnaire98. Participants were asked if they had experienced any unwanted sexual acts against their will before the age of 18 years, including exhibitionism (being forced to look or forced to show genitals), sexual fondling or touching, and completed or attempted sexual intercourse by use of bribes or threats, force, or drugs and/or alcohol. The scores from these scales were combined using confirmatory factor analysis to form an overall Abuse index.

Childhood trajectories of anxiousness and disruptiveness

Items from the teacher-rated Social Behavior Questionnaire99, administered yearly from age 6 to 12, were considered to evaluate trajectories of childhood anxiousness and childhood disruptiveness. The trajectories of childhood anxiousness were assessed using the following items: fearful or afraid of things or new situations; is worried, worries about many things; cries easily; has a tendency to work alone; looks sad, unhappy, tearful; easily distracted (age 6 Cronbach’s α = 0.74). The trajectories of childhood disruptiveness items included: destroys one’s own things or those of others; fights with other kids; is not liked by peers; irritable; disobedient; lies; mistreats, intimidates peers; does not share material used for a particular task; blames others; inconsiderate of others; hits and kicks others; fidgets, squirms, cannot keep still; agitated, always running and jumping, restless (age 6 Cronbach’s α = 0.90). The trajectories (low, average and high) were characterized using semi-parametric group-based modeling100. Children classified in the low trajectory were consistently rated as displaying low levels of anxiousness or disruptiveness, while participants in the high trajectory had elevated ratings of anxiousness or disruptiveness throughout childhood. Participants in the average trajectories varied between high and low yearly ratings of anxiousness or disruptiveness from age 6 to 12.

OXTR Target Sequence Selection

To identify potential regulatory regions of the gene, ChIP-sequencing experiment data available in an open-access database part of the ENCODE histone project were used101. H3K4Me1 and H3K4Me3 signal tracks using PBMC and HepG2 cell lines were obtained through the ENCODE database and visualized via the UCSC genome browser using the genome build hg19102. H3K4Me1 and H3K4Me3 are specific histone modifying proteins that are indicative of active enhancers and promoters respectively. Regions of peak signal intensity are regions within the DNA sequence that bind specifically to H3K4Me1 and H3K4Me3 antibodies with high signal strength, representing robust areas of protein binding when taking into account the background noise inherent in the ChIP-sequencing procedure103. The exact DNA sequence corresponding to a signal region was extracted with full annotation of its CpG dinucleotides. For the promoter, the region length was defined from 700 base pairs upstream of the gene to the transcription start site (TSS). For enhancers, a DNA sequence was extracted 100 base pairs upstream and downstream from the H3K4Me1 signal peak, which corresponds to the region with the highest confidence of protein binding activity101. As depicted in Fig. 1, CpG dinucleotides within 3 genomic regions were selected for the present study.

Figure 1
Schematic representation of the genomic regions of interest within the OXTR gene. The promoter region (chr3: 8811303-8811915) was identified using the H3K4Me3 signals from the ENCODE database. The Intron 1 region (chr3: 8810699-8810875) was identified ...

Sample Preparation

Peripheral blood samples were collected from participants and stored in EDTA coated tubes at 4 °C before extraction. DNA extraction was performed using the QIAamp DNA Mini Kit (Qiagen, #51304) according to the manufacturer’s instructions and stored at −20 °C. Extracted DNA was stored in a −80 °C freezer.


To investigate OXTR DNA methylation, a total of three sets of outside primers and four sets of inside primers were developed to probe all CpG sites within the target regions (promoter, intron 1, enhancers). The nested reverse primers were biotinylated for pyrosequencing (IDT Technologies). 500 ng of DNA was treated with sodium bisulphite (EZ Methylation Gold, Zymo Research) and underwent two rounds of PCR amplification (#1, 95 C × 15 min, [94 × 1 min, Primer TM * 1 min, 72 C * 1 min] for 35 cycles, 72 C * 10 min; #2 95 C × 15 min, [94 × 1 min, Primer TM * 1 min, 72 C * 1 min] for 40 cycles, 72 C * 10 min). The subsequent PCR product then underwent gel electrophoresis to confirm the purity and the success of the amplification protocol. 20 ul of the PCR product was then used to perform pyrosequencing using PyroMarkQ24 (Qiagen) according to the manufacturer protocol. All PCR primers used are listed in Supplementary Table S1. The methylation percentage at each individual CpG site was analyzed and exported using PyroMark Q24 software (Qiagen). Triplicate analyses were performed per sample to assure accuracy. Data are reported as the average of the triplicates. We were unable to quantify Enhancer 2 methylation because of difficulties in designing the sequencing primer that adheres to the target segment with high affinity due to the repetitive composition of the Enhancer Region 2 DNA sequence. A list of all successfully quantified CpG sites with their respective genomic positions is listed in Supplementary Table S2. Site-specific methylation analyses were performed at CFI Imaging and Molecular Biology Platform at McGill University in the Department of Pharmacology and Therapeutics.

Luciferase Reporter Construct

Three artificial constructs of the OXTR promoter were generated (602 base pair whole promoter sense, 602 base pair whole promoter antisense, and 50 base pair promoter) using primers with artificially introduced restriction sites to allow for sense and antisense insertion of DNA fragments. A sense sequence is a DNA sequence in its 5′ to 3′ direction while an antisense sequence is the complement of the sense sequence in the reverse 3′ to 5′ direction. Human genomic DNA from whole blood was subject to two rounds of PCR amplification with HotStar Taq DNA Polymerase (Qiagen). The resulting DNA was digested with BamHI and HindIII and subcloned into the CpG-free pCpGL-basic luciferase reporter plasmid. SssI DNA methyltransferase (New England Biolabs) was then used to methylate the plasmid construct in vitro.

The pCpGL-basic plasmid is a construct with no-inherent CpG sites; methylation of such plasmid occurs exclusively on the introduced gene vector sequence. This strategy ensures that there is no confounding effect of methylation sites along the whole plasmid construct. The resultant constructs were then validated via sequencing (Genome Quebec, Montreal). Both the methylated and non-methylated plasmids were then transfected into HEK2 93 (human embryonic kidney) cells using standard methods. HEK293 cells were cultured in DMEM 1X (GIBCO, Invitrogen) with 10% fetal bovine serum (GIBCO, Invitrogen), plated and transfected using a standard calcium phosphate method104. Cells were lysed and harvested 48 hours post transfection, and luciferase activity was measured via the Luciferase Assay System (Promega). The primers used are listed in Supplementary Table S3.

Statistical Analysis

All data were inspected for normality of the distributions. All variables approximated the normal distribution except for promoter CpG3. The distribution of promoter CpG 3 showed a significant departure from normality driven by one outlier value. However, even when this outlier value was removed from the analysis, the effect of ELA remained statistically significant after a Bonferroni correction. Given the small sample size, we decided to retain all the data in the analyses. Multivariate and univariate general linear models explored the associations among ELA, OXTR methylation, and childhood trajectories of anxiousness and disruptiveness. Weighted posterior probabilities were used to test the associations between OXTR DNA methylation and the childhood trajectories of anxiousness and disruptiveness100. It was hypothesized that OXTR methylation would mediate the relationships between ELA and childhood trajectories of anxiousness among females. Mediation analyses were performed by testing the significance of the indirect effect using bootstrapping resampling methods with the PROCESS module for SPSS105. Using this approach, the indirect effect is considered statistically significant when its confidence interval does not include zero. Logistic regressions were used for the mediation given the categorical dependent variable. This method is superior to traditional tests of mediation because it does not require large sample sizes nor does it assume normal distribution of the indirect effect106. For the luciferase assay, paired two-tailed t-tests were used to examine the signal differences in control and test plasmids. Given the small number of missing data, mean substitution was used to deal with missing data in multivariate tests. All the data available were used for the univariate analyses. Statistical significance was reported at p < 0.05. Family-wise Bonferroni corrections were used to limit type I error due to multiple comparisons. Given the risk of type II error in a small sample, interpretation focused on large effect sizes (Cohen’s d > 0.80). Effect sizes quantified using Cohen’s d were calculated by dividing the mean difference between the two groups by the pooled standard deviation107. Partialη 2 was used as measure of effect size in multivariate analyses. It represents the proportion of the total variance in a dependent variable that is associated with the independent grouping variable after partialing out the effects of other independent variables and interaction effects108. Statistical analyses were conducted using SPSS v. 22.


OXTR DNA Methylation and Early Life Adversity

DNA methylation frequency of three genomic regions within OXTR was quantified in 46 individuals via pyrosequencing. All participants were 27 years of age and of Western European ancestry, with an equal number of males and females (23 males, 23 females) and a similar number of individuals in high (n = 24) and low (n = 22) ELA groups. DNA methylation frequency across the 16 CpG sites that were successfully quantified ranged from 1% to 95%. The average methylation frequencies of all CpG sites as a function of ELA exposure for males and females are presented in Figs 2 and and3,3, respectively.

Figure 2
OXTR Methylation as a Function of ELA Exposure in Males. *Represents statistically significant ELA group differences at p < 0.05; **represents statistically significant difference after a Bonferonni correction where p < 0.003. ...
Figure 3
OXTR Methylation as a Function of ELA Exposure in Females. *Represents statistically significant ELA group differences at p < 0.05; **represents statistically significant difference after a Bonferonni correction where p < 0.003. ...

There was no overall DNA methylation difference between the high and low ELA groups, F(16,29) = 1.03, p = 0.46, η 2 = 0.36. Given that the association between ELA and OXTR methylation varied across genomic regions109, 110, the associations of ELA with individual CpG sites were evaluated. Individuals in the high ELA group had significantly greater promoter CpG 7 methylation, compared to participants in the low ELA group, F (1,44) = 5.3, p = 0.03, d = 0.69, mean methylation difference = 5.9%. However, this difference was no longer statistically significant after a Bonferroni correction, adjusted p = 0.48.

Sex Differences in OXTR Methylation

Given previous evidence showing a sexually dimorphic effect of OT on social behavior111, we examined sex differences in OXTR methylation, regardless of ELA. Three CpG sites exhibited significant sex differences. Females had significantly lower DNA methylation in promoter CpG 7 than males, F(1,44) = 6.0, p = 0.02, d = 0.64. In contrast, females had higher DNA methylation in Intron CpG 2, F(1,44) = 7.49, p = 0.01, d = 0.82, and Intron CpG 3, F(1,44) = 7.11, p = 0.01, d = 0.81, compared to males.

ELA and OXTR Methylation Among Females and Males

The association between ELA and global OXTR methylation was a trend, F(6,16) = 3.06, p = 0.09, partial η 2 = 0.89, among females, while there was no significant association among males, F(6,16) = 0.62, p = 0.79, partial η 2 = 0.62. Tables 1 and and22 provide ELA group differences for each individual CpG site for females and males respectively. For females, the association between ELA and OXTR methylation had a large effect size, i.e., a Cohen’s d > 0.80, for 5 CpG sites. For promoter CpG 3, promoter CpG 7, intron 1 CpG 4, and intron 1 CpG 5, females in the high ELA group had higher OXTR DNA methylation than females in the low ELA group. In contrast, high ELA was associated with lower methylation for enhancer 1 CpG 2. After applying a Bonferroni correction for multiple comparisons, there was a significant group difference only for intron 1 CpG 5. There was no ELA group difference with a large effect size in males.

Table 1
Effect Sizes of the Association Between ELA and OXTR Methylation for the Different CpG Sites among Females.
Table 2
Effect Sizes of the Association between ELA and OXTR Methylation for the Different CpG Sites among Males.

OXTR Methylation and Childhood Trajectories of Anxiousness and Disruptiveness

In this sample selected for high and low ELA exposure, 39.13% of the participants belonged to the low childhood anxiousness trajectory, 47.83% were assigned to average childhood anxiousness, and 13.04% displayed elevated anxiousness throughout childhood. There was an equivalent number of males and females in each trajectory, χ2 (45) = 0.0001, p = 0.99. For the childhood disruptiveness trajectory, 52.17% of the participants were assigned to the low trajectory, 28.26% to the average trajectory, and 19.56% to the high trajectory. There was no significant sex difference in the childhood disruptiveness trajectories, χ2(45) = 3.19, p = 0.20.

A MANOVA tested the associations between childhood trajectories of anxiousness and methylation of the 5 OXTR CpG sites that showed large effect sizes with ELA. OXTR methylation was significantly related to childhood anxiousness among females, F(5,17) = 4.03, p = 0.01, partial η 2 = 0.54, but not among males, F(5,17) = 1.10, p = 0.40, partial η 2 = 0.37. Post-hoc tests indicated that for promoter CpG 3, F(2,20) = 4.45, p = 0.03, partial η 2 = 0.31, and promoter CpG 7, F(2,30) = 8.89, p = 0.002, partial η 2 = 0.47, higher methylation was associated with greater childhood anxiousness. However, only the association of promoter CpG 7 with childhood anxiousness survived correction for multiple comparisons. Intron 1 CpG 4, Intron1 CpG 5, and Enhancer 1 CpG 2 were not significantly related to childhood anxiousness among females, all p-values > 0.07. Moreover, for both females, F(5,17) = 1.84, p = 0.16, partial η 2 = 0.35, and males, F(5,17) = 0.57, p = 0.11, partial η 2 = 0.14, there was no significant association between childhood trajectory of disruptiveness and OXTR methylation. Figure 4 displays the relationship between OXTR methylation and childhood trajectories of anxiousness.

Figure 4
Associations between OXTR methylation and childhood trajectories of anxiousness. Figure 4 depicts the association between promoter CpG 7 methylation and childhood anxiousness. Individuals within the high childhood anxiousness trajectory had significant ...

Exploratory mediation analysis

Prior analyses indicated that there were large Cohen’s d effect sizes for the associations among ELA, OXTR DNA methylation at the promoter CpG 7 site, and childhood trajectories of anxiousness. We conducted an exploratory analysis to examine whether OXTR methylation may act as a mediator of the association between ELA and childhood trajectories of anxiousness. Given the sex-specific associations observed, the mediation models tested whether OXTR methylation mediated the relationship between ELA and childhood trajectories of anxiousness in the female subsample only. The indirect effect tested the extent to which the association between ELA and childhood anxiousness is explained by their common association with OXTR methylation. The indirect effect for promoter CpG 7 was statistically significant, b = 0.35 (SE = 0.17), CI (0.07–0.78). In contrast, the indirect effect of an opposite model whereby promoter CpG 7 methylation predicted early life adversity through increased anxiousness was not statistically significant, b = 0.75, SE = 1.23 CI: −0.03–5.85. However, these results should be interpreted as exploratory given that not all individual paths from the model remained significant after correction for multiple comparisons. The mediation model is illustrated in Fig. 5.

Figure 5
Mediation model. !Indicates the direct effect of ELA on childhood anxiousness (c path); *indicates the direct effect path coefficient after adjusting for all the other effects in the model (c’ path). In this model, the confidence interval of the ...

Functional validation of the differentially methylated region in OXTR

While there were significant correlations among ELA, OXTR DNA methylation within the promoter, and childhood trajectories of anxiousness, it is still unclear whether methylation of CpG sites within the promoter does indeed have an effect on the transcriptional machinery that may alter OXTR expression or other affected downstream genes. To determine the functional activity of the CpG within the promoter, we introduced the differentially methylated regions in OXTR (Ctrl-no insert, promoter sense, promoter antisense, and promoter subregion) to the pCpGL-basic CpG-free luciferase reporter plasmid and performed in vitro methylation with SSSI methyltransferase. Since the plasmid and the reporter do not contain CpG sequences, SSSI methylates the OXTR regions exclusively. We could therefore measure the effects of DNA methylation of this OXTR region without confounding effects of vector methylation. Introduction of the putative promoter region to the reporter vector induced reporter luciferase activity, confirming that the region is indeed an active promoter. We then compared the luciferase activity driven by the unmethylated promoter with the unmethylated promoter sense plasmid construct, the empty vector and methylated and unmethylated anti-sense constructs (Promoter sense unmethylated versus sense methylated, t = 7.38, p = 0.002. promoter sense unmethylated versus antisense unmethylated, t = −4.41, p = 0.01, promoter sense methylated versus antisense unmethylated, t = 2.81, p = 0.04, promoter antisense unmethylated versus promoter antisense methylated, t = 0.40, p = 0.70) (see Fig. 6). These data are consistent with the idea that CpG sites within the promoter region are important for regulating the expression of the OXTR gene.

Figure 6
In vitro methylation regulates gene expression at OXTR promoter. Normalized luciferase activity (RLU/mg protein) in the HEK293 cell line for the antisense, methylated antisense, sense, and methylated sense promoter constructs are shown. Values are expressed ...


We investigated associations among ELA, OXTR DNA methylation, and childhood trajectories of anxiousness and disruptiveness among young adults. Among females, one CpG site within the first intron was significantly associated with ELA and one CpG site within the promoter was significantly related to childhood trajectories of anxiousness. No significant associations were found between ELA, childhood anxiousness or disruptiveness, and OXTR methylation among males. Lastly, the functional significance of promoter CpG methylation was validated using in vitro methylation of a plasmid construct with inserted promoter sequence.

In the present study, an ELA measure encompassing both childhood SES and abuse was associated with large Cohen’s d effect sizes in differences in 5 CpG sites within different regulatory regions of the gene. However, only one site within the first intron remained significantly associated with ELA after correction for multiple comparisons. Although prior studies have reported associations between ELA exposure and OXTR DNA methylation, the specific CpG sites varied across studies. Unternaehrer et al.110 reported an association between low maternal care and OXTR DNA methylation within the 3rd exon, while Smearman et al.90 found that early abuse exposure was associated with greater methylation in two CpG sites within the OXTR promoter. In contrast, Needleman et al.109 reported that low childhood SES was associated with OXTR methylation in non-promoter regions of the genes, but not in the OXTR promoter. The present study investigated 3 genomic regions within OXTR, including an a priori defined region within intron 1 that has been previously associated with autism symptoms, externalizing behaviours, psychological distress, and social perception112, as well as the promoter region upstream of the transcription start site next to exon 1, and a distal enhancer element located within intron 3, identified using ChIP-sequencing experiment data as potentially regulatory regions of the gene. While there was no overall change in methylation in these different genomic regions, specific CpG sites within each region were differentially methylated as a function of ELA exposure among females. Results from past studies are not directly comparable with the present findings because they used different conceptualizations of ELA and they differ in the specific CpG sites assessed, with some studies relaying on commercial arrays and other studies using different a priori defined regions. This highlights the importance of using a similar set of CpG dinucleotides in future research to facilitate comparison across studies.

Among females, high ELA was related to higher mean DNA methylation in intron 1 CpG 5 and trajectories of greater childhood anxiousness were related to higher methylation of promoter CpG 7. The effect sizes of the differences in methylation between ELA and childhood anxiousness groups were equivalent to large effect sizes according to Cohen’s criteria107, with 8.13 and 22.22% mean methylation difference between groups. This is consistent with other studies of early life adversity, DNA methylation, and psychiatric disorders whereby mean DNA methylation differences related to phenotypic or environmental factors ranged from 1% to 5% using salivary or blood samples23, 113115. Greater DNA methylation, especially in the promoter of the gene, is typically associated with lower gene expression116. Results of the in vitro luciferase experiment confirmed that greater methylation within the promoter region resulted in lower gene expression.

The oxytocinergic system is involved in the expression of social behaviors and has been related to risk for anxiety and depression. Oxytocin exerts most of its behavioural effects by binding to receptors located in the brain35. There is evidence that OXTR methylation from peripheral cells correlates with OXTR methylation in brain cells76, 87. Furthermore, a number of studies suggest that OXTR methylation from whole blood is significantly associated with structural and functional neural processes. OXTR methylation was associated with neural activation in response to ambiguous social stimuli81, decreased functioning coupling between the amygdala and the prefrontal cortex in response to angry or fearful faces82, and volumetric differences in temporal-limbic and prefrontal regions involved in social cognition117. In addition, OXTR methylation has been related to social cognition processes and styles such as emotion recognition as well as interpersonal distrust86, 117. OXTR methylation may thus lead to subtle changes in the neural networks supporting social cognition, leading to increased risk for anxiety and depression later in life. An exploratory analysis indicated that there was an indirect effect of the relationship between ELA and childhood anxiousness through OXTR promoter CpG 7 methylation among females. Future studies should attempt to replicate this finding in larger samples.

In this sample, OXTR DNA methylation was not associated with childhood trajectories of disruptiveness. This is in contrast with studies showing an association between OXTR DNA methylation and callous-unemotional traits among children78. However, other studies observed that OXTR methylation was associated with callous-unemotional traits only among adolescents, but not among children77. This suggests that our childhood measure of disruptiveness might have been less sensitive to OXTR DNA methylation than if it had been collected during adolescence. Alternatively, differences between our disruptiveness measure, including aspects of hyperactivity and oppositional behaviors, and the callous-unemotional trait assessments performed in other studies might also explain the discrepant results.

The association between OXTR DNA methylation and ELA was stronger among females than among males. These data parallel findings by Rubin et al.117 who observed sex-specific associations between OXTR methylation and plasma OT levels, performance at an emotion recognition task, and volumetric differences in brain regions involved in fear and social cognition. Notably, OXTR methylation was associated with these behavioral and physiological phenotypes among females, but not among males117. These data are broadly consistent with the sexually dimorphic effects of oxytocin on social behavior in different rodent species111 as well as with the sex differences in the effect of intra-nasal oxytocin administration in some118120, but not all human studies121. While these results need replication, they suggest that women’s greater sensitivity to the impact of ELA on OXTR methylation may underlie some of the sex differences in anxiety and depression later in life122.

The exact mechanism through which early life experiences are translated into DNA methylation changes is not yet clear. It has been postulated that external stimuli may trigger specific signaling pathways, which then recruit DNA methylation-specific proteins to exert their actions on specific genes and gene segments123. For example, it is generally thought that maternal behavior triggers specific serotonin signaling pathways in the brain, followed by release of secondary messenger signals (cAMP), which recruit DNA methylation-specific enzymes that aid in targeting specific genome locations and methylating/demethylating the corresponding positions5. This is thought to lead to active transcription or repression of gene expression. If women are preferentially recruiting the oxytocinergic system during social interactions, this may explain the sexual dimorphism observed in this study.

ELA not only affects DNA methylation within candidate genes, but also has genome-wide and system-wide effects124. Given that genes do not act individually but in clusters of functional circuitries, it is likely that changes in one candidate gene will have ramifications in other genes within the same functional pathway125. Szyf & Bick5 suggest that OXTR DNA methylation may affect other downstream effectors within the oxytocinergic system such as the CD38 or OXT genes as well as genes related to the vasopressinergic system and other genes of the same functional circuitry involved in the regulation of social behavior126. Future studies should assess methylation of genes within the same functional circuitry as the OXTR gene.

Recent studies suggest that DNA methylation of certain CpG islands is influenced by the individual’s genotype127. With regards to the OXTR gene, there is preliminary evidence that OXTR genotype may influence DNA methylation of specific OXTR CpG dinucleotides90. Furthermore, OXTR single nucleotide polymorphisms (SNPs) interacted with OXTR methylation to predict current anxiety and depression symptoms90. These OXTR genotype by methylation interactions predicting current psychological distress have been observed in other studies85, 128, 129. However, interactions of specific SNPs and CpG sites have not been consistent across these studies. Genotype information for the OXTR gene was not available in the current study. Future studies should assess the interaction between OXTR methylation and genotype in the prediction of behavioral phenotypes. Larger samples encompassing a broad representation of the different OXTR genotypes will be necessary to examine this issue.

This study possesses several unique strengths. The study participants were selected from a subsample of a longitudinal cohort that was recruited at the same age and from the same ethnic background, Caucasians of Western European Ancestry, minimizing age- and race-related methylation differences. A second strength is the prospective assessment for early socio-economic status during childhood rather than a retrospective assessment used in most studies. Third, our ELA assessment considered both early socioeconomic status as well as exposure to physical and sexual abuse, providing a more comprehensive assessment of early life adversity. Fourth, childhood trajectories of anxiousness and disruptiveness were teacher-rated, minimizing the risk for gene-environment correlations associated with mother- or participant-reports. Finally, the functional significance of promoter CpG methylation was confirmed using an in vitro luciferase experiment.

One limitation of this study is related to the timing of the DNA methylation assessment. Although OXTR methylation was conceptualized as the mediator of the effects of ELA, its assessment at the age of 27 followed rather than preceded the assessment of the childhood trajectories of anxiousness. This approach is defendable because ELA-related methylation has been found to be relatively stable in other studies and to be more strongly related to early life stress than current stress17, 130, 131. However, given that no information is available on the stability of the methylation signature of this specific gene from this specific tissue132, an alternative explanation is also plausible. Namely, the possibility that childhood anxiousness could also be driving the changes in methylation cannot be ruled out in this study design. Furthermore, given the correlational design of the study, no causal inference can be made from these data. Moreover, DNA methylation frequency was assessed using whole blood. It is possible that the cell type composition of our sample might have influenced the DNA methylation frequency observed133. Furthermore, this study did not include information about current anxiety and depression. These factors may also impact current OXTR methylation83, 85. Also, given that participants were selected for extreme exposure to ELA, it is possible that a different pattern of results may emerge among individuals with less severe exposure to ELA. Other forms of ELA such as child neglect and household dysfunction should also be assessed in other studies97. Furthermore, given that many CpG sites were included in the promoter construct, results of the luciferase experiment indicate that the promoter region influences gene expression, but it does not provide evidence for the role of specific CpG sites. Moreover, this study was based on relatively small number of participants. Despite large effect sizes, most of the significant associations became non-significant after adjustment for multiple comparisons due to the large number of CpG sites assessed. Furthermore, only a sex-specific correction for multiple comparisons was applied in the present study; all significant findings would become non-significant if the number of comparisons across sexes was considered for the Bonferroni correction. Given the goal to identify novel CpG dinucleotides sensitive to ELA, a relatively large number of CpG sites was evaluated in this small study, thereby increasing risk of Type I error. Replications of the present results using a priori hypotheses regarding specific CpG sites in larger samples are thus paramount.

In sum, this study provides preliminary evidence that ELA is associated with lasting changes in OXTR methylation and that these changes are related to phenotypic differences in internalizing symptoms among females. This suggests that the early social environment may shape the epigenetic regulation of the OXTR gene in a way that increases risk for anxiety and depression later in life. By identifying potential CpG sites sensitive to ELA and internalizing symptoms, these preliminary results may help uncover the molecular, physiological, and psychosocial processes underlying the protracted risk of ELA on anxiety and depression later in life among females.

Electronic supplementary material


The authors would like to thank the participants for their continued involvement in the cohort. This study was part of Q.Q. Zhou’s master’s thesis. This study was funded by grants from: Quebec Network on Suicide, Mood Disorders, and Related Disorders (JPG), the Canadian Institutes of Health Research (RT), the Social Sciences and Humanities Research Council of Canada (RT), the Fonds de Recherche du Québec- Santé (LB), the Canada Research Chair Secretariat (JPG), and the CIHR-GSK chair in pharmacology (MS).

Author Contributions

Author Contributions

J.P.G., L.B., M.S. and Q.Q.Z. designed the study, R.T., M.B., F.V., G.T. and S.C. recruited and followed the participants, M.H. and I.O.M. designed the adversity measures, Q.Q.Z. and M.S. conducted the methylation analysis, J.P.G. and Q.Q.Z. analyzed the data and wrote the manuscript. All authors reviewed and approved the manuscript.


Competing Interests

The authors declare that they have no competing interests.


Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-07950-x

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Li M, D’Arcy C, Meng X. Maltreatment in childhood substantially increases the risk of adult depression and anxiety in prospective cohort studies: systematic review, meta-analysis, and proportional attributable fractions. Psychol Med. 2016;46:717–730. doi: 10.1017/S0033291715002743. [PubMed] [Cross Ref]
2. Gilman SE, Kawachi I, Fitzmaurice GM, Buka SL. Socioeconomic status in childhood and the lifetime risk of major depression. Int J Epidemiol. 2002;31:359–367. doi: 10.1093/intjepid/31.2.359. [PubMed] [Cross Ref]
3. Danese A, McEwen BS. Adverse childhood experiences, allostasis, allostatic load, and age-related disease. Physiology & behavior. 2012;106:29–39. doi: 10.1016/j.physbeh.2011.08.019. [PubMed] [Cross Ref]
4. Lutz PE, Turecki G. DNA methylation and childhood maltreatment: from animal models to human studies. Neuroscience. 2014;264:142–156. doi: 10.1016/j.neuroscience.2013.07.069. [PMC free article] [PubMed] [Cross Ref]
5. Szyf M, Bick J. DNA methylation: a mechanism for embedding early life experiences in the genome. Child Dev. 2013;84:49–57. doi: 10.1111/j.1467-8624.2012.01793.x. [PMC free article] [PubMed] [Cross Ref]
6. Razin A. CpG methylation, chromatin structure and gene silencing-a three-way connection. EMBO J. 1998;17:4905–4908. doi: 10.1093/emboj/17.17.4905. [PubMed] [Cross Ref]
7. Hochberg, Z. et al. Child health, developmental plasticity, and epigenetic programming. Endocr Rev. 32, 159–224, doi:er.2009-0039 (2011). [PubMed]
8. Szyf, M. Early life, the epigenome and human health. Acta Paediatr. 98, 1082–1084, doi:APA1382 (2009). [PubMed]
9. Weaver IC, et al. Epigenetic programming by maternal behavior. Nat Neurosci. 2004;7:847–854. doi: 10.1038/nn1276. [PubMed] [Cross Ref]
10. Champagne, F. A. et al. Maternal care associated with methylation of the estrogen receptor-alpha1b promoter and estrogen receptor-alpha expression in the medial preoptic area of female offspring. Endocrinology. 147, 2909–2915, doi:en.2005–1119 (2006). [PubMed]
11. Korosi, A. et al. Early-life experience reduces excitation to stress-responsive hypothalamic neurons and reprograms the expression of corticotropin-releasing hormone. J Neurosci. 30, 703–713, doi:30/2/703 (2010). [PMC free article] [PubMed]
12. Roth, T. L., Lubin, F. D., Funk, A. J. & Sweatt, J. D. Lasting epigenetic influence of early-life adversity on the BDNF gene. Biol Psychiatry. 65, 760–769, doi:S0006-3223(08)01530-8 (2009). [PMC free article] [PubMed]
13. Murgatroyd, C. et al. Dynamic DNA methylation programs persistent adverse effects of early-life stress. Nat Neurosci. 12, 1559–1566, doi:nn.2436 (2009). [PubMed]
14. Provencal N, et al. The signature of maternal rearing in the methylome in rhesus macaque prefrontal cortex and T cells. J Neurosci. 2012;32:15626–15642. doi: 10.1523/JNEUROSCI.1470-12.2012. [PMC free article] [PubMed] [Cross Ref]
15. Naumova OY, et al. Differential patterns of whole-genome DNA methylation in institutionalized children and children raised by their biological parents. Dev Psychopathol. 2012;24:143–155. doi: 10.1017/S0954579411000605. [PMC free article] [PubMed] [Cross Ref]
16. Essex MJ, et al. Epigenetic vestiges of early developmental adversity: childhood stress exposure and DNA methylation in adolescence. Child Dev. 2013;84:58–75. doi: 10.1111/j.1467-8624.2011.01641.x. [PMC free article] [PubMed] [Cross Ref]
17. Borghol N, et al. Associations with early-life socio-economic position in adult DNA methylation. Int J Epidemiol. 2012;41:62–74. doi: 10.1093/ije/dyr147. [PMC free article] [PubMed] [Cross Ref]
18. Labonte B, et al. Genome-wide epigenetic regulation by early-life trauma. Arch Gen Psychiatry. 2012;69:722–731. doi: 10.1001/archgenpsychiatry.2011.2287. [PMC free article] [PubMed] [Cross Ref]
19. McGowan PO, et al. Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat Neurosci. 2009;12:342–348. doi: 10.1038/nn.2270. [PMC free article] [PubMed] [Cross Ref]
20. Essex MJ, et al. Epigenetic Vestiges of Early Developmental Adversity: Childhood Stress Exposure and DNA Methylation in Adolescence. Child Dev. 2011 [PMC free article] [PubMed]
21. Labonte, B. et al. Differential glucocorticoid receptor exon 1(B), 1(C), and 1(H) expression and methylation in suicide completers with a history of childhood abuse. Biol Psychiatry72. 41–48, doi:S0006-3223(12)00132-1 (2012). [PubMed]
22. Tyrka AR, Price LH, Marsit C, Walters OC, Carpenter LL. Childhood adversity and epigenetic modulation of the leukocyte glucocorticoid receptor: preliminary findings in healthy adults. PLoS One. 2012;7:e30148. doi: 10.1371/journal.pone.0030148. [PMC free article] [PubMed] [Cross Ref]
23. Perroud N, et al. Increased methylation of glucocorticoid receptor gene (NR3C1) in adults with a history of childhood maltreatment: a link with the severity and type of trauma. Transl Psychiatry. 2011;1:e59. doi: 10.1038/tp.2011.60. [PMC free article] [PubMed] [Cross Ref]
24. Beach SR, Brody GH, Todorov AA, Gunter TD, Philibert RA. Methylation at SLC6A4 is linked to family history of child abuse: an examination of the Iowa Adoptee sample. Am J Med Genet B Neuropsychiatr Genet. 2010;153B:710–713. doi: 10.1002/ajmg.b.31028. [PMC free article] [PubMed] [Cross Ref]
25. Perroud N, et al. Response to psychotherapy in borderline personality disorder and methylation status of the BDNF gene. Transl Psychiatry. 2013;3:e207. doi: 10.1038/tp.2012.140. [PMC free article] [PubMed] [Cross Ref]
26. Bick J, et al. Childhood adversity and DNA methylation of genes involved in the hypothalamus-pituitary-adrenal axis and immune system: whole-genome and candidate-gene associations. Dev Psychopathol. 2012;24:1417–1425. doi: 10.1017/S0954579412000806. [PMC free article] [PubMed] [Cross Ref]
27. Klengel T, et al. Allele-specific FKBP5 DNA demethylation mediates gene-childhood trauma interactions. Nat Neurosci. 2013;16:33–41. doi: 10.1038/nn.3275. [PMC free article] [PubMed] [Cross Ref]
28. Booij L, et al. DNA methylation of the serotonin transporter gene in peripheral cells and stress-related changes in hippocampal volume: a study in depressed patients and healthy controls. PloS one. 2015;10:e0119061. doi: 10.1371/journal.pone.0119061. [PMC free article] [PubMed] [Cross Ref]
29. Tyrka AR, et al. Childhood maltreatment and methylation of FK506 binding protein 5 gene (FKBP5) Dev Psychopathol. 2015;27:1637–1645. doi: 10.1017/S0954579415000991. [PMC free article] [PubMed] [Cross Ref]
30. Kang HJ, et al. Association of SLC6A4 methylation with early adversity, characteristics and outcomes in depression. Prog Neuropsychopharmacol Biol Psychiatry. 2013;44:23–28. doi: 10.1016/j.pnpbp.2013.01.006. [PubMed] [Cross Ref]
31. McGowan PO, et al. Promoter-wide hypermethylation of the ribosomal RNA gene promoter in the suicide brain. PLoS One. 2008;3:e2085. doi: 10.1371/journal.pone.0002085. [PMC free article] [PubMed] [Cross Ref]
32. Cao-Lei L, et al. DNA methylation mediates the effect of exposure to prenatal maternal stress on cytokine production in children at age 13(1/2) years: Project Ice Storm. Clin Epigenetics. 2016;8:54. doi: 10.1186/s13148-016-0219-0. [PMC free article] [PubMed] [Cross Ref]
33. Cao-Lei L, et al. DNA methylation mediates the impact of exposure to prenatal maternal stress on BMI and central adiposity in children at age 13(1/2) years: Project Ice Storm. Epigenetics. 2015;10:749–761. doi: 10.1080/15592294.2015.1063771. [PMC free article] [PubMed] [Cross Ref]
34. Kumsta R, et al. Severe psychosocial deprivation in early childhood is associated with increased DNA methylation across a region spanning the transcription start site of CYP2E1. Transl Psychiatry. 2016;6:e830. doi: 10.1038/tp.2016.95. [PMC free article] [PubMed] [Cross Ref]
35. Carter CS. Oxytocin pathways and the evolution of human behavior. Annual review of psychology. 2014;65:17–39. doi: 10.1146/annurev-psych-010213-115110. [PubMed] [Cross Ref]
36. Bartz, J. A., Zaki, J., Bolger, N. & Ochsner, K. N. Social effects of oxytocin in humans: context and person matter. Trends Cogn Sci. 15, 301–309, doi:S1364-6613(11)00083-0 (2011). [PubMed]
37. Matsushita H, et al. Oxytocin mediates the antidepressant effects of mating behavior in male mice. Neurosci Res. 2010;68:151–153. doi: 10.1016/j.neures.2010.06.007. [PubMed] [Cross Ref]
38. Waldherr, M. & Neumann, I. D. Centrally released oxytocin mediates mating-induced anxiolysis in male rats. Proc Natl Acad Sci USA104, 16681–16684, doi:0705860104 (2007). [PubMed]
39. Arletti R, Bertolini A. Oxytocin acts as an antidepressant in two animal models of depression. Life Sci. 1987;41:1725–1730. doi: 10.1016/0024-3205(87)90600-X. [PubMed] [Cross Ref]
40. Chaviaras S, Mak P, Ralph D, Krishnan L, Broadbear JH. Assessing the antidepressant-like effects of carbetocin, an oxytocin agonist, using a modification of the forced swimming test. Psychopharmacology (Berl). 2010;210:35–43. doi: 10.1007/s00213-010-1815-x. [PubMed] [Cross Ref]
41. Norman GJ, et al. Social interaction prevents the development of depressive-like behavior post nerve injury in mice: a potential role for oxytocin. Psychosomatic medicine. 2010;72:519–526. doi: 10.1097/PSY.0b013e3181de8678. [PubMed] [Cross Ref]
42. Grippo AJ, et al. Peripheral oxytocin administration buffers autonomic but not behavioral responses to environmental stressors in isolated prairie voles. Stress. 2012;15:149–161. doi: 10.3109/10253890.2011.605486. [PMC free article] [PubMed] [Cross Ref]
43. Parker, K. J. et al. Preliminary evidence that plasma oxytocin levels are elevated in major depression. Psychiatry Res. 178, 359–362, doi:S0165-1781(09)00362-X (2010). [PMC free article] [PubMed]
44. Scantamburlo G, et al. Plasma oxytocin levels and anxiety in patients with major depression. Psychoneuroendocrinology. 2007;32:407–410. doi: 10.1016/j.psyneuen.2007.01.009. [PubMed] [Cross Ref]
45. Yuen KW, et al. Plasma oxytocin concentrations are lower in depressed vs. healthy control women and are independent of cortisol. J Psychiatr Res. 2014;51:30–36. doi: 10.1016/j.jpsychires.2013.12.012. [PubMed] [Cross Ref]
46. Cyranowski, J. M. et al. Evidence of dysregulated peripheral oxytocin release among depressed women. Psychosom Med. 70, 967–975, doi:70/9/967 (2008). [PMC free article] [PubMed]
47. Tabak, B. A., McCullough, M. E., Szeto, A., Mendez, A. J. & McCabe, P. M. Oxytocin indexes relational distress following interpersonal harms in women. Psychoneuroendocrinology. 36, 115–122, doi:S0306-4530(10)00172-1 (2011). [PMC free article] [PubMed]
48. Rodrigues SM, Saslow LR, Garcia N, John OP, Keltner D. Oxytocin receptor genetic variation relates to empathy and stress reactivity in humans. Proc Natl Acad Sci USA. 2009;106:21437–21441. doi: 10.1073/pnas.0909579106. [PubMed] [Cross Ref]
49. Saphire-Bernstein S, Way BM, Kim HS, Sherman DK, Taylor SE. Oxytocin receptor gene (OXTR) is related to psychological resources. Proc Natl Acad Sci USA. 2011;108:15118–15122. doi: 10.1073/pnas.1113137108. [PubMed] [Cross Ref]
50. Kawamura Y, et al. The association between oxytocin receptor gene (OXTR) polymorphisms and affective temperaments, as measured by TEMPS-A. Journal of Affective Disorders. 2010;127:31–37. doi: 10.1016/j.jad.2010.04.014. [PubMed] [Cross Ref]
51. Costa B, et al. Oxytocin receptor polymorphisms and adult attachment style in patients with depression. Psychoneuroendocrinology. 2009;34:1506–1514. doi: 10.1016/j.psyneuen.2009.05.006. [PubMed] [Cross Ref]
52. McInnis OA, McQuaid RJ, Matheson K, Anisman H. The moderating role of an oxytocin receptor gene polymorphism in the relation between unsupportive social interactions and coping profiles: implications for depression. Front Psychol. 2015;6:1133. doi: 10.3389/fpsyg.2015.01133. [PMC free article] [PubMed] [Cross Ref]
53. Bradley B, et al. Association between childhood maltreatment and adult emotional dysregulation in a low-income, urban, African American sample: Moderation by oxytocin receptor gene. Development and psychopathology. 2011;23:439–452. doi: 10.1017/S0954579411000162. [PMC free article] [PubMed] [Cross Ref]
54. McQuaid RJ, McInnis OA, Stead JD, Matheson K, Anisman H. A paradoxical association of an oxytocin receptor gene polymorphism: early-life adversity and vulnerability to depression. Front Neurosci. 2013;7:128. doi: 10.3389/fnins.2013.00128. [PMC free article] [PubMed] [Cross Ref]
55. Thompson RJ, Parker KJ, Hallmayer JF, Waugh CE, Gotlib IH. Oxytocin receptor gene polymorphism (rs2254298) interacts with familial risk for psychopathology to predict symptoms of depression and anxiety in adolescent girls. Psychoneuroendocrinology. 2011;36:144–147. doi: 10.1016/j.psyneuen.2010.07.003. [PMC free article] [PubMed] [Cross Ref]
56. Myers AJ, et al. Variation in the oxytocin receptor gene is associated with increased risk for anxiety, stress and depression in individuals with a history of exposure to early life stress. J Psychiatr Res. 2014;59:93–100. doi: 10.1016/j.jpsychires.2014.08.021. [PMC free article] [PubMed] [Cross Ref]
57. Bales KL, Carter CS. Sex differences and developmental effects of oxytocin on aggression and social behavior in prairie voles (Microtus ochrogaster) Horm Behav. 2003;44:178–184. doi: 10.1016/S0018-506X(03)00154-5. [PubMed] [Cross Ref]
58. Champagne, F. A. & Meaney, M. J. Transgenerational effects of social environment on variations in maternal care and behavioral response to novelty. Behav Neurosci. 121, 1353–1363, doi:2007-18058-021 (2007). [PubMed]
59. Francis, D. D., Champagne, F. C. & Meaney, M. J. Variations in maternal behaviour are associated with differences in oxytocin receptor levels in the rat. J Neuroendocrinol. 12, 1145–1148, doi:599 (2000). [PubMed]
60. Francis, D. D., Young, L. J., Meaney, M. J. & Insel, T. R. Naturally occurring differences in maternal care are associated with the expression of oxytocin and vasopressin (V1a) receptors: gender differences. J Neuroendocrinol. 14, 349–353, doi:776 (2002). [PubMed]
61. Champagne, F. A. & Meaney, M. J. Stress during gestation alters postpartum maternal care and the development of the offspring in a rodent model. Biol Psychiatry. 59, 1227–1235, doi:S0006-3223(05)01370-3 (2006). [PubMed]
62. Curley JP, Davidson S, Bateson P, Champagne FA. Social enrichment during postnatal development induces transgenerational effects on emotional and reproductive behavior in mice. Front Behav Neurosci. 2009;3:25. doi: 10.3389/neuro.08.025.2009. [PMC free article] [PubMed] [Cross Ref]
63. Ahern TH, Young LJ. The impact of early life family structure on adult social attachment, alloparental behavior, and the neuropeptide systems regulating affiliative behaviors in the monogamous prairie vole (microtus ochrogaster) Front Behav Neurosci. 2009;3:17. doi: 10.3389/neuro.08.017.2009. [PMC free article] [PubMed] [Cross Ref]
64. Lukas, M., Bredewold, R., Neumann, I. D. & Veenema, A. H. Maternal separation interferes with developmental changes in brain vasopressin and oxytocin receptor binding in male rats. Neuropharmacology. 58, 78–87, doi:S0028-3908(09)00174-9 (2010). [PubMed]
65. Bales KL, Boone E, Epperson P, Hoffman G, Carter CS. Are behavioral effects of early experience mediated by oxytocin? Front Psychiatry. 2011;2:24. doi: 10.3389/fpsyt.2011.00024. [PMC free article] [PubMed] [Cross Ref]
66. Todeschin, A. S. et al. Effects of neonatal handling on social memory, social interaction, and number of oxytocin and vasopressin neurons in rats. Horm Behav. 56, 93–100, doi:S0018-506X(09)00058-0 (2009). [PubMed]
67. Winkelmann-Duarte, E. C. et al. Plastic changes induced by neonatal handling in the hypothalamus of female rats. Brain Res. 1170, 20–30, doi:S0006-8993(07)01656-3 (2007). [PubMed]
68. Winslow JT. Neuropeptides and non-human primate social deficits associated with pathogenic rearing experience. Int. J. Dev. Neurosci. 2005;23:245–251. doi: 10.1016/j.ijdevneu.2004.03.003. [PubMed] [Cross Ref]
69. Heim, C. et al. Lower CSF oxytocin concentrations in women with a history of childhood abuse. Mol Psychiatry. 14, 954–958, doi:mp2008112 (2009). [PubMed]
70. Meinlschmidt, G. & Heim, C. Sensitivity to intranasal oxytocin in adult men with early parental separation. Biol Psychiatry. 61, 1109–1111, doi:S0006-3223(06)01113-9 (2007). [PubMed]
71. Bakermans-Kranenburg MJ, van Ijzendoorn MH, Riem MM, Tops M, Alink LR. Oxytocin decreases handgrip force in reaction to infant crying in females without harsh parenting experiences. Soc Cogn Affect Neurosci. 2012;7:951–957. doi: 10.1093/scan/nsr067. [PMC free article] [PubMed] [Cross Ref]
72. Grimm S, et al. Early life stress modulates oxytocin effects on limbic system during acute psychosocial stress. Soc Cogn Affect Neurosci. 2014;9:1828–1835. doi: 10.1093/scan/nsu020. [PMC free article] [PubMed] [Cross Ref]
73. Wismer Fries, A. B., Ziegler, T. E., Kurian, J. R., Jacoris, S. & Pollak, S. D. Early experience in humans is associated with changes in neuropeptides critical for regulating social behavior. Proc Natl Acad Sci USA102, 17237–17240, doi:102/47/17237 (2005). [PubMed]
74. Kiss A, Mikkelsen JD. Oxytocin–anatomy and functional assignments: a minireview. Endocr Regul. 2005;39:97–105. [PubMed]
75. Kusui, C. et al. DNA Methylation of the Human Oxytocin Receptor Gene Promoter Regulates Tissue-Specific Gene Suppression. Biochemical and Biophysical Research Communications. 289, 681–686, doi:10.1006/bbrc.2001.6024 (2001). [PubMed]
76. Gregory SG, et al. Genomic and epigenetic evidence for oxytocin receptor deficiency in autism. BMC Med. 2009;7:62. doi: 10.1186/1741-7015-7-62. [PMC free article] [PubMed] [Cross Ref]
77. Dadds MR, et al. Methylation of the oxytocin receptor gene and oxytocin blood levels in the development of psychopathy. Dev Psychopathol. 2014;26:33–40. doi: 10.1017/S0954579413000497. [PubMed] [Cross Ref]
78. Cecil CA, et al. Environmental risk, Oxytocin Receptor Gene (OXTR) methylation and youth callous-unemotional traits: a 13-year longitudinal study. Mol Psychiatry. 2014;19:1071–1077. doi: 10.1038/mp.2014.95. [PMC free article] [PubMed] [Cross Ref]
79. Kim YR, Kim JH, Kim MJ, Treasure J. Differential methylation of the oxytocin receptor gene in patients with anorexia nervosa: a pilot study. PloS one. 2014;9:e88673. doi: 10.1371/journal.pone.0088673. [PMC free article] [PubMed] [Cross Ref]
80. Unternaehrer E, et al. Dynamic changes in DNA methylation of stress-associated genes (OXTR, BDNF) after acute psychosocial stress. Transl Psychiatry. 2012;2:e150. doi: 10.1038/tp.2012.77. [PMC free article] [PubMed] [Cross Ref]
81. Jack A, Connelly JJ, Morris JP. DNA methylation of the oxytocin receptor gene predicts neural response to ambiguous social stimuli. Frontiers in human neuroscience. 2012;6:280. doi: 10.3389/fnhum.2012.00280. [PMC free article] [PubMed] [Cross Ref]
82. Puglia MH, Lillard TS, Morris JP, Connelly JJ. Epigenetic modification of the oxytocin receptor gene influences the perception of anger and fear in the human brain. Proceedings of the National Academy of Sciences of the United States of America. 2015;112:3308–3313. doi: 10.1073/pnas.1422096112. [PubMed] [Cross Ref]
83. Ziegler C, et al. Oxytocin receptor gene methylation: converging multilevel evidence for a role in social anxiety. Neuropsychopharmacology: official publication of the American College of Neuropsychopharmacology. 2015;40:1528–1538. doi: 10.1038/npp.2015.2. [PMC free article] [PubMed] [Cross Ref]
84. Kimmel M, et al. Oxytocin receptor DNA methylation in postpartum depression. Psychoneuroendocrinology. 2016;69:150–160. doi: 10.1016/j.psyneuen.2016.04.008. [PubMed] [Cross Ref]
85. Reiner I, et al. Methylation of the oxytocin receptor gene in clinically depressed patients compared to controls: The role of OXTR rs53576 genotype. J Psychiatr Res. 2015;65:9–15. doi: 10.1016/j.jpsychires.2015.03.012. [PubMed] [Cross Ref]
86. Simons, R. L., Lei, M. K., Beach, S. R., Cutrona, C. E. & Philibert, R. A. Methylation of the oxytocin receptor gene mediates the effect of adversity on negative schemas and depression. Dev Psychopathol. 1–12, doi:10.1017/S0954579416000420 (2016). [PMC free article] [PubMed]
87. Beery AK, McEwen LM, MacIsaac JL, Francis DD, Kobor MS. Natural variation in maternal care and cross-tissue patterns of oxytocin receptor gene methylation in rats. Horm Behav. 2016;77:42–52. doi: 10.1016/j.yhbeh.2015.05.022. [PMC free article] [PubMed] [Cross Ref]
88. Needham, B. L. et al. Life Course Socioeconomic Status and DNA Methylation in Genes Related to Stress Reactivity and Inflammation: The Multi-Ethnic Study of Atherosclerosis. Epigenetics. 00–00, doi:10.1080/15592294.2015.1085139 (2015). [PMC free article] [PubMed]
89. Unternaehrer, E. et al. Childhood maternal care is associated with DNA methylation of the genes for brain-derived neurotrophic factor (BDNF) and oxytocin receptor (OXTR) in peripheral blood cells in adult men and women. Stress. 1–11, doi:10.3109/10253890.2015.1038992 (2015). [PubMed]
90. Smearman EL, et al. Oxytocin Receptor Genetic and Epigenetic Variations: Association With Child Abuse and Adult Psychiatric Symptoms. Child Dev. 2016;87:122–134. doi: 10.1111/cdev.12493. [PMC free article] [PubMed] [Cross Ref]
91. Saxonov S, Berg P, Brutlag DL. A genome-wide analysis of CpG dinucleotides in the human genome distinguishes two distinct classes of promoters. Proc Natl Acad Sci USA. 2006;103:1412–1417. doi: 10.1073/pnas.0510310103. [PubMed] [Cross Ref]
92. Hon GC, Hawkins RD, Ren B. Predictive chromatin signatures in the mammalian genome. Hum Mol Genet. 2009;18:R195–201. doi: 10.1093/hmg/ddp409. [PMC free article] [PubMed] [Cross Ref]
93. Blackwood EM, Kadonaga JT. Going the distance: a current view of enhancer action. Science. 1998;281:60–63. doi: 10.1126/science.281.5373.60. [PubMed] [Cross Ref]
94. Suderman M, et al. Childhood abuse is associated with methylation of multiple loci in adult DNA. BMC Med Genomics. 2014;7:13. doi: 10.1186/1755-8794-7-13. [PMC free article] [PubMed] [Cross Ref]
95. Ouellet-Morin I, et al. Effects of the MAOA gene and levels of exposure to violence on antisocial outcomes. The British journal of psychiatry: the journal of mental science. 2016;208:42–48. doi: 10.1192/bjp.bp.114.162081. [PubMed] [Cross Ref]
96. Straus, M. A., Hamby, S. L., Boney-McCoy, S., Sugarman, D. B. The Revised Conflict Tactics Scales (CTS2). Development and preliminary psychometric data. J Fam Issues. 17 (1996).
97. Felitti VJ, et al. Relationship of childhood abuse and household dysfunction to many of the leading causes of death in adults. The Adverse Childhood Experiences (ACE) Study. Am J Prev Med. 1998;14:245–258. doi: 10.1016/S0749-3797(98)00017-8. [PubMed] [Cross Ref]
98. Finkelhor, D. Sexually victimized children. (Free Press, 1979).
99. Masse LC, Tremblay RE. Behavior of boys in kindergarten and the onset of substance use during adolescence. Arch Gen Psychiatry. 1997;54:62–68. doi: 10.1001/archpsyc.1997.01830130068014. [PubMed] [Cross Ref]
100. Nagin D, Tremblay RE. Trajectories of boys’ physical aggression, opposition, and hyperactivity on the path to physically violent and nonviolent juvenile delinquency. Child Dev. 1999;70:1181–1196. doi: 10.1111/1467-8624.00086. [PubMed] [Cross Ref]
101. Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74. doi: 10.1038/nature11247. [PMC free article] [PubMed] [Cross Ref]
102. Karolchik, D., Hinrichs, A. S. & Kent, W. J. The UCSC Genome Browser. Curr Protoc Bioinformatics. Chapter 1, Unit1 4, doi:10.1002/0471250953.bi0104s40 (2012).
103. Landt SG, et al. ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 2012;22:1813–1831. doi: 10.1101/gr.136184.111. [PubMed] [Cross Ref]
104. Rouleau J, Tanigawa G, Szyf M. The mouse DNA methyltransferase 5′-region. A unique housekeeping gene promoter. J Biol Chem. 1992;267:7368–7377. [PubMed]
105. Preacher KJ, Hayes AF. SPSS and SAS procedures for estimating indirect effects in simple mediation models. Behav Res Methods Instrum Comput. 2004;36:717–731. doi: 10.3758/BF03206553. [PubMed] [Cross Ref]
106. Hayes AF, Preacher KJ. Statistical mediation analysis with a multicategorical independent variable. Br J Math Stat Psychol. 2014;67:451–470. doi: 10.1111/bmsp.12028. [PubMed] [Cross Ref]
107. Cohen J. A power primer. Psychol Bull. 1992;112:155–159. doi: 10.1037/0033-2909.112.1.155. [PubMed] [Cross Ref]
108. Richardson JTE. Eta squared and partial eta squared as measures of effect size in educational research. Educational Research Review. 2011;6:135–147. doi: 10.1016/j.edurev.2010.12.001. [Cross Ref]
109. Needham, B. L. et al. Life Course Socioeconomic Status and DNA Methylation in Genes Related to Stress Reactivity and Inflammation: The Multi-Ethnic Study of Atherosclerosis. Epigenetics. 0, doi:10.1080/15592294.2015.1085139 (2015). [PMC free article] [PubMed]
110. Unternaehrer E, et al. Childhood maternal care is associated with DNA methylation of the genes for brain-derived neurotrophic factor (BDNF) and oxytocin receptor (OXTR) in peripheral blood cells in adult men and women. Stress. 2015;18:451–461. doi: 10.3109/10253890.2015.1038992. [PubMed] [Cross Ref]
111. Uhl-Bronner S, Waltisperger E, Martinez-Lorenzana G, Condes Lara M, Freund-Mercier MJ. Sexually dimorphic expression of oxytocin binding sites in forebrain and spinal cord of the rat. Neuroscience. 2005;135:147–154. doi: 10.1016/j.neuroscience.2005.05.025. [PubMed] [Cross Ref]
112. Kumsta R, Hummel E, Chen FS, Heinrichs M. Epigenetic regulation of the oxytocin receptor gene: implications for behavioral neuroscience. Front Neurosci. 2013;7:83. doi: 10.3389/fnins.2013.00083. [PMC free article] [PubMed] [Cross Ref]
113. Yehuda R, et al. Lower methylation of glucocorticoid receptor gene promoter 1F in peripheral blood of veterans with posttraumatic stress disorder. Biol Psychiatry. 2015;77:356–364. doi: 10.1016/j.biopsych.2014.02.006. [PubMed] [Cross Ref]
114. Melas PA, et al. Genetic and epigenetic associations of MAOA and NR3C1 with depression and childhood adversities. Int J Neuropsychopharmacol. 2013;16:1513–1528. doi: 10.1017/S1461145713000102. [PubMed] [Cross Ref]
115. Martín-Blanco A, et al. Association between methylation of the glucocorticoid receptor gene, childhood maltreatment, and clinical severity in borderline personality disorder. Journal of Psychiatric Research. 2014;57:34–40. doi: 10.1016/j.jpsychires.2014.06.011. [PubMed] [Cross Ref]
116. Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13:484–492. doi: 10.1038/nrg3230. [PubMed] [Cross Ref]
117. Rubin LH, et al. Sex and diagnosis specific associations between DNA methylation of the oxytocin receptor gene with emotion processing and temporal-limbic and prefrontal brain volumes in psychotic disorders. Biol Psychiatry Cogn Neurosci Neuroimaging. 2016;1:141–151. doi: 10.1016/j.bpsc.2015.10.003. [PMC free article] [PubMed] [Cross Ref]
118. Tseng HC, et al. Sex-specific associations between plasma oxytocin levels and schizotypal personality features in healthy individuals. J Psychiatr Res. 2014;51:37–41. doi: 10.1016/j.jpsychires.2013.12.011. [PubMed] [Cross Ref]
119. Fischer-Shofty M, Levkovitz Y, Shamay-Tsoory SG. Oxytocin facilitates accurate perception of competition in men and kinship in women. Soc Cogn Affect Neurosci. 2013;8:313–317. doi: 10.1093/scan/nsr100. [PMC free article] [PubMed] [Cross Ref]
120. Ditzen B, et al. Sex-specific effects of intranasal oxytocin on autonomic nervous system and emotional responses to couple conflict. Soc Cogn Affect Neurosci. 2013;8:897–902. doi: 10.1093/scan/nss083. [PMC free article] [PubMed] [Cross Ref]
121. Cardoso C, Ellenbogen MA, Linnen AM. Acute intranasal oxytocin improves positive self-perceptions of personality. Psychopharmacology. 2012;220:741–749. doi: 10.1007/s00213-011-2527-6. [PubMed] [Cross Ref]
122. Piccinelli M, Wilkinson G. Gender differences in depression. Critical review. The British journal of psychiatry: the journal of mental science. 2000;177:486–492. doi: 10.1192/bjp.177.6.486. [PubMed] [Cross Ref]
123. Szyf M, McGowan P, Meaney MJ. The social environment and the epigenome. Environ Mol Mutagen. 2008;49:46–60. doi: 10.1002/em.20357. [PubMed] [Cross Ref]
124. Suderman, M. et al. Childhood abuse is associated with methylation of multiple loci in adult DNA. Bmc Med Genomics. 7, doi:Artn 1310.1186/1755-8794-7-13 (2014). [PMC free article] [PubMed]
125. van Weerd JH, Koshiba-Takeuchi K, Kwon C, Takeuchi JK. Epigenetic factors and cardiac development. Cardiovasc Res. 2011;91:203–211. doi: 10.1093/cvr/cvr138. [PMC free article] [PubMed] [Cross Ref]
126. Ebstein RP, Knafo A, Mankuta D, Chew SH, Lai PS. The contributions of oxytocin and vasopressin pathway genes to human behavior. Horm Behav. 2012;61:359–379. doi: 10.1016/j.yhbeh.2011.12.014. [PubMed] [Cross Ref]
127. Teh AL, et al. The effect of genotype and in utero environment on interindividual variation in neonate DNA methylomes. Genome Res. 2014;24:1064–1074. doi: 10.1101/gr.171439.113. [PubMed] [Cross Ref]
128. Chagnon YC, Potvin O, Hudon C, Preville M. DNA methylation and single nucleotide variants in the brain-derived neurotrophic factor (BDNF) and oxytocin receptor (OXTR) genes are associated with anxiety/depression in older women. Front Genet. 2015;6:230. doi: 10.3389/fgene.2015.00230. [PMC free article] [PubMed] [Cross Ref]
129. Bell AF, et al. Interaction between oxytocin receptor DNA methylation and genotype is associated with risk of postpartum depression in women without depression in pregnancy. Front Genet. 2015;6:243. doi: 10.3389/fgene.2015.00243. [PMC free article] [PubMed] [Cross Ref]
130. Lam LL, et al. Factors underlying variable DNA methylation in a human community cohort. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(Suppl 2):17253–17260. doi: 10.1073/pnas.1121249109. [PubMed] [Cross Ref]
131. Mehta D, et al. Childhood maltreatment is associated with distinct genomic and epigenetic profiles in posttraumatic stress disorder. Proceedings of the National Academy of Sciences of the United States of America. 2013;110:8302–8307. doi: 10.1073/pnas.1217750110. [PubMed] [Cross Ref]
132. Levesque ML, et al. Genome-wide DNA methylation variability in adolescent monozygotic twins followed since birth. Epigenetics. 2014;9:1410–1421. doi: 10.4161/15592294.2014.970060. [PMC free article] [PubMed] [Cross Ref]
133. Wang D, et al. Peripheral SLC6A4 DNA methylation is associated with in vivo measures of human brain serotonin synthesis and childhood physical aggression. PloS one. 2012;7:e39501. doi: 10.1371/journal.pone.0039501. [PMC free article] [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group