Despite extensive research on DNA methylation (DNAm) signatures associated with alcohol use disorder (AUD), findings are often inconsistent and not replicated. We conducted a large-scale meta-analysis of epigenome-wide association studies (EWAS) to identify reliable, reproducible epigenetic markers of AUD. Seven cohorts, comprising 3,775 individuals (1,325 with AUD), contributed to this meta-analysis within the framework of the Psychiatric Genomics Consortium Substance Use Disorders Epigenetics Working Group. Downstream analyses included the identification of differentially methylated regions, overrepresentation analyses, and the construction of a methylation risk score (MRS). We identified 118 significant CpG sites associated with AUD, with the strongest association found at cg24889777 ( p =5.12×10 -17 ) in the long non-coding RNA LOC100505942. CpG sites were enriched for pathways related to GTPase signaling and transmembrane transporter activity, as well as EWAS signals of alcohol consumption. The MRS explained 10.44% of variance in heavy drinking in an independent cohort (N=2,534, AUC=0.657). This large-scale meta-analysis offers key insights into the epigenetic mechanisms of AUD and lays the groundwork for future research on methylation risk scores for the diagnosis, prognosis, and treatment in AUD.

This study tests the association of whole-blood DNA methylation and antidepressant exposure in 16,531 individuals from Generation Scotland (GS), using self-report and prescription-derived measures. We identify 8 associations and a high concordance of results between self-report and prescription-derived measures. Sex-stratified analyses observe nominally significant increased effect estimates in females for four CpGs. There is observed enrichment for genes expressed in the Amygdala and annotated to synaptic vesicle membrane ontology. Two CpGs (cg15071067; DGUOK-AS1 and cg26277237; KANK1) show correlation between DNA methylation with the time in treatment. There is a significant overlap in the top 1% of CpGs with another independent methylome-wide association study of antidepressant exposure. Finally, a methylation profile score trained on this sample shows a significant association with antidepressant exposure in a meta-analysis of eight independent external datasets. In this large investigation of antidepressant exposure and DNA methylation, we demonstrate robust associations which warrant further investigation to inform on the design of more effective and tolerated treatments for depression.

Background Adverse life events (ALE) represent a major risk factor many psychiatric disorders, most notably major depressive disorder, anxiety, and post-traumatic stress disorder. However, the biological consequences of ALE are not well understood. Therefore, we examined the effects of ALE in a large, diverse population of postpartum women with RNA-sequencing and DNA methylation from whole blood. Methods Participants were recruited at six weeks postpartum from four obstetrical clinics in North Carolina. Postpartum depression status was determined using the MINI-plus. Participants were administered validated, self-report instruments to assess depression, adverse life events, and everyday life stressors. For each participant, genotyping, DNA methylation (Illumina 450k), and RNA-sequencing were performed. Results This population (n = 1,517) was diverse (68 % Black, 13 % Latina, and 18 % European). Participants reported high rates of abuse and trauma with 66.6 % reporting a history of at least one traumatic event. Two DNA methylation sites (CpGs) and no transcriptional changes were observed in whole blood. Cell specific effects were observed with B-cells exhibiting a majority of methylome- (64 CpGs) and transcriptome-wide significant results (412 transcripts representing 409 genes). Results from pathway and QTL analyses are underway and will be presented at the conference. Discussion This work highlights the impact adverse life events can have on an individual even years after their exposure. Furthermore, these events have cell specific effects at multiple levels of regulation (DNA methylation and RNA transcription). Identifying such biological signatures may act as biomarkers that can indicate potential response to treatments as part of trauma informed clinical care.

Adverse neonatal outcomes are a prevailing risk factor for both short- and long-term mortality and morbidity in infants. Given the importance of these outcomes, refining their assessment is paramount for improving prevention and care. Here we aim to enhance the assessment of these often correlated and multifaceted neonatal outcomes. To achieve this, we employ factor analysis to identify common and unique effects and further confirm these effects using criterion-related validity testing. This validation leverages methylome-wide profiles from neonatal blood. Specifically, we investigate nine neonatal health risk variables, including gestational age, Apgar score, three indicators of body size, jaundice, birth diagnosis, maternal preeclampsia, and maternal age. The methylomic profiles used for this research capture data from nearly all 28 million methylation sites in human blood, derived from the blood spot collected from 333 neonates, within 72 h post-birth. Our factor analysis revealed two common factors, size factor, that captured the shared effects of weight, head size, height, and gestational age and disease factor capturing the orthogonal shared effects of gestational age, combined with jaundice and birth diagnosis. To minimize false positives in the validation studies, validation was limited to variables with significant cumulative association as estimated through an in-sample replication procedure. This screening resulted in that the two common factors and the unique effects for gestational age, jaundice and Apgar were further investigated with full-scale cell-type specific methylome-wide association analyses. Highly significant, cell-type specific, associations were detected for both common effect factors and for Apgar. Gene Ontology analyses revealed multiple significant biologically relevant terms for the five fully investigated neonatal health risk variables. Given the established links between adverse neonatal outcomes and both immediate and long-term health, the distinct factor effects (representing the common and unique effects of the risk variables) and their biological profiles confirmed in our work, suggest their potential role as clinical biomarkers for assessing health risks and enhancing personalized care.

Objective: To test whether childhood mental health symptoms, substance use, and early adversity accelerate the rate of DNA methylation (DNAm) aging from adolescence to adulthood. Method: DNAm was assayed from blood samples in 381 participants in both adolescence (mean [SD] age = 13.9 [1.6] years) and adulthood (mean [SD] age = 25.9 [2.7] years). Structured diagnostic interviews were completed with participants and their parents at multiple childhood observations (1,950 total) to assess symptoms of common mental health disorders (attention-deficit/hyperactivity disorder, oppositional defiant disorder, conduct disorder, anxiety, and depression) and common types of substance use (alcohol, cannabis, nicotine) and early adversities. Results: Neither childhood mental health symptoms nor substance use variables were associated with DNAm aging cross-sectionally. In contrast, the following mental health symptoms and substance variables were associated with accelerated DNAm aging from adolescence to adulthood: depressive symptoms (b = 0.314, SE = 0.127, p = .014), internalizing symptoms (b = 0.108, SE = 0.049, p = .029), weekly cannabis use (b =1.665, SE = 0.591, p = .005), and years of weekly cannabis use (b = 0.718, SE = 0.283, p = .012). In models testing all individual variables simultaneously, the combined effect of the variables was equivalent to a potential difference of 3.17 to 3.76 years in DNAm aging. A final model tested a variable assessing cumulative exposure to mental health symptoms, substance use, and early adversities. This cumulative variable was strongly associated with accelerated aging (b = 0.126, SE = 0.044, p = .005). Conclusion: Mental health symptoms and substance use accelerated DNAm aging into adulthood in a manner consistent with a shared risk mechanism. Plain language summary: Using data from 381 participants in the Great Smoky Mountains Study, the authors examined whether childhood mental health symptoms, substance use, and early adversity accelerate biological aging, as measured by DNA methylation age, from adolescence to adulthood. Depressive symptoms and cannabis use were found to significantly accelerate biological aging. Models that tested the combined effect of mental health symptoms, substance use, and early adversity demonstrated that there was a shared effect across these types of childhood problems on accelerated aging.

Background. Cross-sectional studies have identified health risks associated with epigenetic aging. However, it is unclear whether these risks make epigenetic clocks 'tick faster' (i.e. accelerate biological aging). The current study examines concurrent and lagged within-person changes of a variety of health risks associated with epigenetic aging. Methods. Individuals from the Great Smoky Mountains Study were followed from age 9 to 35 years. DNA methylation profiles were assessed from blood, at multiple timepoints (i.e. waves) for each individual. Health risks were psychiatric, lifestyle, and adversity factors. Concurrent (N = 539 individuals; 1029 assessments) and lagged (N = 380 individuals; 760 assessments) analyses were used to determine the link between health risks and epigenetic aging. Results. Concurrent models showed that BMI (r = 0.15, P-FDR < 0.01) was significantly correlated to epigenetic aging at the subject-level but not wave-level. Lagged models demonstrated that depressive symptoms (b = 1.67 months per symptom, P-FDR = 0.02) in adolescence accelerated epigenetic aging in adulthood, also when models were fully adjusted for BMI, smoking, and cannabis and alcohol use. Conclusions. Within-persons, changes in health risks were unaccompanied by concurrent changes in epigenetic aging, suggesting that it is unlikely for risks to immediately 'accelerate' epigenetic aging. However, time lagged analyses indicated that depressive symptoms in childhood/adolescence predicted epigenetic aging in adulthood. Together, findings suggest that age-related biological embedding of depressive symptoms is not instant but provides prognostic opportunities. Repeated measurements and longer follow-up times are needed to examine stable and dynamic contributions of childhood experiences to epigenetic aging across the lifespan.

DNA methylation (DNAm) is one of the most reliable biomarkers of aging across mammalian tissues. While the age-dependent global loss of DNAm has been well characterized, DNAm gain is less characterized. Studies have demonstrated that CpGs which gain methylation with age are enriched in Polycomb Repressive Complex 2 (PRC2) targets. However, whole-genome examination of all PRC2 targets as well as determination of the pan-tissue or tissue-specific nature of these associations is lacking. Here, we show that low-methylated regions (LMRs) which are highly bound by PRC2 in embryonic stem cells (PRC2 LMRs) gain methylation with age in all examined somatic mitotic cells. We estimated that this epigenetic change represents around 90% of the age-dependent DNAm gain genome-wide. Therefore, we propose the "PRC2-AgeIndex," defined as the average DNAm in PRC2 LMRs, as a universal biomarker of cellular aging in somatic cells which can distinguish the effect of different anti-aging interventions.

Although it is widely known that various pharmaceuticals affect the methylome, the knowledge of the effects from anesthesia is limited, and nearly nonexistent regarding the effects of obstetric anesthesia on the newborn child. Using sequencing based-methylation data and a reference-based statistical deconvolution approach we performed methylome-wide association studies (MWAS) of neonatal whole blood, and for each cell-type specifically, to detect methylation variations that are associated with the pain relief administered to the mother during delivery. Significant findings were replicated in a different dataset and followed-up with gene ontology analysis to pinpoint biological functions of potential relevance to these neonatal methylation alterations. The MWAS analyses detected methylome-wide significant (q<0.1) alterations in the newborn for laughing gas in granulocytes (two CpGs, p<5.50x10(-9), q = 0.067), and for pudendal block in monocytes (five CpGs across three loci, p<1.51 x10(-8), q = 0.073). Suggestively significant findings (p<1.00x10(-6)) were detected for both treatments for bulk and all cell-types, and replication analyses showed consistent significant enrichment (odds ratios ranging 3.47-39.02; p<4.00x10(-4)) for each treatment, suggesting our results are robust. In contrast, we did not observe any overlap across treatments, suggesting that the treatments are associated with different alterations of the neonatal blood methylome. Gene ontology analyses of the replicating suggestively significant results indicated functions related to, for example, cell differentiation, intracellular membrane-bound organelles and calcium transport. In conclusion, for the first time, we investigated and detected effect of obstetric pain-relief on the blood methylome in the newborn child. The observed differences suggest that anesthetic treatment, such as laughing gas or pudendal block, may alter the neonatal methylome in a cell-type specific manner. Some of the observed alterations are part of gene ontology terms that previously have been suggested in relation to anesthetic treatment, supporting its potential role also in obstetric anesthesia.

Brain disorders are leading causes of disability worldwide. Gene expression studies provide promising opportunities to better understand their etiology but it is critical that expression is studied on a cell-type level. Cell-type specific association studies can be performed with bulk expression data using statistical methods that capitalize on cell-type proportions estimated with the help of a reference panel. To create a fine-grained reference panel for the human prefrontal cortex, we performed an integrated analysis of the seven largest single nucleus RNA-seq studies. Our panel included 17 cell-types that were robustly detected across all studies, subregions of the prefrontal cortex, and sex and age groups. To estimate the cell-type proportions, we used an empirical Bayes estimator that substantially outperformed three estimators recommended previously after a comprehensive evaluation of methods to estimate cell-type proportions from brain transcriptome data. This is important as being able to precisely estimate the cell-type proportions may avoid unreliable results in downstream analyses particularly for the multiple cell-types that had low abundances. Transcriptome-wide association studies performed with permuted bulk expression data showed that it is possible to perform transcriptome-wide association studies for even the rarest cell-types without an increased risk of false positives.

Gene expression studies offer promising opportunities to better understand the processes underlying alcohol use disorder (AUD). As cell types differ in their function, gene expression profiles will typically vary across cell types. When studying bulk tissue, failure to account for this cellular diversity has a detrimental impact on the ability to detect disease associations. We therefore assayed the transcriptomes of 32,531 individual nuclei extracted from the nucleus accumbens (NAc) of nine donors with AUD and nine controls (72% male). Our study identified 17 clearly delineated cell types. We detected 26 transcriptome-wide significant differentially expressed genes (DEGs) that mainly involved medium spiny neurons with both D1-type and D2-type dopamine receptors, microglia (MGL) and oligodendrocytes. A higher than expected number of DEGs replicated in an existing single nucleus gene expression study of alcohol dependence in the prefrontal cortex (enrichment ratio 1.91, p value 0.019) with two genes remaining significant after a Bonferroni correction. Our most compelling result involved CD53 in MGL that replicated in the same cell type in the prefrontal cortex and was previously implicated in studies of DNA methylation, bulk gene expression and genetic variants. Several DEGs were previously reported to be associated with AUD (e.g., PER1 and MGAT5). The DEGs for MSN.3 seemed involved in neurodegeneration, disruption of circadian rhythms, alterations in glucose metabolism and changes in synaptic plasticity. For MGL, the DEGs implicated neuroinflammation and immune-related processes and for OLI, disruptions in myelination. This identification of the specific cell-types from which the association signals originate will be key for designing proper follow-up experiments and, eventually, novel clinical interventions.

Schizophrenia is a disabling disorder involving genetic predisposition in combination with environmental influences that likely act via dynamic alterations of the epigenome and the transcriptome but its detailed pathophysiology is largely unknown. We performed cell-type specific methylome-wide association study of neonatal blood (N = 333) from individuals who later in life developed schizophrenia and controls. Suggestively significant associations (P < 1.0 × 10−6) were detected in all cell-types and in whole blood with methylome-wide significant associations in monocytes (P = 2.85 × 10−9–4.87 × 10−9), natural killer cells (P = 1.72 × 10−9–7.82 × 10−9) and B cells (P = 3.8 × 10−9). Validation of methylation findings in post-mortem brains (N = 596) from independent schizophrenia cases and controls showed significant enrichment of transcriptional differences (enrichment ratio = 1.98–3.23, P = 2.3 × 10−3–1.0 × 10−5), with specific highly significant differential expression for, for example, BDNF (t = −6.11, P = 1.90 × 10−9). In addition, expression difference in brain significantly predicted schizophrenia (multiple correlation = 0.15–0.22, P = 3.6 × 10−4–4.5 × 10−8). In summary, using a unique design combining pre-disease onset (neonatal) blood methylomic data and post-disease onset (post-mortem) brain transcriptional data, we have identified genes of likely functional relevance that are associated with schizophrenia susceptibility, rather than confounding disease associated artifacts. The identified loci may be of clinical value as a methylation-based biomarker for early detection of increased schizophrenia susceptibility.

Reproducibility is a cornerstone of scientific progress. In epigenome- and transcriptome-wide association studies (E/TWAS) failure to reproduce may be the result of false discoveries. Whereas multiple methods exist to control false discoveries due to sampling error, minimizing false discoveries due to outliers and other data artefacts remains challenging. We propose a robust E/TWAS approach that outperforms alternative methods to improve reproducibility such as split-half replication. Furthermore, robust E/TWAS results in only a minor loss of power if there are no outliers and can in the presence of outliers, likely a more realistic scenario, even be more powerful than regular E/TWAS.

Anxiety Disorders (ANX) such as panic disorder, generalized anxiety disorder, and phobias, are highly prevalent conditions that are moderately heritable. Evidence suggests that DNA methylation may play a role, as it is involved in critical adaptations to changing environments. Applying an enrichment-based sequencing approach covering nearly 28 million autosomal CpG sites, we conducted a methylome-wide association study (MWAS) of lifetime ANX in 1132 participants (618 cases/514 controls) from the Netherlands Study of Depression and Anxiety. Using epigenomic deconvolution, we performed MWAS for the main cell types in blood: granulocytes, T-cells, B-cells and monocytes. Cell-type specific analyses identified 280 and 82 methylome-wide significant associations (q-value < 0.1) in monocytes and granulocytes, respectively. Our top finding in monocytes was located in ZNF823 on chromosome 19 (p = 1.38 × 10−10) previously associated with schizophrenia. We observed significant overlap (p < 1 × 10−06) with the same direction of effect in monocytes (210 sites), T-cells (135 sites), and B-cells (727 sites) between this Discovery MWAS signal and a comparable replication dataset from the Great Smoky Mountains Study (N = 433). Overlapping Discovery-Replication MWAS signal was enriched for findings from published GWAS of ANX, major depression, and post-traumatic stress disorder. In monocytes, two specific sites in the FZR1 gene showed significant replication after Bonferroni correction with an additional 15 nominally replicated sites in monocytes and 4 in T-cells. FZR1 regulates neurogenesis in the hippocampus, and its knockout leads to impairments in associative fear memory and long-term potentiation in mice. In the largest and most extensive methylome-wide study of ANX, we identified replicable methylation sites located in genes of potential relevance for brain mechanisms of psychiatric conditions.

Postpartum depression (PPD) affects 1 in 7 women and has negative mental health consequences for both mother and child. However, the precise biological mechanisms behind the disorder are unknown. Therefore, we performed the largest transcriptome-wide association study (TWAS) for PPD (482 cases, 859 controls) to date using RNA-sequencing in whole blood and deconvoluted cell types. No transcriptional changes were observed in whole blood. B-cells showed a majority of transcriptome-wide significant results (891 transcripts representing 789 genes) with pathway analyses implicating altered B-cell activation and insulin resistance. Integration of other data types revealed cell type-specific DNA methylation loci and disease-associated eQTLs (deQTLs), but not hormones/neuropeptides (estradiol, progesterone, oxytocin, BDNF), serve as regulators for part of the transcriptional differences between cases and controls. Further, deQTLs were enriched for several brain region-specific eQTLs, but no overlap with MDD risk loci was observed. Altogether, our results constitute a convergence of evidence for pathways most affected in PPD with data across different biological mechanisms.

Brain disorders are leading causes of disability worldwide. Gene expression studies provide promising opportunities to better understand their etiology. When studying bulk tissue, cellular diversity may cause many genes that are differentially expressed in cases and controls to remain undetected. Furthermore, identifying the specific cell-types from which association signals originate is key to formulating refined hypotheses of disease etiology, designing proper follow-up experiments and, eventually, developing novel clinical interventions. Cell-type effects can be deconvoluted statistically from bulk expression data using cell-type proportions estimated with the help of a reference panel. To create a fine-grained reference panel for the human prefrontal cortex, we analyzed data from the seven largest single nucleus RNA-seq (s n RNA-seq) studies. Seventeen cell-types were robustly detected across all seven studies. To estimate the cell-type proportions, we proposed an empirical Bayes estimator that is suitable for the new panel that involves multiple low abundant cell-types. Furthermore, to avoid the use of a very large reference panel and prevent challenges with public access of nuclei level data, our estimator uses a panel comprising mean expression levels rather than the nuclei level s n RNA-seq data. Evaluations show that our empirical Bayes estimator produces highly accurate and unbiased cell-type proportion estimates. Transcriptome-wide association studies performed with permuted bulk RNA-seq data showed that it is possible to perform TWASs for even the rarest cell-types without an increased risk of false positives. Furthermore, we determined that for optimal statistical power the best approach is to analyze all cell-types in the panel as opposed to grouping or omitting (rare) cell-types.

Cross-sectional studies have identified a variety of correlates of epigenetic aging. However, it is unclear whether these make epigenetic clocks “tick faster” (i.e., accelerate biological aging). This longitudinal study examined whether a change in psychiatric problems, lifestyle variables or adversities at a particular time point was accompanied by a concurrent change in epigenetic aging or whether exposure to these correlates in childhood and adolescence lead to lagged epigenetic aging in adulthood.

The interplay of genetic factors and environmental adversities in ADHD can be reflected by cell type-specific DNA methylation (DNAm). Previous methylome-wide association studies (MWAS) for ADHD have been performed in bulk tissue only, where meaningful DNAm signal could be left undetected. We aim to identify cell type-specific DNAm profiles associated with childhood ADHD symptoms.

BACKGROUND:Longitudinal studies are needed to clarify whether early adversities are associated with advanced methylation age or if they actually accelerate methylation aging. This study test whether different dimensions of childhood adversity accelerate biological aging from childhood to adulthood, and, if so, via which mechanisms.METHODS:381 participants provided one blood sample in childhood (average age 15.0; SD = 2.3) and another in young adulthood (average age 23.1; SD = 2.8). Participants and their parents provided a median of 6 childhood assessments (total = 1,950 childhood observations), reporting exposures to different types of adversity dimensions (i.e. threat, material deprivation, loss, unpredictability). The blood samples were assayed to estimate DNA methylation age in both childhood and adulthood and also change in methylation age across this period.RESULTS:Cross-sectional associations between the childhood adversity dimensions and childhood measures of methylation age were non-significant. In contrast, multiple adversity dimensions were associated with accelerated within-person change in methylation age from adolescence to young adulthood. These associations attenuated in model testing all dimensions at the same time. Accelerated aging increased with increasing number of childhood adversities: Individuals with highest number of adversities experienced 2+ additional years of methylation aging compared to those with no exposure to childhood adversities. The association between total childhood adversity exposure and accelerated aging was partially explained by childhood depressive symptoms, but not anxiety or behavioral symptoms.CONCLUSIONS:Early adversities accelerate epigenetic aging long after they occur, in proportion to the total number of such experiences, and in a manner consistent with a shared effect that crosses multiple early dimensions of risk.

Childhood trauma is robustly linked to a broad range of adverse outcomes with consequences persisting far into adulthood. We conducted a prospective longitudinal study to predict psychiatric disorders and other adverse outcomes from trauma-related methylation changes 16.9 years after trauma exposure in childhood. Methylation was assayed using a sequencing-based approach that provides near-complete coverage of all 28 million sites in the blood methylome. Methylation data involved 673 assays from 489 participants aged 13.6 years (SD = 1.9) with outcomes measures collected at age 30.4 (SD = 2.26). For a subset of 303 participants we also generated methylation data in adulthood. Trauma-related methylation risk scores (MRSs) significantly predicted adult depression, externalizing problems, nicotine dependence, alcohol use disorder, serious medical problems, social problems and poverty. The predictive power of the MRSs was higher than that of reported trauma and could not be explained by the reported trauma, correlations with demographic variables, or a continuity of the predicted health problems from childhood to adulthood. Rather than measuring the occurrence of traumatic events, the MRSs seemed to capture the subject-specific impact of trauma. The majority of predictive sites did not remain associated with the outcomes suggesting the signatures of trauma do not become biologically embedded in the blood methylome. Instead, the long-term effects of trauma therefore seemed more consistent with a developmental mechanism where the initial subject-specific impacts of trauma are magnified over time. The MRSs have the potential to be a novel clinical biomarker for the assessment of trauma-related health risks.