Neuronal-specific methylome and hydroxymethylome analysis reveal significant loci associated with alcohol use disorder.

Background: Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5 mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5 hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5 mC and 5 hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Methods: Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5 mC and 5 hmC at the genome-wide level. Differential 5 mC and 5 hmC were evaluated using the methylKit R package and significance was set at false discovery rate < 0.05 and differential methylation > 2. Functional enrichment analyses were performed, and gene-level convergence was evaluated in an independent dataset that assessed 5 mC and 5 hmC of AUD in bulk cortical tissue. Results: We identified 417 5 mC and 363 5hmC significant differential CpG sites associated with AUD, with 59% in gene promoters. Some of the identified genes have been previously implicated in alcohol consumption, including SYK, DNMT3A for 5 mC, GAD1, DLX1, DLX2, for 5 hmC and GATA4 in both. Convergence with a previous AUD 5 mC and 5 hmC study was observed for 28 genes. We also identified 5 and 35 differential regions for 5 mC and 5 hmC, respectively. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5 mC genes. Discussion: This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD, identifying both previously reported and potentially novel gene associations with AUD. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

Associations of depression and anxiety symptoms in childhood and adolescence with epigenetic aging.

BACKGROUND: Childhood anxiety and depression symptoms are potential risk factors for accelerated biological aging. In child and adolescent twins, we tested whether these symptoms were associated with DNA methylation (DNAm) aging, a measure of biological aging. METHODS: 276 twins (135 pairs, 6 singletons) had DNAm assayed from saliva in middle childhood (mean = 7.8 years). Residuals of five different DNAm age estimates regressed on chronological age were used to indicate accelerated aging. Anxiety and depression symptoms were assessed in middle childhood and early adolescence using the Child Behavior Checklist. Mixed effect regression was used to examine potential relationships between anxiety or depression symptoms, and accelerated DNAm age. MZ twin difference analysis was then utilized to determine if associations were environmentally-driven or due to genetic or shared-environment confounding. RESULTS: Anxiety and depression symptoms were not associated with accelerated DNAm aging in middle childhood. In early adolescence, only the Wu clock was significant and indicated that each one symptom increase in anxiety symptoms had an associated age acceleration of 0.03 years (~0.4 months; p = 0.019). MZ twin difference analysis revealed non-significant within-pair effects, suggesting genetic and shared environmental influences. LIMITATIONS: Sample is predominantly male and white. Generalizability to other populations may be limited. CONCLUSION: Accelerated DNAm aging of the Wu clock in middle childhood is associated with anxiety, but not depression, symptoms in early adolescence. Further, this association may be the result of shared genetic and environmental influences. Accelerated DNAm aging may serve as an early risk factor or predictor of later anxiety symptoms.

Alcohol Use Disorder-Associated DNA Methylation in the Nucleus Accumbens and Dorsolateral Prefrontal Cortex.

Background: Alcohol use disorder (AUD) has a profound public health impact. However, understanding of the molecular mechanisms underlying the development and progression of AUD remain limited. Here, we interrogate AUD-associated DNA methylation (DNAm) changes within and across addiction-relevant brain regions: the nucleus accumbens (NAc) and dorsolateral prefrontal cortex (DLPFC). Methods: Illumina HumanMethylation EPIC array data from 119 decedents of European ancestry (61 cases, 58 controls) were analyzed using robust linear regression, with adjustment for technical and biological variables. Associations were characterized using integrative analyses of public gene regulatory data and published genetic and epigenetic studies. We additionally tested for brain region-shared and -specific associations using mixed effects modeling and assessed implications of these results using public gene expression data. Results: At a false discovery rate ≤ 0.05, we identified 53 CpGs significantly associated with AUD status for NAc and 31 CpGs for DLPFC. In a meta-analysis across the regions, we identified an additional 21 CpGs associated with AUD, for a total of 105 unique AUD-associated CpGs (120 genes). AUD-associated CpGs were enriched in histone marks that tag active promoters and our strongest signals were specific to a single brain region. Of the 120 genes, 23 overlapped with previous genetic associations for substance use behaviors; all others represent novel associations. Conclusions: Our findings identify AUD-associated methylation signals, the majority of which are specific within NAc or DLPFC. Some signals may constitute predisposing genetic and epigenetic variation, though more work is needed to further disentangle the neurobiological gene regulatory differences associated with AUD.

Neuronal-specific methylome and hydroxymethylome analysis reveal replicated and novel loci associated with alcohol use disorder.

Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5mC and 5hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5mC and 5hmC at the genome-wide level. Differential 5mC and 5hmC were evaluated using the methylKit R package and significance was set at false discovery rate <0.05 and differential methylation >2. Functional enrichment analyses were performed and replication was evaluated replication in an independent dataset that assessed 5mC and 5hmC of AUD in bulk cortical tissue. We identified 417 5mC and 363 5hmC genome-wide significant differential CpG sites associated with AUD, with 59% in gene promoters. We also identified genes previously implicated in alcohol consumption, such as SYK, CHRM2, DNMT3A, and GATA4, for 5mC and GATA4, and GAD1, GATA4, DLX1 for 5hmC. Replication was observed for 28 CpG sites from a previous AUD 5mC and 5hmC study, including FOXP1. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5mC genes. This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD. We replicated previous findings and identified novel associations with AUD for both 5mC and 5hmC marks within the OFC. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

The Impact of Childhood Mental Health and Substance Use on Methylation Aging Into Adulthood.

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.

Characterizing patterns of substance use in trauma exposed youth.

OBJECTIVE: Previous work investigating the impact of childhood trauma on substance use and co-occurring psychiatric disorders has primarily been conducted in adults or on specific trauma types. This limits understanding of traumas impact in childhood and how different types of traumas play a role. We sought to characterize substance use in a sample of trauma-exposed youth in the context of psychiatric comorbidities. METHOD: 1152 youth from the Texas Childhood Trauma Research Network (TX-CTRN) that were exposed to at least one trauma meeting DSM-5 Criterion A were assessed for current substance use and psychiatric diagnoses. Latent class analysis was used to identify patterns of substance use. To characterize these patterns, we examined if demographics, number of trauma types experienced, or childhood psychiatric disorders predicted class membership. RESULTS: We identified four primary patterns of substance use: Non-use (66.1%), predominantly alcohol use (19.7%), predominantly cannabis use (4.5%), and polysubstance use (9.7%). Compared to the non-users, polysubstance users tended to be older, Non-Hispanic White, have experienced more types of trauma. They were also more likely to have fulfilled diagnostic criteria for suicidality and ADHD. Comparisons among the substance using classes were more nuanced. CONCLUSION: The findings highlight the need for universal assessments of trauma, substance misuse, and mental health symptoms in youth as the presence or absence of their co-occurrence has implications for treatment.

Methylome-wide association study of anxiety disorders.

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.

A single-nucleus transcriptomics study of alcohol use disorder in the nucleus accumbens.

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.

Longitudinal patterns of posttraumatic stress disorder symptoms among sexual assault survivors: A latent transition analysis.

OBJECTIVE: The aim of this study was to identify latent classes of posttraumatic stress disorder (PTSD) symptoms in a community sample of sexual assault survivors whose assaults occurred varying lengths of time in the past and to explore patterns of transition between those latent classes over time. METHOD: Latent class analysis was used to identify naturally occurring subgroups of PTSD symptoms in a sample of sexual assault survivors who completed two mailed surveys 1 year apart (N = 1,271). Latent transition analysis was then used to examine individuals' probabilities of transitioning into each latent class at Time 2 based on their latent class membership at Time 1. RESULTS: A four-class model emerged as the best fitting model at both Time 1 and Time 2. Classes demonstrated overall severity and symptom cluster severity differences. Transition into a lower severity class was more common than transition into a higher severity class, though escalation was demonstrated by 6-20% of participants in each latent class. CONCLUSIONS: The substantial heterogeneity in sexual assault survivors' PTSD symptoms highlights the variety of ways that posttraumatic stress may be experienced years after a sexual assault. Future research should explore factors that affect long-term symptoms, including cumulative lifetime trauma and social support. (PsycInfo Database Record (c) 2022 APA, all rights reserved).

Transcriptome-wide association study for postpartum depression implicates altered B-cell activation and insulin resistance.

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.

Dual methylation and hydroxymethylation study of alcohol use disorder.

Using an integrative, multi-tissue design, we sought to characterize methylation and hydroxymethylation changes in blood and brain associated with alcohol use disorder (AUD). First, we used epigenomic deconvolution to perform cell-type-specific methylome-wide association studies within subpopulations of granulocytes/T-cells/B-cells/monocytes in 1132 blood samples. Blood findings were then examined for overlap with AUD-related associations with methylation and hydroxymethylation in 50 human post-mortem brain samples. Follow-up analyses investigated if overlapping findings mediated AUD-associated transcription changes in the same brain samples. Lastly, we replicated our blood findings in an independent sample of 412 individuals and aimed to replicate published alcohol methylation findings using our results. Cell-type-specific analyses in blood identified methylome-wide significant associations in monocytes and T-cells. The monocyte findings were significantly enriched for AUD-related methylation and hydroxymethylation in brain. Hydroxymethylation in specific sites mediated AUD-associated transcription in the same brain samples. As part of the most comprehensive methylation study of AUD to date, this work involved the first cell-type-specific methylation study of AUD conducted in blood, identifying and replicating a finding in DLGAP1 that may be a blood-based biomarker of AUD. In this first study to consider the role of hydroxymethylation in AUD, we found evidence for a novel mechanism for cognitive deficits associated with AUD. Our results suggest promising new avenues for AUD research.

Clarifying the Genetic Influences on Nicotine Dependence and Quantity of Use in Cigarette Smokers.

Nicotine dependence and smoking quantity are both robustly associated with the CHRNA5-A3-B4 gene cluster in the 15q25 region, and SNP rs16969968 in particular. The purpose of this paper is to use structural equation modeling techniques (SEM) to disentangle the complex pattern of relationships between rs16969968, nicotine quantity (as measured by the number of cigarettes an individual smokes per day; CPD) and nicotine dependence (as measured by the Fagerström Test for Nicotine Dependence; FTND). CPD is an indicator, but also a potential cause, of FTND, complicating the interpretation of associations between these constructs and requires a more detailed investigation than standard GWAS or general linear regression models can provide. FTND items and genotypes were collected in four samples, with a combined sample size of 5,373 respondents. A mega-analysis was conducted using a multiple group SEM approach to test competing hypotheses regarding the relationships between the SNP rs16969968, FTND and CPD. In the best fitting model, the FTND items loaded onto two correlated factors. The first, labeled "maintenance," assesses the motivation to maintain constant levels of nicotine through out the day. The second was labeled "urgency" as its items concern the urgency to restore nicotine levels after abstinence. We focus our attention on the "maintenance" factor, of which CPD was an indicator. The best fitting model included a negative feedback loop between the Maintenance factor and CPD. Accordingly, the motivation to maintain higher levels of nicotine increased the quantity of nicotine consumed, which subsequently decreases the maintenance motivation. The fact that the Maintenance-CPD feedback model fits the data best implies that there are at least two biological pathways that lead from rs16969968 to smoking behaviors. The model is consistent with a supply and demand system, which allows individuals to achieve a homeostatic equilibrium for their nicotine concentration.

Methylomic Investigation of Problematic Adolescent Cannabis Use and Its Negative Mental Health Consequences.

OBJECTIVE: The impact of adolescent cannabis use is a pressing public health question owing to the high rates of use and links to negative outcomes. This study considered the association between problematic adolescent cannabis use and methylation. METHOD: Using an enrichment-based sequencing approach, a methylome-wide association study (MWAS) was performed of problematic adolescent cannabis use in 703 adolescent samples from the Great Smoky Mountain Study. Using epigenomic deconvolution, MWASs were performed for the main cell types in blood: granulocytes, T cells, B cells, and monocytes. Enrichment testing was conducted to establish overlap between cannabis-associated methylation differences and variants associated with negative mental health effects of adolescent cannabis use. RESULTS: Whole-blood analyses identified 45 significant CpGs, and cell type-specific analyses yielded 32 additional CpGs not identified in the whole-blood MWAS. Significant overlap was observed between the B-cell MWAS and genetic studies of education attainment and intelligence. Furthermore, the results from both T cells and monocytes overlapped with findings from an MWAS of psychosis conducted in brain tissue. CONCLUSION: In one of the first methylome-wide association studies of adolescent cannabis use, several methylation sites located in genes of importance for potentially relevant brain functions were identified. These findings resulted in several testable hypotheses by which cannabis-associated methylation can impact neurological development and inflammation response as well as potential mechanisms linking cannabis-associated methylation to potential downstream mental health effects.

GW-SEM 2.0: Efficient, Flexible, and Accessible Multivariate GWAS.

Most genome-wide association study (GWAS) analyses test the association between single-nucleotide polymorphisms (SNPs) and a single trait or outcome. While valuable second-step analyses of these associations (e.g., calculating genetic correlations between traits) are common, single-step multivariate analyses of GWAS data are rarely performed. This is unfortunate because multivariate analyses can reveal information which is irrevocably obscured in multi-step analysis. One simple example is the distinction between variance common to a set of measures, and variance specific to each. Neither GWAS of sum- or factor-scores, nor GWAS of the individual measures will deliver a clean picture of loci associated with each measure's specific variance. While multivariate GWAS opens up a broad new landscape of feasible and informative analyses, its adoption has been slow, likely due to the heavy computational demands and difficulties specifying models it requires. Here we describe GW-SEM 2.0, which is designed to simplify model specification and overcome the inherent computational challenges associated with multivariate GWAS. In addition, GW-SEM 2.0 allows users to accurately model ordinal items, which are common in behavioral and psychological research, within a GWAS context. This new release enhances computational efficiency, allows users to select the fit function that is appropriate for their analyses, expands compatibility with standard genomic data formats, and outputs results for seamless reading into other standard post-GWAS processing software. To demonstrate GW-SEM's utility, we conducted (1) a series of GWAS using three substance use frequency items from data in the UK Biobank, (2) a timing study for several predefined GWAS functions, and (3) a Type I Error rate study. Our multivariate GWAS analyses emphasize the utility of GW-SEM for identifying novel patterns of associations that vary considerably between genomic loci for specific substances, highlighting the importance of differentiating between substance-specific use behaviors and polysubstance use. The timing studies demonstrate that the analyses take a reasonable amount of time and show the cost of including additional items. The Type I Error rate study demonstrates that hypothesis tests for genetic associations with latent variable models follow the hypothesized uniform distribution. Taken together, we suggest that GW-SEM may provide substantially deeper insights into the underlying genomic architecture for multivariate behavioral and psychological systems than is currently possible with standard GWAS methods. The current release of GW-SEM 2.0 is available on CRAN (stable release) and GitHub (beta release), and tutorials are available on our github wiki ( https://jpritikin.github.io/gwsem/ ).

Methylome-wide association findings for major depressive disorder overlap in blood and brain and replicate in independent brain samples.

We present the first large-scale methylome-wide association studies (MWAS) for major depressive disorder (MDD) to identify sites of potential importance for MDD etiology. Using a sequencing-based approach that provides near-complete coverage of all 28 million common CpGs in the human genome, we assay methylation in MDD cases and controls from both blood (N = 1132) and postmortem brain tissues (N = 61 samples from Brodmann Area 10, BA10). The MWAS for blood identified several loci with P ranging from 1.91 × 10-8 to 4.39 × 10-8 and a resampling approach showed that the cumulative association was significant (P = 4.03 × 10-10) with the signal coming from the top 25,000 MWAS markers. Furthermore, a permutation-based analysis showed significant overlap (P = 5.4 × 10-3) between the MWAS findings in blood and brain (BA10). This overlap was significantly enriched for a number of features including being in eQTLs in blood and the frontal cortex, CpG islands and shores, and exons. The overlapping sites were also enriched for active chromatin states in brain including genic enhancers and active transcription start sites. Furthermore, three loci located in GABBR2, RUFY3, and in an intergenic region on chromosome 2 replicated with the same direction of effect in the second brain tissue (BA25, N = 60) from the same individuals and in two independent brain collections (BA10, N = 81 and 64). GABBR2 inhibits neuronal activity through G protein-coupled second-messenger systems and RUFY3 is implicated in the establishment of neuronal polarity and axon elongation. In conclusion, we identified and replicated methylated loci associated with MDD that are involved in biological functions of likely importance to MDD etiology.

A methylation study of long-term depression risk.

Recurrent and chronic major depressive disorder (MDD) accounts for a substantial part of the disease burden because this course is most prevalent and typically requires long-term treatment. We associated blood DNA methylation profiles from 581 MDD patients at baseline with MDD status 6 years later. A resampling approach showed a highly significant association between methylation profiles in blood at baseline and future disease status (P = 2.0 × 10-16). Top MWAS results were enriched specific pathways, overlapped with genes found in GWAS of MDD disease status, autoimmune disease and inflammation, and co-localized with eQTLS and (genic enhancers of) of transcription sites in brain and blood. Many of these findings remained significant after correction for multiple testing. The major themes emerging were cellular responses to stress and signaling mechanisms linked to immune cell migration and inflammation. This suggests that an immune signature of treatment-resistant depression is already present at baseline. We also created a methylation risk score (MRS) to predict MDD status 6 years later. The AUC of our MRS was 0.724 and higher than risk scores created using a set of five putative MDD biomarkers, genome-wide SNP data, and 27 clinical, demographic and lifestyle variables. Although further studies are needed to examine the generalizability to different patient populations, these results suggest that methylation profiles in blood may present a promising avenue to support clinical decision making by providing empirical information about the likelihood MDD is chronic or will recur in the future.

A Whole Methylome Study of Ethanol Exposure in Brain and Blood: An Exploration of the Utility of Peripheral Blood as Proxy Tissue for Brain in Alcohol Methylation Studies.

BACKGROUND: Recent reviews have highlighted the potential use of blood-based methylation biomarkers as diagnostic and prognostic tools of current and future alcohol use and addiction. Due to the substantial overlap that often exists between methylation patterns across different tissues, including blood and brain, blood-based methylation may track methylation changes in brain; however, little work has explored the overlap in alcohol-related methylation in these tissues. METHODS: To study the effects of alcohol on the brain methylome and identify possible biomarkers of these changes in blood, we performed a methylome-wide association study in brain and blood from 40 male DBA/2J mice that received either an acute ethanol (EtOH) or saline intraperitoneal injection. To investigate all 22 million CpGs in the mouse genome, we enriched for the methylated genomic fraction using methyl-CpG binding domain (MBD) protein capture followed by next-generation sequencing (MBD-seq). We performed association tests in blood and brain separately followed by enrichment testing to determine whether there was overlapping alcohol-related methylation in the 2 tissues. RESULTS: The top result for brain was a CpG located in an intron of Ttc39b (p = 5.65 × 10-08 ), and for blood, the top result was located in Espnl (p = 5.11 × 10-08 ). Analyses implicated pathways involved in inflammation and neuronal differentiation, such as CXCR4, IL-7, and Wnt signaling. Enrichment tests indicated significant overlap among the top results in brain and blood. Pathway analyses of the overlapping genes converge on MAPKinase signaling (p = 5.6 × 10-05 ) which plays a central role in acute and chronic responses to alcohol and glutamate receptor pathways, which can regulate neuroplastic changes underlying addictive behavior. CONCLUSIONS: Overall, we have shown some methylation changes in brain and blood after acute EtOH administration and that the changes in blood partly mirror the changes in brain suggesting the potential for DNA methylation in blood to be biomarkers of alcohol use.

Convergence of evidence from a methylome-wide CpG-SNP association study and GWAS of major depressive disorder.

DNA methylation is an epigenetic modification that provides stability and diversity to the cellular phenotype. It is influenced by both genetic sequence variation and environmental factors, and can therefore potentially account for variation of heritable phenotypes and disorders. Therefore, methylome-wide association studies (MWAS) are promising complements to genome-wide association studies (GWAS) of genetic variants. Of particular interest are methylation sites (CpGs) that are created or destroyed by the alleles of single-nucleotide polymorphisms (SNPs), as these so-called CpG-SNPs may show variation in methylation levels on top of what can be explained by the sequence variation. Using sequencing-based data from 1132 major depressive disorder (MDD) cases and controls, we performed a MWAS of 970,414 common CpG-SNPs. The analysis identified 27 suggestively significant (P < 1.00 × 10-5) CpG-SNPs associations. Furthermore, the MWAS results were over-represented (odds ratios ranging 1.36-5.00; P ranging 4.9 × 10-3-8.1 × 10-2) among findings from three recent GWAS for MDD-related phenotypes. Overlapping loci included, e.g., ROBO2, ASIC2, and DCC. As the CpG-SNP analysis accounts for the number of alleles that creates CpGs, the methylation differences could not be explained by differences in allele frequencies. Thus, the results show that the MWAS and GWASs provide independent lines of evidence for the involvement of these loci in MDD. In conclusion, our methylation study of MDD contributes novel information about loci of relevance that complements previous findings and generates new hypothesis about MDD etiology, such as that the functional effects of genetic association may be partly mediated and/or enhanced by the methylation status in these loci.

Epigenetic Aging in Major Depressive Disorder.

OBJECTIVE: Major depressive disorder is associated with an increased risk of mortality and aging-related diseases. The authors examined whether major depression is associated with higher epigenetic aging in blood as measured by DNA methylation (DNAm) patterns, whether clinical characteristics of major depression have a further impact on these patterns, and whether the findings replicate in brain tissue. METHOD: DNAm age was estimated using all methylation sites in blood of 811 depressed patients and 319 control subjects with no lifetime psychiatric disorders and low depressive symptoms from the Netherlands Study of Depression and Anxiety. The residuals of the DNAm age estimates regressed on chronological age were calculated to indicate epigenetic aging. Major depression diagnosis and clinical characteristics were assessed with questionnaires and psychiatric interviews. Analyses were adjusted for sociodemographic characteristics, lifestyle, and health status. Postmortem brain samples of 74 depressed patients and 64 control subjects were used for replication. Pathway enrichment analysis was conducted using ConsensusPathDB to gain insight into the biological processes underlying epigenetic aging in blood and brain. RESULTS: Significantly higher epigenetic aging was observed in patients with major depression compared with control subjects (Cohen's d=0.18), with a significant dose effect with increasing symptom severity in the overall sample. In the depression group, epigenetic aging was positively and significantly associated with childhood trauma score. The case-control difference was replicated in an independent data set of postmortem brain samples. The top significantly enriched Gene Ontology terms included neuronal processes. CONCLUSIONS: As compared with control subjects, patients with major depression exhibited higher epigenetic aging in blood and brain tissue, suggesting that they are biologically older than their corresponding chronological age. This effect was even more profound in the presence of childhood trauma.

RaMWAS: fast methylome-wide association study pipeline for enrichment platforms.

Motivation: Enrichment-based technologies can provide measurements of DNA methylation at tens of millions of CpGs for thousands of samples. Existing tools for methylome-wide association studies cannot analyze datasets of this size and lack important features like principal component analysis, combined analysis with SNP data and outcome predictions that are based on all informative methylation sites. Results: We present a Bioconductor R package called RaMWAS with a full set of tools for large-scale methylome-wide association studies. It is free, cross-platform, open source, memory efficient and fast. Availability and implementation: Release version and vignettes with small case study at bioconductor.org/packages/ramwas Development version at github.com/andreyshabalin/ramwas. Supplementary information: Supplementary data are available at Bioinformatics online.