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.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value.

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Trans-ancestry epigenome-wide association meta-analysis of DNA methylation with lifetime cannabis use.

Cannabis is widely used worldwide, yet its links to health outcomes are not fully understood. DNA methylation can serve as a mediator to link environmental exposures to health outcomes. We conducted an epigenome-wide association study (EWAS) of peripheral blood-based DNA methylation and lifetime cannabis use (ever vs. never) in a meta-analysis including 9436 participants (7795 European and 1641 African ancestry) from seven cohorts. Accounting for effects of cigarette smoking, our trans-ancestry EWAS meta-analysis revealed four CpG sites significantly associated with lifetime cannabis use at a false discovery rate of 0.05 [Formula: see text]: cg22572071 near gene ADGRF1, cg15280358 in ADAM12, cg00813162 in ACTN1, and cg01101459 near LINC01132. Additionally, our EWAS analysis in participants who never smoked cigarettes identified another epigenome-wide significant CpG site, cg14237301 annotated to APOBR. We used a leave-one-out approach to evaluate methylation scores constructed as a weighted sum of the significant CpGs. The best model can explain 3.79% of the variance in lifetime cannabis use. These findings unravel the DNA methylation changes associated with lifetime cannabis use that are independent of cigarette smoking and may serve as a starting point for further research on the mechanisms through which cannabis exposure impacts health outcomes.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWAS) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N=52) and nonsmokers (N=171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and using a two-stage multiple testing approach with eigenMT and Bonferroni corrections. We found >2 million significant meQTL variants (padj<0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects; five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTLs for 958 unique eGenes (padj<0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTIN1 and ITIH4 colocalized across all data types (GWAS + meQTL + eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Cis-meQTL for cocaine use-associated DNA methylation in an HIV-positive cohort show pleiotropic effects on multiple traits.

BACKGROUND: Cocaine use (CU) is associated with psychiatric and medical diseases. Little is known about the mechanisms of CU-related comorbidities. Findings from preclinical and clinical studies have suggested that CU is associated with aberrant DNA methylation (DNAm) that may be influenced by genetic variants [i.e., methylation quantitative trait loci (meQTLs)]. In this study, we mapped cis-meQTLs for CU-associated DNAm sites (CpGs) in an HIV-positive cohort (Ntotal = 811) and extended the meQTLs to multiple traits. RESULTS: We conducted cis-meQTL analysis for 224 candidate CpGs selected for their association with CU in blood. We identified 7,101 significant meQTLs [false discovery rate (FDR) < 0.05], which mostly mapped to genes involved in immunological functions and were enriched in immune pathways. We followed up the meQTLs using phenome-wide association study and trait enrichment analyses, which revealed 9 significant traits. We tested for causal effects of CU on these 9 traits using Mendelian Randomization and found evidence that CU plays a causal role in increasing hypertension (p-value = 2.35E-08) and decreasing heel bone mineral density (p-value = 1.92E-19). CONCLUSIONS: These findings suggest that genetic variants for CU-associated DNAm have pleiotropic effects on other relevant traits and provide new insights into the causal relationships between cocaine use and these complex traits.

Africa-specific human genetic variation near CHD1L associates with HIV-1 load.

HIV-1 remains a global health crisis1, highlighting the need to identify new targets for therapies. Here, given the disproportionate HIV-1 burden and marked human genome diversity in Africa2, we assessed the genetic determinants of control of set-point viral load in 3,879 people of African ancestries living with HIV-1 participating in the international collaboration for the genomics of HIV3. We identify a previously undescribed association signal on chromosome 1 where the peak variant associates with an approximately 0.3 log10-transformed copies per ml lower set-point viral load per minor allele copy and is specific to populations of African descent. The top associated variant is intergenic and lies between a long intergenic non-coding RNA (LINC00624) and the coding gene CHD1L, which encodes a helicase that is involved in DNA repair4. Infection assays in iPS cell-derived macrophages and other immortalized cell lines showed increased HIV-1 replication in CHD1L-knockdown and CHD1L-knockout cells. We provide evidence from population genetic studies that Africa-specific genetic variation near CHD1L associates with HIV replication in vivo. Although experimental studies suggest that CHD1L is able to limit HIV infection in some cell types in vitro, further investigation is required to understand the mechanisms underlying our observations, including any potential indirect effects of CHD1L on HIV spread in vivo that our cell-based assays cannot recapitulate.

Long-term air pollution exposure and markers of cardiometabolic health in the National Longitudinal Study of Adolescent to Adult Health (Add Health).

BACKGROUND: Air pollution exposure is associated with cardiovascular morbidity and mortality. Although exposure to air pollution early in life may represent a critical window for development of cardiovascular disease risk factors, few studies have examined associations of long-term air pollution exposure with markers of cardiovascular and metabolic health in young adults. OBJECTIVES: By combining health data from the National Longitudinal Study of Adolescent to Adult Health (Add Health) with air pollution data from the Fused Air Quality Surface using Downscaling (FAQSD) archive, we: (1) calculated multi-year estimates of exposure to ozone (O3) and particulate matter with an aerodynamic diameter ≤ 2.5 µm (PM2.5) for Add Health participants; and (2) estimated associations between air pollution exposures and multiple markers of cardiometabolic health. METHODS: Add Health is a nationally representative longitudinal cohort study of over 20,000 adolescents aged 12-19 in the United States (US) in 1994-95 (Wave I). Participants have been followed through adolescence and into adulthood with five in-home interviews. Estimated daily concentrations of O3 and PM2.5 at census tracts were obtained from the FAQSD archive and used to generate tract-level annual averages of O3 and PM2.5 concentrations. We estimated associations between average O3 and PM2.5 exposures from 2002 to 2007 and markers of cardiometabolic health measured at Wave IV (2008-09), including hypertension, hyperlipidemia, body mass index (BMI), diabetes, C-reactive protein, and metabolic syndrome. RESULTS: The final sample size was 11,259 individual participants. The average age of participants at Wave IV was 28.4 years (range: 24-34 years). In models adjusting for age, race/ethnicity, and sex, long-term O3 exposure (2002-07) was associated with elevated odds of hypertension, with an odds ratio (OR) of 1.015 (95% confidence interval [CI]: 1.011, 1.029); obesity (1.022 [1.004, 1.040]); diabetes (1.032 [1.009,1.054]); and metabolic syndrome (1.028 [1.014, 1.041]); PM2.5 exposure (2002-07) was associated with elevated odds of hypertension (1.022 [1.001, 1.045]). CONCLUSION: Findings suggest that long-term ambient air pollution exposure, particularly O3 exposure, is associated with cardiometabolic health in early adulthood.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing.

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Polygenic risk scores and the need for pharmacotherapy in neonatal abstinence syndrome.

BACKGROUND: The aim of this study was to identify genetic variants associated with NAS through a genome-wide association study (GWAS) and estimate a Polygenic Risk Score (PRS) model for NAS. METHODS: A prospective case-control study included 476 in utero opioid-exposed term neonates. A GWAS of 1000 genomes-imputed genotypes was performed to identify variants associated with need for pharmacotherapy for NAS. PRS models for estimating genetic predisposition were generated via a nested cross-validation approach using 382 neonates of European ancestry. PRS predictive ability, discrimination, and calibration were assessed. RESULTS: Cross-ancestry GWAS identified one intergenic locus on chromosome 7 downstream of SNX13 exhibiting genome-wide association with need for pharmacotherapy. PRS models derived from the GWAS for a subset of the European ancestry neonates reliably discriminated between need for pharmacotherapy using cis variant effect sizes within validation sets of European and African American ancestry neonates. PRS were less effective when applying variant effect sizes across datasets and in calibration analyses. CONCLUSIONS: GWAS has the potential to identify genetic loci associated with need for pharmacotherapy for NAS and enable development of clinically predictive PRS models. Larger GWAS with additional ancestries are needed to confirm the observed SNX13 association and the accuracy of PRS in NAS risk prediction models. IMPACT: Genetic associations appear to be important in neonatal abstinence syndrome. This is the first genome-wide association in neonates with neonatal abstinence syndrome. Polygenic risk scores can be developed examining single-nucleotide polymorphisms across the entire genome. Polygenic risk scores were higher in neonates receiving pharmacotherapy for treatment of their neonatal abstinence syndrome. Future studies with larger cohorts are needed to better delineate these genetic associations.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing.

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.

Chromatin architecture in addiction circuitry identifies risk genes and potential biological mechanisms underlying cigarette smoking and alcohol use traits.

Cigarette smoking and alcohol use are among the most prevalent substances used worldwide and account for a substantial proportion of preventable morbidity and mortality, underscoring the public health significance of understanding their etiology. Genome-wide association studies (GWAS) have successfully identified genetic variants associated with cigarette smoking and alcohol use traits. However, the vast majority of risk variants reside in non-coding regions of the genome, and their target genes and neurobiological mechanisms are unknown. Chromosomal conformation mappings can address this knowledge gap by charting the interaction profiles of risk-associated regulatory variants with target genes. To investigate the functional impact of common variants associated with cigarette smoking and alcohol use traits, we applied Hi-C coupled MAGMA (H-MAGMA) built upon cortical and newly generated midbrain dopaminergic neuronal Hi-C datasets to GWAS summary statistics of nicotine dependence, cigarettes per day, problematic alcohol use, and drinks per week. The identified risk genes mapped to key pathways associated with cigarette smoking and alcohol use traits, including drug metabolic processes and neuronal apoptosis. Risk genes were highly expressed in cortical glutamatergic, midbrain dopaminergic, GABAergic, and serotonergic neurons, suggesting them as relevant cell types in understanding the mechanisms by which genetic risk factors influence cigarette smoking and alcohol use. Lastly, we identified pleiotropic genes between cigarette smoking and alcohol use traits under the assumption that they may reveal substance-agnostic, shared neurobiological mechanisms of addiction. The number of pleiotropic genes was ~26-fold higher in dopaminergic neurons than in cortical neurons, emphasizing the critical role of ascending dopaminergic pathways in mediating general addiction phenotypes. Collectively, brain region- and neuronal subtype-specific 3D genome architecture helps refine neurobiological hypotheses for smoking, alcohol, and general addiction phenotypes by linking genetic risk factors to their target genes.

Convergence of case-specific epigenetic alterations identify a confluence of genetic vulnerabilities tied to opioid overdose.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQß1.

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQß1Leu26 (p valueMeta = 1.24 × 10-14) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQß1Pro55 (p valueMeta = 8.23 × 10-11) located in the peptide binding region was positively associated, independently of HLA-DQß1Leu26. These two amino acids are not in linkage disequilibrium (r2 < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1∗03:01, and HLA-DRB1∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQß1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC50 nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10-38). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQß1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.

Studying the Utility of Using Genetics to Predict Smoking-Related Outcomes in a Population-Based Study and a Selected Cohort.

INTRODUCTION: The purpose of this study is to examine the predictive utility of polygenic risk scores (PRSs) for smoking behaviors. AIMS AND METHODS: Using summary statistics from the Sequencing Consortium of Alcohol and Nicotine use consortium, we generated PRSs of ever smoking, age of smoking initiation, cigarettes smoked per day, and smoking cessation for participants in the population-based Atherosclerosis Risk in Communities (ARIC) study (N = 8638), and the Collaborative Genetic Study of Nicotine Dependence (COGEND) (N = 1935). The outcomes were ever smoking, age of smoking initiation, heaviness of smoking, and smoking cessation. RESULTS: In the European ancestry cohorts, each PRS was significantly associated with the corresponding smoking behavior outcome. In the ARIC cohort, the PRS z-score for ever smoking predicted smoking (odds ratio [OR]: 1.37; 95% confidence interval [CI]: 1.31, 1.43); the PRS z-score for age of smoking initiation was associated with age of smoking initiation (OR: 0.87; 95% CI: 0.82, 0.92); the PRS z-score for cigarettes per day was associated with heavier smoking (OR: 1.17; 95% CI: 1.11, 1.25); and the PRS z-score for smoking cessation predicted successful cessation (OR: 1.24; 95% CI: 1.17, 1.32). In the African ancestry cohort, the PRSs did not predict smoking behaviors. CONCLUSIONS: Smoking-related PRSs were associated with smoking-related behaviors in European ancestry populations. This improvement in prediction is greatest in the lowest and highest genetic risk categories. The lack of prediction in African ancestry populations highlights the urgent need to increase diversity in research so that scientific advances can be applied to populations other than those of European ancestry. IMPLICATIONS: This study shows that including both genetic ancestry and PRSs in a single model increases the ability to predict smoking behaviors compared with the model including only demographic characteristics. This finding is observed for every smoking-related outcome. Even though adding genetics is more predictive, the demographics alone confer substantial and meaningful predictive power. However, with increasing work in PRSs, the predictive ability will continue to improve.

Integration of evidence across human and model organism studies: A meeting report.

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Alcohol and cigarette smoking consumption as genetic proxies for alcohol misuse and nicotine dependence.

PURPOSE: To investigate the role of consumption phenotypes as genetic proxies for alcohol misuse and nicotine dependence. METHODS: We leveraged GWAS data from well-powered studies of consumption, alcohol misuse, and nicotine dependence phenotypes measured in individuals of European ancestry from the UK Biobank (UKB) and other population-based cohorts (largest total N = 263,954), and performed genetic correlations within a medical-center cohort, BioVU (N = 66,915). For alcohol, we used quantitative measures of consumption and misuse via AUDIT from UKB. For smoking, we used cigarettes per day from UKB and non-UKB cohorts comprising the GSCAN consortium, and nicotine dependence via ICD codes from UKB and Fagerström Test for Nicotine Dependence from non-UKB cohorts. RESULTS: In a large phenome-wide association study, we show that smoking consumption and dependence phenotypes show similar strongly negatively associations with a plethora of diseases, whereas alcohol consumption shows patterns of genetic association that diverge from those of alcohol misuse. CONCLUSIONS: Our study suggests that cigarette smoking consumption, which can be easily measured in the general population, may be good a genetic proxy for nicotine dependence, whereas alcohol consumption is not a direct genetic proxy of alcohol misuse.

Epigenome-wide analysis uncovers a blood-based DNA methylation biomarker of lifetime cannabis use.

Cannabis use is highly prevalent and is associated with adverse and beneficial effects. To better understand the full spectrum of health consequences, biomarkers that accurately classify cannabis use are needed. DNA methylation (DNAm) is an excellent candidate, yet no blood-based epigenome-wide association studies (EWAS) in humans exist. We conducted an EWAS of lifetime cannabis use (ever vs. never) using blood-based DNAm data from a case-cohort study within Sister Study, a prospective cohort of women at risk of developing breast cancer (Discovery N = 1,730 [855 ever users]; Replication N = 853 [392 ever users]). We identified and replicated an association with lifetime cannabis use at cg15973234 (CEMIP): combined p = 3.3 × 10-8 . We found no overlap between published blood-based cis-meQTLs of cg15973234 and reported lifetime cannabis use-associated single nucleotide polymorphism (SNPs; p < .05), suggesting that the observed DNAm difference was driven by cannabis exposure. We also developed a multi-CpG classifier of lifetime cannabis use using penalized regression of top EWAS CpGs. The resulting 50-CpG classifier produced an area under the curve (AUC) = 0.74 (95% CI [0.72, 0.76], p = 2.00 × 10-5 ) in the discovery sample and AUC = 0.54 ([0.51, 0.57], p = 2.87 × 10-2 ) in the replication sample. Our EWAS findings provide evidence that blood-based DNAm is associated with lifetime cannabis use.

Genome-wide DNA methylation differences in nucleus accumbens of smokers vs. nonsmokers.

Numerous DNA methylation (DNAm) biomarkers of cigarette smoking have been identified in peripheral blood studies, but because of tissue specificity, blood-based studies may not detect brain-specific smoking-related DNAm differences that may provide greater insight as neurobiological indicators of smoking and its exposure effects. We report the first epigenome-wide association study (EWAS) of smoking in human postmortem brain, focusing on nucleus accumbens (NAc) as a key brain region in developing and reinforcing addiction. Illumina HumanMethylation EPIC array data from 221 decedents (120 European American [23% current smokers], 101 African American [26% current smokers]) were analyzed. DNAm by smoking (current vs. nonsmoking) was tested within each ancestry group using robust linear regression models adjusted for age, sex, cell-type proportion, DNAm-derived negative control principal components (PCs), and genotype-derived PCs. The resulting ancestry-specific results were combined via meta-analysis. We extended our NAc findings, using published smoking EWAS results in blood, to identify DNAm smoking effects that are unique (tissue-specific) vs. shared between tissues (tissue-shared). We identified seven CpGs (false discovery rate < 0.05), of which three CpGs are located near genes previously indicated with blood-based smoking DNAm biomarkers: ZIC1, ZCCHC24, and PRKDC. The other four CpGs are novel for smoking-related DNAm changes: ABLIM3, APCDD1L, MTMR6, and CTCF. None of the seven smoking-related CpGs in NAc are driven by genetic variants that share association signals with predisposing genetic risk variants for smoking, suggesting that the DNAm changes reflect consequences of smoking. Our results provide the first evidence for smoking-related DNAm changes in human NAc, highlighting CpGs that were undetected as peripheral biomarkers and may reflect brain-specific responses to smoking exposure.