Genome-wide DNA methylation investigation of glucocorticoid exposure within buccal samples

Epigenetic signals associated with delirium replicated across four independent cohorts

Delirium is risky and indicates poor outcomes for patients. Therefore, it is crucial to create an effective delirium detection method. However, the epigenetic pathophysiology of delirium remains largely unknown. We aimed to discover reliable and replicable epigenetic (DNA methylation: DNAm) markers that are associated with delirium including post-operative delirium (POD) in blood obtained from patients among four independent cohorts. Blood DNA from four independent cohorts (two inpatient cohorts and two surgery cohorts; 16 to 88 patients each) were analyzed using the Illumina EPIC array platform for genome-wide DNAm analysis. We examined DNAm differences in blood between patients with and without delirium including POD. When we compared top CpG sites previously identified from the initial inpatient cohort with three additional cohorts (one inpatient and two surgery cohorts), 11 of the top 13 CpG sites showed statistically significant differences in DNAm values between the delirium group and non-delirium group in the same directions as found in the initial cohort. This study demonstrated the potential value of epigenetic biomarkers as future diagnostic tools. Furthermore, our findings provide additional evidence of the potential role of epigenetics in the pathophysiology of delirium including POD.

Epigenetic associations in HPA axis genes related to bronchopulmonary dysplasia and antenatal steroids

BACKGROUND

Bronchopulmonary dysplasia (BPD), a common morbidity among very preterm infants, is associated with chronic disease and neurodevelopmental impairments. A hypothesized mechanism for these outcomes lies in altered glucocorticoid (GC) activity. We hypothesized that BPD and its treatments may result in epigenetic differences in the hypothalamic-pituitary-adrenal (HPA) axis, which is modulated by GC, and could be ascertained using an established GC risk score and DNA methylation (DNAm) of HPA axis genes.

METHODS

DNAm was quantified from buccal tissue (ECHO-NOVI) and from neonatal blood spots (ELGAN ECHO) via the EPIC microarray. Prenatal maternal characteristics, pregnancy complication, and neonatal medical complication data were collected from medical record review and maternal interviews.

RESULTS

The GC score was not associated with steroid exposure or BPD. However, six HPA genes involved in stress response regulation demonstrated differential methylation with antenatal steroid exposure; two CpGs within FKBP5 and POMC were differentially methylated with BPD severity. These findings were sex-specific in both cohorts; males had greater magnitude of differential methylation within these genes.

CONCLUSIONS

These findings suggest that BPD severity and antenatal steroids are associated with DNAm at some HPA genes in very preterm infants and the effects appear to be sex-, tissue-, and age-specific.

IMPACT

This study addresses bronchopulmonary dysplasia (BPD), an important health outcome among preterm neonates, and interrogates a commonly studied pathway, the hypothalamic-pituitary-adrenal (HPA) axis. The combination of BPD, the HPA axis, and epigenetic markers has not been previously reported. In this study, we found that BPD itself was not associated with epigenetic responses in the HPA axis in infants born very preterm; however, antenatal treatment with steroids was associated with epigenetic responses.

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

BACKGROUND

Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms.

METHODS

We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts.

FINDINGS

We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro.

INTERPRETATION

Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting.

FUNDING

This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.

Integrative multi-omics analysis of genomic, epigenomic, and metabolomics data leads to new insights for Attention-Deficit/Hyperactivity Disorder

The evolving field of multi-omics combines data and provides methods for simultaneous analysis across several omics levels. Here, we integrated genomics (transmitted and non-transmitted polygenic scores [PGSs]), epigenomics, and metabolomics data in a multi-omics framework to identify biomarkers for Attention-Deficit/Hyperactivity Disorder (ADHD) and investigated the connections among the three omics levels. We first trained single- and next multi-omics models to differentiate between cases and controls in 596 twins (cases = 14.8%) from the Netherlands Twin Register (NTR) demonstrating reasonable in-sample prediction through cross-validation. The multi-omics model selected 30 PGSs, 143 CpGs, and 90 metabolites. We confirmed previous associations of ADHD with glucocorticoid exposure and the transmembrane protein family TMEM, show that the DNA methylation of the MAD1L1 gene associated with ADHD has a relation with parental smoking behavior, and present novel findings including associations between indirect genetic effects and CpGs of the STAP2 gene. However, out-of-sample prediction in NTR participants (N = 258, cases = 14.3%) and in a clinical sample (N = 145, cases = 51%) did not perform well (range misclassification was [0.40, 0.57]). The results highlighted connections between omics levels, with the strongest connections between non-transmitted PGSs, CpGs, and amino acid levels and show that multi-omics designs considering interrelated omics levels can help unravel the complex biology underlying ADHD.

A multi-omics data analysis workflow packaged as a FAIR Digital Object

BACKGROUND

Applying good data management and FAIR (Findable, Accessible, Interoperable, and Reusable) data principles in research projects can help disentangle knowledge discovery, study result reproducibility, and data reuse in future studies. Based on the concepts of the original FAIR principles for research data, FAIR principles for research software were recently proposed. FAIR Digital Objects enable discovery and reuse of Research Objects, including computational workflows for both humans and machines. Practical examples can help promote the adoption of FAIR practices for computational workflows in the research community. We developed a multi-omics data analysis workflow implementing FAIR practices to share it as a FAIR Digital Object.

FINDINGS

We conducted a case study investigating shared patterns between multi-omics data and childhood externalizing behavior. The analysis workflow was implemented as a modular pipeline in the workflow manager Nextflow, including containers with software dependencies. We adhered to software development practices like version control, documentation, and licensing. Finally, the workflow was described with rich semantic metadata, packaged as a Research Object Crate, and shared via WorkflowHub.

CONCLUSIONS

Along with the packaged multi-omics data analysis workflow, we share our experiences adopting various FAIR practices and creating a FAIR Digital Object. We hope our experiences can help other researchers who develop omics data analysis workflows to turn FAIR principles into practice.

Aerobic exercise alters DNA hydroxymethylation levels in an experimental rodent model of temporal lobe epilepsy

The therapeutic potential of aerobic exercise in mitigating seizures and cognitive issues in temporal lobe epilepsy (TLE) is recognized, yet the underlying mechanisms are not well understood. Using a rodent TLE model induced by Kainic acid (KA), we investigated the impact of a single bout of exercise (i.e., acute) or 4 weeks of aerobic exercise (i.e., chronic). Blood was processed for epilepsy-associated serum markers, and DNA methylation (DNAme), and hippocampal area CA3 was assessed for gene expression levels for DNAme-associated enzymes. While acute aerobic exercise did not alter serum Brain-Derived Neurotrophic Factor (BDNF) or Interleukin-6 (IL-6), chronic exercise resulted in an exercise-specific decrease in serum BDNF and an increase in serum IL-6 levels in epileptic rats. Additionally, whole blood DNAme levels, specifically 5-hydroxymethylcytosine (5-hmC), decreased in epileptic animals following chronic exercise. Hippocampal CA3 5-hmC levels and ten-eleven translocation protein (TET1) expression mirrored these changes. Furthermore, immunohistochemistry analysis revealed that most 5-hmC changes in response to chronic exercise were neuron-specific within area CA3 of the hippocampus. Together, these findings suggest that DNAme mechanisms in the rodent model of TLE are responsive to chronic aerobic exercise, with emphasis on neuronal 5-hmC DNAme in the epileptic hippocampus.

Epigenetic differences in stress response gene FKBP5 among children with abusive vs accidental injuries

BACKGROUND

Survivors of child abuse experience high rates of adverse physical and mental health outcomes. Epigenetic alterations in the stress response system, the FKBP5 gene specifically, have been implicated as one mechanism that may link abuse to lifelong health issues. Prior studies primarily included older individuals with a remote history of maltreatment; our objective was to test for differential methylation of FKBP5 in children with abusive vs accidental injuries at the time of diagnosis.

METHODS

We conducted a cross-sectional pilot study of acutely injured children <4 years old at two children's hospitals (n = 82). Research personnel collected injury histories, buccal swabs (n = 65), and blood samples (n = 25) to measure DNA methylation. An expert panel classified the injuries as abusive, accidental, or indeterminate.

RESULTS

Children with abusive as compared to accidental injuries had lower methylation of the FKBP5 promoter in buccal and blood cells, even after controlling for injury severity, socioeconomic status, and psychosocial risk factors.

CONCLUSION

These findings suggest that epigenetic variation in FKBP5 may occur at the earliest indication of abuse and may be associated with delayed resolution of the HPA axis stress response. Additional testing for epigenetic differences in larger sample sizes is needed to further verify these findings.

IMPACT

Children (<4 years old) with abusive compared to accidental injuries showed lower methylation of the FKBP5 promoter in buccal and blood cells at the time of initial diagnosis even after controlling for injury severity, socioeconomic status, and psychosocial risk factors. Early childhood physical abuse may impact the epigenetic regulation of the stress response system, including demethylation within promoters and enhancers of the FKBP5 gene, even at the earliest indication of abuse. The findings are important because unmitigated stress is associated with adverse health outcomes throughout the life-course.

Integrative Multi-omics Analysis of Childhood Aggressive Behavior

This study introduces and illustrates the potential of an integrated multi-omics approach in investigating the underlying biology of complex traits such as childhood aggressive behavior. In 645 twins (cases = 42%), we trained single- and integrative multi-omics models to identify biomarkers for subclinical aggression and investigated the connections among these biomarkers. Our data comprised transmitted and two non-transmitted polygenic scores (PGSs) for 15 traits, 78,772 CpGs, and 90 metabolites. The single-omics models selected 31 PGSs, 1614 CpGs, and 90 metabolites, and the multi-omics model comprised 44 PGSs, 746 CpGs, and 90 metabolites. The predictive accuracy for these models in the test (N = 277, cases = 42%) and independent clinical data (N = 142, cases = 45%) ranged from 43 to 57%. We observed strong connections between DNA methylation, amino acids, and parental non-transmitted PGSs for ADHD, Autism Spectrum Disorder, intelligence, smoking initiation, and self-reported health. Aggression-related omics traits link to known and novel risk factors, including inflammation, carcinogens, and smoking.

The genome-wide DNA methylation profiles among neurosurgery patients with and without post-operative delirium

AIMS

There is no previous study demonstrating the differences of genome-wide DNA methylation (DNAm) profiles between patients with and without postoperative delirium (POD). We aimed to discover epigenetic (DNAm) markers that are associated with POD in blood obtained from patients before and after neurosurgery.

METHODS

Pre- and post-surgical blood DNA samples from 37 patients, including 10 POD cases, were analyzed using the Illumina EPIC array genome-wide platform. We examined DNAm differences in blood from patients with and without POD. Enrichment analysis with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes terms were also conducted.

RESULTS

When POD cases were tested for DNAm change before and after surgery, enrichment analyses showed many relevant signals with statistical significance in immune response related-pathways and inflammatory cytokine related-pathways such as "cellular response to cytokine stimulus", "regulation of immune system process", "regulation of cell activation", and "regulation of cytokine production". Furthermore, after excluding the potential effect of common factors related to surgery and anesthesia between POD cases and non-POD controls, the enrichment analyses showed significant signals such as "immune response" and "T cell activation", which are same pathways previously identified from an independent non-surgical inpatient cohort.

CONCLUSIONS

Our first genome-wide DNAm investigation of POD showed promising signals related to immune response, inflammatory response and other relevant signals considered to be associated with delirium pathophysiology. Our data supports the hypothesis that epigenetics play an important role in the pathophysiological mechanism of delirium and suggest the potential usefulness of an epigenetics-based biomarker of POD.

An explorative epigenome-wide association study of plasma renin and aldosterone concentration in a Ghanaian population: the RODAM study

BACKGROUND

The epigenetic regulation of the renin-angiotensin-aldosterone system (RAAS) potentially plays a role in the pathophysiology underlying the high burden of hypertension in sub-Saharan Africans (SSA). Here we report the first epigenome-wide association study (EWAS) of plasma renin and aldosterone concentrations and the aldosterone-to-renin ratio (ARR).

METHODS

Epigenome-wide DNA methylation was measured using the Illumina 450K array on whole blood samples of 68 Ghanaians. Differentially methylated positions (DMPs) were assessed for plasma renin concentration, aldosterone, and ARR using linear regression models adjusted for age, sex, body mass index, diabetes mellitus, hypertension, and technical covariates. Additionally, we extracted methylation loci previously associated with hypertension, kidney function, or that were annotated to RAAS-related genes and associated these with renin and aldosterone concentration.

RESULTS

We identified one DMP for renin, ten DMPs for aldosterone, and one DMP associated with ARR. Top DMPs were annotated to the PTPRN2, SKIL, and KCNT1 genes, which have been reported in relation to cardiometabolic risk factors, atherosclerosis, and sodium-potassium handling. Moreover, EWAS loci previously associated with hypertension, kidney function, or RAAS-related genes were also associated with renin, aldosterone, and ARR.

CONCLUSION

In this first EWAS on RAAS hormones, we identified DMPs associated with renin, aldosterone, and ARR in a SSA population. These findings are a first step in understanding the role of DNA methylation in regulation of the RAAS in general and in a SSA population specifically. Replication and translational studies are needed to establish the role of these DMPs in the hypertension burden in SSA populations.

Neural correlates of attentional control in social anxiety disorder: the impact of early-life adversity and DNA methylation

BACKGROUND

Social anxiety disorder is characterized by intense fear and avoidance of social interactions and scrutiny by others. Although alterations in attentional control seem to play a central role in the psychopathology of social anxiety disorder, the neural underpinnings in prefrontal brain regions have not yet been fully clarified.

METHODS

The present study used functional MRI in participants (age 18-50 yr) with social anxiety disorder (n = 42, 31 female) and without (n = 58, 33 female). It investigated the interrelation of the effects of social anxiety disorder and early-life adversity (a main environmental risk factor of social anxiety disorder) on brain activity during an attentional control task. We applied DNA methylation analysis to determine whether epigenetic modulation in the gene encoding the glucocorticoid receptor, NR3C1, might play a mediating role in this process.

RESULTS

We identified 2 brain regions in the left and medial prefrontal cortex that exhibited an interaction effect of social anxiety disorder and early-life adversity. In participants with low levels of early-life adversity, neural activity in response to disorder-related stimuli was increased in association with social anxiety disorder. In participants with high levels of early-life adversity, neural activity was increased only in participants without social anxiety disorder. NR3C1 DNA methylation partly mediated the effect of social anxiety disorder on brain activity as a function of early-life adversity.

LIMITATIONS

The absence of behavioural correlates associated with social anxiety disorder limited functional interpretation of the results.

CONCLUSION

These findings demonstrate that the neurobiological processes that underlie social anxiety disorder might be fundamentally different depending on experiences of early-life adversity. Long-lasting effects of early-life adversity might be encoded in NR3C1 DNA methylation and entail alterations in social anxiety disorder-related activity patterns in the neural network of attentional control.

DNA methylation in the inflammatory genes after neurosurgery and diagnostic ability of post-operative delirium

The pathophysiological mechanisms of postoperative delirium (POD) are still not clear, and no reliable biomarker is available to differentiate those with and without POD. Pre- and post-surgery blood from epilepsy subjects undergoing neurosurgery were collected. DNA methylation (DNAm) levels of the TNF gene, IL1B gene, and IL6 gene by the Illumina EPIC array method, and DNAm levels of the TNF gene by pyrosequencing, were analyzed. Blood from 37 subjects were analyzed by the EPIC array method, and blood from 27 subjects were analyzed by pyrosequencing. Several CpGs in the TNF gene in preoperative blood showed a negative correlation between their DNAm and age both in the POD group and in the non-POD group. However, these negative correlations were observed only in the POD group after neurosurgery. Neurosurgery significantly altered DNAm levels at 17 out of 24 CpG sites on the TNF gene, 8 out of 14 CpG sites on the IL1B gene, and 4 out of 14 CpG sites on the IL6 gene. Furthermore, it was found that the Inflammatory Methylation Index (IMI), which was based on the post-surgery DNAm levels at the selected five CpG sites, can be a potential detection tool for delirium with moderate accuracy; area under the curve (AUC) value was 0.84. The moderate accuracy of this IMI was replicated using another cohort from our previous study, in which the AUC was 0.79. Our findings provide further evidence of the potential role of epigenetics and inflammation in the pathophysiology of delirium.

Epigenetic alterations associated with dexamethasone sodium phosphate through DNMT and TET in RPE cells

PURPOSE

To elucidate the mechanism behind epigenetic alteration associated with dexamethasone (DEX) sodium phosphate treatment.

METHODS

We performed enzyme-linked immunosorbent assay to quantify changes in global DNA methylation and hydroxymethylation, quantitative real-time PCR (qRT-PCR) of the DNA methylation- and hydroxymethylation-related gene, in vitro DNA methyltransferase (DNMT) enzymatic activity assays with purified DNMTs, and DNA hydroxymethylation pattern with super-resolution imaging.

RESULTS

We identified global DNA hypomethylation and hyper-hydroxymethylation upon DEX treatment, associated with aberrant mRNA expression levels of DNMT and ten-eleven translocation (TET) proteins. Additionally, DEX exposure could directly hinder DNMT activities.

CONCLUSIONS

We showed that DEX-induced epigenetic alterations are linked to aberrant DNMT and TET expression, potentially through an essential role of DNMT.

Dynamic DNA Methylation Changes in the COMT Gene Promoter Region in Response to Mental Stress and Its Modulation by Transcranial Direct Current Stimulation

Changes in epigenetic modifications present a mechanism how environmental factors, such as the experience of stress, can alter gene regulation. While stress-related disorders have consistently been associated with differential DNA methylation, little is known about the time scale in which these alterations emerge. We investigated dynamic DNA methylation changes in whole blood of 42 healthy male individuals in response to a stressful cognitive task, its association with concentration changes in cortisol, and its modulation by transcranial direct current stimulation (tDCS). We observed a continuous increase in COMT promotor DNA methylation which correlated with higher saliva cortisol levels and was still detectable one week later. However, this lasting effect was suppressed by concurrent activity-enhancing anodal tDCS to the dorsolateral prefrontal cortex. Our findings support the significance of gene-specific DNA methylation in whole blood as potential biomarkers for stress-related effects. Moreover, they suggest alternative molecular mechanisms possibly involved in lasting behavioral effects of tDCS.

Crabp1 Modulates HPA Axis Homeostasis and Anxiety-like Behaviors by Altering FKBP5 Expression

Retinoic acid (RA), the principal active metabolite of vitamin A, is known to be involved in stress-related disorders. However, its mechanism of action in this regard remains unclear. This study reports that, in mice, endogenous cellular RA binding protein 1 (Crabp1) is highly expressed in the hypothalamus and pituitary glands. Crabp1 knockout (CKO) mice exhibit reduced anxiety-like behaviors accompanied by a lowered stress induced-corticosterone level. Furthermore, CRH/DEX tests show an increased sensitivity (hypersensitivity) of their feedback inhibition in the hypothalamic-pituitary-adrenal (HPA) axis. Gene expression studies show reduced FKBP5 expression in CKO mice; this would decrease the suppression of glucocorticoid receptor (GR) signaling thereby enhancing their feedback inhibition, consistent with their dampened corticosterone level and anxiety-like behaviors upon stress induction. In AtT20, a pituitary gland adenoma cell line elevating or reducing Crabp1 level correspondingly increases or decreases FKBP5 expression, and its endogenous Crabp1 level is elevated by GR agonist dexamethasone or RA treatment. This study shows, for the first time, that Crabp1 regulates feedback inhibition of the the HPA axis by modulating FKBP5 expression. Furthermore, RA and stress can increase Crabp1 level, which would up-regulate FKBP5 thereby de-sensitizing feedback inhibition of HPA axis (by decreasing GR signaling) and increasing the risk of stress-related disorders.

DNA methylation in the TNF-alpha gene decreases along with aging among delirium inpatients

It has been suggested that aging and inflammation play key roles in the development of delirium. In the present study, we investigated the differences of the DNAm patterns in the TNF gene between patients with delirium and without. The data and samples derived from previous and ongoing cohort studies were analyzed. DNAm levels of the TNF gene were analyzed using the Illumina EPIC array genome-wide method and pyrosequencing method. Correlations between age and DNAm levels of each CpG were calculated. Several CpG in the TNF gene in blood showed negative correlation between their DNAm and age in delirium cases both with the EPIC array and by the pyrosequencing method. However, there was no CpG that had significant correlation between their DNAm and age regardless of delirium status among buccal samples. On the other hand, among peripheral blood mononuclear cells samples, it was found that several CpG showed negative correlation between their DNAm and age in delirium cases. The evidence of DNAm change in the TNF gene among delirious subjects was demonstrated.

Bound Together: How Psychoanalysis Diminishes Inter-generational DNA Trauma

The concept of intergenerational transmission of trauma plays a fundamental role in psychoanalysis. While it is known that intergenerational trauma can be transmitted through attachment relationships, a new branch of genetics (epigenetics) has emerged to study the interaction between human behavior and changes in DNA expression. Therefore, psychoanalysis, which has proven to reduce the intergenerational transmission of trauma from a behavioral perspective, can play a positive role in regulating DNA changes caused by environmental stress. The present paper focuses on recent research suggesting a direct correlation between psychological trauma and DNA modifications. In particular, DNA changes caused by psychological trauma can be transmitted from generation to generation, validating the psychoanalytic concept of intergenerational transmission of trauma. This evidence not only supports the essential role psychoanalysis has in influencing human behavior, but also suggests that it affects not only the individuals who undergo it but their offspring, as well, via the epigenetic passage of DNA.

Differential DNA methylation in high-grade serous ovarian cancer (HGSOC) is associated with tumor behavior

The epigenome offers an additional facet of cancer that can help categorize patients into those at risk of disease, recurrence, or treatment failure. We conducted a retrospective, nested, case-control study of advanced and recurrent high-grade serous ovarian cancer (HGSOC) patients in which we assessed epigenome-wide association using Illumina methylationEPIC arrays to characterize DNA methylation status and RNAseq to evaluate gene expression. Comparing HGSOC tumors with normal fallopian tube tissues we observe global hypomethylation but with skewing towards hypermethylation when interrogating gene promoters. In total, 5,852 gene interrogating probes revealed significantly different methylation. Within HGSOC, 57 probes highlighting 17 genes displayed significant differential DNA methylation between primary and recurrent disease. Between optimal vs suboptimal surgical outcomes 99 probes displayed significantly different methylation but only 29 genes showed an inverse correlation between methylation status and gene expression. Overall, differentially methylated genes point to several pathways including RAS as well as hippo signaling in normal vs primary HGSOC; valine, leucine, and isoleucine degradation and endocytosis in primary vs recurrent HGSOC; and pathways containing immune driver genes in optimal vs suboptimal surgical outcomes. Thus, differential DNA methylation identified numerous genes that could serve as potential biomarkers and/or therapeutic targets in HGSOC.

DNA methylation is associated with inhaled corticosteroid response in persistent childhood asthmatics

BACKGROUND

Response to inhaled corticosteroids is highly variable, and the association between DNA methylation and treatment response is not known.

OBJECTIVE

To examine the association between peripheral blood DNA methylation and inhaled corticosteroid response in children with persistent asthma.

METHODS

Epigenome-wide DNA methylation was analysed in individuals on inhaled corticosteroids in three independent and ethnically diverse cohorts-Childhood Asthma Management Program (CAMP); Children, Allergy, Milieu, Stockholm, Epidemiology (BAMSE); and Genetic Epidemiology of Asthma in Costa Rica Study (GACRS). Treatment response was evaluated using two definitions, the absence of emergency department visits and/or hospitalizations and the absence oral corticosteroid use while on inhaled corticosteroid therapy. CpG sites meeting nominal significance (P < 0.05) for each outcome were combined in a three-cohort meta-analysis with adjustment for multiple testing. DNA methylation was correlated with gene expression using Pearson and partial correlations.

RESULTS

In 154 subjects from CAMP, 72 from BAMSE, and 168 from GACRS, relative hypomethylation of cg00066816 (171 bases upstream of IL12B) was associated with the absence of emergency department visits and/or hospitalizations (Q = 0.03) in all cohorts and lower IL12B expression (ρ = 0.34, P = 0.01) in BAMSE. Relative hypermethylation of cg04256470 (688 bases upstream of CORT) was associated with the absence of oral corticosteroid use (Q = 0.04) in all cohorts and higher CORT expression (ρ = 0.20, P = 0.045) in CAMP.

CONCLUSION AND CLINICAL RELEVANCE

Differential DNA methylation of IL12B and CORT are associated with inhaled corticosteroid treatment response in persistent childhood asthmatics. Pharmaco-methylation can identify novel markers of treatment sensitivity in asthma.