Front. Neurosci. Frontiers in Neuroscience Front. Neurosci. 1662-453X Frontiers Media S.A. 10.3389/fnins.2021.678350 Neuroscience Original Research Using Monozygotic Twins to Dissect Common Genes in Posttraumatic Stress Disorder and Migraine Bainomugisa Charlotte K. 1 Sutherland Heidi G. 1 2 Parker Richard 3 Mcrae Allan F. 4 Haupt Larisa M. 1 2 Griffiths Lyn R. 1 2 Heath Andrew 5 Nelson Elliot C. 5 Wright Margaret J. 6 7 Hickie Ian B. 8 Martin Nicholas G. 3 Nyholt Dale R. 1 Mehta Divya 1 * 1Centre for Genomics and Personalised Health, School of Biomedical Science, Faculty of Health, Queensland University of Technology, Kelvin Grove, QLD, Australia 2Centre for Genomics and Personalised Health, Genomics Research Centre, School of Biomedical Sciences, Institute of Health and Biomedical Innovation, Kelvin Grove, QLD, Australia 3QIMR Berghofer Medical Research Institute, Royal Brisbane Hospital, Herston, QLD, Australia 4Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD, Australia 5Department of Psychiatry, Washington University School of Medicine in St. Louis, St. Louis, MO, United States 6Queensland Brain Institute, The University of Queensland, Brisbane, QLD, Australia 7Centre for Advanced Imaging, The University of Queensland, Brisbane, QLD, Australia 8Brain and Mind Centre, The University of Sydney, Sydney, NSW, Australia

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

*Correspondence: Divya Mehta, divya.mehta@qut.edu.au

This article was submitted to Neurogenomics, a section of the journal Frontiers in Neuroscience

22 06 2021 2021 15 678350 09 03 2021 31 05 2021 Copyright © 2021 Bainomugisa, Sutherland, Parker, Mcrae, Haupt, Griffiths, Heath, Nelson, Wright, Hickie, Martin, Nyholt and Mehta. 2021 Bainomugisa, Sutherland, Parker, Mcrae, Haupt, Griffiths, Heath, Nelson, Wright, Hickie, Martin, Nyholt and Mehta

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 DOCK2, DICER1, and ADCYAP1. In the epigenome-wide scan, seven novel CpGs were significantly associated with PTSD within/near IL37, WNT3, ADNP2, HTT, SLFN11, and NQO2, with all CpGs except the IL37 CpG hypermethylated in PTSD. These results were significantly enriched for genes whose DNA methylation was previously associated with migraine (p-value = 0.036). At 10% FDR, 132 CpGs in 99 genes associated with PTSD were also associated with migraine in the migraine twin samples. Genes associated with PTSD were overrepresented in vascular smooth muscle, axon guidance and oxytocin signaling pathways, while genes associated with both PTSD and migraine were enriched for AMPK signaling and longevity regulating pathways. In conclusion, these results suggest that common genes and pathways are likely involved in PTSD and migraine, explaining at least in part the co-morbidity between the two disorders.

DNA methylation twins posttraumatic stress disorder migraine genes

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      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 (Matosin et al., 2017). It develops as a result of failure to contain the normal stress response hence dysregulation of the hypothalamic–pituitary–adrenal (HPA) axis which is one of the body’s major stress response systems (Wong, 2002; Yehuda and LeDoux, 2007). Although nearly 90% of individuals are exposed to a traumatic event during their lifetime (Breslau et al., 1998), only a small proportion go on to develop PTSD (Matosin et al., 2017). There is significant variation in the risk of PTSD among individuals experiencing the same trauma; this is likely determined by genetic predisposition and epigenetic mechanisms (Kessler et al., 2017). One of the best understood epigenetic mechanisms is DNA methylation through the addition of a methyl chemical bond to the cytosine C5 in the cytosine-phosphate-guanine (CpG) dinucleotides of the DNA, influencing gene activity and considered to be a stable epigenetic mark in post-mitotic cells (Guo et al., 2011). Generally, DNA methylation, especially if occurring close to the promoter region of the gene is known to inhibit gene transcription.

      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 (Almli et al., 2014). Most EWAS examining epigenetic effects of trauma exposure use DNA methylation profiles in peripheral tissues like blood and saliva, which are useful for biomarker identification (Almli et al., 2014). Several genes previously associated with stress or epigenetic regulation of neuronal function have been identified through EWAS of PTSD (Melroy-Greif et al., 2017; Nievergelt et al., 2018; Howie et al., 2019). Additionally, multiple genome-wide studies have also uncovered genetic contributions to PTSD risk and symptomatology through candidate genes which have encoded for varied proteins (Gormley et al., 2016; Kilaru et al., 2016; Blacker et al., 2019). EWAS show that there are significant epigenetic differences between individuals with PTSD compared to healthy controls (Uddin et al., 2010; Smith et al., 2011), particularly in genes involved in inflammation, immune and nervous system function (Uddin et al., 2010; Smith et al., 2011; Rusiecki et al., 2012). Psychosocial stress may alter global and gene-specific DNA methylation patterns potentially associated with peripheral immune dysregulation (Uddin et al., 2010; Smith et al., 2011; Rusiecki et al., 2012).

      A larger number of studies have investigated DNA methylation of specific candidate genes, including genes involved in the regulation of the HPA axis (Almli et al., 2014). Epigenetic DNA methylation changes may accompany lifetime experiences and alter gene expression profiles (Yehuda et al., 2015). Significant methylation changes in early life, specifically in genes implicated in developing severe psychiatric conditions including DLG4, DRD2, NOS1, NRXN1, and SOX10 also indicate vulnerability to the effects of stress and psychiatric disorders via epigenetic mechanisms (Numata et al., 2012). Large longitudinal twin and molecular genetic cohort studies suggest that the impact of adverse life events is probably moderated by genetic variants through genetic and environmental interactions (Meaney and Szyf, 2005; Probst et al., 2009).

      Post-traumatic stress disorder is associated with the occurrence of multiple comorbidities including depression and coronary heart disease which have been well studied (Sareen, 2014). PTSD is also highly comorbid with chronic pain conditions that often co-occur such as migraine headaches, tension headaches, temporomandibular disorder, irritable bowel syndrome, fibromyalgia, chronic fatigue syndrome and chronic prostatitis/chronic pelvic pain syndrome (Gasperi et al., 2021). Patients with PTSD have a risk of developing pain disorders, which may produce long-lasting changes in the threshold for migraine attacks by inducing epigenetic modifications throughout the brain (Eising et al., 2013). Although several epidemiological studies have reported that PTSD is a predictor of migraine and is much more prevalent in patients with migraine than in the general population (Peterlin et al., 2008; Smitherman et al., 2009; Minen et al., 2016; Zarei et al., 2016), migraine is an understudied comorbidity of PTSD.

      Migraine may be aggravated by stress, exercise, sleep deficiency, hormonal changes, head traumas, major depression, PTSD and environmental cues (Theeler et al., 2012; Guglielmetti et al., 2020). In addition, migraine-related pain may cause sensitisation of certain pain pathways via inflammation-induced changes in epigenetic gene regulation (Eising et al., 2013). A recent GWAS reported that migraine showed a higher genetic correlation with psychiatric disorders when compared to other neurological disorders, suggesting common genetic basis or pathways (Anttila et al., 2018). Little is understood about the link between PTSD and migraine disorders, however, it has been suggested that several systems such as the immune system are likely to be involved in the co-occurrence of these disorders (Peterlin et al., 2011). Biological, environmental and genetic risk factors may converge to produce a brain state which predisposes an individual to both PTSD and migraine (Antonaci et al., 2011). Causal pathways shared between migraine and its comorbid disorders may be modulated by epigenetic mechanisms (Eising et al., 2013); these have been suggested to play a role in development of both disorders (McGowan, 2013).

      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 (Blacker et al., 2019). In this study we used the disease-discordant monozygotic (MZ) twin design to investigate PTSD-migraine comorbidity. This is a powerful strategy in genetic and epigenetic epidemiology as participants are genetically identical and well-matched by age, sex, maternal environment, population cohort effects and exposure to many shared environmental factors. Recent studies have used this design and uncovered considerable epigenetic (methylation) variation between MZ twins for several complex phenotypic traits which are detectable in blood DNA samples (Bell and Spector, 2011; Bolund et al., 2017; Kaut et al., 2017; Gerring et al., 2018; Peng et al., 2018; Gasperi et al., 2021).

      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.

      Materials and Methods Participants and Samples

      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 (Wright and Martin, 2004). All twins were of Caucasian descent. Informed written consent was obtained from each participant. All questions were administered using a computer-administered telephone interview (for PTSD). In addition, respondents were requested to complete a brief questionnaire either online or in person which included questions on physical health, personality and other measures (for PTSD and migraine studies). Ethical clearance for this study was obtained through QIMR Berghofer Medical Research Institute along with Queensland University of Technology (QUT) Human Research Ethics Committee approval.

      Assessment of PTSD

      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 (Bell, 1994). All twins had experienced a PTSD-qualifying potentially traumatic event as per the DSM-IV criteria. PTSD was assessed via the DSM-IV criteria through structured interview questions that were asked over the phone by an experienced interviewer. PTSD diagnosis was based on the self-reported DSM-IV criteria.

      Assessment of Migraine

      The MZ twins discordant for migraine were part of the 25-UP (Mitchell et al., 2019) and Memory attention and problem solving (MAPS) study (Wright and Martin, 2004) at QIMR. Migraine were assessed using the International Headache Society (IHS) diagnostic criteria (the International Classification of Headache Disorders, ICHD-3) together with a diagnosis of migraine with or without aura [Headache Classification Committee of the International Headache Society (IHS), 2013]. For the collection of detailed ICHD-3 diagnostic criteria, participants answering “yes” to ever having “migraine or recurrent attacks of headache” (screening positive), then answered a number of questions relating to their symptoms. Diagnoses were determined for the two major varieties of migraine: migraine without aura and migraine with aura (Launer et al., 1999). MZ twins discordant for migraine with aura were selected for the study.

      DNA Methylation Microarray Analyses

      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.

      Statistical Analysis

      Raw scan data from the Illumina EPIC arrays were exported into R (V4.0.2) for statistical analysis. Samples with probe detection call rates <95% and those with an average intensity value of either <50% of the experiment-wide sample mean or <2000 arbitrary units (AU) were excluded from further analysis. The raw DNA methylation beta values were background and control-normalized using the Bioconductor MINFI package (1.4.0) (Aryee et al., 2014). Cell counts were analyzed using the Houseman method (Houseman et al., 2012). All samples were run in a single batch and were of Caucasian ethnicity. Epigenome-wide differential methylation analysis between case and control groups was performed using linear mixed effects models in R (lmer) to account for the twin pairs and adjusting for cell-counts, age and sex. To correct for multiple testing and identify significant CpGs, we calculated a stringent Bonferroni threshold for significance (P = 5.77 × 10–8) and a less stringent threshold of 10% false discovery rate (FDR).

      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 (Graw et al., 2019), with a classical empirical power of 66 and 83% to detect differentially methylated DNA methylation, respectively.

      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 p-value of enrichment.

      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 (Wang et al., 2013) to identify enriched pathways using a hypergeometric test for enrichment evaluation analysis and significance level of 10% FDR.

      Results Demographics

      The study comprized of a total of 21 pairs of MZ twins (n = 42). Of these, six pairs of MZ twins were discordant for PTSD while 15 pairs of MZ twins were discordant for migraine. Of the six pairs of MZ twins discordant for PTSD, five pairs were male while one pair was female, with a mean age of 43.3 years [SD = 8.58] across the pairs. Of these MZ twins, 10 of them were male while two were female. Of the 15 pairs of MZ twins discordant for migraine, seven pairs were male while eight pairs were female, with a mean age of 23 years [SD = 9.95] across the pairs. Table 1 shows details of the study participants.

      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; n
      Male Female

      PTSD (n = 6) 43.3 5 1
      Control (n = 6) 43.3 5 1
      Total 43.3 (±8.58) 10 2
      Migraine (n = 15) 23.4 7 8
      Control (n = 15) 22.6 7 8
      Total 23 ± 9.95 14 16
      PTSD Candidate Genes

      We investigated 60 PTSD candidate genes that had been previously reported to be associated with PTSD in different studies (Cornelis et al., 2010; Almli et al., 2014; Mehta et al., 2019; Bian et al., 2020; Stein et al., 2021), including genes uncovered in a GWAS of over 250,000 participants in the Million Veteran Program (Stein et al., 2021). For these 60 candidate genes, there were 4411 CpGs detected in the current study on the EPIC microarray that could be tested. We tested for association between DNA methylation and PTSD status in the twin pairs using linear mixed effects (lmer) models with age, sex and cell counts as covariates.

      Of the 4411 CpGs tested, 440 CpGs from 54 genes were significantly associated with PTSD (p < 0.05). There were 14 genes with at least one CpG surviving the locus-specific Bonferroni correction, including DOCK2, SLC6A3, DICER1, DRD2, ADCYAP1, ADCYAP1R1, SKA2, OXTR, STMN1, SLC6A4, DBH, ZNF626, TRAIP, TSNAIRE1, and IMMP2L. Details of the top candidate genes including the number of CpGs significant within each gene are shown in Table 2. Based on 1,000 permutations, this indicates significant overlap above what would be expected by chance (enrichment p-value = 0.0098). Full results are shown in Supplementary Table 1. Figure 1 illustrates DNA methylation differences between PTSD and non-PTSD for three of the top candidate genes.

      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 p ≦ 0.05 (No of CpGs with p ≦ 0.05) 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.

      Genome Wide Analysis for PTSD Status and DNA Methylation

      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 p < 5.77 × 10–8 (Table 3). Full results across all the CpG sites are shown in Figure 2 and Supplementary Table 2, a total of 159 CpGs (115 genes) were significant at 5% FDR. The strongest association was detected for chromosomes 1, 2, 4, 6, 17, and 18. The top CpGs were included IL37 cg26483669, CSF1 cg26433527, ADNP2 cg06405715, WNT3 cg26575738, NQO2 cg11037719, HTT cg11432275, and SLFN11 cg13341380 (Figure 3). The CpG in IL37 showed hypomethylation in the PTSD twins while the other CpGs in CSF1, ADNP2, WNT3, NQO2, HTT, and SLFN11 were hypermethylated in PTSD. As DNA methylation signals across neighboring CpGs can be highly correlated, we used a less conservative threshold of significance (10% FDR) to identify other genes of interest in PTSD. This included 1453 CpG sites in 1036 genes that were significant at 10% FDR. Several of these genes such as HDAC4 and NRG1 were known to be affected by a known drug (n = 84 genes, clinically actionable∗∗) or had genomes that could be used to build new drugs (n = 248 genes, druggable genomes) as per the Drug Gene Interaction Database (Cotto et al., 2017; Supplementary Table 2).

      List of CpG sites significantly associated with PTSD at Bonferroni threshold for significance (p < 5.77 × 10–8).

      Cpg P-value 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
      aPosition in Human GRCh37/hg19 Assembly. TSS, distance from transcription start site.

      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 (p-value < 5.77 × 10– 8) in the MZ twins is depicted. The mean DNA methylation is higher in PTSD than in controls for all the genome wide significant genes except for the IL37 locus probe.

      Overlapping Genes and Pathways in PTSD and Migraine

      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.

      Known PTSD Candidate Genes Associated With Migraine in the MZ Twins

      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 (P < 0.05). Of these genes, 11 genes survived locus specific Bonferroni corrections for multiple testing including ADCYAP1, AIM2, CRHR1, DBH, DOCK2, FKBP5, HTR3A, OXTR, RORA, WWC1, and TSNAIRE1 (Table 4 and Supplementary Table 3).

      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 p ≦ 0.05 (No of CpGs) 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
      Epigenome-Wide Overlap of Genes Associated With PTSD Also Associated With Migraine in the MZ Twins

      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 (p < 0.05, Figure 4). At 10% FDR, 47 CpGs were down-regulated (hypo-methylated) in PTSD of which 25 were also downregulated in migraine while 85 CpGs were up-regulated (hypermethylated) in PTSD of which 53 were also up-regulated in migraine (Supplementary Table 4). The overlapping genes included RERE, MEG8, SPOPL, C1orf187 (DRAXIN), DAPK2, and TM6SF2 (Figure 5).

      Overlap of PTSD and migraine: Venn diagram of CpGs and genes significantly associated with PTSD (10% FDR) and also associated with migraine (p < 0.05).

      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 SPOPL and C1orf187.

      Genes Previously Associated With Migraine DNA Methylation

      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 Gerring et al., 2018. The study had identified a total of 62 genes associated with migraine (Gerring et al., 2018). Using these genes (62 genes across 2,351 CpGs) we tested how many CpGs within these genes were also associated with migraine and PTSD using the MZ twin sample.

      In the migraine MZ twins, for 46 genes out of 62, there was at least one CpG (p < 0.05) significantly associated with migraine in the current study. Of these, for 6 genes (KCNG2, DGKG, SND1, LHX6, ADIRF, and RPTOR), the results survived multiple testing correction at 10% FDR.

      In the PTSD MZ twins, we identified that 216 out of 2351 tested CpGs (4 genes i.e., RPTOR, NUFIP1, SLC38A4, and KCP) were also significantly associated with PTSD (10% FDR) in the current study. Based on 1000 permutations, this overlap is significantly higher than expected by chance alone (enrichment p-value = 0.036), providing support for overlapping genes between migraine and PTSD.

      Functional Annotation of Overlapping Genes

      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 (Liao et al., 2019). For genes associated only with PTSD, vascular smooth muscle, axon guidance and oxytocin signaling pathways were overrepresented, all these are well-known pathways for PTSD. For genes associated with both PTSD and migraine, the overrepresented pathways were AMPK signaling (Enrichment Ratio = 9.73, p-value = 1.40E-04) and longevity regulating (Enrichment Ratio = 10.49; p-value = 5.25E-04). Detailed results are shown in Table 5.

      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 P-value 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
      PTSD only
      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
      Discussion

      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 (Bell and Spector, 2011; Tan et al., 2015).

      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 IL37. This gene has not been directly implicated in PTSD and is an anti-inflammatory cytokine that has been found to be increased in the amygdala and dorsolateral prefrontal cortex of children with autism spectrum disorder (Tsilioni et al., 2019). The cg26433527 site is upstream of the CSF1 gene, a cytokine that controls the production, differentiation, and function of macrophages.

      Other genome-wide significant genes included WNT3, ADNP2, SLFN11 and HTT. WNT3 has been associated with stress-induced depression-like behaviors (Zhou et al., 2016) while ADNP2 has been suggested to cause changes in cellular viability under oxidative stress (Kushnir et al., 2008). SLFN11 acts as a global regulator of chromatin structure with the potential to engage the innate immune activation in response to replicative stress (Murai et al., 2020). The HTT is the Huntingtin gene involved in Huntington’s disease, a neurodegenerative disorder characterized by loss of striatal neurons. This is in line with our previous finding of an involvement of DOCK2 in PTSD, a gene which has also been implicated in the formation of amyloid plaques in the brain in Alzheimer’s disease (Mehta et al., 2017), suggesting the role of genes common to both PTSD and neurodegenerative disorders.

      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 Table 2). The candidate genes included DOCK2, SLC6A3, DICER1, DRD2, ADCYAP1, ADCYAP1R1, SKA2, OXTR, STMN1, SLC6A4, DBH, ZNF626, TRAIP, TSNAIRE1, and IMMP2L. Our team has previously found DOCK2, an amyloid-plaque associated gene in Alzheimer’s, to be associated with PTSD in an epigenome-wide study in Australian veterans (Mehta et al., 2017) and here have validated this gene in the MZ PTSD twins. Another PTSD candidate gene SLC6A3 encodes the dopamine transporter, individuals with PTSD carrying the SLC6A3 9-repeat allele were found to be at higher risk for PTSD when also having higher methylation in the SLC6A3 promoter locus (Chang et al., 2012). SLC6A4 encodes the serotonin transporter. A study reported increased methylation of SLC6A4 in bullied twins at age 10 compared to their non-bullied monozygotic twins (Ouellet-Morin et al., 2013). In genome wide studies, reduced SLC6A4 methylation levels were associated with more traumatic events and increased risk for PTSD in individuals carrying a specific SLC6A4 risk allele genotype while higher SLC6A4 methylation appeared protective against the development of PTSD (Koenen et al., 2011). The contrary results of SLC6A4 are likely driven by a gene-by-environment interaction whereby the particular genotype interacts with environment (trauma) to induce DNA methylation changes and confer susceptibility towards PTSD. Other stress-related candidate genes include SKA2 and ADCYAP1R1; a study showed associations between the methylation status, a polymorphic site in the 3′UTR of the SKA2 gene—involved in mitosis—with reduced thickness of several cortical areas and symptom severity in PTSD (Daskalakis and Yehuda, 2014). Methylation in peripheral blood samples of ADCYAP1 (PACAP protein) and its receptor ADCYAP1R1 (PAC1 protein), genes involved in regulating the cellular stress response, were associated with PTSD diagnosis and symptom severity, specifically in females (Ressler et al., 2011). There is emerging evidence of a link between PACAP and migraine, with clinical trials targeting PACAP as a likely prevention for migraine (Sundrum and Walker, 2018; Rustichelli et al., 2020).

      Posttraumatic stress disorder and chronic pain conditions often occur together, and the underlying mechanisms are diverse and multifactorial (Smitherman et al., 2009; Minen et al., 2016; Gasperi et al., 2021). To investigate genes involved in both PTSD and migraine, we compared genes associated with PTSD in the current study and tested whether the same genes were also associated with migraine in an independent sample of 15 MZ twins discordant for migraine. We identified DAPK2 and TM6SF2 as two of the top overlapping genes between the two disorders. DAPK2 is a calmodulin-regulated protein kinase, it has been implicated in the intracellular degradation process essential for adaptation to metabolic stress (autophagy) (Shiloh et al., 2018). TM6SF2 is associated with cardiovascular disease and plays a role in oxidative stress (Taliento et al., 2019). These findings suggest that epigenetic changes in response to different types of stress may “mediate” stress phenotypes. Stress is a possible underlying cause of both PTSD and migraine and can impact the epigenome (Fuchikami et al., 2009; Schouten et al., 2013; Chakravarty et al., 2014; Kenworthy et al., 2014) and damage the brain in severe cases (McEwen, 2006, 2013). For instance, a single immobilization stress alters hippocampal brain-derived neurotrophic factor (BDNF) gene expression and histone acetylation at BDNF gene promoters (Fuchikami et al., 2009). BDNF is involved in the neural plasticity underlying the extinction of fear (Chhatwal et al., 2006; Heldt et al., 2007; Soliman et al., 2010) and recovery from stress. The BDNF gene has been studied in relation to anxiety disorders such as PTSD (Ressler et al., 2011; Andero and Ressler, 2012). Increased methylation of a CpG site in BDNF and peripheral blood of adults with PTSD have been reported (Smith et al., 2011). Studies have shown that BDNF is increased during migraine attacks, and in cluster headache (Fischer et al., 2012), and there is some genetic evidence suggesting a role of BDNF in migraine (Sutherland et al., 2014).

      Another well-studied PTSD candidate is Catechol-O-methyltransferase (COMT), which encodes enzymes that degrade neurotransmitters such as dopamine and the serotonin transporter SERT. Increased COMT promoter methylation was associated with impaired fear inhibition in individuals with PTSD carrying the COMT met/met genotype (Meyer-Lindenberg and Weinberger, 2006). Polymorphisms in the BDNF, COMT, and SERT lead to disturbances in the normal brain pathway neurotransmitters making patients more susceptible to several neuropsychiatric disorders (Andersen and Skorpen, 2009; Fischer and Jan, 2020). Interestingly, COMT has been associated with headache medication overuse (Fischer and Jan, 2020). These studies further supporting the involvement of similar stress-related genes in the pathophysiology of both PTSD and migraine. Other proposed neurobiological mechanisms underlying PTSD-migraine comorbidity includes dysfunction of the autonomic nervous system, Hypothalamic Pituitary Adrenal (HPA) axis and brain serotonergic dysfunction (Smitherman et al., 2009; Minen et al., 2016).

      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 (Mason and Russo, 2018; van den Maagdenberg et al., 2019; Guo et al., 2020). For genes associated with both migraine and PTSD, AMPK signaling and longevity regulating pathways were overrepresented. Adenosine monophosphate-activated protein kinase (AMPK) pathway activation might be a therapeutic target for PTSD (Wang et al., 2017) and AMPK activators might also be effective for treatment of chronic pain disorders by inhibiting signaling pathways that promote changes in the function of peripheral nociceptive neurons (Asiedu et al., 2016). There is not much known about the role of longevity pathways in PTSD or migraine.

      Epigenetic regulation of gene expression is a dynamic and reversible process and therefore a good pathway target with drugs (Eising et al., 2013). In our study we used the Drug Gene Interaction Database (DGIdb) for all genes associated with PTSD to check if they were known to be affected by a known drug (clinically actionable) or are genes or encode gene products that are known or predicted to be targets for new drugs (i.e., a druggable genome). Of genes significantly associated with PTSD at 10% FDR, we found that 7.4% of these genes were known and clinically actionable targets for drugs while a further 22% of genes had a druggable genome. It has been suggested that defining the exact nature of this association and the pathophysiological mechanisms underlying the comorbidity are relevant in clinical practice as it might influence both the response to treatment and likelihood to achieve remission (Dresler et al., 2019; Seng and Seng, 2016). Current acute and prophylactic treatments are effective in less than half of the patients (Goadsby et al., 2002), indicating the need for more effective drugs. There are currently no epigenetic drugs that have gone into clinical trials for the treatment of neuropsychiatric disorders or migraine with the exception of valproate (Peedicayil and Kumar, 2012; Eising et al., 2013). Identifying factors that predispose to migraine attacks is therefore crucial to provide specific molecular targets to design novel migraine drugs. Understanding which genes are involved in co-morbid disorders such as PTSD and migraine will support the development of more advanced psychological and/or pharmacological interventions that target key pathways underlying both disorders.

      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 (Ligthart et al., 2016). The modest overlap in genes and pathways between PTSD and migraine might be due to several factors including limited power and differences across the two MZ twin samples compared. Given the cross-sectional nature of the DNA methylation measurements, it is difficult to disentangle cause from effect. It is recommended that future studies should use large cohorts and at various time points to identify biological underpinnings of PTSD and migraine.

      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 (Tan et al., 2015). Twin studies provide a useful reference for hypotheses to be replicated and validated to understand epigenetics of complex diseases (Bell and Spector, 2011).

      There are several proposed mechanisms underlying the psychiatric comorbidities and migraine relationship, but their exact etiology and biological mechanisms are not entirely known (Smitherman et al., 2009; Eising et al., 2013). Although there have been promising developments in genome-wide studies on psychopathologies in recent years, data on the biological basis of PTSD and migraine are limited. Few studies with small populations have been conducted to examine genetic and epigenetic (methylation) associations in PTSD (Uddin et al., 2010; Mehta et al., 2018). There is need for comprehensive evaluation and integrated model of care of psychiatric disorders in migraine (Dresler et al., 2019). More research is required to identify epigenetic targets that affect migraine pathophysiology and drugs that specifically act to modulate chromatin structure at migraine pathways (Eising et al., 2013; Minen et al., 2016).

      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.

      Data Availability Statement

      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 GSE172464.

      Ethics Statement

      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.

      Author Contributions

      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.

      Conflict of Interest

      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.

      Funding. We acknowledge Queensland University of Technology for funding the epigenetics study. This QIMR twin study was primarily supported by the National Institute of Health grants to ECN (AA13446; AA011998_5978). Additional support includes grants to ACH (AA10249, AA07728, AA11998, and AA13321) and NM (AA13326). The 25-UP study was funded by NHMRC Project Grant APP1069141 to IH and NM.

      The authors would like to thank all the participants involved in this study.

      Supplementary Material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fnins.2021.678350/full#supplementary-material

      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.

      References Almli L. M. Fani N. Smith A. K. Ressler K. J. (2014). Genetic approaches to understanding post-traumatic stress disorder. Int. J. Neuropsychopharmacol. 17 355370. 10.1017/s1461145713001090 24103155 Andero R. Ressler K. J. (2012). Fear extinction and BDNF: translating animal models of PTSD to the clinic. Genes Brain Behav. 11 503512. 10.1111/j.1601-183x.2012.00801.x 22530815 Andersen S. Skorpen F. (2009). Variation in the COMT gene: implications for pain perception and pain treatment. Pharmacogenomics 10 669684. 10.2217/pgs.09.13 19374521 Antonaci F. Nappi G. Galli F. Manzoni G. C. Calabresi P. Costa A. (2011). Migraine and psychiatric comorbidity: a review of clinical findings. J. Headache Pain 12 115125. Anttila V. Bulik-Sullivan B. Finucane H. K. Walters R. K. Bras J. Duncan L. (2018). Analysis of shared heritability in common disorders of the brain. Science 360:eaa8757. Aryee M. J. Jaffe A. E. Corrada-Bravo H. Ladd-Acosta C. Feinberg A. P. Hansen K. D. (2014). Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays. Bioinformatics 30 13631369. 10.1093/bioinformatics/btu049 24478339 Asiedu M. N. Dussor G. Price T. J. (2016). Targeting AMPK for the alleviation of pathological pain. Exp. Suppl. 107 257285. 10.1007/978-3-319-43589-3_11 Bell C. C. (1994). DSM-IV diagnostic and statistical manual of mental disorders. JAMA 272 828829. Bell J. T. Spector T. D. (2011). A twin approach to unraveling epigenetics. Trends Genet. 27 116125. 10.1016/j.tig.2010.12.005 21257220 Bian Y. Y. Yang L. L. Zhang B. Li W. Li Z. J. Li W. L. (2020). Identification of key genes involved in post-traumatic stress disorder: evidence from bioinformatics analysis. World J. Psychiatry 10 286298. 10.5498/wjp.v10.i12.286 33392005 Blacker C. J. Frye M. A. Morava E. Kozicz T. Veldic M. A. (2019). Review of epigenetics of PTSD in comorbid psychiatric conditions. Genes 10:140. 10.3390/genes10020140 30781888 Bolund A. C. S. Starnawska A. Miller M. R. Schlünssen V. Backer V. Børglum A. D. (2017). Lung function discordance in monozygotic twins and associated differences in blood DNA methylation. Clin. Epigenetics 9:132. Breslau N. Kessler R. C. Chilcoat H. D. Schultz L. R. Davis G. C. Andreski P. (1998). Trauma and posttraumatic stress disorder in the community: the 1996 Detroit Area Survey of Trauma. Arch. Gen. Psychiatry 55 626632. 10.1001/archpsyc.55.7.626 9672053 Chakravarty S. Pathak S. S. Maitra S. Khandelwal N. Chandra K. B. Kumar A. (2014). Epigenetic regulatory mechanisms in stress-induced behavior. Int. Rev. Neurobiol. 115 117154. 10.1016/b978-0-12-801311-3.00004-4 25131544 Chang S. C. Koenen K. C. Galea S. Aiello A. E. Soliven R. Wildman D. E. (2012). Molecular variation at the SLC6A3 locus predicts lifetime risk of PTSD in the Detroit neighborhood health study. PLoS One 7:e39184. 10.1371/journal.pone.0039184 22745713 Chhatwal J. P. Stanek-Rattiner L. Davis M. Ressler K. J. (2006). Amygdala BDNF signaling is required for consolidation but not encoding of extinction. Nat. Neurosci. 9 870872. 10.1038/nn1718 16783370 Cornelis M. C. Nugent N. R. Amstadter A. B. Koenen K. C. (2010). Genetics of post-traumatic stress disorder: review and recommendations for genome-wide association studies. Curr. Psychiatry Rep. 12 313326. 10.1007/s11920-010-0126-6 20549395 Cotto K. C. Wagner A. H. Feng Y. Y. Kiwala S. Coffman A. C. Spies G. (2017). DGIdb 3.0: a redesign and expansion of the drug–gene interaction database. Nucleic Acids Res. 46 D1068D1073. Daskalakis N. P. Yehuda R. (2014). Site-specific methylation changes in the glucocorticoid receptor exon 1F promoter in relation to life adversity: systematic review of contributing factors. Front. Neurosci. 8:369. 10.3389/fnins.2014.00369 25484853 Dresler T. Caratozzolo S. Guldolf K. Huhn J. I. Loiacono C. Niiberg-Pikksööt T. (2019). Understanding the nature of psychiatric comorbidity in migraine: a systematic review focused on interactions and treatment implications. J. Headache Pain 20:51. Eising E. Datson N. A. van den Maagdenberg A. M. Ferrari M. D. (2013). Epigenetic mechanisms in migraine: a promising avenue? BMC Med. 11:26. 10.1186/1741-7015-11-26 23379668 Fischer M. A. Jan A. (2020). Medication-Overuse Headache. Treasure Island, FL: StatPearls [Internet]. Fischer M. Wille G. Klien S. Shanib H. Holle D. Gaul C. (2012). Brain-derived neurotrophic factor in primary headaches. J. Headache Pain 13 469475. 10.1007/s10194-012-0454-5 22584531 Fuchikami M. Morinobu S. Kurata A. Yamamoto S. Yamawaki S. (2009). Single immobilization stress differentially alters the expression profile of transcripts of the brain-derived neurotrophic factor (BDNF) gene and histone acetylation at its promoters in the rat hippocampus. Int. J. Neuropsychopharmacol. 12 7382. 10.1017/s1461145708008997 18544182 Gasperi M. Panizzon M. Goldberg J. Buchwald D. Afari N. (2021). Posttraumatic stress disorder and chronic pain conditions in men: a twin study. Psychosom. Med. 83 109117. 10.1097/psy.0000000000000899 33337593 Gerring Z. F. McRae A. F. Montgomery G. W. Nyholt D. R. (2018). Genome-wide DNA methylation profiling in whole blood reveals epigenetic signatures associated with migraine. BMC Genomics 19:69. 10.1186/s12864-018-4450-2 29357833 Goadsby P. J. Lipton R. B. Ferrari M. D. (2002). Migraine–current understanding and treatment. N. Engl. J. Med. 346 257270. 10.1056/nejmra010917 11807151 Gormley P. Anttila V. Winsvold B. S. Palta P. Esko T. Pers T. H. (2016). Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine. Nat. Genet. 48 856866. Graw S. Henn R. Thompson J. A. Koestler D. C. (2019). pwrEWAS: a user-friendly tool for comprehensive power estimation for epigenome wide association studies (EWAS). BMC Bioinformatics 20:218. 10.1186/s12859-019-2804-7 31035919 Guglielmetti M. Serafini G. Amore M. Martelletti P. (2020). The relation between persistent post-traumatic headache and PTSD: similarities and possible differences. Int. J. Environ. Res. Public Health 17:4024. 10.3390/ijerph17114024 32516965 Guo J. U. Ma D. K. Mo H. Ball M. P. Jang M. H. Bonaguidi M. A. (2011). Neuronal activity modifies the DNA methylation landscape in the adult brain. Nat. Neurosci. 14 13451351. 10.1038/nn.2900 21874013 Guo Y. Rist P. M. Daghlas I. Giulianini F. Kurth T. Chasman D. I. (2020). A genome-wide cross-phenotype meta-analysis of the association of blood pressure with migraine. Nat. Commun. 11:3368. Headache Classification Committee of the International Headache Society (IHS) (2013). The international classification of headache disorders, 3rd edition (beta version). Cephalalgia 33 629808. 10.1177/0333102413485658 23771276 Heldt S. A. Stanek L. Chhatwal J. P. Ressler K. J. (2007). Hippocampus-specific deletion of BDNF in adult mice impairs spatial memory and extinction of aversive memories. Mol. Psychiatry 12 656670. 10.1038/sj.mp.4001957 17264839 Houseman E. A. Accomando W. P. Koestler D. C. Christensen B. C. Marsit C. J. Nelson H. H. (2012). DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 13:86. 10.1186/1471-2105-13-86 22568884 Howie H. Rijal C. M. Ressler K. J. (2019). A review of epigenetic contributions to post-traumatic stress disorder. Dialogues Clin. Neurosci. 21 417428. 10.31887/dcns.2019.21.4/kressler 31949409 Kaut O. Schmitt I. Tost J. Busato F. Liu Y. Hofmann P. (2017). Epigenome-wide DNA methylation analysis in siblings and monozygotic twins discordant for sporadic Parkinson’s disease revealed different epigenetic patterns in peripheral blood mononuclear cells. Neurogenetics 18 722. 10.1007/s10048-016-0497-x 27709425 Kenworthy C. A. Sengupta A. Luz S. M. Ver Hoeve E. S. Meda K. Bhatnagar S. (2014). Social defeat induces changes in histone acetylation and expression of histone modifying enzymes in the ventral hippocampus, prefrontal cortex, and dorsal raphe nucleus. Neuroscience 264 8898. 10.1016/j.neuroscience.2013.01.024 23370319 Kessler R. C. Aguilar-Gaxiola S. Alonso J. Benjet C. Bromet E. J. Cardoso G. (2017). Trauma and PTSD in the WHO world mental health surveys. Eur. J. Psychotraumatol. 8(Suppl. 5):1353383. Kilaru V. Iyer S. V. Almli L. M. Stevens J. S. Lori A. Jovanovic T. (2016). Genome-wide gene-based analysis suggests an association between Neuroligin 1 (NLGN1) and post-traumatic stress disorder. Transl. Psychiatry 6:e820. 10.1038/tp.2016.69 27219346 Koenen K. C. Uddin M. Chang S. C. Aiello A. E. Wildman D. E. Goldmann E. (2011). SLC6A4 methylation modifies the effect of the number of traumatic events on risk for posttraumatic stress disorder. Depress Anxiety 28 639647. 10.1002/da.20825 21608084 Kushnir M. Dresner E. Mandel S. Gozes I. (2008). Silencing of the ADNP-family member, ADNP2, results in changes in cellular viability under oxidative stress. J. Neurochem. 105 537545. 10.1111/j.1471-4159.2007.05173.x 18179478 Launer L. J. Terwindt G. M. Ferrari M. D. (1999). The prevalence and characteristics of migraine in a population-based cohort: the GEM study. Neurology 53 537542. 10.1212/wnl.53.3.537 10449117 Liao A. Wang J. Jaehnig E. J. Shi Z. Zhang B. (2019). WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 47 W199W205. Ligthart S. Marzi C. Aslibekyan S. Mendelson M. M. Conneely K. N. Tanaka T. (2016). DNA methylation signatures of chronic low-grade inflammation are associated with complex diseases. Genome Biol. 17:255. Mason B. N. Russo A. F. (2018). Vascular contributions to migraine: time to revisit? Front. Cell. Neurosci. 12:233. 10.3389/fncel.2018.00233 30127722 Matosin N. Cruceanu C. Binder E. B. (2017). Preclinical and clinical evidence of DNA methylation changes in response to trauma and chronic stress. Chronic Stress 1:2470547017710764 McEwen B. S. (2006). Protective and damaging effects of stress mediators: central role of the brain. Dialogues Clin. Neurosci. 8 367381. 10.31887/dcns.2006.8.4/bmcewen McEwen B. S. (2013). The brain on stress: toward an integrative approach to brain, body, and behavior. Perspect. Psychol. Sci. 8 673675. 10.1177/1745691613506907 25221612 McGowan P. O. (2013). Epigenomic mechanisms of early adversity and HPA dysfunction: considerations for PTSD research. Front. Psychiatry 4:110. 10.3389/fpsyt.2013.00110 24133457 Meaney M. J. Szyf M. (2005). Environmental programming of stress responses through DNA methylation: life at the interface between a dynamic environment and a fixed genome. Dialogues Clin. Neurosci. 7 103123. 10.31887/dcns.2005.7.2/mmeaney Mehta D. Bruenig D. Carrillo-Roa T. Lawford B. Harvey W. Morris C. P. (2017). Genomewide DNA methylation analysis in combat veterans reveals a novel locus for PTSD. Acta Psychiatr. Scand. 136 493505. 10.1111/acps.12778 28795405 Mehta D. Bruenig D. Lawford B. Harvey W. Carrillo-Roa T. Morris C. P. (2018). Accelerated DNA methylation aging and increased resilience in veterans: the biological cost for soldiering on. Neurobiol. Stress 8 112119. 10.1016/j.ynstr.2018.04.001 29888306 Mehta D. Pelzer E. S. Bruenig D. Lawford B. McLeay S. Morris C. P. (2019). DNA methylation from germline cells in veterans with PTSD. J. Psychiatr. Res. 116 4250. 10.1016/j.jpsychires.2019.06.001 31195163 Melroy-Greif W. E. Wilhelmsen K. C. Yehuda R. Ehlers C. L. (2017). Genome-wide association study of post-traumatic stress disorder in two high-risk populations. Twin Res. Hum. Genet. 20 197207. 10.1017/thg.2017.12 28262088 Meyer-Lindenberg A. Weinberger D. R. (2006). Intermediate phenotypes and genetic mechanisms of psychiatric disorders. Nat. Rev. Neurosci. 7 818827. 10.1038/nrn1993 16988657 Minen M. T. Begasse De Dhaem O. Kroon Van Diest A. Powers S. Schwedt T. J. Lipton R. (2016). Migraine and its psychiatric comorbidities. J. Neurol. Neurosurg. Psychiatry 87 741749. Mitchell B. L. Campos A. I. Renteria M. E. Parker R. Sullivan L. McAloney K. (2019). Twenty-five and up (25up) study: a new wave of the Brisbane longitudinal twin study. Twin Res. Hum. Genet. 22 154163. 10.1017/thg.2019.27 31198126 Murai J. Zhang H. Pongor L. Tang S. W. Jo U. Moribe F. (2020). Chromatin remodeling and immediate early gene activation by SLFN11 in response to replication stress. Cell Rep. 30 41374151.e6. Nievergelt C. M. Ashley-Koch A. E. Dalvie S. Hauser M. A. Morey R. A. Smith A. K. (2018). Genomic approaches to posttraumatic stress disorder: the psychiatric genomic consortium initiative. Biol. Psychiatry 83 831839. 10.1016/j.biopsych.2018.01.020 29555185 Numata S. Ye T. Hyde T. M. Guitart-Navarro X. Tao R. Wininger M. (2012). DNA methylation signatures in development and aging of the human prefrontal cortex. Am. J. Hum. Genet. 90 260272. 10.1016/j.ajhg.2011.12.020 22305529 Ouellet-Morin I. Wong C. C. Danese A. Pariante C. M. Papadopoulos A. S. Mill J. (2013). Increased serotonin transporter gene (SERT) DNA methylation is associated with bullying victimization and blunted cortisol response to stress in childhood: a longitudinal study of discordant monozygotic twins. Psychol. Med. 43 18131823. 10.1017/s0033291712002784 23217646 Peedicayil J. Kumar A. (2012). Time for clinical trials of epigenetic drugs in psychiatric disorders? Br. J. Clin. Pharmacol. 73 309310. 10.1111/j.1365-2125.2011.04100.x 22023150 Peng H. Zhu Y. Strachan E. Fowler E. Bacus T. Roy-Byrne P. (2018). Childhood trauma, DNA methylation of stress-related genes, and depression: findings from two monozygotic twin studies. Psychosom. Med. 80 599608. 10.1097/psy.0000000000000604 29781947 Peterlin B. L. Rosso A. L. Sheftell F. D. Libon D. J. Mossey J. M. Merikangas K. R. (2011). Post-traumatic stress disorder, drug abuse and migraine: new findings from the National Comorbidity Survey Replication (NCS-R). Cephalalgia 31 235244. 10.1177/0333102410378051 20813779 Peterlin B. L. Tietjen G. Meng S. Lidicker J. Bigal M. (2008). Post-traumatic stress disorder in episodic and chronic migraine. Headache 48 517522. 10.1111/j.1526-4610.2008.00917.x 18377377 Probst A. V. Dunleavy E. Almouzni G. (2009). Epigenetic inheritance during the cell cycle. Nat. Rev. Mol. Cell Biol. 10 192206. 10.1038/nrm2640 19234478 Ressler K. J. Mercer K. B. Bradley B. Jovanovic T. Mahan A. Kerley K. (2011). Post-traumatic stress disorder is associated with PACAP and the PAC1 receptor. Nature 470 492497. Rusiecki J. A. Chen L. Srikantan V. Zhang L. Yan L. Polin M. L. (2012). DNA methylation in repetitive elements and post-traumatic stress disorder: a case-control study of US military service members. Epigenomics 4 2940. 10.2217/epi.11.116 22332656 Rustichelli C. Lo Castro F. Baraldi C. Ferrari A. (2020). Targeting pituitary adenylate cyclase-activating polypeptide (PACAP) with monoclonal antibodies in migraine prevention: a brief review. Expert Opin. Investig. Drugs 29 12691275. 10.1080/13543784.2020.1811966 32877252 Sareen J. (2014). Posttraumatic stress disorder in adults: impact, comorbidity, risk factors, and treatment. Can. J. Psychiatry 59 460467. 10.1177/070674371405900902 25565692 Schouten M. Aschrafi A. Bielefeld P. Doxakis E. Fitzsimons C. P. (2013). microRNAs and the regulation of neuronal plasticity under stress conditions. Neuroscience 241 188205. 10.1016/j.neuroscience.2013.02.065 23500144 Seng E. K. Seng C. D. (2016). Understanding migraine and psychiatric comorbidity. Curr. Opin. Neurol. 29 309313. 10.1097/wco.0000000000000309 26886355 Shiloh R. Gilad Y. Ber Y. Eisenstein M. Aweida D. Bialik S. (2018). Non-canonical activation of DAPK2 by AMPK constitutes a new pathway linking metabolic stress to autophagy. Nat. Commun. 9:1759. Smith A. K. Conneely K. N. Kilaru V. Mercer K. B. Weiss T. E. Bradley B. (2011). Differential immune system DNA methylation and cytokine regulation in post-traumatic stress disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 156b 700708. 10.1002/ajmg.b.31212 21714072 Smitherman T. A. Rains J. C. Penzien D. B. (2009). Psychiatric comorbidities and migraine chronification. Curr. Pain Headache Rep. 13 326331. 10.1007/s11916-009-0052-7 19586598 Soliman F. Glatt C. E. Bath K. G. Levita L. Jones R. M. Pattwell S. S. (2010). A genetic variant BDNF polymorphism alters extinction learning in both mouse and human. Science 327 863866. 10.1126/science.1181886 20075215 Stein M. B. Levey D. F. Cheng Z. Wendt F. R. Harrington K. Pathak G. A. (2021). Genome-wide association analyses of post-traumatic stress disorder and its symptom subdomains in the Million Veteran Program. Nat. Genet. 53 174184. 10.1038/s41588-020-00767-x 33510476 Sundrum T. Walker C. S. (2018). Pituitary adenylate cyclase-activating polypeptide receptors in the trigeminovascular system: implications for migraine. Br. J. Pharmacol. 175 41094120. 10.1111/bph.14053 28977676 Sutherland H. Maher B. Rodriguez Acevedo A. J. Haupt L. Griffiths L. (2014). Investigation of brain-derived neurotrophic factor (BDNF) gene variants in migraine. Headache 54 11841193. 10.1111/head.12351 24708359 Taliento A. E. Dallio M. Federico A. Prati D. Valenti L. (2019). Novel insights into the genetic landscape of nonalcoholic fatty liver disease. Int. J. Environ. Res. Public Health 16:2755. 10.3390/ijerph16152755 31375010 Tan Q. Christiansen L. von Bornemann Hjelmborg J. Christensen K. (2015). Twin methodology in epigenetic studies. J. Exp. Biol. 218:134. 10.1242/jeb.107151 25568460 Theeler B. J. Flynn F. G. Erickson J. C. (2012). Chronic daily headache in U.S. soldiers after concussion. Headache 52 732738. 10.1111/j.1526-4610.2012.02112.x 22404747 Tsilioni I. Patel A. B. Pantazopoulos H. Berretta S. Conti P. Leeman S. E. (2019). IL-37 is increased in brains of children with autism spectrum disorder and inhibits human microglia stimulated by neurotensin. Proc. Natl. Acad. Sci. U.S.A. 116 2165921665. 10.1073/pnas.1906817116 31591201 Uddin M. Aiello A. E. Wildman D. E. Koenen K. C. Pawelec G. de Los Santos R. (2010). Epigenetic and immune function profiles associated with posttraumatic stress disorder. Proc. Natl. Acad. Sci. U.S.A. 107 94709475. 10.1073/pnas.0910794107 20439746 van den Maagdenberg A. Nyholt D. R. Anttila V. (2019). Novel hypotheses emerging from GWAS in migraine? J. Headache Pain 20:5. Wang J. Duncan D. Shi Z. Zhang B. (2013). WEB-based GEne SeT analysis toolkit (WebGestalt): update 2013. Nucleic Acids Res. 41 W77W83. Wang J. Xiao B. Han F. Shi Y. Ibata K. (2017). Metformin alleviated the neuronal oxidative stress in hippocampus of rats under single prolonged stress. J. Mol. Neurosci. 63 2835. 10.1007/s12031-017-0953-6 28752417 Wong C. M. (2002). Post-traumatic stress disorder: advances in psychoneuroimmunology. Psychiatr. Clin. North Am. 25 369383. Wright M. J. Martin N. G. (2004). Brisbane adolescent twin study: outline of study methods and research projects. Aust. J. Psychol. 56 6578. 10.1080/00049530410001734865 Yehuda R. LeDoux J. (2007). Response variation following trauma: a translational neuroscience approach to understanding PTSD. Neuron 56 1932. 10.1016/j.neuron.2007.09.006 17920012 Yehuda R. Hoge C. W. McFarlane A. C. Vermetten E. Lanius R. A. Nievergelt C. M. (2015). Post-traumatic stress disorder. Nat. Rev. Dis. Primers 1:15057. Zarei M. R. Shabani M. Chamani G. Abareghi F. Razavinasab M. Nazeri M. (2016). Migraine patients have a higher prevalence of PTSD symptoms in comparison to chronic tension-type headache and healthy subjects: a case-control study. Acta Odontol. Scand. 74 633635. 10.1080/00016357.2016.1232435 27644346 Zhou W. J. Xu N. Kong L. Sun S. C. Xu X. F. Jia M. Z. (2016). The antidepressant roles of Wnt2 and Wnt3 in stress-induced depression-like behaviors. Transl. Psychiatry 6:e892. 10.1038/tp.2016.122 27622936
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016www.imersia.org.cn
      www.fdyxsd.com.cn
      i-33.com.cn
      eydfpx.com.cn
      www.lbirti.com.cn
      lykxgm.org.cn
      keaibo.com.cn
      www.mmcpew.com.cn
      www.qxkpoo.com.cn
      pocketgf.org.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p