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.

Common genes associated with antidepressant response in mouse and man identify key role of glucocorticoid receptor sensitivity.

Response to antidepressant treatment in major depressive disorder (MDD) cannot be predicted currently, leading to uncertainty in medication selection, increasing costs, and prolonged suffering for many patients. Despite tremendous efforts in identifying response-associated genes in large genome-wide association studies, the results have been fairly modest, underlining the need to establish conceptually novel strategies. For the identification of transcriptome signatures that can distinguish between treatment responders and nonresponders, we herein submit a novel animal experimental approach focusing on extreme phenotypes. We utilized the large variance in response to antidepressant treatment occurring in DBA/2J mice, enabling sample stratification into subpopulations of good and poor treatment responders to delineate response-associated signature transcript profiles in peripheral blood samples. As a proof of concept, we translated our murine data to the transcriptome data of a clinically relevant human cohort. A cluster of 259 differentially regulated genes was identified when peripheral transcriptome profiles of good and poor treatment responders were compared in the murine model. Differences in expression profiles from baseline to week 12 of the human orthologues selected on the basis of the murine transcript signature allowed prediction of response status with an accuracy of 76% in the patient population. Finally, we show that glucocorticoid receptor (GR)-regulated genes are significantly enriched in this cluster of antidepressant-response genes. Our findings point to the involvement of GR sensitivity as a potential key mechanism shaping response to antidepressant treatment and support the hypothesis that antidepressants could stimulate resilience-promoting molecular mechanisms. Our data highlight the suitability of an appropriate animal experimental approach for the discovery of treatment response-associated pathways across species.

The Atherosclerosis Risk Variant rs2107595 Mediates Allele-Specific Transcriptional Regulation of HDAC9 via E2F3 and Rb1.

Background and Purpose- Genome-wide association studies have identified the HDAC9 (histone deacetylase 9) gene region as a major risk locus for atherosclerotic stroke and coronary artery disease in humans. Previous results suggest a role of altered HDAC9 expression levels as the underlying disease mechanism. rs2107595, the lead single nucleotide polymorphism for stroke and coronary artery disease resides in noncoding DNA and colocalizes with histone modification marks suggestive of enhancer elements. Methods- To determine the mechanisms by which genetic variation at rs2107595 regulates HDAC9 expression and thus vascular risk we employed targeted resequencing, proteome-wide search for allele-specific nuclear binding partners, chromatin immunoprecipitation, genome-editing, reporter assays, circularized chromosome conformation capture, and gain- and loss-of-function experiments in cultured human cell lines and primary immune cells. Results- Targeted resequencing of the HDAC9 locus in patients with atherosclerotic stroke and controls supported candidacy of rs2107595 as the causative single nucleotide polymorphism. A proteomic search for nuclear binding partners revealed preferential binding of the E2F3/TFDP1/Rb1 complex (E2F transcription factor 3/transcription factor Dp-1/Retinoblastoma 1) to the rs2107595 common allele, consistent with the disruption of an E2F3 consensus site by the risk allele. Gain- and loss-of-function studies showed a regulatory effect of E2F/Rb proteins on HDAC9 expression. Compared with the common allele, the rs2107595 risk allele exhibited higher transcriptional capacity in luciferase assays and was associated with higher HDAC9 mRNA levels in primary macrophages and genome-edited Jurkat cells. Circularized chromosome conformation capture revealed a genomic interaction of the rs2107595 region with the HDAC9 promoter, which was stronger for the common allele as was the in vivo interaction with E2F3 and Rb1 determined by chromatin immunoprecipitation. Gain-of-function experiments in isogenic Jurkat cells demonstrated a key role of E2F3 in mediating rs2107595-dependent transcriptional regulation of HDAC9. Conclusions- Collectively, our findings imply allele-specific transcriptional regulation of HDAC9 via E2F3 and Rb1 as a major mechanism mediating vascular risk at rs2107595.

Eosinophilia and reduced STAT3 signaling affect neutrophil cell death in autosomal-dominant Hyper-IgE syndrome.

The autosomal-dominant hyper-IgE syndrome (HIES), caused by mutations in STAT3, is a rare primary immunodeficiency that predisposes to mucocutaneous candidiasis and staphylococcal skin and lung infections. This infection phenotype is suggestive of defects in neutrophils, but data on neutrophil functions in HIES are inconsistent. This study was undertaken to functionally characterize neutrophils in STAT3-deficient HIES patients and to analyze whether the patients` eosinophilia affects the neutrophil phenotype in S. aureus infection. Neutrophil functions and cell death kinetics were studied in eight STAT3-deficient patients. Moreover, the response of STAT3-deficient neutrophils to S. aureus and the impact of autologous eosinophils on pathogen-induced cell death were analyzed. No specific aberrations in neutrophil functions were detected within this cohort. However, the half-life of STAT3-deficient neutrophils ex vivo was reduced, which was partially attributable to the presence of eosinophils. Increased S. aureus-induced cell lysis, dependent on the staphylococcal virulence controlling accessory gene regulator (agr)-locus, was observed in STAT3-deficient neutrophils and upon addition of eosinophils. Accelerated neutrophil cell death kinetics may underlie the reported variability in neutrophil function testing in HIES. Increased S. aureus-induced lysis of STAT3-deficient neutrophils might affect pathogen control and contribute to tissue destruction during staphylococcal infections in HIES.

Evidence for oestrogen sensitivity in perinatal depression: pharmacological sex hormone manipulation study.

BACKGROUND: Enhanced sensitivity to oestrogen signalling may drive increased risk for depressive symptoms when exposed to peripartum sex-steroid hormone fluctuations. AIM: Testing if 116 pre-identified sex steroid-responsive transcripts that predicted perinatal depression (PND) translates to a pharmacological model of hormone-induced mood changes. METHOD: We generated longitudinal, genome-wide gene-expression and DNA-methylation data from 60 women exposed to a gonadotrophin-releasing hormone agonist (GnRHa) or placebo. We used linear mixed-effect models to assess differences between baseline and follow-up for gene expression and DNA methylation in the biphasic ovarian response to GnRHa. RESULTS: Of the 116 PND-predictive transcripts, a significant (19%) overlap was observed with those differentially expressed post-GnRHa at both early and later follow-up, indicating sustained effects. Similarly, 49% of tested genes were differentially methylated post-GnRHa at the late follow-up. Within the GnRHa group, a large proportion of PND genes were significantly associated (gene expression; DNA methylation) with changes in depressive symptoms (28%; 66%), oestradiol levels (49%; 66%) and neocortex serotonin transporter binding (8%; 45%) between baseline and follow-up. CONCLUSIONS: Our data bridge clinical PND biomarkers with a pharmacological model of sex hormone-induced mood changes and directly relate oestrogen-induced biological changes with depressive symptoms and associated serotonin-signalling changes. Our data highlight that individual variations in molecular sensitivity to oestrogen associate with susceptibility to hormone-induced mood changes and hold promise for candidate biomarkers. DECLARATION OF INTEREST: V.G.F. received honorarium for being a speaker for H. Lundbeck A/S. E.B.B. receives research funding from Böhringer Ingelheim to investigate FKBP5 as a potential drug target for depression.

Childhood maltreatment is associated with distinct genomic and epigenetic profiles in posttraumatic stress disorder.

Childhood maltreatment is likely to influence fundamental biological processes and engrave long-lasting epigenetic marks, leading to adverse health outcomes in adulthood. We aimed to elucidate the impact of different early environment on disease-related genome-wide gene expression and DNA methylation in peripheral blood cells in patients with posttraumatic stress disorder (PTSD). Compared with the same trauma-exposed controls (n = 108), gene-expression profiles of PTSD patients with similar clinical symptoms and matched adult trauma exposure but different childhood adverse events (n = 32 and 29) were almost completely nonoverlapping (98%). These differences on the level of individual transcripts were paralleled by the enrichment of several distinct biological networks between the groups. Moreover, these gene-expression changes were accompanied and likely mediated by changes in DNA methylation in the same loci to a much larger proportion in the childhood abuse (69%) vs. the non-child abuse-only group (34%). This study is unique in providing genome-wide evidence of distinct biological modifications in PTSD in the presence or absence of exposure to childhood abuse. The findings that nonoverlapping biological pathways seem to be affected in the two PTSD groups and that changes in DNA methylation appear to have a much greater impact in the childhood-abuse group might reflect differences in the pathophysiology of PTSD, in dependence of exposure to childhood maltreatment. These results contribute to a better understanding of the extent of influence of differences in trauma exposure on pathophysiological processes in stress-related psychiatric disorders and may have implications for personalized medicine.

Cohort profile: BioMD-Y (biopsychosocial factors of major depression in youth) - a biobank study on the molecular genetics and environmental factors of depression in children and adolescents in Munich.

PURPOSE: BioMD-Y is a comprehensive biobank study of children and adolescents with major depression (MD) and their healthy peers in Germany, collecting a host of both biological and psychosocial information from the participants and their parents with the aim of exploring genetic and environmental risk and protective factors for MD in children and adolescents. PARTICIPANTS: Children and adolescents aged 8-18 years are recruited to either the clinical case group (MD, diagnosis of MD disorder) or the typically developing control group (absence of any psychiatric condition). FINDINGS TO DATE: To date, four publications on both genetic and environmental risk and resilience factors (including FKBP5, glucocorticoid receptor activation, polygenic risk scores, psychosocial and sociodemographic risk and resilience factors) have been published based on the BioMD-Y sample. FUTURE PLANS: Data collection is currently scheduled to continue into 2026. Research questions will be further addressed using available measures.

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.

Polygenic risk for schizophrenia converges on alternative polyadenylation as molecular mechanism underlying synaptic impairment.

Schizophrenia (SCZ) is a genetically heterogenous psychiatric disorder of highly polygenic nature. Correlative evidence from genetic studies indicate that the aggregated effects of distinct genetic risk factor combinations found in each patient converge onto common molecular mechanisms. To prove this on a functional level, we employed a reductionistic cellular model system for polygenic risk by differentiating induced pluripotent stem cells (iPSCs) from 104 individuals with high polygenic risk load and controls into cortical glutamatergic neurons (iNs). Multi-omics profiling identified widespread differences in alternative polyadenylation (APA) in the 3' untranslated region of many synaptic transcripts between iNs from SCZ patients and healthy donors. On the cellular level, 3'APA was associated with a reduction in synaptic density of iNs. Importantly, differential APA was largely conserved between postmortem human prefrontal cortex from SCZ patients and healthy donors, and strongly enriched for transcripts related to synapse biology. 3'APA was highly correlated with SCZ polygenic risk and affected genes were significantly enriched for SCZ associated common genetic variation. Integrative functional genomic analysis identified the RNA binding protein and SCZ GWAS risk gene PTBP2 as a critical trans-acting factor mediating 3'APA of synaptic genes in SCZ subjects. Functional characterization of PTBP2 in iNs confirmed its key role in 3'APA of synaptic transcripts and regulation of synapse density. Jointly, our findings show that the aggregated effects of polygenic risk converge on 3'APA as one common molecular mechanism that underlies synaptic impairments in SCZ.

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.

vipR: variant identification in pooled DNA using R.

MOTIVATION: High-throughput-sequencing (HTS) technologies are the method of choice for screening the human genome for rare sequence variants causing susceptibility to complex diseases. Unfortunately, preparation of samples for a large number of individuals is still very cost- and labor intensive. Thus, recently, screens for rare sequence variants were carried out in samples of pooled DNA, in which equimolar amounts of DNA from multiple individuals are mixed prior to sequencing with HTS. The resulting sequence data, however, poses a bioinformatics challenge: the discrimination of sequencing errors from real sequence variants present at a low frequency in the DNA pool. RESULTS: Our method vipR uses data from multiple DNA pools in order to compensate for differences in sequencing error rates along the sequenced region. More precisely, instead of aiming at discriminating sequence variants from sequencing errors, vipR identifies sequence positions that exhibit significantly different minor allele frequencies in at least two DNA pools using the Skellam distribution. The performance of vipR was compared with three other models on data from a targeted resequencing study of the TMEM132D locus in 600 individuals distributed over four DNA pools. Performance of the methods was computed on SNPs that were also genotyped individually using a MALDI-TOF technique. On a set of 82 sequence variants, vipR achieved an average sensitivity of 0.80 at an average specificity of 0.92, thus outperforming the reference methods by at least 0.17 in specificity at comparable sensitivity. AVAILABILITY: The code of vipR is freely available via: http://sourceforge.net/projects/htsvipr/ CONTACT: altmann@mpipsykl.mpg.de.

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.

Cohort profile: InTraUterine sampling in early pregnancy (ITU), a prospective pregnancy cohort study in Finland: study design and baseline characteristics.

PURPOSE: The InTraUterine sampling in early pregnancy (ITU) is a prospective pregnancy cohort study. The overarching aim of ITU is to unravel genomic, epigenomic, transcriptomic, endocrine, inflammatory and metabolic maternal-placental-fetal mechanisms involved in the programming of health and disease after exposure to prenatal environmental adversity, such as maternal malnutrition, cardiometabolic disorders, infections, medical interventions, mental disorders and psychosocial stress. This paper describes the study protocol, design and baseline characteristics of the cohort. PARTICIPANTS: We included 944 pregnant Finnish women, their partners and children born alive between April 2012 and December 2017. The women were recruited through the national, voluntary trisomy 21 screening between 9+0 and 21+6 gestational weeks. Of the participating women, 543 were screen positive and underwent fetal chromosomal testing. Test result of these women suggested no fetal chromosomal abnormality. Further, we recruited 401 women who were screen negative and who did not undergo fetal chromosomal testing. FINDINGS TO DATE: We have collected chorionic villi and amniotic fluid from the screen-positive women; blood, urine, buccal swabs and diurnal salivary samples from all women; blood and buccal swabs from all partners; and placenta, cord blood and buccal swabs from all newborns for analyses of the genome, epigenome, transcriptome, and endocrine, inflammatory and metabolic markers. These data are coupled with comprehensive phenotypes, including questions on demographic characteristics, health and well-being of the women and their partners during pregnancy and of the women and their children at the child's age of 1.7 and 3 years. Data also come from patient records and nationwide registers covering health, lifestyle and medication data. FUTURE PLANS: Multiple layers of ITU data allow integrative data analyses, which translate to biomarker identification and allow risk stratification and understanding of the biological mechanisms involved in prenatal programming of health and disease.

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.

The pediatric buccal epigenetic clock identifies significant ageing acceleration in children with internalizing disorder and maltreatment exposure.

BACKGROUND: Studies reporting accelerated ageing in children with affective disorders or maltreatment exposure have relied on algorithms for estimating epigenetic age derived from adult samples. These algorithms have limited validity for epigenetic age estimation during early development. We here use a pediatric buccal epigenetic (PedBE) clock to predict DNA methylation-based ageing deviation in children with and without internalizing disorder and assess the moderating effect of maltreatment exposure. We further conduct a gene set enrichment analysis to assess the contribution of glucocorticoid signaling to PedBE clock-based results. METHOD: DNA was isolated from saliva of 158 children [73 girls, 85 boys; mean age (SD) = 4.25 (0.8) years] including children with internalizing disorder and maltreatment exposure. Epigenetic age was estimated based on DNA methylation across 94 CpGs of the PedBE clock. Residuals of epigenetic age regressed against chronological age were contrasted between children with and without internalizing disorder. Maltreatment was coded in 3 severity levels and entered in a moderation model. Genome-wide dexamethasone-responsive CpGs were derived from an independent sample and enrichment of these CpGs within the PedBE clock was identified. RESULTS: Children with internalizing disorder exhibited significant acceleration of epigenetic ageing as compared to children without internalizing disorder (F1,147 = 6.67, p = .011). This association was significantly moderated by maltreatment severity (b = 0.49, 95% CI [0.073, 0.909], t = 2.322, p = .022). Children with internalizing disorder who had experienced maltreatment exhibited ageing acceleration relative to children with no internalizing disorder (1-2 categories: b = 0.50, 95% CI [0.170, 0.821], t = 3.008, p = .003; 3 or more categories: b = 0.99, 95% CI [0.380, 1.593], t = 3.215, p = .002). Children with internalizing disorder who were not exposed to maltreatment did not show epigenetic ageing acceleration. There was significant enrichment of dexamethasone-responsive CpGs within the PedBE clock (OR = 4.36, p = 1.65*10-6). Among the 94 CpGs of the PedBE clock, 18 (19%) were responsive to dexamethasone. CONCLUSION: Using the novel PedBE clock, we show that internalizing disorder is associated with accelerated epigenetic ageing in early childhood. This association is moderated by maltreatment severity and may, in part, be driven by glucocorticoids. Identifying developmental drivers of accelerated epigenetic ageing after maltreatment will be critical to devise early targeted interventions.

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.