Edited by: Adriana Lori, Emory University, United States
Reviewed by: Andrew F. Russo, The University of Iowa, United States; Gita A. Pathak, Yale University, United States
This article was submitted to Neurogenomics, a section of the journal Frontiers in Neuroscience
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Epigenetic mechanisms have been associated with genes involved in Posttraumatic stress disorder (PTSD). PTSD often co-occurs with other health conditions such as depression, cardiovascular disorder and respiratory illnesses. PTSD and migraine have previously been reported to be symptomatically positively correlated with each other, but little is known about the genes involved. The aim of this study was to understand the comorbidity between PTSD and migraine using a monozygotic twin disease discordant study design in six pairs of monozygotic twins discordant for PTSD and 15 pairs of monozygotic twins discordant for migraine. DNA from peripheral blood was run on Illumina EPIC arrays and analyzed. Multiple testing correction was performed using the Bonferroni method and 10% false discovery rate (FDR). We validated 11 candidate genes previously associated with PTSD including
香京julia种子在线播放
Post-traumatic stress disorder (PTSD) is a debilitating, stress-related psychiatric condition, which occurs among persons exposed to traumatic events involving life threats, serious injury, or death (
Epigenome-wide Association Studies (EWAS) provide an unbiased approach to identify DNA methylation loci associated with a disease independently of what is known about the pathophysiology (
A larger number of studies have investigated DNA methylation of specific candidate genes, including genes involved in the regulation of the HPA axis (
Post-traumatic stress disorder is associated with the occurrence of multiple comorbidities including depression and coronary heart disease which have been well studied (
Migraine may be aggravated by stress, exercise, sleep deficiency, hormonal changes, head traumas, major depression, PTSD and environmental cues (
It is important to understand the role genetic, environmental and epigenetic factors play in determining an individual’s susceptibility to PTSD and other co-occurring symptoms (
The study aimed to identify PTSD associated genes overlapping with migraine and evaluate similarity in biological pathways between these disorders. To the best of our knowledge this is the first study to use the monozygotic twin design to determine DNA methylation differences between twins discordant for PTSD and migraine.
The study was based on a subset of six pairs of monozygotic (MZ) twins discordant for PTSD and 15 pairs of MZ twins discordant for migraine that were part of a larger cohort of twins recruited by the QIMR Berghofer Medical Research Institute (
The MZ twins discordant for PTSD were part of a larger Missouri Alcoholism Research Center Project 7 (MARC7) study investigating the effects of early experiences and alcohol use in twins, siblings and their spouses in 2010–2013. e-Trauma is a semi-structured interview, which capitalizes on prior research in psychiatric epidemiology and is based on items previously validated by other research interviews, including SCID (Structured Clinical Interview for DSM Disorders) and DSM-IV (
The MZ twins discordant for migraine were part of the 25-UP (
Blood samples were collected from all participants. Purified DNA was quantified on a Qubit Fluorometer (Thermo Fisher Scientific, United States) and for each sample 500 ng was bisulphite-converted using EZ DNA Methylation Kits (Zymo Research, United States). The samples were then assayed for genome-wide DNA methylation levels using Illumina EPIC DNA methylation arrays that offer a high coverage of CpGs >850,000 CpG sites at single-nucleotide resolution, covering all known genes (96% Refseq genes). All procedures were performed according to the manufacturer’s protocol, and arrays were scanned on an Illumina HiScan (Illumina, United States) at the Genomics Research Centre, QUT.
Raw scan data from the Illumina EPIC arrays were exported into R (V
The power of a sample was determined using a traditional power calculator based on 50 simulations and expected target delta ranging from 0.20–0.50 at 5% FDR (
To test whether the overlap of genes between two analyses was more than expected by chance, enrichment testing was performed using 1,000 permutations (using random sets) and applying a two-sided Binominal test in R to give a
To assess the biological and molecular mechanisms in PTSD and migraine, pathway and gene set analyses were performed using the KEGG pathway analysis tool via the Webgestalt interface (
The study comprized of a total of 21 pairs of MZ twins (
Demographics of the study participants including six pairs of MZ twins discordant for PTSD and 15 pairs of MZ twins discordant for migraine.
Status | Age; mean (SD) | Sex; |
|
43.3 | 5 | 1 | |
43.3 | 5 | 1 | |
23.4 | 7 | 8 | |
22.6 | 7 | 8 | |
We investigated 60 PTSD candidate genes that had been previously reported to be associated with PTSD in different studies (
Of the 4411 CpGs tested, 440 CpGs from 54 genes were significantly associated with PTSD (
PTSD Candidate genes with at least one CpG significant at gene-wise Bonferroni threshold for significance.
Gene symbol | Number of CpGs tested | ≧1 CpG with |
Survive bonferroni |
DOCK2 | 96 | YES (6) | YES |
SLC6A3 | 81 | YES (8) | YES |
DICER1 | 49 | YES (5) | YES |
DRD2 | 41 | YES (4) | YES |
ADCYAP1 | 40 | YES (3) | YES |
STMN1 | 38 | YES (1) | YES |
ADCYAP1R1 | 36 | YES (6) | YES |
SLC6A4 | 31 | YES (2) | YES |
SKA2 | 25 | YES (4) | YES |
OXTR | 22 | YES (3) | YES |
DBH | 18 | YES (3) | YES |
ZNF626 | 9 | YES (1) | YES |
TRAIP | 22 | YES (2) | YES |
TSNARE1 | 215 | YES (27) | YES |
IMMP2L | 96 | YES (10) | YES |
Boxplots for significant PTSD candidate genes: Box plots of DNA methylation level beta-values at significant candidate genes for PTSD versus non-PTSD (control) MZ twins that were significant at gene-wise Bonferroni correction. The mean DNA methylation is higher in PTSD than in controls for two of the three CpGs for DICER and DOCK2, while ADCYAP1 has only one CpG with a higher mean DNA methylation in individuals with PTSD compared to control twins. The remainder of the CpGs have a higher mean DNA methylation in controls than with PTSD.
Next, we performed a hypothesis-free epigenome-wide association analysis to test for DNA methylation differences between the MZ twins discordant for PTSD across all 866K CpG sites on the array. Association testing was performed using lmer models and adjusting for cell counts, age and sex in the model as possible confounders.
A total of seven CpGs were significant at the Bonferroni level of significance of
List of CpG sites significantly associated with PTSD at Bonferroni threshold for significance (
Cpg | Chromosome | Basepaira | Gene symbol (closest gene) | CpG location to gene | Direction PTSD | |
cg26483669 | 2.47E-10 | 2 | 113717037 | IL37 | Intergenic | Downregulated |
cg06405715 | 9.90E-09 | 18 | 77895056 | ADNP2 | Gene body | Upregulated |
cg26433527 | 8.45E-09 | 1 | 110370709 | CSF1 | Intergenic | Upregulated |
cg26575738 | 6.11E-09 | 17 | 44896168 | WNT3 | TSS200 | Upregulated |
cgl 1037719 | 4.72E-08 | 6 | 2999695 | NQ02 | TSS1500 | Upregulated |
cgl 1432275 | 3.46E-08 | 4 | 3239323 | HTT | Gene body | Upregulated |
cgl3341380 | 4.15E-08 | 17 | 33701529 | SLFN11 | TSS 1500 | Upregulated |
Manhattan plot of CpGs associated with PTSD: the blue line indicates Bonferroni threshold of significance (<5.77 × 10– 8) and the red line indicates 10% FDR. A total of 7 CpGs were significant at Bonferroni threshold and 1,585 CpGs were significant at 10% FDR.
Box plots of most significant genes: boxplots of genes significantly associated with PTSD (
Posttraumatic stress disorder often occurs with other comorbidities, including migraine. As little is known about the basis for this comorbidity, we aimed to assess the occurrence of DNA methylation variation on PTSD candidate genes in association with migraine. Using the same study design as the PTSD sample, we investigated 15 pairs of MZ twins discordant for migraine. We tested for association of methylation at genetic loci that overlap between PTSD and migraine using several different analyses.
First, we tested 2,569 CpGs for the 60 PTSD candidate genes identified from earlier studies among the twins with migraine. There were 368 CpGs across 51 genes that had at least one significant CpG associated with migraine (
PTSD Candidate genes also associated with migraine in MZ twins with at least one Bonferroni significant CpG site.
Gene symbol | Number of CpGs tested | ≧1 CpG with |
Survive bonferroni |
ADCYAP1 | 40 | YES (4) | YES |
AIM2 | 16 | YES (1) | YES |
CRHR1 | 70 | YES (7) | YES |
DBH | 18 | YES (2) | YES |
DOCK2 | 96 | YES (8) | YES |
FKBP5 | 53 | YES (9) | YES |
HTR3A | 20 | YES (4) | YES |
OXTR | 22 | YES (3) | YES |
RORA | 238 | YES (21) | YES |
WWC1 | 77 | YES (5) | YES |
TSNARE1 | 215 | YES (13) | YES |
At the epigenome-wide level, we assessed how many of the 1036 genes (1,453 CpGs) associated with PTSD at 5 and 10% FDR overlapped with those significantly associated with migraine in the discordant migraine MZ twins. At 5% FDR, DNA methylation of 13 CpGs (six genes) and at 10% FDR DNA methylation of 99 genes (132 CpGs loci) associated with PTSD was also associated with migraine (
Overlap of PTSD and migraine: Venn diagram of CpGs and genes significantly associated with PTSD (10% FDR) and also associated with migraine (
Boxplots of overlapping genes: box plots of PTSD associated genes overlapping in migraine. The mean DNA methylation is higher in PTSD than in controls for all the overlapping genes in migraine except for
To the best of our knowledge there has only been one study so far investigating comprehensive epigenome-wide DNA methylation changes in migraine performed by
In the migraine MZ twins, for 46 genes out of 62, there was at least one CpG (
In the PTSD MZ twins, we identified that 216 out of 2351 tested CpGs (4 genes i.e.,
Pathway analysis was performed to determine the biological and molecular function associated with PTSD only and both migraine and PTSD using the KEGG pathway through an online interface (
Biological pathways overrepresented among genes associated among genes associated with PTSD (only) and genes associated with both PTSD and migraine.
PTSD Migraine overlap |
|||||
Pathway | Number of genes | Enrichment ratio | FDR | Genes within pathway | |
AMPK signaling pathway | 5 | 9.7252 | 1.40E-04 | 0.0433 | AKT3; CREB5; IRS1; MLYCD; TSC2 |
Longevity regulating pathway | 4 | 10.4901 | 5.25E-04 | 0.0811 | AKT3; CREB5; IRS1; TSC2 |
Vascular smooth muscle contraction | 17 | 3.4978 | 5.86E-06 | 0.0017 | ACTA2; ADCY2; ADCY9; ARHGEF12; CACNA1C; CACNA1S; CALD1; GNA13; KCNMA1; MRVI1; MYH11; MYL6B; MYLK2; MYLK3; NPR2; PLA2G6; PPP1R12A |
Axon guidance | 19 | 2.7030 | 7.05E-05 | 0.0107 | ABL1; ABLIM2; ARHGEF12; CAMK2B; CXCL12; EFNA5; ENAH; EPHA4; EPHA5; NGEF; NTN1; PAK6; PARD3; PLCG1; PLCG2; PPP3CB; PRKCZ; ROBO2; TRPC1 |
Oxytocin signaling pathway | 17 | 2.7844 | 1.18E-04 | 0.0119 | ADCY2; ADCY9; CACNA1C; CACNA1S; CACNB2; CAMK1D; CAMK2B; KCNJ5; CACNA2D3; MAP2K5; MYL6B; MYLK2; MYLK3; NPR2; PPP1R12A; PPP3CB; RYR3 |
Leukocyte transendothelial migration | 13 | 2.8897 | 5.18E-04 | 0.0396 | ACTN1; CLDN1; CTNNA1; CTNNA2; CTNND1; CXCL12; JAM2; MYL12B; PLCG1; PLCG2; RAPGEF4; SIPA1; THY1 |
Phospholipase D signaling pathway | 15 | 2.5578 | 7.34E-04 | 0.0449 | ADCY2; ADCY9; AGPAT1; AGPAT3; AGPAT4; ARF1; DGKH; GNA13; PDGFC; PLCG1; PLCG2; RALB; RAPGEF4; SHC1; TSC2 |
PTSD often co-occurs with other disorders and migraine being one of the less-studied comorbidities of PTSD. In the current study, we analyzed epigenome-wide DNA methylation data in MZ twins discordant for PTSD and migraine to identify common genes and pathways in PTSD and migraine. DNA methylation can be affected by both genetic and environmental factors. Therefore, the MZ twin’s disease-discordant design is a powerful approach to dissect overlap between disorders as the participants are genetically and demographically matched (
Using this unique study design, we first investigated six pairs of MZ twins that were all exposed to stress but were discordant for PTSD. Genome-wide analyses of DNA methylation differences across the twins identified 7 CpG loci that were significantly associated with PTSD even after stringent Bonferroni correction for multiple testing. The top CpG was in an intergenic region on chromosome 2, located near gene
Other genome-wide significant genes included
When assessing 45 candidate genes known to be involved in PTSD through other studies, we found that DNA methylation of these genes were significantly more likely to be associated with PTSD in the current study than expected by random chance. Of the PTSD candidate genes, we found 11 genes were also associated with PTSD in the current twin study after gene-wise Bonferroni correction (see
Posttraumatic stress disorder and chronic pain conditions often occur together, and the underlying mechanisms are diverse and multifactorial (
Another well-studied PTSD candidate is Catechol-O-methyltransferase (
When assessing biological pathways implicated in genes associated with PTSD and migraine, we found that several different pathways were involved in genes implicated in PTSD only versus those implicated in PTSD and migraine. For instance, vascular smooth muscle, axon guidance and oxytocin signaling pathways were overrepresented for genes associated with PTSD; all these are plausible and well-known pathways for PTSD. Interestingly, migraine GWAS studies have also identified an abundance of vascular genes associated with migraine (
Epigenetic regulation of gene expression is a dynamic and reversible process and therefore a good pathway target with drugs (
There are several limitations of this study including small sample size (42 samples) which has limited power to detect small DNA methylation changes. Given the unique study design of genetically and demographically matched samples, this is a small yet homogenous sample to detect changes in disease. Another limitation of the study is the use of peripheral blood to assess DNA methylation in PTSD and migraine and while this has been a topic for discussion for many years, it is now evident that whole blood acts as a relevant and easily accessible surrogate to investigate brain-disorders. We also used an online tool by Hannon et al. to investigate the correlation of DNA methylation in blood for the significant genes with four brain regions (
Nevertheless, we were able to validate and replicate several previously reported findings from PTSD in this study and this is the first study of its kind to use a MZ-twin design to disentangle the genes overlapping in PTSD and migraine. Using the disease discordant MZ twin design is a powerful approach in EWAS as the participants are genetically matched with similar environmental exposure and lifestyle especially in the earlier years (
There are several proposed mechanisms underlying the psychiatric comorbidities and migraine relationship, but their exact etiology and biological mechanisms are not entirely known (
To the best of our knowledge, this is the first study to examine overlapping genes and biological processes in PTSD and migraine comorbidity using the monozygotic co-twin design. These results are important and suggest that common genes and pathways might be associated with PTSD and migraine, with implications for diagnosis of comorbidities in PTSD and common treatments for these co-morbid disorders.
The datasets presented in this study can be found in online repositories. The names of the repository and accession number can be found below: Gene Expression Omnibus
The studies involving human participants were reviewed and approved by QIMR Berghofer Medical Research Institute and Queensland University of Technology (QUT) Human Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.
CB analyzed the data and wrote the first draft of the manuscript with DM support and guidance. NM, EN, AH, IH, and MW led studies and provided data from QIMR Berghofer Medical Research Institute. HS and RP conducted the laboratory sample analysis. LH and LG contributed to head of the Genomics Research Centre at QUT and supervised all experimental protocols. AM provided statistical advice and input. DM and DN conceived the study and supervised the project. All authors approved the manuscript and provided critical comments on the manuscript.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank all the participants involved in this study.
The Supplementary Material for this article can be found online at:
Correlation of DNA methylation in blood and brain for MEG8 gene.
A list of all PTSD Candidate genes tested in the MZ twins.
List of genes significantly associated with PTSD (10% FDR) and overlapping with migraine in the MZ twins.
A list of all PTSD Candidate genes associated with migraine in the MZ migraine twins.