Adrian Boström

Epigenomics of Total Acute Sleep Deprivation in Relation to Genome-Wide DNA Methylation Profiles and RNA Expression.

Abstract Despite an established link between sleep deprivation and epigenetic processes in humans, it remains unclear to what extent sleep deprivation modulates DNA methylation. We performed a within-subject randomized blinded study with 16 healthy subjects to examine the effect of one night of total sleep deprivation (TSD) on the genome-wide methylation profile in blood compared with that in normal sleep. Genome-wide differences in methylation between both conditions were assessed by applying a paired regression model that corrected for monocyte subpopulations. In addition, the correlations between the methylation of genes detected to be modulated by TSD and gene expression were examined in a separate, publicly available cohort of 10 healthy male donors (E-GEOD-49065). Sleep deprivation significantly affected the DNA methylation profile both independently and in dependency of shifts in monocyte composition. Our study detected differential methylation of 269 probes. Notably, one CpG site was located 69 bp upstream of ING5, which has been shown to be differentially expressed after sleep deprivation. Gene set enrichment analysis detected the Notch and Wnt signaling pathways to be enriched among the differentially methylated genes. These results provide evidence that total acute sleep deprivation alters the methylation profile in healthy human subjects. This is, to our knowledge, the first study that systematically investigated the impact of total acute sleep deprivation on genome-wide DNA methylation profiles in blood and related the epigenomic findings to the expression data.

Lower inflammatory markers in women with antenatal depression brings the M1/M2 balance into focus from a new direction

BACKGROUND Antenatal depression and use of serotonin reuptake inhibitors (SSRI) in pregnancy have both been associated with an increased risk of poor pregnancy outcomes, such as preterm birth and impaired fetal growth. While the underlying biological pathways for these complications are poorly understood, it has been hypothesized that inflammation may be a common physiological pathway. The aim of the present study was to assess peripheral inflammatory markers in healthy women, women with antenatal depression, and in women using SSRI during pregnancy. METHODS 160 healthy pregnant controls, 59 women with antenatal depression and 39 women on treatment with SSRIs were included. The relative levels of 92 inflammatory proteins were analyzed by proximity extension assay technology. RESULTS Overall, 23 of the inflammatory markers were significantly lower in women with antenatal depression and in women on treatment with SSRIs in comparison with the healthy controls. No difference in any of the inflammatory markers was observed between women with antenatal depression and those who were using SSRI. Top three inflammatory markers that were down-regulated in women with antenatal depression were TNF-related apoptosis-inducing ligand (TRAIL), p=0.000001, macrophage colony-stimulating factor 1 (CSF-1), p=0.000004, and fractalkine (CX3CL1), p=0.000005. Corresponding inflammatory markers in SSRI users were CSF-1, p=0.000011, vascular endothelial growth factor A (VEGF-A), p=0.000016, and IL-15 receptor subunit alpha (IL-15RA), p=0.000027. The inflammatory markers were negatively correlated with cortisone serum concentrations in controls, but not in the cases. Differential DNA methylation of was found for seven of these inflammatory markers in an independent epigenetics cohort. CONCLUSION Women with antenatal depression or on SSRI treatment have lower levels of a number of peripheral inflammatory markers than healthy pregnant controls. Hypothetically, this could be due to dysregulated switch to the pro-M2 milieu that characterizes normal third trimester pregnancy. However, longitudinal blood sampling is needed to elucidate whether the presumably dysregulated M2 shift is driving the development of antenatal depression or is a result of the depression.

Inflammatory markers in late pregnancy in association with postpartum depression—A nested case-control study

Recent studies indicate that the immune system adaptation during pregnancy could play a significant role in the pathophysiology of perinatal depression. The aim of this study was to investigate if inflammation markers in a late pregnancy plasma sample can predict the presence of depressive symptoms at eight weeks postpartum. Blood samples from 291 pregnant women (median and IQR for days to delivery, 13 and 7-23days respectively) comprising 63 individuals with postpartum depressive symptoms, as assessed by the Edinburgh postnatal depression scale (EPDS≥12) and/or the Mini International Neuropsychiatric Interview (M.I.N.I.) and 228 controls were analyzed with an inflammation protein panel using multiplex proximity extension assay technology, comprising of 92 inflammation-associated markers. A summary inflammation variable was also calculated. Logistic regression, LASSO and Elastic net analyses were implemented. Forty markers were lower in late pregnancy among women with depressive symptoms postpartum. The difference remained statistically significant for STAM-BP (or otherwise AMSH), AXIN-1, ADA, ST1A1 and IL-10, after Bonferroni correction. The summary inflammation variable was ranked as the second best variable, following personal history of depression, in predicting depressive symptoms postpartum. The protein-level findings for STAM-BP and ST1A1 were validated in relation to methylation status of loci in the respective genes in a different population, using openly available data. This explorative approach revealed differences in late pregnancy levels of inflammation markers between women presenting with depressive symptoms postpartum and controls, previously not described in the literature. Despite the fact that the results do not support the use of a single inflammation marker in late pregnancy for assessing risk of postpartum depression, the use of STAM-BP or the novel notion of a summary inflammation variable developed in this work might be used in combination with other biological markers in the future.

Methylation of HPA axis related genes in men with hypersexual disorder

Hypersexual Disorder (HD) defined as non-paraphilic sexual desire disorder with components of compulsivity, impulsivity and behavioral addiction, and proposed as a diagnosis in the DSM 5, shares some overlapping features with substance use disorder including common neurotransmitter systems and dysregulated hypothalamic-pituitary-adrenal (HPA) axis function. In this study, comprising 67 HD male patients and 39 male healthy volunteers, we aimed to identify HPA-axis coupled CpG-sites, in which modifications of the epigenetic profile are associated with hypersexuality. The genome-wide methylation pattern was measured in whole blood using the Illumina Infinium Methylation EPIC BeadChip, measuring the methylation state of over 850K CpG sites. Prior to analysis, the global DNA methylation pattern was pre-processed according to standard protocols and adjusted for white blood cell type heterogeneity. We included CpG sites located within 2000bp of the transcriptional start site of the following HPA-axis coupled genes: Corticotropin releasing hormone (CRH), corticotropin releasing hormone binding protein (CRHBP), corticotropin releasing hormone receptor 1 (CRHR1), corticotropin releasing hormone receptor 2 (CRHR2), FKBP5 and the glucocorticoid receptor (NR3C1). We performed multiple linear regression models of methylation M-values to a categorical variable of hypersexuality, adjusting for depression, dexamethasone non-suppression status, Childhood Trauma Questionnaire total score and plasma levels of TNF-alpha and IL-6. Of 76 tested individual CpG sites, four were nominally significant (p<0.05), associated with the genes CRH, CRHR2 and NR3C1. Cg23409074-located 48bp upstream of the transcription start site of the CRH gene - was significantly hypomethylated in hypersexual patients after corrections for multiple testing using the FDR-method. Methylation levels of cg23409074 were positively correlated with gene expression of the CRH gene in an independent cohort of 11 healthy male subjects. The methylation levels at the identified CRH site, cg23409074, were significantly correlated between blood and four different brain regions. CRH is an important integrator of neuroendocrine stress responses in the brain, with a key role in the addiction processes. Our results show epigenetic changes in the CRH gene related to hypersexual disorder in men.

Plasma 1-deoxysphingolipids are early predictors of incident type 2 diabetes mellitus

1-Deoxysphingolipids (1-deoxySLs) are atypical sphingolipids, which are formed in a side reaction during sphingolipid de-novo synthesis. Recently, we demonstrated that 1-deoxySLs are biomarkers for the prediction of T2DM in obese, non-diabetic patients. Here we investigated the relevance of 1-deoxySLs as long-term predictive biomarkers for the incidence of T2DM in an asymptomatic population. Here, we analyzed the plasma sphingoid base profile in a nested group of non-diabetic individuals (N = 605) selected from a population-based study including 5 year follow-up data (CoLaus study). 1-DeoxySLs at baseline were significantly elevated in individuals who developed T2DM during the follow-up (p<0.001), together with increased glucose (p<5.11E-14), triglycerides (p<0.001) and HOMA-IR indices (p<0.001). 1-Deoxy-sphinganine (1-deoxySA) and 1-deoxy-sphingosine (1-deoxySO) were predictive for T2DM, even after adjusting for fasting glucose levels in the binary regression analyses. The predictive value of the combined markers 1-deoxySA+glucose were superior to glucose alone in normal-weight subjects (p<0.001) but decreased substantially with increasing BMI. Instead, plasma adiponectin and waist-to-hip ratio appeared to be better risk predictors for obese individuals (BMI>30kg/m2). In conclusion, elevated plasma 1-deoxySL levels are strong and independent risk predictors of future T2DM, especially for non-obese individuals in the general population.

Chenodeoxycholic acid significantly impacts the expression of miRNAs and genes involved in lipid, bile acid and drug metabolism in human hepatocytes.

AIMS Bile acids (BAs) are important gut signaling hormones, influencing lipid, glucose, and energy homeostasis. The exact mechanisms behind these effects are not yet fully understood. Lately, they have come to the fore as putative therapeutics in metabolic diseases, such as e.g. nonalcoholic fatty liver disease (NAFLD). We elucidate to what extent BAs impacts on the mRNAome and microRNAome in hepatocytes to gather novel insights into the mechanisms behind metabolic and toxicologic effects of bile acids. MAIN METHODS Five batches of primary human hepatocytes were treated with 50μmol/l chenodeoxycholic acid (CDCA) for 24 or 48h. Total RNA was extracted, size fractionated and subjected to Next Generation Sequencing to generate mRNA and miRNA profiles. KEY FINDINGS Expression of 738 genes and 52 miRNAs were CDCA dependently decreased, whereas 1566 genes and 29 miRNAs were significantly increased in hepatocytes. Distinct gene clusters controlling BA and lipid homeostasis (FGF(R), APO and FABP family members, HMGCS2) and drug metabolism (CYP, UGT and SULT family members) were significantly modulated by CDCA. Importantly, CDCA affected distinct microRNAs, including miR-34a, -505, -885, -1260 and -552 that systematically correlated in expression with gene clusters responsible for bile acid, lipid and drug homeostasis incorporating genes, such as e.g. SLCO1B1, SLC22A7, FGF19, CYP2E1, CYP1A2, APO family members and FOXO3. SIGNIFICANCE Bile acids significantly modulate metabolic and drug associated gene networks that are connected to distinct shifts in the microRNAome These findings give novel insights on how BA enfold metabolic and system toxic effects.

A methylome-wide mQTL analysis reveals associations of methylation sites with GAD1 and HDAC3 SNPs and a general psychiatric risk score

Genome-wide association studies have identified a number of single-nucleotide polymorphisms (SNPs) that are associated with psychiatric diseases. Increasing body of evidence suggests a complex connection of SNPs and the transcriptional and epigenetic regulation of gene expression, which is poorly understood. In the current study, we investigated the interplay between genetic risk variants, shifts in methylation and mRNA levels in whole blood from 223 adolescents distinguished by a risk for developing psychiatric disorders. We analyzed 37 SNPs previously associated with psychiatric diseases in relation to genome-wide DNA methylation levels using linear models, with Bonferroni correction and adjusting for cell-type composition. Associations between DNA methylation, mRNA levels and psychiatric disease risk evaluated by the Development and Well-Being Assessment (DAWBA) score were identified by robust linear models, Pearson’s correlations and binary regression models. We detected five SNPs (in HCRTR1, GAD1, HADC3 and FKBP5) that were associated with eight CpG sites, validating five of these SNP–CpG pairs. Three of these CpG sites, that is, cg01089319 (GAD1), cg01089249 (GAD1) and cg24137543 (DIAPH1), manifest in significant gene expression changes and overlap with active regulatory regions in chromatin states of brain tissues. Importantly, methylation levels at cg01089319 were associated with the DAWBA score in the discovery group. These results show how distinct SNPs linked with psychiatric diseases are associated with epigenetic shifts with relevance for gene expression. Our findings give a novel insight on how genetic variants may modulate risks for the development of psychiatric diseases.

A targeted analysis reveals relevant shifts in the methylation and transcription of genes responsible for bile acid homeostasis and drug metabolism in non-alcoholic fatty liver disease

BackgroundNon-alcoholic fatty liver disease (NAFLD) is associated with a high risk for liver cirrhosis and cancer. Recent studies demonstrate that NAFLD significantly impacts on the genome wide methylation and expression reporting top hit genes to be associated with e.g. diabetes mellitus. In a targeted analysis we specifically investigate to what extent NAFLD is associated with methylation and transcriptional changes in gene networks responsible for drug metabolism (DM) and bile acid (BA) homeostasis, which may trigger liver and system toxic events.MethodsWe performed a systematic analysis of 73 genes responsible for BA homeostasis and DM based on liver derived methylation and expression data from three cohort studies including 103 NAFLD and 75 non-NAFLD patients. Using multiple linear regression models, we detected methylation differences in proximity to the transcriptional start site of these genes in two NAFLD cohorts and correlated the methylation of significantly changed CpG sites to transcriptional expression in a third cohort using robust multiple linear regression approaches.ResultsWe detected 64 genes involved in BA homeostasis and DM to be significantly differentially methylated. In 26 of these genes, methylation significantly correlated with RNA expression, detecting i.e. genes such as CYP27A1, OSTɑ, and SLC27A5 (BA homeostasis), and SLCO2B1, SLC47A1, and several UGT and CYP genes (DM) to be NAFLD dependently modulated.ConclusionsNAFLD is associated with significant shifts in the methylation of key genes responsible for BA and DM that are associated with transcriptional modulations. These findings have implications for BA composition, BA regulated metabolic pathways and for drug safety and efficacy.

A MIR4646 associated methylation locus is hypomethylated in adolescent depression

BACKGROUND Studies of epigenetics and transcriptional activity in adolescents may provide knowledge about possible preventive strategies of depression. METHODS We present a methylome-wide association study (MWAS) and cohort validation analysis of depression in adolescents, in two separate cohorts: discovery (n=93) and validation data set 1 (n=78). The genome-wide methylation pattern was measured from whole blood using the Illumina 450K array. A second validation cohort, validation data set 2, consists of post-mortem brain biopsies from depressed adults (n=58). We performed a MWAS by robust multiple linear regressions of methylation to a modified risk-score assessment of depression. Methylation levels of candidate CpG sites were correlated with expression levels of the associated gene in an independent cohort of 11 healthy volunteers. RESULTS The methylation state of two CpG sites reliably predicted ratings of depression in adolescents (cg13227623 and cg04102384) (p<10E-06). Cohort validation analysis confirmed cg04102384 - located in the promoter region of microRNA 4646 (MIR4646) - to be hypomethylated in both validation data set 1 and validation data set 2 (p<0.05). Cg04102384 was inversely correlated to expression levels of MIR4646-3p in healthy controls (p<0.05). LIMITATIONS MIR4646 was not differentially expressed in a subset of samples with adolescent depression measured by qRT-PCR measurements. CONCLUSION We identify a specific MIR4646 associated epigenetic risk site to be associated with depression in adolescents. Cg04102384 putatively regulates gene expression of MIR4646-3p. Target gene prediction and gene set overrepresentation analysis revealed involvement of this miRNA in fatty acid elongation, a process related to omega-3 fatty acids, previously associated with depression.

Epigenetic changes in the CRH gene are related to severity of suicide attempt and a general psychiatric risk score in adolescents

The aim of this study, comprising 88 suicide attempters, was to identify hypothalamic-pituitary-adrenal (HPA) -axis coupled CpG-sites showing methylation shifts linked to severity of the suicide attempt. Candidate methylation loci were further investigated as risk loci for a general psychiatric risk score in two cohorts of adolescents (cohort 1 and 2). The genome-wide methylation pattern was measured in whole blood using the Illumina Infinium Methylation EPIC BeadChip. Subjects were stratified into high-risk and low-risk groups based on the severity of the suicidal behavior. We included CpG sites located within 2000 basepairs away from transcriptional start site of the following HPA-axis coupled genes: corticotropin releasing hormone (CRH), corticotropin releasing hormone binding protein (CRHBP), corticotropin releasing hormone receptor 1 (CRHR1), corticotropin releasing hormone receptor 2 (CRHR2), FK506-binding protein 51 (FKBP5) and the glucocorticoid receptor (NR3C1). The methylation state of two corticotropin releasing hormone (CRH)-associated CpG sites were significantly hypomethylated in the high-risk group of suicide attempters (n = 31) (cg19035496 and cg23409074) (p < 0.001). Adolescent cohort 1 and 2 consisted of 129 and 93 subjects, respectively, and were stratified by the in silico generated DAWBA measurements of a general psychiatric risk score into high-risk group (>~50% risk) or controls. In adolescent cohort 2, cg19035496 was hypermethylated in subjects with a high general psychiatric risk score. Our results show epigenetic changes in the CRH gene related to severity of suicide attempt in adults and a general psychiatric risk score in adolescents.