High-throughput screening of glucocorticoid-induced enhancer activity reveals mechanisms of stress-related psychiatric disorders.

Exposure to stressful life events increases the risk for psychiatric disorders. Mechanistic insight into the genetic factors moderating the impact of stress can increase our understanding of disease processes. Here, we test 3,662 single nucleotide polymorphisms (SNPs) from preselected expression quantitative trait loci in massively parallel reporter assays to identify genetic variants that modulate the activity of regulatory elements sensitive to glucocorticoids, important mediators of the stress response. Of the tested SNP sequences, 547 were located in glucocorticoid-responsive regulatory elements of which 233 showed allele-dependent activity. Transcripts regulated by these functional variants were enriched for those differentially expressed in psychiatric disorders in the postmortem brain. Phenome-wide Mendelian randomization analysis in 4,439 phenotypes revealed potentially causal associations specifically in neurobehavioral traits, including major depression and other psychiatric disorders. Finally, a functional gene score derived from these variants was significantly associated with differences in the physiological stress response, suggesting that these variants may alter disease risk by moderating the individual set point of the stress response.

DNA methylation patterns of FKBP5 regulatory regions in brain and blood of humanized mice and humans.

Humanized mouse models can be used to explore human gene regulatory elements (REs), which frequently lie in non-coding and less conserved genomic regions. Epigenetic modifications of gene REs, also in the context of gene x environment interactions, have not yet been explored in humanized mouse models. We applied high-accuracy measurement of DNA methylation (DNAm) via targeted bisulfite sequencing (HAM-TBS) to investigate DNAm in three tissues/brain regions (blood, prefrontal cortex and hippocampus) of mice carrying the human FK506-binding protein 5 (FKBP5) gene, an important candidate gene associated with stress-related psychiatric disorders. We explored DNAm in three functional intronic glucocorticoid-responsive elements (at introns 2, 5, and 7) of FKBP5 at baseline, in cases of differing genotype (rs1360780 single nucleotide polymorphism), and following application of the synthetic glucocorticoid dexamethasone. We compared DNAm patterns in the humanized mouse (N = 58) to those in human peripheral blood (N = 447 and N = 89) and human postmortem brain prefrontal cortex (N = 86). Overall, DNAm patterns in the humanized mouse model seem to recapitulate DNAm patterns observed in human tissue. At baseline, this was to a higher extent in brain tissue. The animal model also recapitulated effects of dexamethasone on DNAm, especially in peripheral blood and to a lesser extent effects of genotype on DNAm. The humanized mouse model could thus assist in reverse translation of human findings in psychiatry that involve genetic and epigenetic regulation in non-coding elements.

Dissecting depression symptoms: Multi-omics clustering uncovers immune-related subgroups and cell-type specific dysregulation.

In a subset of patients with mental disorders, such as depression, low-grade inflammation and altered immune marker concentrations are observed. However, these immune alterations are often assessed by only one data type and small marker panels. Here, we used a transdiagnostic approach and combined data from two cohorts to define subgroups of depression symptoms across the diagnostic spectrum through a large-scale multi-omics clustering approach in 237 individuals. The method incorporated age, body mass index (BMI), 43 plasma immune markers and RNA-seq data from peripheral mononuclear blood cells (PBMCs). Our initial clustering revealed four clusters, including two immune-related depression symptom clusters characterized by elevated BMI, higher depression severity and elevated levels of immune markers such as interleukin-1 receptor antagonist (IL-1RA), C-reactive protein (CRP) and C-C motif chemokine 2 (CCL2 or MCP-1). In contrast, the RNA-seq data mostly differentiated a cluster with low depression severity, enriched in brain related gene sets. This cluster was also distinguished by electrocardiography data, while structural imaging data revealed differences in ventricle volumes across the clusters. Incorporating predicted cell type proportions into the clustering resulted in three clusters, with one showing elevated immune marker concentrations. The cell type proportion and genes related to cell types were most pronounced in an intermediate depression symptoms cluster, suggesting that RNA-seq and immune markers measure different aspects of immune dysregulation. Lastly, we found a dysregulation of the SERPINF1/VEGF-A pathway that was specific to dendritic cells by integrating immune marker and RNA-seq data. This shows the advantages of combining different data modalities and highlights possible markers for further stratification research of depression symptoms.

Methylation Patterns of the FKBP5 Gene in Association with Childhood Maltreatment and Depressive Disorders.

Childhood maltreatment is an important risk factor for adult depression and has been associated with changes in the hypothalamic pituitary adrenal (HPA) axis, including cortisol secretion and methylation of the FKBP5 gene. Furthermore, associations between depression and HPA changes have been reported. This study investigated the associations of whole-blood FKBP5 mRNA levels, serum cortisol levels, childhood maltreatment, and depressive symptoms with the whole-blood methylation status (assessed via target bisulfite sequencing) of 105 CpGs at the FKBP5 locus using data from the general population-based Study of Health in Pomerania (SHIP) (N = 203). Both direct and interaction effects with the rs1360780 single-nucleotide polymorphism were investigated. Nominally significant associations of main effects on methylation of a single CpG site were observed at intron 3, intron 7, and the 3'-end of the gene. Additionally, methylation at two clusters at the 3'-end and intron 7 were nominally associated with childhood maltreatment × rs1360780 and depressive symptoms × rs1360780, respectively. The results add to the understanding of molecular mechanisms underlying the emergence of depression and could aid the development of personalised depression therapy and drug development.

Extensive evaluation of DNA methylation of functional elements in the murine Fkbp5 locus using high-accuracy DNA methylation measurement via targeted bisulfite sequencing.

FKBP5 is an important stress-regulatory gene implicated in stress-related psychiatric diseases. Single nucleotide polymorphisms of the FKBP5 gene were shown to interact with early life stress to alter the glucocorticoid-related stress response and moderate disease risk. Demethylation of cytosine-phosphate-guanine-dinucleotides (CpGs) in regulatory glucocorticoid-responsive elements was suggested to be the mediating epigenetic mechanism for long-term stress effects, but studies on Fkbp5 DNA methylation (DNAm) in rodents are so far limited. We evaluated the applicability of high-accuracy DNA methylation measurement via targeted bisulfite sequencing (HAM-TBS), a next-generation sequencing-based technology, to allow a more in-depth characterisation of the DNA methylation of the murine Fkbp5 locus in three different tissues (blood, frontal cortex and hippocampus). In this study, we not only increased the number of evaluated sites in previously described regulatory regions (in introns 1 and 5), but also extended the evaluation to novel, possibly relevant regulatory regions of the gene (in intron 8, the transcriptional start site, the proximal enhancer and CTCF-binding sites within the 5'UTR). We here describe the assessment of HAM-TBS assays for a panel of 157 CpGs with possible functional relevance in the murine Fkbp5 gene. DNAm profiles were tissue-specific, with lesser differences between the two brain regions than between the brain and blood. Moreover, we identified DNAm changes in the Fkbp5 locus after early life stress exposure in the frontal cortex and blood. Our findings indicate that HAM-TBS is a valuable tool for broader exploration of the DNAm of the murine Fkbp5 locus and its involvement in the stress response.

Genome-Wide Copy Number Variant and High-Throughput Transcriptomics Analyses of Placental Tissues Underscore Persisting Child Susceptibility in At-Risk Pregnancies Cleared in Standard Genetic Testing.

Several studies have shown that children from pregnancies with estimated first-trimester risk based on fetal nuchal translucency thickness and abnormal maternal serum pregnancy protein and hormone levels maintain a higher likelihood of adverse outcomes, even if initial testing for known genetic conditions is negative. We used the Finnish InTraUterine cohort (ITU), which is a comprehensively characterized perinatal cohort consisting of 943 mothers and their babies followed throughout pregnancy and 18 months postnatally, including mothers shortlisted for prenatal genetic testing but cleared for major aneuploidies (cases: n = 544, 57.7%) and control pregnancies (n = 399, 42.3%). Using genome-wide genotyping and RNA sequencing of first-trimester and term placental tissue, combined with medical information from registry data and maternal self-report data, we investigated potential negative medical outcomes and genetic susceptibility to disease and their correlates in placenta gene expression. Case mothers did not present with higher levels of depression, perceived stress, or anxiety during pregnancy. Case children were significantly diagnosed more often with congenital malformations of the circulatory system (4.12 (95% CI [1.22−13.93]) higher hazard) and presented with significantly more copy number duplications as compared to controls (burden analysis, based on all copy number variants (CNVs) with at most 10% frequency, 823 called duplications in 297 cases versus 626 called duplications in 277 controls, p = 0.01). Fifteen genes showed differential gene expression (FDR < 0.1) in association with congenital malformations in first-trimester but not term placenta. These were significantly enriched for genes associated with placental dysfunction. In spite of normal routine follow-up prenatal testing results in early pregnancy, case children presented with an increased likelihood of negative outcomes, which should prompt vigilance in follow-up during pregnancy and after birth.

Single-nucleus transcriptomic profiling of human orbitofrontal cortex reveals convergent effects of aging and psychiatric disease.

Aging is a complex biological process and represents the largest risk factor for neurodegenerative disorders. The risk for neurodegenerative disorders is also increased in individuals with psychiatric disorders. Here, we characterized age-related transcriptomic changes in the brain by profiling ~800,000 nuclei from the orbitofrontal cortex from 87 individuals with and without psychiatric diagnoses and replicated findings in an independent cohort with 32 individuals. Aging affects all cell types, with LAMP5+LHX6+ interneurons, a cell-type abundant in primates, by far the most affected. Disrupted synaptic transmission emerged as a convergently affected pathway in aged tissue. Age-related transcriptomic changes overlapped with changes observed in Alzheimer's disease across multiple cell types. We find evidence for accelerated transcriptomic aging in individuals with psychiatric disorders and demonstrate a converging signature of aging and psychopathology across multiple cell types. Our findings shed light on cell-type-specific effects and biological pathways underlying age-related changes and their convergence with effects driven by psychiatric diagnosis.

Human cortical neurogenesis is altered via glucocorticoid-mediated regulation of ZBTB16 expression.

Glucocorticoids are important for proper organ maturation, and their levels are tightly regulated during development. Here, we use human cerebral organoids and mice to study the cell-type-specific effects of glucocorticoids on neurogenesis. We show that glucocorticoids increase a specific type of basal progenitors (co-expressing PAX6 and EOMES) that has been shown to contribute to cortical expansion in gyrified species. This effect is mediated via the transcription factor ZBTB16 and leads to increased production of neurons. A phenome-wide Mendelian randomization analysis of an enhancer variant that moderates glucocorticoid-induced ZBTB16 levels reveals causal relationships with higher educational attainment and altered brain structure. The relationship with postnatal cognition is also supported by data from a prospective pregnancy cohort study. This work provides a cellular and molecular pathway for the effects of glucocorticoids on human neurogenesis that relates to lasting postnatal phenotypes.

Transdiagnostic evaluation of epigenetic age acceleration and burden of psychiatric disorders.

Different psychiatric disorders as well as exposure to adverse life events have individually been associated with multiple age-related diseases and mortality. Age acceleration in different epigenetic clocks can serve as biomarker for such risk and could help to disentangle the interplay of psychiatric comorbidity and early adversity on age-related diseases and mortality. We evaluated five epigenetic clocks (Horvath, Hannum, PhenoAge, GrimAge and DunedinPoAm) in a transdiagnostic psychiatric sample using epigenome-wide DNA methylation data from peripheral blood of 429 subjects from two studies at the Max Planck Institute of Psychiatry. Burden of psychiatric disease, represented by a weighted score, was significantly associated with biological age acceleration as measured by GrimAge and DunedinPoAm (R2-adj. 0.22 and 0.33 for GrimAge and DunedinPoAm, respectively), but not the other investigated clocks. The relation of burden of psychiatric disease appeared independent of differences in socioeconomic status and medication. Our findings indicate that increased burden of psychiatric disease may associate with accelerated biological aging. This highlights the importance of medical management of patients with multiple psychiatric comorbidities and the potential usefulness of specific epigenetic clocks for early detection of risk and targeted intervention to reduce mortality in psychiatric patients.

FKBP5 intron 7 methylation is associated with higher anxiety proneness and smaller right thalamus volume in adolescents.

Dysregulation of stress response systems may mediate the detrimental effects of childhood trauma (CT) on mental health. FKBP5 regulates glucocorticoid receptor sensitivity and exerts pleiotropic effects on intracellular signaling, neurobiology and behavior. We investigated whether CT, alone and in combination with rs1360780 genotype, is associated with altered FKBP5 methylation and whether CT-associated methylation profiles are associated with anxiety proneness (AP) and structural brain volumes. Ninety-four adolescents completed the Childhood Trauma Questionnaire, and a composite AP score was generated from the Childhood Anxiety Sensitivity Index and the State-Trait Anxiety Inventory-Trait measure. Mean methylation values for 12 regulatory regions and 25 individual CpG sites were determined using high-accuracy measurement via targeted bisulfite sequencing. FKBP5 rs1360780 genotype and structural MRI data were available for a subset of participants (n = 71 and n = 75, respectively). Regression models revealed an inverse association between methylation of three intron 7 CpG sites (35558438, 35558566 and 35558710) and right thalamus volume. CpG35558438 methylation was positively associated with AP scores. Our data indicate that an intron 7 methylation profile, consistent with lower FKBP5 expression and elevated high sensitivity glucocorticoid receptor levels, is associated with higher AP and smaller right thalamus volume. Research into the mechanisms underlying the intron 7 methylation-thalamus volume relationship, and whether it confers increased risk for long-term psychopathology by altering the regulatory threshold of stress responding, is required.

DiffBrainNet: Differential analyses add new insights into the response to glucocorticoids at the level of genes, networks and brain regions.

Genome-wide gene expression analyses are invaluable tools for studying biological and disease processes, allowing a hypothesis-free comparison of expression profiles. Traditionally, transcriptomic analysis has focused on gene-level effects found by differential expression. In recent years, network analysis has emerged as an important additional level of investigation, providing information on molecular connectivity, especially for diseases associated with a large number of linked effects of smaller magnitude, like neuropsychiatric disorders. Here, we describe how combined differential expression and prior-knowledge-based differential network analysis can be used to explore complex datasets. As an example, we analyze the transcriptional responses following administration of the glucocorticoid/stress receptor agonist dexamethasone in 8 mouse brain regions important for stress processing. By applying a combination of differential network- and expression-analyses, we find that these explain distinct but complementary biological mechanisms of the glucocorticoid responses. Additionally, network analysis identifies new differentially connected partners of risk genes and can be used to generate hypotheses on molecular pathways affected. With DiffBrainNet (http://diffbrainnet.psych.mpg.de), we provide an analysis framework and a publicly available resource for the study of the transcriptional landscape of the mouse brain which can identify molecular pathways important for basic functioning and response to glucocorticoids in a brain-region specific manner.

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.

Childhood adversity correlates with stable changes in DNA methylation trajectories in children and converges with epigenetic signatures of prenatal stress.

Childhood maltreatment (CM) is an established major risk factor for a number of negative health outcomes later in life. While epigenetic mechanisms, such as DNA methylation (DNAm), have been proposed as a means of embedding this environmental risk factor, little is known about its timing and trajectory, especially in very young children. It is also not clear whether additional environmental adversities, often experienced by these children, converge on similar DNAm changes. Here, we calculated a cumulative adversity score, which additionally to CM includes socioeconomic status (SES), other life events, parental psychopathology and epigenetic biomarkers of prenatal smoking and alcohol consumption. We investigated the effects of CM alone as well as the adversity score on longitudinal DNAm trajectories in the Berlin Longitudinal Child Study. This is a cohort of 173 children aged 3-5 years at baseline of whom 86 were exposed to CM. These children were followed-up for 2 years with extensive psychometric and biological assessments as well as saliva collection at 5 time points providing genome-wide DNAm levels. Overall, only a few DNAm patterns were stable over this timeframe, but less than 10 DNAm regions showed significant changes. At baseline, neither CM nor the adversity score associated with DNAm changes. However, in 6 differentially methylated regions (DMRs), CM and the adversity score significantly moderated DNAm trajectories over time. A number of these DMRs have previously been associated with adverse prenatal exposures. In our study, children exposed to CM also presented with epigenetic signatures indicative of increased prenatal exposure to tobacco and alcohol, as compared to non-CM exposed children. These epigenetic signatures of prenatal exposure strongly correlate with DNAm regions associated with CM and the adversity score. Finally, weighted correlation network analysis revealed a module of CpGs exclusively associated with CM. While our study identifies DNAm loci specifically associated with CM, especially within long non-coding RNAs, the majority of associations were found with the adversity score with convergent association with indicators of adverse prenatal exposures. This study highlights the importance of mapping not only of the epigenome but also the exposome and extending the observational timeframe to well before birth.

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.

HAM-TBS: high-accuracy methylation measurements via targeted bisulfite sequencing.

BACKGROUND: The ability to accurately and efficiently measure DNA methylation is critical to advance the understanding of this epigenetic mechanism and its contribution to common diseases. Here, we present a highly accurate method to measure methylation using bisulfite sequencing (termed HAM-TBS). This novel method is able to assess DNA methylation in multiple samples with high accuracy in a cost-effective manner. We developed this assay for the FKBP5 locus, an important gene in the regulation of the stress system and previously linked to stress-related disorders, but the method is applicable to any locus of interest. RESULTS: HAM-TBS enables multiplexed analyses of up to 96 samples and regions spanning 10 kb using the Illumina MiSeq. It incorporates a triplicate bisulfite conversion step, pooled target enrichment via PCR, PCR-free library preparation and a minimum coverage of 1000×. TBS was able to resolve DNA methylation levels with a mean accuracy of 0.72%. Using this method, we designed and validated a targeted panel to specifically assess regulatory regions within the FKBP5 locus that are not covered in commercially available DNA methylation arrays. CONCLUSIONS: HAM-TBS represents a highly accurate, medium-throughput sequencing approach for robust detection of DNA methylation changes in specific target regions.