A Glucocorticoid-Sensitive Hippocampal Gene Network Moderates the Impact of Early-Life Adversity on Mental Health Outcomes

BACKGROUND

Early stress increases the risk for psychiatric disorders. Glucocorticoids are stress mediators that regulate transcriptional activity and morphology in the hippocampus, which is implicated in the pathophysiology of multiple psychiatric conditions. We aimed to establish the relevance of hippocampal glucocorticoid-induced transcriptional activity as a mediator of the effects of early life on later psychopathology in humans.

METHODS

RNA sequencing was performed with anterior and posterior hippocampal dentate gyrus from adult female macaques (n = 12/group) that were chronically treated with betamethasone (glucocorticoid receptor agonist) or vehicle. Coexpression network analysis identified a preserved gene network in the posterior hippocampal dentate gyrus that was strongly associated with glucocorticoid exposure. The single nucleotide polymorphisms in the genes in this network were used to create an expression-based polygenic score in humans.

RESULTS

The expression-based polygenic score significantly moderated the association between early adversity and psychotic disorders in adulthood (UK Biobank, women, n = 44,519) and on child peer relations (ALSPAC [Avon Longitudinal Study of Parents and Children], girls, n = 1666 for 9-year-olds and n = 1594 for 11-year-olds), an endophenotype for later psychosis. Analyses revealed that this network was enriched for glucocorticoid-induced epigenetic remodeling in human hippocampal cells. We also found a significant association between single nucleotide polymorphisms from the expression-based polygenic score and adult brain gray matter density.

CONCLUSIONS

We provide an approach for the use of transcriptomic data from animal models together with human data to study the impact of environmental influences on mental health. The results are consistent with the hypothesis that hippocampal glucocorticoid-related transcriptional activity mediates the effects of early adversity on neural mechanisms implicated in psychiatric disorders.

Associations of psychiatric disease and ageing with FKBP5 expression converge on superficial layer neurons of the neocortex

Identification and characterisation of novel targets for treatment is a priority in the field of psychiatry. FKBP5 is a gene with decades of evidence suggesting its pathogenic role in a subset of psychiatric patients, with potential to be leveraged as a therapeutic target for these individuals. While it is widely reported that FKBP5/FKBP51 mRNA/protein (FKBP5/1) expression is impacted by psychiatric disease state, risk genotype and age, it is not known in which cell types and sub-anatomical areas of the human brain this occurs. This knowledge is critical to propel FKBP5/1-targeted treatment development. Here, we performed an extensive, large-scale postmortem study (n = 1024) of FKBP5/1, examining neocortical areas (BA9, BA11 and ventral BA24/BA24a) derived from subjects that lived with schizophrenia, major depression or bipolar disorder. With an extensive battery of RNA (bulk RNA sequencing, single-nucleus RNA sequencing, microarray, qPCR, RNAscope) and protein (immunoblot, immunohistochemistry) analysis approaches, we thoroughly investigated the effects of disease state, ageing and genotype on cortical FKBP5/1 expression including in a cell type-specific manner. We identified consistently heightened FKBP5/1 levels in psychopathology and with age, but not genotype, with these effects strongest in schizophrenia. Using single-nucleus RNA sequencing (snRNAseq; BA9 and BA11) and targeted histology (BA9, BA24a), we established that these disease and ageing effects on FKBP5/1 expression were most pronounced in excitatory superficial layer neurons of the neocortex, and this effect appeared to be consistent in both the granular and agranular areas examined. We then found that this increase in FKBP5 levels may impact on synaptic plasticity, as FKBP5 gex levels strongly and inversely correlated with dendritic mushroom spine density and brain-derived neurotrophic factor (BDNF) levels in superficial layer neurons in BA11. These findings pinpoint a novel cellular and molecular mechanism that has potential to open a new avenue of FKBP51 drug development to treat cognitive symptoms in psychiatric disorders.

Cell-Type-Specific Impact of Glucocorticoid Receptor Activation on the Developing Brain: A Cerebral Organoid Study

OBJECTIVE

A fine-tuned balance of glucocorticoid receptor (GR) activation is essential for organ formation, with disturbances influencing many health outcomes. In utero, glucocorticoids have been linked to brain-related negative outcomes, with unclear underlying mechanisms, especially regarding cell-type-specific effects. An in vitro model of fetal human brain development, induced human pluripotent stem cell (hiPSC)-derived cerebral organoids, was used to test whether cerebral organoids are suitable for studying the impact of prenatal glucocorticoid exposure on the developing brain.

METHODS

The GR was activated with the synthetic glucocorticoid dexamethasone, and the effects were mapped using single-cell transcriptomics across development.

RESULTS

The GR was expressed in all cell types, with increasing expression levels through development. Not only did its activation elicit translocation to the nucleus and the expected effects on known GR-regulated pathways, but also neurons and progenitor cells showed targeted regulation of differentiation- and maturation-related transcripts. Uniquely in neurons, differentially expressed transcripts were significantly enriched for genes associated with behavior-related phenotypes and disorders. This human neuronal glucocorticoid response profile was validated across organoids from three independent hiPSC lines reprogrammed from different source tissues from both male and female donors.

CONCLUSIONS

These findings suggest that excessive glucocorticoid exposure could interfere with neuronal maturation in utero, leading to increased disease susceptibility through neurodevelopmental processes at the interface of genetic susceptibility and environmental exposure. Cerebral organoids are a valuable translational resource for exploring the effects of glucocorticoids on early human brain development.

Effects of stressful life-events on DNA methylation in panic disorder and major depressive disorder

BACKGROUND

Panic disorder (PD) is characterized by recurrent panic attacks and higher affection of women as compared to men. The lifetime prevalence of PD is about 2-3% in the general population leading to tremendous distress and disability. Etiologically, genetic and environmental factors, such as stress, contribute to the onset and relapse of PD. In the present study, we investigated epigenome-wide DNA methylation (DNAm) in respond to a cumulative, stress-weighted life events score (wLE) in patients with PD and its boundary to major depressive disorder (MDD), frequently co-occurring with symptoms of PD.

METHODS

DNAm was assessed by the Illumina HumanMethylation450 BeadChip. In a meta-analytic approach, epigenome-wide DNAm changes in association with wLE were first analyzed in two PD cohorts (with a total sample size of 183 PD patients and 85 healthy controls) and lastly in 102 patients with MDD to identify possible overlapping and opposing effects of wLE on DNAm. Additionally, analysis of differentially methylated regions (DMRs) was conducted to identify regional clusters of association.

RESULTS

Two CpG-sites presented with p-values below 1 × 10^-05 in PD: cg09738429 (p = 6.40 × 10^-06, located in an intergenic shore region in next proximity of PYROXD1) and cg03341655 (p = 8.14 × 10^-06, located in the exonic region of GFOD2). The association of DNAm at cg03341655 and wLE could be replicated in the independent MDD case sample indicating a diagnosis independent effect. Genes mapping to the top hits were significantly upregulated in brain and top hits have been implicated in the metabolic system. Additionally, two significant DMRs were identified for PD only on chromosome 10 and 18, including CpG-sites which have been reported to be associated with anxiety and other psychiatric phenotypes.

CONCLUSION

This first DNAm analysis in PD reveals first evidence of small but significant DNAm changes in PD in association with cumulative stress-weighted life events. Most of the top associated CpG-sites are located in genes implicated in metabolic processes supporting the hypothesis that environmental stress contributes to health damaging changes by affecting a broad spectrum of systems in the body.

Polygenic risk for immuno-metabolic markers and specific depressive symptoms: A multi-sample network analysis study

BACKGROUND

About every fourth patient with major depressive disorder (MDD) shows evidence of systemic inflammation. Previous studies have shown inflammation-depression associations of multiple serum inflammatory markers and multiple specific depressive symptoms. It remains unclear, however, if these associations extend to genetic/lifetime predisposition to higher inflammatory marker levels and what role metabolic factors such as Body Mass Index (BMI) play. It is also unclear whether inflammation-symptom associations reflect direct or indirect associations, which can be disentangled using network analysis.

METHODS

This study examined associations of polygenic risk scores (PRSs) for immuno-metabolic markers (C-reactive protein [CRP], interleukin [IL]-6, IL-10, tumour necrosis factor [TNF]-α, BMI) with seven depressive symptoms in one general population sample, the UK Biobank study (n = 110,010), and two patient samples, the Munich Antidepressant Response Signature (MARS, n = 1058) and Sequenced Treatment Alternatives to Relieve Depression (STAR*D, n = 1143) studies. Network analysis was applied jointly for these samples using fused graphical least absolute shrinkage and selection operator (FGL) estimation as primary analysis and, individually, using unregularized model search estimation. Stability of results was assessed using bootstrapping and three consistency criteria were defined to appraise robustness and replicability of results across estimation methods, network bootstrapping, and samples.

RESULTS

Network analysis results displayed to-be-expected PRS-PRS and symptom-symptom associations (termed edges), respectively, that were mostly positive. Using FGL estimation, results further suggested 28, 29, and six PRS-symptom edges in MARS, STAR*D, and UK Biobank samples, respectively. Unregularized model search estimation suggested three PRS-symptom edges in the UK Biobank sample. Applying our consistency criteria to these associations indicated that only the association of higher CRP PRS with greater changes in appetite fulfilled all three criteria. Four additional associations fulfilled at least two consistency criteria; specifically, higher CRP PRS was associated with greater fatigue and reduced anhedonia, higher TNF-α PRS was associated with greater fatigue, and higher BMI PRS with greater changes in appetite and anhedonia. Associations of the BMI PRS with anhedonia, however, showed an inconsistent valence across estimation methods.

CONCLUSIONS

Genetic predisposition to higher systemic inflammatory markers are primarily associated with somatic/neurovegetative symptoms of depression such as changes in appetite and fatigue, consistent with previous studies based on circulating levels of inflammatory markers. We extend these findings by providing evidence that associations are direct (using network analysis) and extend to genetic predisposition to immuno-metabolic markers (using PRSs). Our findings can inform selection of patients with inflammation-related symptoms into clinical trials of immune-modulating drugs for MDD.

Versatile knowledge guided network inference method for prioritizing key regulatory factors in multi-omics data

Constantly decreasing costs of high-throughput profiling on many molecular levels generate vast amounts of multi-omics data. Studying one biomedical question on two or more omic levels provides deeper insights into underlying molecular processes or disease pathophysiology. For the majority of multi-omics data projects, the data analysis is performed level-wise, followed by a combined interpretation of results. Hence the full potential of integrated data analysis is not leveraged yet, presumably due to the complexity of the data and the lacking toolsets. We propose a versatile approach, to perform a multi-level fully integrated analysis: The Knowledge guIded Multi-Omics Network inference approach, KiMONo (https://github.com/cellmapslab/kimono). KiMONo performs network inference by using statistical models for combining omics measurements coupled to a powerful knowledge-guided strategy exploiting prior information from existing biological sources. Within the resulting multimodal network, nodes represent features of all input types e.g. variants and genes while edges refer to knowledge-supported and statistically derived associations. In a comprehensive evaluation, we show that our method is robust to noise and exemplify the general applicability to the full spectrum of multi-omics data, demonstrating that KiMONo is a powerful approach towards leveraging the full potential of data sets for detecting biomarker candidates.

Dissecting the Association Between Inflammation, Metabolic Dysregulation, and Specific Depressive Symptoms: A Genetic Correlation and 2-Sample Mendelian Randomization Study

IMPORTANCE

Observational studies highlight associations of C-reactive protein (CRP), a general marker of inflammation, and interleukin 6 (IL-6), a cytokine-stimulating CRP production, with individual depressive symptoms. However, it is unclear whether inflammatory activity is associated with individual depressive symptoms and to what extent metabolic dysregulation underlies the reported associations.

OBJECTIVE

To explore the genetic overlap and associations between inflammatory activity, metabolic dysregulation, and individual depressive symptoms.

GWAS DATA SOURCES

Genome-wide association study (GWAS) summary data of European individuals, including the following: CRP levels (204 402 individuals); 9 individual depressive symptoms (3 of which did not differentiate between underlying diametrically opposite symptoms [eg, insomnia and hypersomnia]) as measured with the Patient Health Questionnaire 9 (up to 117 907 individuals); summary statistics for major depression, including and excluding UK Biobank participants, resulting in sample sizes of 500 199 and up to 230 214 individuals, respectively; insomnia (up to 386 533 individuals); body mass index (BMI) (up to 322 154 individuals); and height (up to 253 280 individuals).

DESIGN

In this genetic correlation and 2-sample mendelian randomization (MR) study, linkage disequilibrium score (LDSC) regression was applied to infer single-nucleotide variant-based heritability and genetic correlation estimates. Two-sample MR tested potential causal associations of genetic variants associated with CRP levels, IL-6 signaling, and BMI with depressive symptoms. The study dates were November 2019 to April 2020.

RESULTS

Based on large GWAS data sources, genetic correlation analyses revealed consistent false discovery rate (FDR)-controlled associations (genetic correlation range, 0.152-0.362; FDR P = .006 to P < .001) between CRP levels and depressive symptoms that were similar in size to genetic correlations of BMI with depressive symptoms. Two-sample MR analyses suggested that genetic upregulation of IL-6 signaling was associated with suicidality (estimate [SE], 0.035 [0.010]; FDR plus Bonferroni correction P = .01), a finding that remained stable across statistical models and sensitivity analyses using alternative instrument selection strategies. Mendelian randomization analyses did not consistently show associations of higher CRP levels or IL-6 signaling with other depressive symptoms, but higher BMI was associated with anhedonia, tiredness, changes in appetite, and feelings of inadequacy.

CONCLUSIONS AND RELEVANCE

This study reports coheritability between CRP levels and individual depressive symptoms, which may result from the potentially causal association of metabolic dysregulation with anhedonia, tiredness, changes in appetite, and feelings of inadequacy. The study also found that IL-6 signaling is associated with suicidality. These findings may have clinical implications, highlighting the potential of anti-inflammatory approaches, especially IL-6 blockade, as a putative strategy for suicide prevention.

A polyepigenetic glucocorticoid exposure score at birth and childhood mental and behavioral disorders

BACKGROUND

Maternal depression and anxiety during pregnancy may enhance fetal exposure to glucocorticoids (GCs) and harm neurodevelopment. We tested whether a novel cross-tissue polyepigenetic biomarker indicative of in utero exposure to GC is associated with mental and behavioral disorders and their severity in children, possibly mediating the associations between maternal prenatal depressive and anxiety symptoms and these child outcomes.

METHODS

Children (n = 814) from the Prediction and Prevention of Preeclampsia and Intrauterine Growth Restriction (PREDO) study were followed-up from birth to age 7.1-10.7 years. A weighted polyepigenetic GC exposure score was calculated based on the methylation profile of 24 CpGs from umbilical cord blood. Child diagnosis of mental and behavioral disorder (n = 99) and its severity, defined as the number of days the child had received treatment (all 99 had received outpatient treatment and 8 had been additionally in inpatient treatment) for mental or behavioral disorder as the primary diagnosis, came from the Care Register for Health Care. Mothers (n = 408) reported on child total behavior problems at child's age of 2.3-5.8 years and their own depressive and anxiety symptoms during pregnancy (n = 583).

RESULTS

The fetal polyepigenetic GC exposure score at birth was not associated with child hazard of mental and behavioral disorder (HR = 0.82, 95% CI 0.54; 1.24, p = 0.35) or total behavior problems (unstandardized beta = -0.10, 95% CI -0.31; 0.10, p = 0.33). However, for one standard deviation decrease in the polyepigenetic score, the child had spent 2.94 (95%CI 1.59; 5.45, p < 0.001) more days in inpatient or outpatient treatment with any mental and behavioral disorder as the primary diagnosis. Criteria for mediation tests were not met.

CONCLUSIONS

These findings suggest that fetal polyepigenetic GC exposure score at birth was not associated with any mental or behavioral disorder diagnosis or mother-rated total behavior problems, but it may contribute to identifying children at birth who are at risk for more severe mental or behavioral disorders.

DeepWAS: Multivariate genotype-phenotype associations by directly integrating regulatory information using deep learning

Genome-wide association studies (GWAS) identify genetic variants associated with traits or diseases. GWAS never directly link variants to regulatory mechanisms. Instead, the functional annotation of variants is typically inferred by post hoc analyses. A specific class of deep learning-based methods allows for the prediction of regulatory effects per variant on several cell type-specific chromatin features. We here describe "DeepWAS", a new approach that integrates these regulatory effect predictions of single variants into a multivariate GWAS setting. Thereby, single variants associated with a trait or disease are directly coupled to their impact on a chromatin feature in a cell type. Up to 61 regulatory SNPs, called dSNPs, were associated with multiple sclerosis (MS, 4,888 cases and 10,395 controls), major depressive disorder (MDD, 1,475 cases and 2,144 controls), and height (5,974 individuals). These variants were mainly non-coding and reached at least nominal significance in classical GWAS. The prediction accuracy was higher for DeepWAS than for classical GWAS models for 91% of the genome-wide significant, MS-specific dSNPs. DSNPs were enriched in public or cohort-matched expression and methylation quantitative trait loci and we demonstrated the potential of DeepWAS to generate testable functional hypotheses based on genotype data alone. DeepWAS is available at https://github.com/cellmapslab/DeepWAS.

The Role of m6A/m-RNA Methylation in Stress Response Regulation

N⁶-methyladenosine (m⁶A) and N⁶,2′-O-dimethyladenosine (m⁶Am) are abundant mRNA modifications that regulate transcript processing and translation. The role of both, here termed m⁶A/m, in the stress response in the adult brain in vivo is currently unknown. Here, we provide a detailed analysis of the stress epitranscriptome using m⁶A/m-seq, global and gene-specific m⁶A/m measurements. We show that stress exposure and glucocorticoids region and time specifically alter m⁶A/m and its regulatory network. We demonstrate that deletion of the methyltransferase Mettl3 or the demethylase Fto in adult neurons alters the m⁶A/m epitranscriptome, increases fear memory, and changes the transcriptome response to fear and synaptic plasticity. Moreover, we report that regulation of m⁶A/m is impaired in major depressive disorder patients following glucocorticoid stimulation. Our findings indicate that brain m⁶A/m represents a novel layer of complexity in gene expression regulation after stress and that dysregulation of the m⁶A/m response may contribute to the pathophysiology of stress-related psychiatric disorders.

•

m⁶A/m mRNA methylation in the adult mouse brain is regulated by stress

•

m⁶A/m mRNA regulation is brain region, time, and gene specific

•

Mettl3 and Fto cKO alter m⁶A/m, fear memory, expression, and synaptic plasticity

•

The m⁶A/m glucocorticoid response is impaired in major depressive disorder patients

Identification of dynamic glucocorticoid-induced methylation changes at the FKBP5 locus

Background

Epigenetic mechanisms may play a major role in the biological embedding of early-life stress (ELS). One proposed mechanism is that glucocorticoid (GC) release following ELS exposure induces long-lasting alterations in DNA methylation (DNAm) of important regulatory genes of the stress response. Here, we investigate the dynamics of GC-dependent methylation changes in key regulatory regions of the FKBP5 locus in which ELS-associated DNAm changes have been reported.

Results

We repeatedly measured DNAm in human peripheral blood samples from 2 independent cohorts exposed to the GC agonist dexamethasone (DEX) using a targeted bisulfite sequencing approach, complemented by data from Illumina 450K arrays. We detected differentially methylated CpGs in enhancers co-localizing with GC receptor binding sites after acute DEX treatment (1 h, 3 h, 6 h), which returned to baseline levels within 23 h. These changes withstood correction for immune cell count differences. While we observed main effects of sex, age, body mass index, smoking, and depression symptoms on FKBP5 methylation levels, only the functional FKBP5 SNP (rs1360780) moderated the dynamic changes following DEX. This genotype effect was observed in both cohorts and included sites previously shown to be associated with ELS.

Conclusion

Our study highlights that DNAm levels within regulatory regions of the FKBP5 locus show dynamic changes following a GC challenge and suggest that factors influencing the dynamics of this regulation may contribute to the previously reported alterations in DNAm associated with current and past ELS exposure.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0682-5) contains supplementary material, which is available to authorized users.

Stress dynamically regulates co-expression networks of glucocorticoid receptor-dependent MDD and SCZ risk genes

Early-life adversity is an important risk factor for major depressive disorder (MDD) and schizophrenia (SCZ) that interacts with genetic factors to confer disease risk through mechanisms that are still insufficiently understood. One downstream effect of early-life adversity is the activation of glucocorticoid receptor (GR)-dependent gene networks that drive acute and long-term adaptive behavioral and cellular responses to stress. We have previously shown that genetic variants that moderate GR-induced gene transcription (GR-response eSNPs) are significantly enriched among risk variants from genome-wide association studies (GWASs) for MDD and SCZ. Here, we show that the 63 transcripts regulated by these disease-associated functional genetic variants form a tight glucocorticoid-responsive co-expression network (termed GCN). We hypothesized that changes in the correlation structure of this GCN may contribute to early-life adversity-associated disease risk. Therefore, we analyzed the effects of different qualities of social support and stress throughout life on GCN formation across distinct brain regions using a translational mouse model. We observed that different qualities of social experience substantially affect GCN structure in a highly brain region-specific manner. GCN changes were predominantly found in two functionally interconnected regions, the ventral hippocampus and the hypothalamus, two brain regions previously shown to be of relevance for the stress response, as well as psychiatric disorders. Overall, our results support the hypothesis that a subset of genetic variants may contribute to risk for MDD and SCZ by altering circuit-level effects of early and adult social experiences on GCN formation and structure.

Impact on birth weight of maternal smoking throughout pregnancy mediated by DNA methylation

Background

Cigarette smoking has severe adverse health consequences in adults and in the offspring of mothers who smoke during pregnancy. One of the most widely reported effects of smoking during pregnancy is reduced birth weight which is in turn associated with chronic disease in adulthood. Epigenome-wide association studies have revealed that smokers show a characteristic “smoking methylation pattern”, and recent authors have proposed that DNA methylation mediates the impact of maternal smoking on birth weight. The aims of the present study were to replicate previous reports that methylation mediates the effect of maternal smoking on birth weight, and for the first time to investigate whether the observed mediation effects are sex-specific in order to account for known sex-specific differences in methylation levels.

Methods

Methylation levels in the cord blood of 313 newborns were determined using the Illumina HumanMethylation450K Beadchip. A total of 5,527 CpG sites selected on the basis of evidence from the literature were tested. To determine whether the observed association between maternal smoking and birth weight was attributable to methylation, mediation analyses were performed for significant CpG sites. Separate analyses were then performed in males and females.

Results

Following quality control, 282 newborns eventually remained in the analysis. A total of 25 mothers had smoked consistently throughout the pregnancy. The birthweigt of newborns whose mothers had smoked throughout pregnancy was reduced by >200g. After correction for multiple testing, 30 CpGs showed differential methylation in the maternal smoking subgroup including top “smoking methylation pattern” genes AHRR, MYO1G, GFI1, CYP1A1, and CNTNAP2. The effect of maternal smoking on birth weight was partly mediated by the methylation of cg25325512 (PIM1); cg25949550 (CNTNAP2); and cg08699196 (ITGB7). Sex-specific analyses revealed a mediating effect for cg25949550 (CNTNAP2) in male newborns.

Conclusion

The present data replicate previous findings that methylation can mediate the effect of maternal smoking on birth weight. The analysis of sex-dependent mediation effects suggests that the sex of the newborn may have an influence. Larger studies are warranted to investigate the role of both the identified differentially methylated loci and the sex of the newborn in mediating the association between maternal smoking during pregnancy and birth weight.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4652-7) contains supplementary material, which is available to authorized users.

Anxiety Associated Increased CpG Methylation in the Promoter of Asb1: A Translational Approach Evidenced by Epidemiological and Clinical Studies and a Murine Model

Epigenetic regulation in anxiety is suggested, but evidence from large studies is needed. We conducted an epigenome-wide association study (EWAS) on anxiety in a population-based cohort and validated our finding in a clinical cohort as well as a murine model. In the KORA cohort, participants (n=1522, age 32-72 years) were administered the Generalized Anxiety Disorder (GAD-7) instrument, whole blood DNA methylation was measured (Illumina 450K BeadChip), and circulating levels of hs-CRP and IL-18 were assessed in the association between anxiety and methylation. DNA methylation was measured using the same instrument in a study of patients with anxiety disorders recruited at the Max Planck Institute of Psychiatry (MPIP, 131 non-medicated cases and 169 controls). To expand our mechanistic understanding, these findings were reverse translated in a mouse model of acute social defeat stress. In the KORA study, participants were classified according to mild, moderate, or severe levels of anxiety (29.4%/6.0%/1.5%, respectively). Severe anxiety was associated with 48.5% increased methylation at a single CpG site (cg12701571) located in the promoter of the gene encoding Asb1 (β-coefficient=0.56 standard error (SE)=0.10, p (Bonferroni)=0.005), a protein hypothetically involved in regulation of cytokine signaling. An interaction between IL-18 and severe anxiety with methylation of this CpG cite showed a tendency towards significance in the total population (p=0.083) and a significant interaction among women (p=0.014). Methylation of the same CpG was positively associated with Panic and Agoraphobia scale (PAS) scores (β=0.005, SE=0.002, p=0.021, n=131) among cases in the MPIP study. In a murine model of acute social defeat stress, Asb1 gene expression was significantly upregulated in a tissue-specific manner (p=0.006), which correlated with upregulation of the neuroimmunomodulating cytokine interleukin 1 beta. Our findings suggest epigenetic regulation of the stress-responsive Asb1 gene in anxiety-related phenotypes. Further studies are necessary to elucidate the causal direction of this association and the potential role of Asb1-mediated immune dysregulation in anxiety disorders.

"DNA Methylation signatures in panic disorder"

Panic disorder (PD) affects about four million Europeans, with women affected twice as likely as men, causing substantial suffering and high economic costs. The etiopathogenesis of PD remains largely unknown, but both genetic and environmental factors contribute to risk. An epigenome-wide association study (EWAS) was conducted to compare medication-free PD patients (n = 89) with healthy controls (n = 76) stratified by gender. Replication was sought in an independent sample (131 cases, 169 controls) and functional analyses were conducted in a third sample (N = 71). DNA methylation was assessed in whole blood using the Infinium HumanMethylation450 BeadChip. One genome-wide association surviving FDR of 5% (cg07308824, P = 1.094 × 10-7, P-adj = 0.046) was identified in female PD patients (N = 49) compared to controls (N = 48). The same locus, located in an enhancer region of the HECA gene, was also hypermethylated in female PD patients in the replication sample (P = 0.035) and the significance of the association improved in the meta-analysis (P-adj = 0.004). Methylation at this CpG site was associated with HECA mRNA expression in another independent female sample (N = 71) both at baseline (P = 0.046) and after induction by dexamethasone (P = 0.029). Of 15 candidates, 5 previously reported as associated with PD or anxiety traits also showed differences in DNA methylation after gene-wise correction and included SGK1, FHIT, ADCYAP1, HTR1A, HTR2A. Our study examines epigenome-wide differences in peripheral blood for PD patients. Our results point to possible sex-specific methylation changes in the HECA gene for PD but overall highlight that this disorder is not associated with extensive changes in DNA methylation in peripheral blood.

Rapporteur summaries of plenary, symposia, and oral sessions from the XXIIIrd World Congress of Psychiatric Genetics Meeting in Toronto, Canada, 16-20 October 2015

The XXIIIrd World Congress of Psychiatric Genetics meeting, sponsored by the International Society of Psychiatric Genetics, was held in Toronto, ON, Canada, on 16-20 October 2015. Approximately 700 participants attended to discuss the latest state-of-the-art findings in this rapidly advancing and evolving field. The following report was written by trainee travel awardees. Each was assigned one session as a rapporteur. This manuscript represents the highlights and topics that were covered in the plenary sessions, symposia, and oral sessions during the conference, and contains major notable and new findings.

Time-dependent effects of dexamethasone plasma concentrations on glucocorticoid receptor challenge tests

Glucocorticoid challenge tests such as the dexamethasone suppression test (DST) and the combined dexamethasone/corticotropin-releasing hormone (dex-CRH) test are considered to be able to sensitively measure hypothalamic-pituitary-adrenal (HPA) axis activity in stress-related psychiatric and endocrine disorders. We used mass-spectrometry to assess the relationship of plasma dexamethasone concentrations and the outcome of these tests in two independent cohorts. Dexamethasone concentrations were measured after oral ingestion of 1.5mg dexamethasone in two cohorts that underwent a standard (dexamethasone at 23:00h) as well as modified (18:00h) DST and dex-CRH test. The first study population was a case/control cohort of 105 depressed patients and 133 controls in which peripheral blood mRNA expression was also measured. The second was a cohort of 261 depressed patients that underwent a standard dex-CRH test at baseline and after 12 weeks' treatment with cognitive-behavioral therapy or antidepressants. Dexamethasone concentrations explained significant proportions of the variance in the DST in both the first (24.6%) and the second (5.2%) cohort. Dexamethasone concentrations explained a higher proportion of the variance in the dex-CRH test readouts, with 41.9% of the cortisol area under the curve (AUC) in the first sample and 24.7% in the second sample. In contrast to these strong effects at later time points, dexamethasone concentrations did not impact cortisol or ACTH concentrations or mRNA expression 3hours after ingestion. In the second sample, dexamethasone concentrations at baseline and week 12 were highly correlated, independent of treatment type and response status. Importantly, a case/control effect in the Dex-CRH test was only apparent when controlling for dexamethasone concentrations. Our results suggest that the incorporation of plasma dexamethasone concentration or measures of earlier endocrine read-outs may help to improve the assessment of endocrine dysfunction in depression.

Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis.

Lifetime stress accelerates epigenetic aging in an urban, African American cohort: relevance of glucocorticoid signaling

Background

Chronic psychological stress is associated with accelerated aging and increased risk for aging-related diseases, but the underlying molecular mechanisms are unclear.

Results

We examined the effect of lifetime stressors on a DNA methylation-based age predictor, epigenetic clock. After controlling for blood cell-type composition and lifestyle parameters, cumulative lifetime stress, but not childhood maltreatment or current stress alone, predicted accelerated epigenetic aging in an urban, African American cohort (n = 392). This effect was primarily driven by personal life stressors, was more pronounced with advancing age, and was blunted in individuals with higher childhood abuse exposure. Hypothesizing that these epigenetic effects could be mediated by glucocorticoid signaling, we found that a high number (n = 85) of epigenetic clock CpG sites were located within glucocorticoid response elements. We further examined the functional effects of glucocorticoids on epigenetic clock CpGs in an independent sample with genome-wide DNA methylation (n = 124) and gene expression data (n = 297) before and after exposure to the glucocorticoid receptor agonist dexamethasone. Dexamethasone induced dynamic changes in methylation in 31.2 % (110/353) of these CpGs and transcription in 81.7 % (139/170) of genes neighboring epigenetic clock CpGs. Disease enrichment analysis of these dexamethasone-regulated genes showed enriched association for aging-related diseases, including coronary artery disease, arteriosclerosis, and leukemias.

Conclusions

Cumulative lifetime stress may accelerate epigenetic aging, an effect that could be driven by glucocorticoid-induced epigenetic changes. These findings contribute to our understanding of mechanisms linking chronic stress with accelerated aging and heightened disease risk.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0828-5) contains supplementary material, which is available to authorized users.

Re-Annotator: Annotation Pipeline for Microarray Probe Sequences

Microarray technologies are established approaches for high throughput gene expression, methylation and genotyping analysis. An accurate mapping of the array probes is essential to generate reliable biological findings. However, manufacturers of the microarray platforms typically provide incomplete and outdated annotation tables, which often rely on older genome and transcriptome versions that differ substantially from up-to-date sequence databases. Here, we present the Re-Annotator, a re-annotation pipeline for microarray probe sequences. It is primarily designed for gene expression microarrays but can also be adapted to other types of microarrays. The Re-Annotator uses a custom-built mRNA reference database to identify the positions of gene expression array probe sequences. We applied Re-Annotator to the Illumina Human-HT12 v4 microarray platform and found that about one quarter (25%) of the probes differed from the manufacturer's annotation. In further computational experiments on experimental gene expression data, we compared Re-Annotator to another probe re-annotation tool, ReMOAT, and found that Re-Annotator provided an improved re-annotation of microarray probes. A thorough re-annotation of probe information is crucial to any microarray analysis. The Re-Annotator pipeline is freely available at http://sourceforge.net/projects/reannotator along with re-annotated files for Illumina microarrays HumanHT-12 v3/v4 and MouseRef-8 v2.

Genetic Differences in the Immediate Transcriptome Response to Stress Predict Risk-Related Brain Function and Psychiatric Disorders

Depression risk is exacerbated by genetic factors and stress exposure; however, the biological mechanisms through which these factors interact to confer depression risk are poorly understood. One putative biological mechanism implicates variability in the ability of cortisol, released in response to stress, to trigger a cascade of adaptive genomic and non-genomic processes through glucocorticoid receptor (GR) activation. Here, we demonstrate that common genetic variants in long-range enhancer elements modulate the immediate transcriptional response to GR activation in human blood cells. These functional genetic variants increase risk for depression and co-heritable psychiatric disorders. Moreover, these risk variants are associated with inappropriate amygdala reactivity, a transdiagnostic psychiatric endophenotype and an important stress hormone response trigger. Network modeling and animal experiments suggest that these genetic differences in GR-induced transcriptional activation may mediate the risk for depression and other psychiatric disorders by altering a network of functionally related stress-sensitive genes in blood and brain.

An analysis of gene expression in PTSD implicates genes involved in the glucocorticoid receptor pathway and neural responses to stress

We examined the association between posttraumatic stress disorder (PTSD) and gene expression using whole blood samples from a cohort of trauma-exposed white non-Hispanic male veterans (115 cases and 28 controls). 10,264 probes of genes and gene transcripts were analyzed. We found 41 that were differentially expressed in PTSD cases versus controls (multiple-testing corrected p<0.05). The most significant was DSCAM, a neurological gene expressed widely in the developing brain and in the amygdala and hippocampus of the adult brain. We then examined the 41 differentially expressed genes in a meta-analysis using two replication cohorts and found significant associations with PTSD for 7 of the 41 (p<0.05), one of which (ATP6AP1L) survived multiple-testing correction. There was also broad evidence of overlap across the discovery and replication samples for the entire set of genes implicated in the discovery data based on the direction of effect and an enrichment of p<0.05 significant probes beyond what would be expected under the null. Finally, we found that the set of differentially expressed genes from the discovery sample was enriched for genes responsive to glucocorticoid signaling with most showing reduced expression in PTSD cases compared to controls.

A polymorphism in the Crhr1 gene determines stress vulnerability in male mice

Chronic stress is a risk factor for psychiatric disorders but does not necessarily lead to uniform long-term effects on mental health, suggesting modulating factors such as genetic predispositions. Here we address the question whether natural genetic variations in the mouse CRH receptor 1 (Crhr1) locus modulate the effects of adolescent chronic social stress (ACSS) on long-term stress hormone dysregulation in outbred CD1 mice, which allows a better understanding of the currently reported genes × environment interactions of early trauma and CRHR1 in humans. We identified 2 main haplotype variants in the mouse Crhr1 locus that modulate the long-term effects of ACSS on basal hypothalamic-pituitary-adrenal axis activity. This effect is likely mediated by higher levels of CRHR1, because Crhr1 mRNA expression and CRHR1 binding were enhanced in risk haplotype carriers. Furthermore, a CRHR1 receptor antagonist normalized these long-term effects. Deep sequencing of the Crhr1 locus in CD1 mice revealed a large number of linked single-nucleotide polymorphisms with some located in important regulatory regions, similar to the location of human CRHR1 variants implicated in modulating gene × stress exposure interactions. Our data support that the described gene × stress exposure interaction in this animal model is based on naturally occurring genetic variations in the Crhr1 gene associated with enhanced CRHR1-mediated signaling. Our results suggest that patients with a specific genetic predisposition in the CRHR1 gene together with an exposure to chronic stress may benefit from a treatment selectively antagonizing CRHR1 hyperactivity.

Dexamethasone stimulated gene expression in peripheral blood indicates glucocorticoid-receptor hypersensitivity in job-related exhaustion

Work-related stress can lead to various health problems ranging from job-related exhaustion to psychiatric and somatic diseases. Biomarkers of job-related exhaustion could help to improve our understanding of the biological mechanisms and might be useful to guide prevention and treatment strategies. The present study included 12 male cases suffering from job-related exhaustion and 12 matched healthy controls. Severity of exhaustion was assessed with the Maslach Burnout Inventory (MBI) and the Shirom-Melamed Burnout Measure (SMBM). Whole genome expression profiles derived from whole blood cells (baseline and following glucocorticoid-receptor (GR) stimulation with 1.5mg dexamethasone p.o.) and corresponding plasma cortisol levels were analyzed. All cases participated in regular aerobic exercise for 12 consecutive weeks and were then re-assessed at follow-up for exhaustion symptoms as well as for cortisol levels and gene expression profiles. At baseline, we found increased basal cortisol levels and an enhanced suppression of plasma cortisol concentrations following dexamethasone in cases suffering from job-related exhaustion. Gene expression analysis revealed that 1.6-fold more transcripts were significantly regulated by dexamethasone in cases as compared to controls. At follow-up after 12 weeks of regular exercise training which was accompanied by significantly improved exhaustion severity scores, cortisol levels and gene expression profiles of cases normalized to the levels observed in controls. In conclusion, we detected GR-induced neuroendocrine and gene expression changes in cases suffering from job-related exhaustion which are in line with an increased sensitivity of GR function. This GR dysregulation normalized with symptom recovery.

Functional coding variants in SLC6A15, a possible risk gene for major depression

SLC6A15 is a neuron-specific neutral amino acid transporter that belongs to the solute carrier 6 gene family. This gene family is responsible for presynaptic re-uptake of the majority of neurotransmitters. Convergent data from human studies, animal models and pharmacological investigations suggest a possible role of SLC6A15 in major depressive disorder. In this work, we explored potential functional variants in this gene that could influence the activity of the amino acid transporter and thus downstream neuronal function and possibly the risk for stress-related psychiatric disorders. DNA from 400 depressed patients and 400 controls was screened for genetic variants using a pooled targeted re-sequencing approach. Results were verified by individual re-genotyping and validated non-synonymous coding variants were tested in an independent sample (N = 1934). Nine variants altering the amino acid sequence were then assessed for their functional effects by measuring SLC6A15 transporter activity in a cellular uptake assay. In total, we identified 405 genetic variants, including twelve non-synonymous variants. While none of the non-synonymous coding variants showed significant differences in case-control associations, two rare non-synonymous variants were associated with a significantly increased maximal (3)H proline uptake as compared to the wildtype sequence. Our data suggest that genetic variants in the SLC6A15 locus change the activity of the amino acid transporter and might thus influence its neuronal function and the risk for stress-related psychiatric disorders. As statistically significant association for rare variants might only be achieved in extremely large samples (N >70,000) functional exploration may shed light on putatively disease-relevant variants.

Rare variants in TMEM132D in a case-control sample for panic disorder

Genome-wide association studies have identified common variants associated with common diseases. Most variants, however, explain only a small proportion of the estimated heritability, suggesting that rare variants might contribute to a larger extent to common diseases than assumed to date. Here, we use next-generation sequencing to test whether such variants contribute to the risk for anxiety disorders by re-sequencing 40 kb including all exons of the TMEM132D locus which we have previously shown to be associated with panic disorder and anxiety severity measures. DNA from 300 patients suffering from anxiety disorders, mostly panic disorder (84.7%), and 300 healthy controls was screened for the presence of genetic variants using next-generation re-sequencing in a pooled approach. Results were verified by individual re-genotyping. We identified 371 variants of which 247 had not been reported before, including 15 novel non-synonymous variants. The majority, 76% of these variants had a minor allele frequency less than 5%. While we did not identify additional common variants in TMEM132D associated with panic disorders, we observed an overrepresentation of presumably functional coding variants in healthy controls as compared to cases as well as a higher rate of private coding variants in cases, with one non-synonymous coding variant present in four patients but not in any of the matched controls nor in over 5,500 individuals of different ethnic origins from publicly available re-sequencing datasets. Our data suggest that not only common but also putatively functional and/or rare variants within TMEM132D might contribute to the risk to develop anxiety disorders.