Critical evaluation of copy number variant calling methods using DNA methylation.

Recent technological and methodological developments have enabled the use of array-based DNA methylation data to call copy number variants (CNVs). ChAMP, Conumee, and cnAnalysis450k are popular methods currently used to call CNVs using methylation data. However, so far, no studies have analyzed the reliability of these methods using real samples. Data from a cohort of individuals with genotype and DNA methylation data generated using the HumanMethylation450 and MethylationEPIC BeadChips were used to assess the consistency between the CNV calls generated by methylation and genotype data. We also took advantage of repeated measures of methylation data collected from the same individuals to compare the reliability of CNVs called by ChAMP, Conumee, and cnAnalysis450k for both the methylation arrays. ChAMP identified more CNVs than Conumee and cnAnalysis450k for both the arrays and, as a consequence, had a higher overlap (~62%) with the calls from the genotype data. However, all methods had relatively low reliability. For the MethylationEPIC array, Conumee had the highest reliability (57.6%), whereas for the HumanMethylation450 array, cnAnalysis450k had the highest reliability (43.0%). Overall, the MethylationEPIC array provided significant gains in reliability for CNV calling over the HumanMethylation450 array but not for overlap with CNVs called using genotype data.

PTSD is associated with increased DNA methylation across regions of HLA-DPB1 and SPATC1L.

Posttraumatic stress disorder (PTSD) is characterized by intrusive thoughts, avoidance, negative alterations in cognitions and mood, and arousal symptoms that adversely affect mental and physical health. Recent evidence links changes in DNA methylation of CpG cites to PTSD. Since clusters of proximal CpGs share similar methylation signatures, identification of PTSD-associated differentially methylated regions (DMRs) may elucidate the pathways defining differential risk and resilience of PTSD. Here we aimed to identify epigenetic differences associated with PTSD. DNA methylation data profiled from blood samples using the MethylationEPIC BeadChip were used to perform a DMR analysis in 187 PTSD cases and 367 trauma-exposed controls from the Grady Trauma Project (GTP). DMRs were assessed with R package bumphunter. We identified two regions that associate with PTSD after multiple test correction. These regions were in the gene body of HLA-DPB1 and in the promoter of SPATC1L. The DMR in HLA-DPB1 was associated with PTSD in an independent cohort. Both DMRs included CpGs whose methylation associated with nearby sequence variation (meQTL) and that associated with expression of their respective genes (eQTM). This study supports an emerging literature linking PTSD risk to genetic and epigenetic variation in the HLA region.

Blood-Derived DNA Methylation Signatures of Crohn's Disease and Severity of Intestinal Inflammation.

BACKGROUND & AIMS: Crohn's disease is a relapsing and remitting inflammatory disorder with a variable clinical course. Although most patients present with an inflammatory phenotype (B1), approximately 20% of patients rapidly progress to complicated disease, which includes stricturing (B2), within 5 years. We analyzed DNA methylation patterns in blood samples of pediatric patients with Crohn's disease at diagnosis and later time points to identify changes that associate with and might contribute to disease development and progression. METHODS: We obtained blood samples from 164 pediatric patients (1-17 years old) with Crohn's disease (B1 or B2) who participated in a North American study and were followed for 5 years. Participants without intestinal inflammation or symptoms served as controls (n = 74). DNA methylation patterns were analyzed in samples collected at time of diagnosis and 1-3 years later at approximately 850,000 sites. We used genetic association and the concept of Mendelian randomization to identify changes in DNA methylation patterns that might contribute to the development of or result from Crohn's disease. RESULTS: We identified 1189 5'-cytosine-phosphate-guanosine-3' (CpG) sites that were differentially methylated between patients with Crohn's disease (at diagnosis) and controls. Methylation changes at these sites correlated with plasma levels of C-reactive protein. A comparison of methylation profiles of DNA collected at diagnosis of Crohn's disease vs during the follow-up period showed that, during treatment, alterations identified in methylation profiles at the time of diagnosis of Crohn's disease more closely resembled patterns observed in controls, irrespective of disease progression to B2. We identified methylation changes at 3 CpG sites that might contribute to the development of Crohn's disease. Most CpG methylation changes associated with Crohn's disease disappeared with treatment of inflammation and might be a result of Crohn's disease. CONCLUSIONS: Methylation patterns observed in blood samples from patients with Crohn's disease accompany acute inflammation; with treatment, these change to resemble methylation patterns observed in patients without intestinal inflammation. These findings indicate that Crohn's disease-associated patterns of DNA methylation observed in blood samples are a result of the inflammatory features of the disease and are less likely to contribute to disease development or progression.

Association of HLA locus alleles with posttraumatic stress disorder.

BACKGROUND: Immune dysregulation has been widely observed in those with posttraumatic stress disorder (PTSD). An individual's immune response is shaped, in part, by the highly polymorphic Human Leukocyte Antigen (HLA) locus that is associated with major psychiatric disorders such as schizophrenia, major depression and bipolar disorder. The aim of the current study was to investigate the association between common HLA alleles and PTSD. METHODS: Genome-wide association data was used to predict alleles of 7 classical polymorphic HLA genes (A, B, C, DRB1, DQA1, DQB1, DPB1) in 403 lifetime PTSD cases and 369 trauma exposed controls of African ancestry. Association of HLA allelic variations with lifetime PTSD was analyzed using logistic regression, controlling for ancestry, sex and multiple comparisons. The effect of HLA alleles on gene expression was assessed by weighted correlation network analysis (WGCNA), using 353 subjects with available expression data. Enrichment analysis was performed using anRichment to identify associated pathways of each module. RESULTS: HLA-B*58:01 (p = 0.035), HLA-C*07:01 (p = 0.035), HLA-DQA1*01:01 (p = 0.003), HLA-DQB1*05:01 (p = 0.009) and HLA-DPB1*17:01 (p = 0.017) were more common in PTSD cases, while HLA-A*02:01 (p = 0.026), HLA-DQA1*05:05 (p = 0.011) and HLA-DRB1*11:01 (p < 0.001) were more frequent in controls. WGCNA was used to explore expression patterns of the PTSD related alleles. Gene expression modules of PTSD-related HLA alleles were enriched in various pathways, including pathways related to immune and neural activity. CONCLUSIONS: To the best of our knowledge, this is the first study to report an association of HLA alleles with PTSD. Altogether, our results support the link between the immune system, brain and PTSD.

Methylation differences reveal heterogeneity in preterm pathophysiology: results from bipartite network analyses.

BACKGROUND: Recent studies have shown that epigenetic differences can increase the risk of spontaneous preterm birth (PTB). However, little is known about heterogeneity underlying such epigenetic differences, which could lead to hypotheses for biological pathways in specific patient subgroups, and corresponding targeted interventions critical for precision medicine. Using bipartite network analysis of fetal DNA methylation data we demonstrate a novel method for classification of PTB. METHODS: The data consisted of DNA methylation across the genome (HumanMethylation450 BeadChip) in cord blood from 50 African-American subjects consisting of 22 cases of early spontaneous PTB (24-34 weeks of gestation) and 28 controls (>39 weeks of gestation). These data were analyzed using a combination of (1) a supervised method to select the top 10 significant methylation sites, (2) unsupervised "subject-variable" bipartite networks to visualize and quantitatively analyze how those 10 methylation sites co-occurred across all the subjects, and across only the cases with the goal of analyzing subgroups and their underlying pathways, and (3) a simple linear regression to test whether there was an association between the total methylation in the cases, and gestational age. RESULTS: The bipartite network analysis of all subjects and significant methylation sites revealed statistically significant clustering consisting of an inverse symmetrical relationship in the methylation profiles between a case-enriched subgroup and a control-enriched subgroup: the former was predominantly hypermethylated across seven methylation sites, and hypomethylated across three methylation sites, whereas the latter was predominantly hypomethylated across the above seven methylation sites and hypermethylated across the three methylation sites. Furthermore, the analysis of only cases revealed one subgroup that was predominantly hypomethylated across seven methylation sites, and another subgroup that was hypomethylated across all methylation sites suggesting the presence of heterogeneity in PTB pathophysiology. Finally, the analysis found a strong inverse linear relationship between total methylation and gestational age suggesting that methylation differences could be used as predictive markers for gestational length. CONCLUSIONS: The results demonstrate that unsupervised bipartite networks helped to identify a complex but comprehensible data-driven hypotheses related to patient subgroups and inferences about their underlying pathways, and therefore were an effective complement to supervised approaches currently used.

Post-traumatic stress disorder is associated with PACAP and the PAC1 receptor.

Pituitary adenylate cyclase-activating polypeptide (PACAP) is known to broadly regulate the cellular stress response. In contrast, it is unclear if the PACAP-PAC1 receptor pathway has a role in human psychological stress responses, such as post-traumatic stress disorder (PTSD). Here we find, in heavily traumatized subjects, a sex-specific association of PACAP blood levels with fear physiology, PTSD diagnosis and symptoms in females. We examined 44 single nucleotide polymorphisms (SNPs) spanning the PACAP (encoded by ADCYAP1) and PAC1 (encoded by ADCYAP1R1) genes, demonstrating a sex-specific association with PTSD. A single SNP in a putative oestrogen response element within ADCYAP1R1, rs2267735, predicts PTSD diagnosis and symptoms in females only. This SNP also associates with fear discrimination and with ADCYAP1R1 messenger RNA expression in human brain. Methylation of ADCYAP1R1 in peripheral blood is also associated with PTSD. Complementing these human data, ADCYAP1R1 mRNA is induced with fear conditioning or oestrogen replacement in rodent models. These data suggest that perturbations in the PACAP-PAC1 pathway are involved in abnormal stress responses underlying PTSD. These sex-specific effects may occur via oestrogen regulation of ADCYAP1R1. PACAP levels and ADCYAP1R1 SNPs may serve as useful biomarkers to further our mechanistic understanding of PTSD.

DNA Methylation Changes in Whole Blood and CD16+ Neutrophils in Response to Chronic Folic Acid Supplementation in Women of Childbearing Age.

Folate, a water-soluble vitamin, is a key source of one-carbon groups for DNA methylation, but studies of the DNA methylation response to supplemental folic acid yield inconsistent results. These studies are commonly conducted using whole blood, which contains a mixed population of white blood cells that have been shown to confound results. The objective of this study was to determine if CD16+ neutrophils may provide more specific data than whole blood for identifying DNA methylation response to chronic folic acid supplementation. The study was performed in normal weight (BMI 18.5 - 24.9 kg/m2) women (18 - 35 y; n = 12), with blood samples taken before and after 8 weeks of folic acid supplementation at 800 µg/day. DNA methylation patterns from whole blood and isolated CD16+ neutrophils were measured across >485,000 CpG sites throughout the genome using the Infinium HumanMethylation450 BeadChip. Over the course of the 8-week supplementation, 6746 and 7513 CpG sites changed (p < 0.05) in whole blood and CD16+ neutrophils, respectively. DNA methylation decreased in 68.4% (whole blood) and 71.8% (CD16+ neutrophils) of these sites. There were only 182 CpG sites that changed in both the whole blood and CD16+ neutrophils, 139 of which changed in the same direction. These results suggest that the genome-wide DNA methylation response to chronic folic acid supplementation is different between whole blood and CD16+ neutrophils and that a single white blood cell type may function as a more specific epigenetic reporter of folate status than whole blood.

Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline.

Compelling evidence suggests that epigenetic mechanisms such as DNA methylation play a role in stress regulation and in the etiologic basis of stress related disorders such as Post traumatic Stress Disorder (PTSD). Here we describe the purpose and methods of an international consortium that was developed to study the role of epigenetics in PTSD. Inspired by the approach used in the Psychiatric Genomics Consortium, we brought together investigators representing seven cohorts with a collective sample size of N = 1147 that included detailed information on trauma exposure, PTSD symptoms, and genome-wide DNA methylation data. The objective of this consortium is to increase the analytical sample size by pooling data and combining expertise so that DNA methylation patterns associated with PTSD can be identified. Several quality control and analytical pipelines were evaluated for their control of genomic inflation and technical artifacts with a joint analysis procedure established to derive comparable data over the cohorts for meta-analysis. We propose methods to deal with ancestry population stratification and type I error inflation and discuss the advantages and disadvantages of applying robust error estimates. To evaluate our pipeline, we report results from an epigenome-wide association study (EWAS) of age, which is a well-characterized phenotype with known epigenetic associations. Overall, while EWAS are highly complex and subject to similar challenges as genome-wide association studies (GWAS), we demonstrate that an epigenetic meta-analysis with a relatively modest sample size can be well-powered to identify epigenetic associations. Our pipeline can be used as a framework for consortium efforts for EWAS.

Accounting for population stratification in DNA methylation studies.

DNA methylation is an important epigenetic mechanism that has been linked to complex diseases and is of great interest to researchers as a potential link between genome, environment, and disease. As the scale of DNA methylation association studies approaches that of genome-wide association studies, issues such as population stratification will need to be addressed. It is well-documented that failure to adjust for population stratification can lead to false positives in genetic association studies, but population stratification is often unaccounted for in DNA methylation studies. Here, we propose several approaches to correct for population stratification using principal components (PCs) from different subsets of genome-wide methylation data. We first illustrate the potential for confounding due to population stratification by demonstrating widespread associations between DNA methylation and race in 388 individuals (365 African American and 23 Caucasian). We subsequently evaluate the performance of our PC-based approaches and other methods in adjusting for confounding due to population stratification. Our simulations show that (1) all of the methods considered are effective at removing inflation due to population stratification, and (2) maximum power can be obtained with single-nucleotide polymorphism (SNP)-based PCs, followed by methylation-based PCs, which outperform both surrogate variable analysis and genomic control. Among our different approaches to computing methylation-based PCs, we find that PCs based on CpG sites chosen for their potential to proxy nearby SNPs can provide a powerful and computationally efficient approach to adjust for population stratification in DNA methylation studies when genome-wide SNP data are unavailable.

DNA extracted from saliva for methylation studies of psychiatric traits: evidence tissue specificity and relatedness to brain.

DNA methylation has become increasingly recognized in the etiology of psychiatric disorders. Because brain tissue is not accessible in living humans, epigenetic studies are most often conducted in blood. Saliva is often collected for genotyping studies but is rarely used to examine DNA methylation because the proportion of epithelial cells and leukocytes varies extensively between individuals. The goal of this study was to evaluate whether saliva DNA is informative for studies of psychiatric disorders. DNA methylation (HumanMethylation450 BeadChip) was assessed in saliva and blood samples from 64 adult African Americans. Analyses were conducted using linear regression adjusted for appropriate covariates, including estimated cellular proportions. DNA methylation from brain tissues (cerebellum, frontal cortex, entorhinal cortex, and superior temporal gyrus) was obtained from a publically available dataset. Saliva and blood methylation was clearly distinguishable though there was positive correlation overall. There was little correlation in CpG sites within relevant candidate genes. Correlated CpG sites were more likely to occur in areas of low CpG density (i.e., CpG shores and open seas). There was more variability in CpG sites from saliva than blood, which may reflect its heterogeneity. Finally, DNA methylation in saliva appeared more similar to patterns from each of the brain regions examined overall than methylation in blood. Thus, this study provides a framework for using DNA methylation from saliva and suggests that DNA methylation of saliva may offer distinct opportunities for epidemiological and longitudinal studies of psychiatric traits.

A genome-wide identified risk variant for PTSD is a methylation quantitative trait locus and confers decreased cortical activation to fearful faces.

Genetic factors appear to be highly relevant to predicting differential risk for the development of post-traumatic stress disorder (PTSD). In a discovery sample, we conducted a genome-wide association study (GWAS) for PTSD using a small military cohort (Systems Biology PTSD Biomarkers Consortium; SBPBC, N = 147) that was designed as a case-controlled sample of highly exposed, recently returning veterans with and without combat-related PTSD. A genome-wide significant single nucleotide polymorphism (SNP), rs717947, at chromosome 4p15 (N = 147, ß = 31.34, P = 1.28 × 10(-8) ) was found to associate with the gold-standard diagnostic measure for PTSD (the Clinician Administered PTSD Scale). We conducted replication and follow-up studies in an external sample, a larger urban community cohort (Grady Trauma Project, GTP, N = 2006), to determine the robustness and putative functionality of this risk variant. In the GTP replication sample, SNP rs717947 associated with PTSD diagnosis in females (N = 2006, P = 0.005), but not males. SNP rs717947 was also found to be a methylation quantitative trait locus (meQTL) in the GTP replication sample (N = 157, P = 0.002). Further, the risk allele of rs717947 was associated with decreased medial and dorsolateral cortical activation to fearful faces (N = 53, P < 0.05) in the GTP replication sample. These data identify a genome-wide significant polymorphism conferring risk for PTSD, which was associated with differential epigenetic regulation and with differential cortical responses to fear in a replication sample. These results may provide new insight into understanding genetic and epigenetic regulation of PTSD and intermediate phenotypes that contribute to this disorder.

Methylation quantitative trait loci (meQTLs) are consistently detected across ancestry, developmental stage, and tissue type.

BACKGROUND: Individual genotypes at specific loci can result in different patterns of DNA methylation. These methylation quantitative trait loci (meQTLs) influence methylation across extended genomic regions and may underlie direct SNP associations or gene-environment interactions. We hypothesized that the detection of meQTLs varies with ancestral population, developmental stage, and tissue type. We explored this by analyzing seven datasets that varied by ancestry (African American vs. Caucasian), developmental stage (neonate vs. adult), and tissue type (blood vs. four regions of postmortem brain) with genome-wide DNA methylation and SNP data. We tested for meQTLs by constructing linear regression models of methylation levels at each CpG site on SNP genotypes within 50 kb under an additive model controlling for multiple tests. RESULTS: Most meQTLs mapped to intronic regions, although a limited number appeared to occur in synonymous or nonsynonymous coding SNPs. We saw significant overlap of meQTLs between ancestral groups, developmental stages, and tissue types, with the highest rates of overlap within the four brain regions. Compared with a random group of SNPs with comparable frequencies, meQTLs were more likely to be 1) represented among the most associated SNPs in the WTCCC bipolar disorder results and 2) located in microRNA binding sites. CONCLUSIONS: These data give us insight into how SNPs impact gene regulation and support the notion that peripheral blood may be a reliable correlate of physiological processes in other tissues.

Epigenetic changes associated with inflammation in breast cancer patients treated with chemotherapy.

Inflammation has been associated with fatigue during and after various types of breast cancer treatments. We examined whether prior chemotherapy was associated with DNA methylation patterns that could explain persisting inflammation and/or fatigue in women treated for breast cancer. Prior to breast radiation therapy, DNA was extracted from peripheral blood mononuclear cells (PBMCs) of 61 Stage 0-IIIA breast cancer patients who had received partial mastectomy with or without chemotherapy. DNA methylation was assessed at >485,000 CpG sites across the genome along with fatigue and plasma inflammatory markers previously associated with fatigue. Compared to non-chemotherapy-treated, women who had received chemotherapy exhibited significantly decreased methylation at eight CpG sites (p<1.03×10(-7)) including four in exon 11 of transmembrane protein 49 (TMEM49), which demonstrated the largest decreases in methylation. Lower methylation at each identified CpG site was associated with increased plasma soluble tumor necrosis factor receptor 2 (sTNFR2) and interleukin (IL)-6 and mediated the relationship between chemotherapy and increases in these inflammatory biomarkers adjusting for multiple clinical and treatment characteristics. sTNFR2, but not CpG methylation status, was correlated with fatigue. Six months after breast radiation therapy, DNA methylation, inflammatory biomarkers and fatigue assessments were repeated in a subset of subjects (N=39). Reduced methylation in 4 of the 8 identified CpG sites was still observed in chemotherapy versus non-chemotherapy-treated patients, albeit with some decay indicating the dynamic and potentially reversible nature of the changes. Reduced methylation in these 4 CpG sites also continued to correlate with either increased sTNFR2 or IL-6, but not fatigue. In conclusion, prior chemotherapy treatment was associated with decreased methylation of CpG sites in DNA from PBMCs of breast cancer patients, which correlated with increased inflammatory markers prior to and 6months after radiation therapy. Persisting epigenetic changes secondary to chemotherapy may be one factor that contributes to inflammation and its consequences including cancer-related fatigue in vulnerable breast cancer patients.

CpGassoc: an R function for analysis of DNA methylation microarray data.

SUMMARY: With the increasing availability of high-density methylation microarrays, there has been growing interest in analysis of DNA methylation data. We have developed CpGassoc, an R package that can efficiently perform the statistical analysis needed for increasingly large methylation datasets. CpGassoc is a modular, expandable package with functions to perform rapid analyses of DNA methylation data via fixed or mixed effects models, to perform basic quality control, to carry out permutation tests, and to display results via an array of publication-quality plots. AVAILABILITY AND IMPLEMENTATION: CpGassoc is implemented in R and is freely available at http://genetics.emory.edu/conneely; we are in the process of submitting it to CRAN.

Longitudinal DNA methylation profiling of the rectal mucosa identifies cell-specific signatures of disease status, severity and clinical outcomes in ulcerative colitis cell-specific DNA methylation signatures of UC.

BACKGROUND: In peripheral blood, DNA methylation (DNAm) patterns in inflammatory bowel disease patients reflect inflammatory status rather than disease status. Here, we examined DNAm in diseased rectal mucosa from ulcerative colitis (UC) patients, focusing on constituent cell types with the goal of identifying therapeutic targets for UC other than the immune system. We profiled DNAm of rectal mucosal biopsies of pediatric UC at diagnosis (n = 211) and non-IBD control (n = 85) patients and performed epigenome-wide association studies (EWAS) of specific cell types to understand DNAm changes in epithelial, immune and fibroblast cells across disease states, course, and clinical outcomes. We also examined longitudinal analysis on follow-up samples (n = 73), and comparisons were made among patients with clinical outcomes including those undergoing colectomy versus those who did not. Additionally, we included RNA-seq from the same subjects to assess the impact of CpG sites on the transcription of nearby genes during the disease course. RESULTS: At diagnosis, UC rectal mucosa exhibited a lower proportion of epithelial cells and fibroblasts, and higher proportion of immune cells, in conjunction with variation in the DNAm pattern. While treatment had significant effects on the methylation signature of immune cells, its effects on fibroblasts and epithelial cells were attenuated. Individuals who required colectomy exhibited cell composition and DNAm patterns at follow-up more similar to disease onset than patients who did not require colectomy. Combining these results with gene expression profiles, we identify CpG sites whose methylation patterns are most consistent with a contribution to poor disease outcomes and could thus be potential therapeutic targets. CONCLUSIONS: Cell-specific epigenetic changes in the rectal mucosa in UC are associated with disease severity and outcome. Current therapeutics may more effectively target the immune than the epithelial and fibroblast compartments.

Methylation quantitative trait loci are largely consistent across disease states in Crohn's disease.

Recently, we identified 1,189 CpG sites whose DNA methylation level in blood associated with Crohn's disease. Here, we examined associations between DNA methylation and genetic variants to identify methylation quantitative trait loci across disease states in (1) 402 blood samples from 164 newly diagnosed pediatric Crohn's disease patients taken at 2 time points (diagnosis and follow-up), and 74 non-inflammatory bowel disease controls, (2) 780 blood samples from a non-Crohn's disease adult population, and (3) 40 ileal biopsies (17 Crohn's disease cases and 23 non-inflammatory bowel disease controls) from group (1). Genome-wide DNAm profiling and genotyping were performed using the Illumina MethylationEPIC and Illumina Multi-Ethnic arrays. SNP-CpG associations were identified via linear models adjusted for age, sex, disease status, disease subtype, estimated cell proportions, and genotype-based principal components. In total, we observed 535,448 SNP-CpG associations between 287,881 SNPs and 12,843 CpG sites (P < 8.21 × 10-14). Associations were highly consistent across different ages, races, disease states, and tissue types, suggesting that the majority of these methylation quantitative trait loci participate in common gene regulation. However, genes near CpGs associated with inflammatory bowel disease SNPs were enriched for 18 KEGG pathways relevant to inflammatory bowel disease-linked immune function and inflammatory responses. We observed suggestive evidence for a small number of tissue-specific associations and disease-specific associations in ileum, though larger studies will be needed to confirm these results. Our study concludes that the vast majority of blood-derived methylation quantitative trait loci are common across individuals, though a subset may be involved in processes related to Crohn's disease. Independent cohort studies will be required to validate these findings.

Differential immune system DNA methylation and cytokine regulation in post-traumatic stress disorder.

DNA methylation may mediate persistent changes in gene function following chronic stress. To examine this hypothesis, we evaluated African American subjects matched by age and sex, and stratified into four groups by post-traumatic stress disorder (PTSD) diagnosis and history of child abuse. Total Life Stress (TLS) was also assessed in all subjects. We evaluated DNA extracted from peripheral blood using the HumanMethylation27 BeadChip and analyzed both global and site-specific methylation. Methylation levels were examined for association with PTSD, child abuse history, and TLS using a linear mixed model adjusted for age, sex, and chip effects. Global methylation was increased in subjects with PTSD. CpG sites in five genes (TPR, CLEC9A, APC5, ANXA2, and TLR8) were differentially methylated in subjects with PTSD. Additionally, a CpG site in NPFFR2 was associated with TLS after adjustment for multiple testing. Notably, many of these genes have been previously associated with inflammation. Given these results and reports of immune dysregulation associated with trauma history, we compared plasma cytokine levels in these subjects and found IL4, IL2, and TNFα levels associated with PTSD, child abuse, and TLS. Together, these results suggest that psychosocial stress may alter global and gene-specific DNA methylation patterns potentially associated with peripheral immune dysregulation. Our results suggest the need for further research on the role of DNA methylation in stress-related illnesses.

Epigenetic prediction of 17ß-estradiol and relationship to trauma-related outcomes in women.

17ß-estradiol (E2) levels in women correlate with multiple neuropsychiatric symptoms, including those that are stress-related. Furthermore, prior work from our group has demonstrated that E2 status influences DNA methylation (DNAm) across the genome. We developed and validated a DNAm-based predictor of E2 (one of four naturally occurring estrogens) using a training set of 183 females and a test set of 79 females from the same traumatized cohort. We showed that predicted E2 levels were highly correlated with measured E2 concentrations in our testing set (r â€‹= â€‹0.75, p â€‹= â€‹1.8e-15). We further demonstrated that predicted E2 concentrations, in combination with measured values, negatively correlated with current post-traumatic stress disorder (PTSD) (ߠ​= â€‹-0.38, p â€‹= â€‹0.01) and major depressive disorder (MDD) diagnoses (ߠ​= â€‹-0.45, p â€‹= â€‹0.02), as well as a continuous measure of PTSD symptom severity (ߠ​= â€‹-2.3, p â€‹= â€‹0.007) in females. Finally, we tested our predictor in an independent data set (n â€‹= â€‹85) also comprised of recently traumatized female subjects to determine if the predictor would generalize to a different population than the one on which it was developed. We found that the correlation between predicted and actual E2 concentrations in the external validation data set was also high (r â€‹= â€‹0.48, p â€‹= â€‹3.0e-6). While further validation is warranted, a DNAm predictor of E2 concentrations will advance our understanding of hormone-epigenetic interactions. Furthermore, such a DNAm predictor may serve as an epigenetic proxy for E2 concentrations and thus provide an important biomarker to better evaluate the contribution of E2 to current and potentially future psychiatric symptoms in samples for which E2 is not measured.

Genome-wide DNA methylation differences and polychlorinated biphenyl (PCB) exposure in a US population.

Exposure to polychlorinated biphenyls (PCBs), an endocrine-disrupting compound, is ubiquitous despite decades-old bans on the manufacture and use of PCBs. Increased exposure to PCBs is associated with adverse health consequences throughout life, including type 2 diabetes and cancer. PCB exposure is also associated with alterations in epigenetic marks and gene transcription, which could lead to adverse health outcomes, but many of these are population-specific. To further investigate the association between PCB and epigenetic marks, DNA methylation was measured at 787,684 CpG sites in 641 peripheral blood samples from the Michigan Polybrominated Biphenyl (PBB) Registry. 1345 CpGs were associated with increased total PCB level after controlling for age, sex, and 24 surrogate variables (FDR < 0.05). These CpGs were enriched in active promoter and transcription associated regions (p < 0.05), and in regions around the binding sites for transcription factors involved in xenobiotic metabolism and immune function (FDR < 0.05). PCB exposure also associated with proportions of CD4T, NK, and granulocyte cell types, and with the neutrophil to lymphocyte ratio (NLR) (p < 0.05), and the estimated effect sizes of PCB on the epigenome were correlated with the effect sizes previously reported in an epigenome-wide study of C-reactive protein (r = 0.29; p = 2.22e-5), supporting previous studies on the association between PCB and immune dysfunction. These results indicate that PCB exposure is associated with differences in epigenetic marks in active regions of the genome, and future work should investigate whether these may mediate the association between PCB and health consequences.

Sex-specific DNA methylation differences in people exposed to polybrominated biphenyl.

Aim: Michigan residents were exposed to polybrominated biphenyls (PBBs) when it was accidentally added to the food supply. Highly exposed individuals report sex-specific health problems, but the underlying biological mechanism behind these different health risks is not known. Materials and methods: DNA methylation in blood from 381 women and 277 men with PBB exposure was analyzed with the MethylationEPIC BeadChip. Results: 675 CpGs were associated with PBBs levels in males, while only 17 CpGs were associated in females (false discovery rate <0.05). No CpGs were associated in both sexes. These CpGs were enriched in different functional regions and transcription factor binding sites in each sex. Conclusion: Exposure to PBBs may have sex-specific effects on the epigenome that may underlie sex-specific adverse health outcomes.