Antidepressant Exposure and DNA Methylation: Insights from a Methylome-Wide Association Study

Importance

Understanding antidepressant mechanisms could help design more effective and tolerated treatments.

Objective

Identify DNA methylation (DNAm) changes associated with antidepressant exposure.

Design

Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS.

Setting

Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR.

Participants

Participants with DNAm data and self-report/prescription derived antidepressant exposure.

Main Outcomes and Measures

Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions.

Results

The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ).

Conclusions and Relevance

Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments.

3 Key Points

Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

Genome-wide gene-environment interaction in depression: A systematic evaluation of candidate genes: The childhood trauma working-group of PGC-MDD

Gene by environment (GxE) interaction studies have investigated the influence of a number of candidate genes and variants for major depressive disorder (MDD) on the association between childhood trauma and MDD. Most of these studies are hypothesis driven and investigate only a limited number of SNPs in relevant pathways using differing methodological approaches. Here (1) we identified 27 genes and 268 SNPs previously associated with MDD or with GxE interaction in MDD and (2) analyzed their impact on GxE in MDD using a common approach in 3944 subjects of European ancestry from the Psychiatric Genomics Consortium who had completed the Childhood Trauma Questionnaire. (3) We subsequently used the genome-wide SNP data for a genome-wide case-control GxE model and GxE case-only analyses testing for an enrichment of associated SNPs. No genome-wide significant hits and no consistency among the signals of the different analytic approaches could be observed. This is the largest study for systematic GxE interaction analysis in MDD in subjects of European ancestry to date. Most of the known candidate genes/variants could not be supported. Thus, their impact on GxE interaction in MDD may be questionable. Our results underscore the need for larger samples, more extensive assessment of environmental exposures, and greater efforts to investigate new methodological approaches in GxE models for MDD.

GWAS Identifies New Loci for Painful Temporomandibular Disorder: Hispanic Community Health Study/Study of Latinos

Temporomandibular disorder (TMD) is a musculoskeletal condition characterized by pain and reduced function in the temporomandibular joint and/or associated masticatory musculature. Prevalence in the United States is 5% and twice as high among women as men. We conducted a discovery genome-wide association study (GWAS) of TMD in 10,153 participants (769 cases, 9,384 controls) of the US Hispanic Community Health Study/Study of Latinos (HCHS/SOL). The most promising single-nucleotide polymorphisms (SNPs) were tested in meta-analysis of 4 independent cohorts. One replication cohort was from the United States, and the others were from Germany, Finland, and Brazil, totaling 1,911 TMD cases and 6,903 controls. A locus near the sarcoglycan alpha ( SGCA), rs4794106, was suggestive in the discovery analysis ( P = 2.6 × 10⁶) and replicated (i.e., 1-tailed P = 0.016) in the Brazilian cohort. In the discovery cohort, sex-stratified analysis identified 2 additional genome-wide significant loci in females. One lying upstream of the relaxin/insulin-like family peptide receptor 2 ( RXP2) (chromosome 13, rs60249166, odds ratio [OR] = 0.65, P = 3.6 × 10^-8) was replicated among females in the meta-analysis (1-tailed P = 0.052). The other (chromosome 17, rs1531554, OR = 0.68, P = 2.9 × 10^-8) was replicated among females (1-tailed P = 0.002), as well as replicated in meta-analysis of both sexes (1-tailed P = 0.021). A novel locus at genome-wide level of significance (rs73460075, OR = 0.56, P = 3.8 × 10^-8) in the intron of the dystrophin gene DMD (X chromosome), and a suggestive locus on chromosome 7 (rs73271865, P = 2.9 × 10^-7) upstream of the Sp4 Transcription Factor ( SP4) gene were identified in the discovery cohort, but neither of these was replicated. The SGCA gene encodes SGCA, which is involved in the cellular structure of muscle fibers and, along with DMD, forms part of the dystrophin-glycoprotein complex. Functional annotation suggested that several of these variants reside in loci that regulate processes relevant to TMD pathobiologic processes.

Exome Variant Analysis of Chronic Periodontitis in 2 Large Cohort Studies

Periodontitis is characterized by inflammation of the gingival tissue. The main risk factors are socioeconomic factors, sex, age, smoking, and diabetes, but periodontal disease has also a genetic background. Previous genome-wide association studies failed to reveal genome-wide significant associations of single common single-nucleotide polymorphisms with chronic periodontitis. Using the Illumina ExomeChip data of 6,576 participants of the German population-based cohort studies Study of Health in Pomerania (SHIP) and SHIP-Trend, the authors performed single variant and also gene-based association studies of rare and common exonic variations on different periodontal case definitions. Although our study comprised the largest sample size to date to assess genetic predisposition for chronic periodontitis, the authors found no significant association. This study emphasizes that for chronic periodontitis, large sample sizes will be necessary to find genetic associations, even when examining rare genetic variants.

GWAS for executive function and processing speed suggests involvement of the CADM2 gene

To identify common variants contributing to normal variation in two specific domains of cognitive functioning, we conducted a genome-wide association study (GWAS) of executive functioning and information processing speed in non-demented older adults from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) consortium. Neuropsychological testing was available for 5429-32,070 subjects of European ancestry aged 45 years or older, free of dementia and clinical stroke at the time of cognitive testing from 20 cohorts in the discovery phase. We analyzed performance on the Trail Making Test parts A and B, the Letter Digit Substitution Test (LDST), the Digit Symbol Substitution Task (DSST), semantic and phonemic fluency tests, and the Stroop Color and Word Test. Replication was sought in 1311-21860 subjects from 20 independent cohorts. A significant association was observed in the discovery cohorts for the single-nucleotide polymorphism (SNP) rs17518584 (discovery P-value=3.12 × 10(-8)) and in the joint discovery and replication meta-analysis (P-value=3.28 × 10(-9) after adjustment for age, gender and education) in an intron of the gene cell adhesion molecule 2 (CADM2) for performance on the LDST/DSST. Rs17518584 is located about 170 kb upstream of the transcription start site of the major transcript for the CADM2 gene, but is within an intron of a variant transcript that includes an alternative first exon. The variant is associated with expression of CADM2 in the cingulate cortex (P-value=4 × 10(-4)). The protein encoded by CADM2 is involved in glutamate signaling (P-value=7.22 × 10(-15)), gamma-aminobutyric acid (GABA) transport (P-value=1.36 × 10(-11)) and neuron cell-cell adhesion (P-value=1.48 × 10(-13)). Our findings suggest that genetic variation in the CADM2 gene is associated with individual differences in information processing speed.

Gene transcripts associated with muscle strength: a CHARGE meta-analysis of 7,781 persons

Lower muscle strength in midlife predicts disability and mortality in later life. Blood-borne factors, including growth differentiation factor 11 (GDF11), have been linked to muscle regeneration in animal models. We aimed to identify gene transcripts associated with muscle strength in adults. Meta-analysis of whole blood gene expression (overall 17,534 unique genes measured by microarray) and hand-grip strength in four independent cohorts (n = 7,781, ages: 20-104 yr, weighted mean = 56), adjusted for age, sex, height, weight, and leukocyte subtypes. Separate analyses were performed in subsets (older/younger than 60, men/women). Expression levels of 221 genes were associated with strength after adjustment for cofactors and for multiple statistical testing, including ALAS2 (rate-limiting enzyme in heme synthesis), PRF1 (perforin, a cytotoxic protein associated with inflammation), IGF1R, and IGF2BP2 (both insulin like growth factor related). We identified statistical enrichment for hemoglobin biosynthesis, innate immune activation, and the stress response. Ten genes were associated only in younger individuals, four in men only and one in women only. For example, PIK3R2 (a negative regulator of PI3K/AKT growth pathway) was negatively associated with muscle strength in younger (<60 yr) individuals but not older (≥ 60 yr). We also show that 115 genes (52%) have not previously been linked to muscle in NCBI PubMed abstracts. This first large-scale transcriptome study of muscle strength in human adults confirmed associations with known pathways and provides new evidence for over half of the genes identified. There may be age- and sex-specific gene expression signatures in blood for muscle strength.

A common atopy-associated variant in the Th2 cytokine locus control region impacts transcriptional regulation and alters SMAD3 and SP1 binding

BACKGROUND

Type 2 immune responses directed by Th2 cells and characterized by the signature cytokines IL4, IL5, and IL13 play major pathogenic roles in atopic diseases. Single nucleotide polymorphisms in the human Th2 cytokine locus in particular in a locus control region within the DNA repair gene RAD50, containing several RAD50 DNase1-hypersensitive sites (RHS), have been robustly associated with atopic traits in genome-wide association studies (GWAS). Functional variants in IL13 have been intensely studied, whereas no causative variants for the IL13-independent RAD50 signal have been identified yet. This study aimed to characterize the functional impact of the atopy-associated polymorphism rs2240032 located in the human RHS7 on cis-regulatory activity and differential binding of transcription factors.

METHODS

Differential transcription factor binding was analyzed by electrophoretic mobility shift assays (EMSAs) with Jurkat T-cell nuclear extracts. Identification of differentially binding factors was performed using mass spectrometry (LC-MS/MS). Reporter vector constructs carrying either the major or minor allele of rs2240032 were tested for regulating transcriptional activity in Jurkat and HeLa cells.

RESULTS

The variant rs2240032 impacts transcriptional activity and allele-specific binding of SMAD3, SP1, and additional putative protein complex partners. We further demonstrate that rs2240032 is located in an RHS7 subunit which itself encompasses repressor activity and might be important for the fine-tuning of transcription regulation within this region.

CONCLUSION

The human RHS7 critically contributes to the regulation of gene transcription, and the common atopy-associated polymorphism rs2240032 impacts transcriptional activity and transcription factor binding.

Genome-wide association analysis of coffee drinking suggests association with CYP1A1/CYP1A2 and NRCAM

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).

Mendelian randomization suggests non-causal associations of testosterone with cardiometabolic risk factors and mortality

Prospective studies showed that low serum testosterone concentrations are associated with various cardiometabolic risk factors and mortality. However, the causal nature of these associations is controversial. We studied 1 882 men aged 20-79 years with serum testosterone concentrations and genotyping data from the longitudinal population-based Study of Health in Pomerania. Testosterone concentrations were cross-sectionally associated with cardiometabolic risk factors, including anthropometric, lipid, blood pressure and glycaemic parameters; and prospectively with all-cause mortality (277 deaths, 14.7%) during the 10-year follow-up. To overcome problems of residual confounding, reverse causation, or regression dilution bias in the investigated testosterone-outcome associations, we used two-stage least square regression models with previously identified polymorphisms at the SHBG gene (rs12150660) and X chromosome (rs5934505) as multiple genetic instruments in an instrumental variable (IV) approach, also known as Mendelian randomization. In standard regression analyses, testosterone was robustly associated with a wide range of cardiometabolic risk factors. In subsequent IV analyses, no such significant associations were observed. Similarly, prospective analyses showed a consistent association of low testosterone concentrations with increased all-cause mortality risk, which was not apparent in subsequent IV analyses. The present Mendelian randomization analyses did not detect any evidence for causal associations of testosterone concentrations with cardiometabolic risk factors and mortality, suggesting that previously reported associations might largely result from residual confounding or reverse causation. Although testosterone assessment might improve risk prediction, implementation of testosterone replacement therapy requires further evidence of a direct effect on cardiometabolic outcomes from double-blinded randomized controlled trials and large-scale Mendelian randomization meta-analyses.

The effect of catechol-O-methyltransferase polymorphisms on pain is modified by depressive symptoms

BACKGROUND

Variations within the catechol-O-methyltransferase (COMT) gene have been associated with pain severity in temporomandibular disorders (TMDs). Psychological factors such as personal conflicts, life stress and depression, are well known to be associated with onset, severity and chronicity of pain disorders.

AIM

We hypothesized that the relationship between the COMT gene and TMD pain is modified by depressive symptoms.

METHODS

Cross-sectional data from the population-based Study of Health in Pomerania (SHIP) in Germany were used to estimate additive interactions between depressive symptoms and 22 single-nucleotide polymorphisms (SNPs) of the COMT gene and the neighbouring thioredoxin reductase 2 (TXNRD2) gene on TMD pain. All participants were Caucasian subjects from a rural area in Northeast Germany. After exclusion of 79 subjects with antidepressant medication, 29.9% of the remaining 3904 subjects reported lifetime depressive symptoms. TMD pain was assessed by a standardized clinical examination. Among various TMD signs, only those that assessed muscle or joint pain on palpation were used as recommended.

RESULTS

Six SNPs from the first of three COMT/TXNRD2 haploblocks interacted with depressive symptoms on TMD pain (smallest p-value: 2.7 × 10(-10) ). In subjects without depressive symptoms, rs5993882 was identified as the SNP most likely to be related to TMD pain. In subjects with symptoms of depression, rs1544325 was the corresponding top COMT SNP.

CONCLUSIONS

Our results indicate that variants within the COMT gene are associated with pain perception. However, this association is highly moderated by the absence or presence of lifetime depressive symptoms.