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 (GWASs) 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 correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants 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 MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and 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.

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.

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.

Developmental effects of maternal smoking during pregnancy on the human frontal cortex transcriptome.

Cigarette smoking during pregnancy is a major public health concern. While there are well-described consequences in early child development, there is very little known about the effects of maternal smoking on human cortical biology during prenatal life. We therefore performed a genome-wide differential gene expression analysis using RNA sequencing (RNA-seq) on prenatal (N = 33; 16 smoking-exposed) as well as adult (N = 207; 57 active smokers) human postmortem prefrontal cortices. Smoking exposure during the prenatal period was directly associated with differential expression of 14 genes; in contrast, during adulthood, despite a much larger sample size, only two genes showed significant differential expression (FDR < 10%). Moreover, 1,315 genes showed significantly different exposure effects between maternal smoking during pregnancy and direct exposure in adulthood (FDR < 10%)-these differences were largely driven by prenatal differences that were enriched for pathways previously implicated in addiction and synaptic function. Furthermore, prenatal and age-dependent differentially expressed genes were enriched for genes implicated in non-syndromic autism spectrum disorder (ASD) and were differentially expressed as a set between patients with ASD and controls in postmortem cortical regions. These results underscore the enhanced sensitivity to the biological effect of smoking exposure in the developing brain and offer insight into how maternal smoking during pregnancy affects gene expression in the prenatal human cortex. They also begin to address the relationship between in utero exposure to smoking and the heightened risks for the subsequent development of neuropsychiatric disorders.

Age-related changes in the expression of schizophrenia susceptibility genes in the human prefrontal cortex.

The molecular basis of complex neuropsychiatric disorders most likely involves many genes. In recent years, specific genetic variations influencing risk for schizophrenia and other neuropsychiatric disorders have been reported. We have used custom DNA microarrays and qPCR to investigate the expression of putative schizophrenia susceptibility genes and related genes of interest in the normal human brain. Expression of 31 genes was measured in Brodmann's area 10 (BA10) in the prefrontal cortex of 72 postmortem brain samples spanning half a century of human aging (18-67 years), each without history of neuropsychiatric illness, neurological disease, or drug abuse. Examination of expression across age allowed the identification of genes whose expression patterns correlate with age, as well as genes that share common expression patterns and that possibly participate in common cellular mechanisms related to the emergence of schizophrenia in early adult life. The expression of GRM3 and RGS4 decreased across the entire age range surveyed, while that of PRODH and DARPP-32 was shown to increase with age. NRG1, ERBB3, and NGFR show expression changes during the years of greatest risk for the development of schizophrenia. Expression of FEZ1, GAD1, and RGS4 showed especially high correlation with one another, in addition to the strongest mean levels of absolute correlation with all other genes studied here. All microarray data are available at NCBI's Gene Expression Omnibus: GEO Series accession number GSE11546 (http://www.ncbi.nlm.nih.gov/geo) [corrected]

Critical factors in gene expression in postmortem human brain: Focus on studies in schizophrenia.

BACKGROUND: Studies of postmortem human brain are important for investigating underlying pathogenic molecular mechanisms of neuropsychiatric disorders. They are, however, confounded by pre- and postmortem factors. The purpose of this study was to identify sources of variation that will enable a better design of gene expression studies and higher reliability of gene expression data. METHODS: We assessed the contribution of multiple variables to messenger RNA (mRNA) expression of reference (housekeeping) genes measured by reverse transcriptase-polymerase chain reaction (RT-PCR) by multiple regression analysis in a large number (N = 143) of autopsy samples from the hippocampus and white and grey matter of the dorsolateral prefrontal cortex (DLPFC) of patients with schizophrenia and normal control subjects. RESULTS: The strongest predictor of gene expression was total RNA quality. Other significant factors included pH, postmortem interval, age and the duration of the agonal state, but the importance of these factors depended on transcript measured, brain region analyzed, and diagnosis. The quality of RNA obtained from the DLPFC white matter was also adversely affected by smoking. CONCLUSIONS: Our results show that normalization of expression data of target genes with a geometric mean of multiple housekeeping genes should be used to control for differences in RNA quality between samples. The results also suggest that accurate assessment of other confounding factors and their inclusion as regressors in the analysis is critical for obtaining reliable and accurate quantification of mRNA expression.