Varun Kilaru

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.

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.

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.

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.

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.

Neonatal DNA methylation patterns associate with gestational age

Risk for adverse neonatal outcome increases with declining gestational age (GA), and changes in DNA methylation may contribute to the relationship between GA and adverse health outcomes in offspring. To test this hypothesis, we evaluated the association between GA and more than 27,000 CpG sites in neonatal DNA extracted from umbilical cord blood from two prospectively-characterized cohorts: (1) a discovery cohort consisting of 259 neonates from women with a history of neuropsychiatric disorders and (2) a replication cohort consisting of 194 neonates of uncomplicated mothers. GA was determined by obstetrician report and maternal last menstrual period. The associations between proportion of DNA methylated and GA were evaluated by fitting a separate linear mixed effects model for each CpG site, adjusting for relevant covariates including neonatal sex, race, parity, birth weight percentile and chip effects. CpG sites in 39 genes were associated with GA (false discovery rate < 0.05) in the discovery cohort. The same CpG sites in 25 of these genes replicated in the replication cohort, with each association replicating in the same direction. Notably, these CpG sites were located in genes previously implicated in labor and delivery (e.g., AVP, OXT, CRHBP and ESR1) or that may influence the risk for adverse health outcomes later in life (e.g., DUOX2, TMEM176A and CASP8). All associations were independent of method of delivery or induction of labor. These results suggest neonatal DNA methylation varies with GA even in term deliveries. The potential contribution of these changes to clinically significant postnatal outcomes warrants further investigation.

DNA methylation in neonates born to women receiving psychiatric care

Prenatal exposure both to maternal psychiatric illness and psychiatric medication has been linked with adverse child outcomes that affect physiological, emotional and psychiatric development. Studies suggest that epigenetic mechanisms, such as DNA methylation, may facilitate these effects. In this report, we explore the association between maternal psychiatric illness and treatment during pregnancy and neonatal DNA methylation patterns in a prospectively-characterized clinical cohort of 201 dyads. Associations between the percent of umbilical cord blood DNA methylated at 27,578 CpG sites and maternal psychiatric diagnosis, symptoms and antidepressant use were evaluated by fitting a separate linear mixed effects model for each CpG site. There were no significant changes in neonatal DNA methylation attributable to maternal psychiatric diagnosis or depressive symptoms during pregnancy. Exposure to an antidepressant medication was associated with differential methylation of CpG sites in TNFRSF21 and CHRNA2 (false discovery rate < 0.05), but the average difference in methylation for both CpG sites was less than 3% between each group. The results were not specific to type of antidepressant or duration of the exposure. This study suggests that there are no large effects of maternal psychiatric illness, depressive symptoms or prenatal exposure to antidepressants on neonatal DNA methylation. Delineation of the influence of maternal psychiatric illness and pharmacological exposures on the developing fetuses has critical implications for clinical care during pregnancy.

Fetal DNA Methylation Associates with Early Spontaneous Preterm Birth and Gestational Age

Spontaneous preterm birth (PTB, <37 weeks gestation) is a major public health concern, and children born preterm have a higher risk of morbidity and mortality throughout their lives. Recent studies suggest that fetal DNA methylation of several genes varies across a range of gestational ages (GA), but it is not yet clear if fetal epigenetic changes associate with PTB. The objective of this study is to interrogate methylation patterns across the genome in fetal leukocyte DNA from African Americans with early PTB (241/7–340/7 weeks; N = 22) or term births (390/7–406/7weeks; N = 28) and to evaluate the association of each CpG site with PTB and GA. DNA methylation was assessed across the genome with the HumanMethylation450 BeadChip. For each individual sample and CpG site, the proportion of DNA methylation was estimated. The associations between methylation and PTB or GA were evaluated by fitting a separate linear model for each CpG site, adjusting for relevant covariates. Overall, 29 CpG sites associated with PTB (FDR<.05; 5.7×10−10<p<2.9×10−6) independent of GA. Also, 9637 sites associated with GA (FDR<.05; 9.5×10−16<p<1.0×10−3), with 61.8% decreasing in methylation with shorter GA. GA-associated CpG sites were depleted in the CpG islands of their respective genes (p<2.2×10−16). Gene set enrichment analysis (GSEA) supported enrichment of GA-associated CpG sites in genes that play a role in embryonic development as well as the extracellular matrix. Additionally, this study replicated the association of several CpG sites associated with gestational age in other studies (CRHBP, PIK3CD and AVP). Dramatic differences in fetal DNA methylation are evident in fetuses born preterm versus at term, and the patterns established at birth may provide insight into the long-term consequences associated with PTB.

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.

An epigenetic clock for gestational age at birth based on blood methylation data

BackgroundGestational age is often used as a proxy for developmental maturity by clinicians and researchers alike. DNA methylation has previously been shown to be associated with age and has been used to accurately estimate chronological age in children and adults. In the current study, we examine whether DNA methylation in cord blood can be used to estimate gestational age at birth.ResultsWe find that gestational age can be accurately estimated from DNA methylation of neonatal cord blood and blood spot samples. We calculate a DNA methylation gestational age using 148 CpG sites selected through elastic net regression in six training datasets. We evaluate predictive accuracy in nine testing datasets and find that the accuracy of the DNA methylation gestational age is consistent with that of gestational age estimates based on established methods, such as ultrasound. We also find that an increased DNA methylation gestational age relative to clinical gestational age is associated with birthweight independent of gestational age, sex, and ancestry.ConclusionsDNA methylation can be used to accurately estimate gestational age at or near birth and may provide additional information relevant to developmental stage. Further studies of this predictor are warranted to determine its utility in clinical settings and for research purposes. When clinical estimates are available this measure may increase accuracy in the testing of hypotheses related to developmental age and other early life circumstances.