Targeted bisulfite sequencing: A novel tool for the assessment of DNA methylation with high sensitivity and increased coverage

A varying-coefficient model for the analysis of methylation sequencing data

DNA methylation is an important epigenetic modification involved in gene regulation. Advances in the next generation sequencing technology have enabled the retrieval of DNA methylation information at single-base-resolution. However, due to the sequencing process and the limited amount of isolated DNA, DNA-methylation-data are often noisy and sparse, which complicates the identification of differentially methylated regions (DMRs), especially when few replicates are available. We present a varying-coefficient model for detecting DMRs by using single-base-resolved methylation information. The model simultaneously smooths the methylation profiles and allows detection of DMRs, while accounting for additional covariates. The proposed model takes into account possible overdispersion by using a beta-binomial distribution. The overdispersion itself can be modeled as a function of the genomic region and explanatory variables. We illustrate the properties of the proposed model by applying it to two real-life case studies.

Differential methylation of OPRK1 in borderline personality disorder is associated with childhood trauma

According to a growing body of neurobiological evidence, the core symptoms of borderline personality disorder (BPD) may be linked to an opioidergic imbalance between the hedonic and stimulatory activity of mu opioid receptors (MOR) and the reward system inhibiting effects of kappa opioid receptors (KOR). Childhood trauma (CT), which is etiologically relevant to BPD, is also likely to lead to epigenetic and neurobiological adaptations by extensive activation of the stress and endogenous opioid systems. In this study, we investigated the methylation differences in the promoter of the KOR gene (OPRK1) in subjects with BPD (N = 47) and healthy controls (N = 48). Comparing the average methylation rates of regulatorily relevant subregions (specified regions CGI-1, CGI-2, EH1), we found no differences between BPD and HC. Analyzing individual CG nucleotides (N = 175), we found eight differentially methylated CG sites, all of which were less methylated in BPD, with five showing highly interrelated methylation rates. This differentially methylated region (DMR) was found on the falling slope (5') of the promoter methylation gap, whose effect is enhanced by the DMR hypomethylation in BPD. A dimensional assessment of the correlation between disease severity and DMR methylation rate revealed DMR hypomethylation to be negatively associated with BPD symptom severity (measured by BSL-23). Finally, analyzing the influence of CT on DMR methylation, we found DMR hypomethylation to correlate with physical and emotional neglect in childhood (quantified by CTQ). Thus, the newly identified DMR may be a biomarker of the risks caused by CT, which likely epigenetically contribute to the development of BPD.

DNA methylation levels of the serotonin transporter gene are not associated with the outcome of highly standardized one-session exposure-based fear treatment

Epigenetic alterations are regarded as a potential mechanism mediating the effects of environmental risk factors on vulnerability for a range of mental health problems. Recent studies have addressed the question whether DNA methylation patterns predict the outcome of psychological interventions and whether treatment effects might be associated with changes of DNA methylation. We assessed phobic fear symptoms, treatment-relevant traits and treatment response in 308 adults free of psychotropic medication - highly fearful of either spiders, blood-injury-injections, dental-treatments or heights - all subjected to highly standardized exposure-based one-session fear treatment. DNA methylation level of the promotor region of the serotonin transporter gene (SLC6A4) was assessed in either saliva samples (spider and dental treatment fear cohorts) or oral mucosa (BII, heights) to check whether possible effects are independent of the surrogate tissue examined. Moreover, in order to examine possible DNA methylation by genotype effects, patients were assessed for genetic variation of the serotonin transporter-linked polymorphic region (5-HTTLPR). DNA methylation levels were neither associated with pre-treatment fear levels, treatment relevant traits or treatment outcome data even when allelic variation of the 5HTTLPR was considered. Overall DNA methylation levels were higher in saliva samples compared to buccal samples. In saliva samples there was a small pre- to post-treatment increase in DNA methylation, which, however, was also not associated with the investigated phenotypes. We conclude that DNA methylation of SLC6A4 is no suitable biomarker for response efficacy to highly standardized one-session exposure-based fear treatments.

Triggering pyroptosis enhances the antitumor efficacy of PARP inhibitors in prostate cancer

PURPOSE

PARP inhibitors have revolutionized the treatment landscape for advanced prostate cancer (PCa) patients who harboring mutations in homologous recombination repair (HRR) genes. However, the molecular mechanisms underlying PARP inhibitors function beyond DNA damage repair pathways remain elusive, and identifying novel predictive targets that favorably respond to PARP inhibitors in PCa is an active area of research.

METHODS

The expression of GSDME in PCa cell lines and human PCa samples was determined by western blotting. Targeted bisulfite sequencing, gene enrichment analysis (GSEA), clone formation, construction of the stably transfected cell lines, lactate dehydrogenase (LDH) assay, western blotting as well as a mouse model of subcutaneous xenografts were used to investigate the role of GSDME in PCa. The combinational therapeutic effect of olaparib and decitabine was determined using both in vitro and in vivo experiments.

RESULTS

We have found low expression of GSDME in PCa. Interestingly, we demonstrated that GSDME activity is robustly induced in olaparib-treated cells undergoing pyroptosis, and that high methylation of the GSDME promoter dampens its activity in PCa cells. Intriguingly, genetically overexpressing GSDME does not inhibit tumor cell proliferation but instead confers sensitivity to olaparib. Furthermore, pharmacological treatment with the combination of olaparib and decitabine synergistically induces GSDME expression and cleavage through caspase-3 activation, thus promoting pyroptosis and enhancing anti-tumor response, ultimately resulting in tumor remission.

CONCLUSION

Our findings highlight a novel therapeutic strategy for enhancing the long-term response to olaparib beyond HRR-deficient tumors in PCa, underscoring the critical role of GSDME in regulating tumorigenesis.

Introduction of a multiplex amplicon sequencing assay to quantify DNA methylation in target cytosine markers underlying four selected epigenetic clocks

BACKGROUND

DNA methylation analysis has proven to be a powerful tool for age assessment. However, the implementation of epigenetic age prediction in diagnostics or routine forensic casework requires appropriate laboratory methods. In this study, we aimed to compare the performance of large-scale DNA methylation analysis protocols that show promise in terms of accuracy, throughput, multiplexing capacity, and high sensitivity.

RESULTS

The protocols were designed to target a predefined panel of 161 genomic CG/CA sites from four known estimators of epigenetic age-related parameters, optimized and validated using artificially methylated controls or blood samples. We successfully targeted 96% of these loci using two enrichment protocols: Ion AmpliSeq™, an amplicon-based method integrated with Ion Torrent S5, and SureSelectXT Methyl-Seq, a hybridization-based method followed by MiSeq FGx sequencing. Both protocols demonstrated high accuracy and robustness. Although hybridization assays have greater multiplexing capabilities, the best overall performance was observed for the amplicon-based protocol with the lowest variability in DNA methylation at 25 ng of starting DNA, mean observed marker coverage of ~ 6.7 k reads, and accuracy of methylation quantification with a mean absolute difference between observed and expected methylation beta value of 0.054. The Ion AmpliSeq method correlated strongly with genome-scale EPIC microarray data (R = 0.91) and showed superiority in terms of methylation measurement accuracy. Method-to-method bias was accounted for by the use of linear transformation, which provided a highly accurate prediction of calendar age with a mean absolute error of less than 5 years for the VISAGE and Hannum age clocks used. The pace of aging (PoAm) and the mortality risk score (MRS) estimators included in our panel represent next-generation clocks, were found to have low to moderate correlations with the VISAGE and Hannum models (R < 0.75), and thus may capture different aspects of epigenetic aging.

CONCLUSIONS

We propose a laboratory tool that allows the quantification of DNA methylation in cytosines underlying four different clocks, thus providing broad information on epigenetic aging while maintaining a reasonable number of CpG markers, opening the way to a wide range of applications in forensics, medicine, and healthcare.

The DNA methylation landscape of the human oxytocin receptor gene (OXTR): data-driven clusters and their relation to gene expression and childhood adversity

The oxytocin receptor gene (OXTR) is of interest when investigating the effects of early adversity on DNA methylation. However, there is heterogeneity regarding the selection of the most promising CpG sites to target for analyses. The goal of this study was to determine functionally relevant clusters of CpG sites within the OXTR CpG island in 113 mother-infant dyads, with 58 of the mothers reporting childhood maltreatment (CM). OXTR DNA methylation was analyzed in peripheral/umbilical blood mononuclear cells. Different complexity reduction approaches were used to reduce the 188 CpG sites into clusters of co-methylated sites. Furthermore, associations between OXTR DNA methylation (cluster- and site-specific level) and OXTR gene expression and CM were investigated in mothers. Results showed that, first, CpG sections differed strongly regarding their statistical utility for research of individual differences in DNA methylation. Second, cluster analyses and Partial Least Squares (PLS) suggested two clusters consisting of intron1/exon2 and the protein-coding region of exon3, respectively, as most strongly associated with outcome measures. Third, cross-validated PLS regression explained 7% of variance in CM, with low cross-validated variance explained for the prediction of gene expression. Fourth, substantial mother-child correspondence was observed in correlation patterns within the identified clusters, but only modest correspondence outside these clusters. This study makes an important contribution to the mapping of the DNA methylation landscape of the OXTR CpG island by highlighting clusters of CpG sites that show desirable statistical properties and predictive value. We provide a Companion Web Application to facilitate the choice of CpG sites.

Genome-Wide DNA Methylation and Transcriptome Integration Associates DNA Methylation Changes with Bovine Subclinical Mastitis Caused by Staphylococcus chromogenes

Staphylococcus chromogenes (SC) is a common coagulase-negative staphylococcus described as an emerging mastitis pathogen and commonly found in dairy farms. This study investigated the potential involvement of DNA methylation in subclinical mastitis caused by SC. The whole-genome DNA methylation patterns and transcriptome profiles of milk somatic cells from four cows with naturally occurring SC subclinical mastitis (SCM) and four healthy cows were characterized by next-generation sequencing, bioinformatics, and integration analyses. Comparisons revealed abundant DNA methylation changes related to SCM, including differentially methylated cytosine sites (DMCs, n = 2,163,976), regions (DMRs, n = 58,965), and methylation haplotype blocks (dMHBs, n = 53,098). Integration of methylome and transcriptome data indicated a negative global association between DNA methylation at regulatory regions (promoters, first exons, and first introns) and gene expression. A total of 1486 genes with significant changes in the methylation levels of their regulatory regions and corresponding gene expression showed significant enrichment in biological processes and pathways related to immune functions. Sixteen dMHBs were identified as candidate discriminant signatures, and validation of two signatures in more samples further revealed the association of dMHBs with mammary gland health and production. This study demonstrated abundant DNA methylation changes with possible involvement in regulating host responses and potential as biomarkers for SCM.

Choline Regulates SOX4 through miR-129-5p and Modifies H3K27me3 in the Developing Cortex

Choline availability regulates neural progenitor cell proliferation and differentiation in the developing cerebral cortex. Here, we investigated the molecular mechanism underlying this process and demonstrated that choline regulates the transcription factor SOX4 in neural progenitor cells. Specifically, we found that low choline intake during neurogenesis reduces SOX4 protein levels, causing the downregulation of EZH2, a histone methyltransferase. Importantly, we demonstrate that low choline is not involved in SOX4 protein degradation rate and established that protein reduction is caused by aberrant expression of a microRNA (miR-129-5p). To confirm the role of miR-129-5p, we conducted gain-of-function and loss-of-function assays in neural progenitor cells and demonstrated that directly altering miR-129-5p levels could affect SOX4 protein levels. We also observed that the reduction in SOX4 and EZH2 led to decreased global levels of H3K27me3 in the developing cortex, contributing to reduced proliferation and precocious differentiation. For the first time, to our knowledge, we demonstrate that a nutrient, choline, regulates a master transcription factor and its downstream targets, providing a novel insight into the role of choline in brain development.

Epigenetic biomarkers for smoking cessation

Cigarette smoking has been associated with epigenetic alterations that may be reversible upon cessation. As the most-studied epigenetic modification, DNA methylation is strongly associated with smoking exposure, providing a potential mechanism that links smoking to adverse health outcomes. Here, we reviewed the reversibility of DNA methylation in accessible peripheral tissues, mainly blood, in relation to cigarette smoking cessation and the utility of DNA methylation as a biomarker signature to differentiate current, former, and never smokers and to quantify time since cessation. We summarized thousands of differentially methylated Cytosine-Guanine (CpG) dinucleotides and regions associated with smoking cessation from candidate gene and epigenome-wide association studies, as well as the prediction accuracy of the multi-CpG predictors for smoking status. Overall, there is robust evidence for DNA methylation signature of cigarette smoking cessation. However, there are still gaps to fill, including (1) cell-type heterogeneity in measuring blood DNA methylation; (2) underrepresentation of non-European ancestry populations; (3) limited longitudinal data to quantitatively measure DNA methylation after smoking cessation over time; and (4) limited data to study the impact of smoking cessation on other epigenetic features, noncoding RNAs, and histone modifications. Epigenetic machinery provides promising biomarkers that can improve success in smoking cessation in the clinical setting. To achieve this goal, larger and more-diverse samples with longitudinal measures of a broader spectrum of epigenetic marks will be essential to developing a robust DNA methylation biomarker assay, followed by meeting validation requirements for the assay before being implemented as a clinically useful tool.

Mitochondrial DNA as a marker for treatment-response in post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) is a serious mental health condition thought to be mediated by a dysregulated stress response system. Stress, especially chronic stress, affects mitochondrial activity and their efficiency in duplicating their genomes. Human cells contain numerous mitochondria that harbor multiple copies of their own genome, which consist of a mixture of wild type and variant mtDNA - a condition known as mitochondrial heteroplasmy. Number of mitochondrial genomes in a cell and the degree of heteroplasmy may serve as an indicator of mitochondrial allostatic load. Changes in mtDNA copy number and the proportion of variant mtDNA may be related to mental disorders and symptom severity, suggesting an involvement of mitochondrial dysfunction also in PTSD. Therefore, we examined number and composition of mitochondrial DNA before and after six weeks of inpatient psychotherapy treatment in a cohort of 60 female PTSD patients. We extracted DNA from isolated monocytes before and after inpatient treatment and quantified cellular mtDNA using multiplex qPCR. We hypothesized that treatment would lead to changes in cellular mtDNA levels and that change in mtDNA level would be associated with PTSD symptom severity and treatment response. It could be shown that mtDNA copy number and the ratio of variant mtDNA decreased during therapy, however, this change did not correlate with treatment response. Our results suggest that inpatient treatment can reduce signs of mitochondrial allostatic load, which could have beneficial effects on mental health. The quantification of mtDNA and the determination of cellular heteroplasmy could represent valuable biomarkers for the molecular characterization of mental disorders in the future.

Epigenetics and the Exposome: DNA Methylation as a Proxy for Health Impacts of Prenatal Environmental Exposures

The accumulation of every day exposures can impact health across the life course, but our understanding of such exposures is impeded by our ability to delineate the relationship between an individual's early life exposome and later life health effects. Measuring the exposome is challenging. Exposure assessed at a given time point captures a snapshot of the exposome but does not represent the full spectrum of exposures across the life course. In addition, the assessment of early life exposures and their effects is often further challenged by lack of relevant samples and the time gap between exposures and related health outcomes in later life. Epigenetics, specifically DNA methylation, has the potential to overcome these barriers as environmental epigenetic perturbances can be retained through time. In this review, we describe how DNA methylation can be framed in the world of the exposome. We offer three compelling examples of common environmental exposures, including cigarette smoke, the endocrine active compound bisphenol A (BPA), and the metal lead (Pb), to illustrate the application of DNA methylation as a proxy to measure the exposome. We discuss areas for future explorations and current limitations of this approach. Epigenetic profiling is a promising and rapidly developing tool and field of study, offering us a unique and powerful way to assess the early life exposome and its effects across different life stages.

No evidence for intervention-associated DNA methylation changes in monocytes of patients with posttraumatic stress disorder

DNA methylation patterns can be responsive to environmental influences. This observation has sparked interest in the potential for psychological interventions to influence epigenetic processes. Recent studies have observed correlations between DNA methylation changes and therapy outcome. However, most did not control for changes in cell composition. This study had two aims: first, we sought to replicate therapy-associated changes in DNA methylation of commonly assessed candidate genes in isolated monocytes from 60 female patients with post-traumatic stress disorder (PTSD). Our second, exploratory goal was to identify novel genomic regions with substantial pre-to-post intervention DNA methylation changes by performing whole-genome bisulfite sequencing (WGBS) in two patients with PTSD. Equivalence testing and Bayesian analyses provided evidence against physiologically meaningful intervention-associated DNA methylation changes in monocytes of PTSD patients in commonly investigated target genes (NR3C1, FKBP5, SLC6A4, OXTR). Furthermore, WGBS yielded only a limited set of candidate regions with suggestive evidence of differential DNA methylation pre- to post-therapy. These differential DNA methylation patterns did not prove replicable when investigated in the entire cohort. We conclude that there is no evidence for major, recurrent intervention-associated DNA methylation changes in the investigated genes in monocytes of patients with PTSD.

No long-term effects of antenatal synthetic glucocorticoid exposure on epigenetic regulation of stress-related genes

Antenatal synthetic glucocorticoid (sGC) treatment is a potent modifier of the hypothalamic-pituitary-adrenal (HPA) axis. In this context, epigenetic modifications are discussed as potential regulators explaining how prenatal exposure to GCs might translate into persistent changes of HPA axis "functioning". The purpose of this study was to investigate whether DNA methylation and gene expression profiles of stress-associated genes (NR3C1; FKBP5; SLC6A4) may mediate the persistent effects of sGC on cortisol stress reactivity that have been previously observed. In addition, hair cortisol concentrations (hairC) were investigated as a valid biomarker of long-term HPA axis activity. This cross-sectional study comprised 108 term-born children and adolescents, including individuals with antenatal GC treatment and controls. From whole blood, DNA methylation was analyzed by targeted deep bisulfite sequencing. Relative mRNA expression was determined by RT-qPCR experiments and qBase analysis. Acute stress reactivity was assessed by the Trier Social Stress Test (TSST) measuring salivary cortisol by ELISA and hairC concentrations were determined from hair samples by liquid chromatography coupled with tandem mass spectrometry. First, no differences in DNA methylation and mRNA expression levels of the stress-associated genes between individuals treated with antenatal sGC compared to controls were found. Second, DNA methylation and mRNA expression levels were neither associated with cortisol stress reactivity nor with hairC. These findings do not corroborate the belief that DNA methylation and mRNA expression profiles of stress-associated genes (NR3C1; FKBP5; SLC6A4) play a key mediating role of the persistent effects of sGC on HPA axis functioning.

Association between childhood maltreatment, psychopathology and DNA methylation of genes involved in stress regulation: Evidence from a study in Borderline Personality Disorder

Previous research suggests that childhood maltreatment is associated with epigenetic modification of genes involved in hypothalamic-pituitary-adrenal (HPA) functioning, which could cause dysregulation of the stress response system. If pervasive, this may be associated with the development of stress-related disorder in adults, including affective disorders, anxiety disorders, post-traumatic stress disorder (PTSD) or borderline-personality disorder (BPD). The majority of studies have focused on DNA methylation of the glucocorticoid receptor gene (NR3C1) and the FKBP5 encoding gene, which regulates the sensitivity of the glucocorticoid receptor (GR). How methylation of NR3C1 and FKBP5 interferes with childhood adversity and psychopathology as well as empathy is an under-researched issue. Here, we sought to investigate the association of childhood maltreatment in a sample of 89 individuals (44 healthy participants and 45 patients diagnosed with BPD) with the methylation of the 1_F promoter region of NR3C1 and the intron 7 of FKBP5 as well as with different measures of psychopathology and empathy. Methylation of FKBP5 (bin 2) correlated with anxiety (SCL-90-R) and the global psychopathological symptom load index (GSI), as well as with lower empathic perspective-taking abilities. Psychopathology and empathy impairments correlated with the level of childhood maltreatment. No difference in FKBP5 methylation was observed between the clinical and the non-clinical group. Methylation of NR3C1 was lower in BPD patients compared to controls, yet with small differences. The results are discussed regarding their biological relevance, including possible evolutionary explanations. In short, the regulation of the GR sensitivity by methylation of FKBP5 correlated with psychopathology and empathy scores, while no correlation emerged with the severity of childhood adversity.