pwrEWAS: a user-friendly tool for comprehensive power estimation for epigenome wide association studies (EWAS)

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC v1.0 BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC v1.0 arrays. We conducted a comprehensive assessment of the EPIC v1.0 array probe reliability using 69 blood DNA samples, each measured twice, generated by the Alzheimer's Disease Neuroimaging Initiative study. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliability information for probes on the EPIC v1.0 array, will serve as a valuable resource for future DNAm studies.

Genetic variation and pesticide exposure influence blood DNA methylation signatures in females with early-stage Parkinson's disease

Although sex, genetics, and exposures can individually influence risk for sporadic Parkinson's disease (PD), the joint contributions of these factors to the epigenetic etiology of PD have not been comprehensively assessed. Here, we profiled sex-stratified genome-wide blood DNAm patterns, SNP genotype, and pesticide exposure in agricultural workers (71 early-stage PD cases, 147 controls) and explored replication in three independent samples of varying demographics (n = 218, 222, and 872). Using a region-based approach, we found more associations of blood DNAm with PD in females (69 regions) than in males (2 regions, Δβadj| ≥0.03, padj ≤ 0.05). For 48 regions in females, models including genotype or genotype and pesticide exposure substantially improved in explaining interindividual variation in DNAm (padj ≤ 0.05), and accounting for these variables decreased the estimated effect of PD on DNAm. The results suggested that genotype, and to a lesser degree, genotype-exposure interactions contributed to variation in PD-associated DNAm. Our findings should be further explored in larger study populations and in experimental systems, preferably with precise measures of exposure.

The race-based stress reduction intervention (RiSE) study on African American women in NYC and Chicago: Design and methods for complex genomic analysis

RiSE study aims to evaluate a race-based stress-reduction intervention as an effective strategy to improve coping and decrease stress-related symptoms, inflammatory burden, and modify DNA methylation of stress response-related genes in older AA women. This article will describe genomic analytic methods to be utilized in this longitudinal, randomized clinical trial of older adult AA women in Chicago and NYC that examines the effect of the RiSE intervention on DNAm pre- and post-intervention, and its overall influence on inflammatory burden. Salivary DNAm will be measured at baseline and 6 months following the intervention, using the Oragene-DNA kit. Measures of perceived stress, depressive symptoms, fatigue, sleep, inflammatory burden, and coping strategies will be assessed at 4 time points including at baseline, 4 weeks, 8 weeks, and 6 months. Genomic data analysis will include the use of pre-processed and quality-controlled methylation data expressed as beta (β) values. Association analyses will be performed to detect differentially methylated sites on the targeted candidate genes between the intervention and non-intervention groups using the Δβ (changes in methylation) with adjustment for age, health behaviors, early life adversity, hybridization batch, and top principal components of the probes as covariates. To account for multiple testing, we will use FDR adjustment with a corrected p-value of <0.05 regarded as statistically significant. To assess the relationship between inflammatory burden and Δβ among the study samples, we will repeat association analyses with the inclusion of individual inflammation protein measures. ANCOVA will be used because it is more statistically powerful to detect differences.

Cord blood DNA methylation signatures associated with preeclampsia are enriched for cardiovascular pathways: insights from the VDAART trial

BACKGROUND

Preeclampsia has been associated with maternal epigenetic changes, in particular DNA methylation changes in the placenta. It has been suggested that preeclampsia could also cause DNA methylation changes in the neonate. We examined DNA methylation in relation to gene expression in the cord blood of offspring born to mothers with preeclampsia.

METHODS

This study included 128 mother-child pairs who participated in the Vitamin D Antenatal Asthma Reduction Trial (VDAART), where assessment of preeclampsia served as secondary outcome. We performed an epigenome-wide association study of preeclampsia and cord blood DNA methylation (Illumina 450 K chip). We then examined gene expression of the same subjects for validation and replicated the gene signatures in independent DNA methylation datasets. Lastly, we applied functional enrichment and network analyses to identify biological pathways that could potentially be involved in preeclampsia.

FINDINGS

In the cord blood samples (n = 128), 263 CpGs were differentially methylated (FDR <0.10) in preeclampsia (n = 16), of which 217 were annotated. Top pathways in the functional enrichment analysis included apelin signaling pathway and other endothelial and cardiovascular pathways. Of the 217 genes, 13 showed differential expression (p's < 0.001) in preeclampsia and 11 had been previously related to preeclampsia (p's < 0.0001). These genes were linked to apelin, cGMP and Notch signaling pathways, all having a role in angiogenic process and cardiovascular function.

INTERPRETATION

Preeclampsia is related to differential cord blood DNA methylation signatures of cardiovascular pathways, including the apelin signaling pathway. The association of these cord blood DNA methylation signatures with offspring's long-term morbidities due to preeclampsia should be further investigated.

FUNDING

VDAART is funded by National Heart, Lung, and Blood Institute grants of R01HL091528 and UH3OD023268. HMK is supported by Jane and Aatos Erkko Foundation, Paulo Foundation, and the Pediatric Research Foundation. HM is supported by K01 award from NHLBI (1K01HL146977-01A1). PK is supported by K99HL159234 from NIH/NHLBI.

Epigenome-wide association study of peripheral immune cell populations in Parkinson's disease

Understanding the contribution of immune mechanisms to Parkinson's disease pathogenesis is an important challenge, potentially of major therapeutic implications. To further elucidate the involvement of peripheral immune cells, we studied epigenome-wide DNA methylation in isolated populations of CD14+ monocytes, CD19+ B cells, CD4+ T cells, and CD8+ T cells from Parkinson's disease patients and healthy control participants. We included 25 patients with a maximum five years of disease duration and 25 controls, and isolated four immune cell populations from each fresh blood sample. Epigenome-wide DNA methylation profiles were generated from 186 samples using the Illumina MethylationEpic array and association with disease status was tested using linear regression models. We identified six differentially methylated CpGs in CD14+ monocytes and one in CD8 + T cells. Four differentially methylated regions were identified in monocytes, including a region upstream of RAB32, a gene that has been linked to LRRK2. Methylation upstream of RAB32 correlated negatively with mRNA expression, and RAB32 expression was upregulated in Parkinson's disease both in our samples and in summary statistics from a previous study. Our epigenome-wide association study of early Parkinson's disease provides evidence for methylation changes across different peripheral immune cell types, highlighting monocytes and the RAB32 locus. The findings were predominantly cell-type-specific, demonstrating the value of isolating purified cell populations for genomic studies.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Linking environmental risk factors with epigenetic mechanisms in Parkinson's disease

Sporadic Parkinson's disease (PD) is a progressive neurodegenerative disease, with a complex risk structure thought to be influenced by interactions between genetic variants and environmental exposures, although the full aetiology is unknown. Environmental factors, including pesticides, have been reported to increase the risk of developing the disease. Growing evidence suggests epigenetic changes are key mechanisms by which these environmental factors act upon gene regulation, in disease-relevant cell types. We present a systematic review critically appraising and summarising the current body of evidence of the relationship between epigenetic mechanisms and environmental risk factors in PD to inform future research in this area. Epigenetic studies of relevant environmental risk factors in animal and cell models have yielded promising results, however, research in humans is just emerging. While published studies in humans are currently relatively limited, the importance of the field for the elucidation of molecular mechanisms of pathogenesis opens clear and promising avenues for the future of PD research. Carefully designed epidemiological studies carried out in PD patients hold great potential to uncover disease-relevant gene regulatory mechanisms. Therefore, to advance this burgeoning field, we recommend broadening the scope of investigations to include more environmental exposures, increasing sample sizes, focusing on disease-relevant cell types, and recruiting more diverse cohorts.

Epigenome-wide association study of diabetic chronic kidney disease progression in the Korean population: the KNOW-CKD study

Since the etiology of diabetic chronic kidney disease (CKD) is multifactorial, studies on DNA methylation for kidney function deterioration have rarely been performed despite the need for an epigenetic approach. Therefore, this study aimed to identify epigenetic markers associated with CKD progression based on the decline in the estimated glomerular filtration rate in diabetic CKD in Korea. An epigenome-wide association study was performed using whole blood samples from 180 CKD recruited from the KNOW-CKD cohort. Pyrosequencing was also performed on 133 CKD participants as an external replication analysis. Functional analyses, including the analysis of disease-gene networks, reactome pathways, and protein-protein interaction networks, were conducted to identify the biological mechanisms of CpG sites. A phenome-wide association study was performed to determine the associations between CpG sites and other phenotypes. Two epigenetic markers, cg10297223 on AGTR1 and cg02990553 on KRT28 indicated a potential association with diabetic CKD progression. Based on the functional analyses, other phenotypes (blood pressure and cardiac arrhythmia for AGTR1) and biological pathways (keratinization and cornified envelope for KRT28) related to CKD were also identified. This study suggests a potential association between the cg10297223 and cg02990553 and the progression of diabetic CKD in Koreans. Nevertheless, further validation is needed through additional studies.

Effects of prenatal exposure to (es)citalopram and maternal depression during pregnancy on DNA methylation and child neurodevelopment

Studies assessing associations between prenatal exposure to antidepressants, maternal depression, and offspring DNA methylation (DNAm) have been inconsistent. Here, we investigated whether prenatal exposure to citalopram or escitalopram ((es)citalopram) and maternal depression is associated with differences in DNAm. Then, we examined if there is an interaction effect of (es)citalopram exposure and DNAm on offspring neurodevelopmental outcomes. Finally, we investigated whether DNAm at birth correlates with neurodevelopmental trajectories in childhood. We analyzed DNAm in cord blood from the Norwegian Mother, Father and Child Cohort Study (MoBa) biobank. MoBa contains questionnaire data on maternal (es)citalopram use and depression during pregnancy and information about child neurodevelopmental outcomes assessed by internationally recognized psychometric tests. In addition, we retrieved ADHD diagnoses from the Norwegian Patient Registry and information on pregnancies from the Medical Birth Registry of Norway. In total, 958 newborn cord blood samples were divided into three groups: (1) prenatal (es)citalopram exposed (n = 306), (2) prenatal maternal depression exposed (n = 308), and (3) propensity score-selected controls (n = 344). Among children exposed to (es)citalopram, there were more ADHD diagnoses and symptoms and delayed communication and psychomotor development. We did not identify differential DNAm associated with (es)citalopram or depression, nor any interaction effects on neurodevelopmental outcomes throughout childhood. Trajectory modeling identified subgroups of children following similar developmental patterns. Some of these subgroups were enriched for children exposed to maternal depression, and some subgroups were associated with differences in DNAm at birth. Interestingly, several of the differentially methylated genes are involved in neuronal processes and development. These results suggest DNAm as a potential predictive molecular marker of later abnormal neurodevelopmental outcomes, but we cannot conclude whether DNAm links prenatal (es)citalopram exposure or maternal depression with child neurodevelopmental outcomes.

DNA Methylation Analysis in Monozygotic Twins Discordant for ALS in Blood Cells

ALS is a fatal motor neuron disease that displays a broad variety of phenotypes ranging from early fatal courses to slowly progressing and rather benign courses. Such divergence can also be seen in genetic ALS cases with varying phenotypes bearing specific mutations, suggesting epigenetic mechanisms like DNA methylation act as disease modifiers. However, the epigenotype dictated by, in addition to other mechanisms, DNA methylation is also strongly influenced by the individual's genotype. Hence, we performed a DNA methylation study using EPIC arrays on 7 monozygotic (MZ) twin pairs discordant for ALS in whole blood, which serves as an ideal model for eliminating the effects of the genetic-epigenetic interplay to a large extent. We found one CpG site showing intra-pair hypermethylation in the affected co-twins, which maps to the Glutamate Ionotropic Receptor Kainate Type Subunit 1 gene (GRIK1). Additionally, we found 4 DMPs which were subsequently confirmed using 2 different statistical approaches. Differentially methylated regions or blocks could not be detected within the scope of this work. In conclusion, we revealed that despite a low sample size, monozygotic twin studies discordant for the disease can bring new insights into epigenetic processes in ALS, pointing to new target loci for further investigations.

DNA-methylation and immunological response in medication overuse headache

OBJECTIVE

To investigate whether medication-overuse headache patients have differential DNA-methylation pattern.

METHODS

We collected blood samples from 120 medication-overuse headache-patients, 57 controls (29 episodic migraine patients and 28 healthy controls) in a hypothesis-generating cross-sectional case-control pilot study; 100 of the medication-overuse headache-patients were followed for six months and samples were collected at two and six months for the longitudinal methylation analyses. Blood cell proportions of leucocytes (neutrophils, NK-cells, monocytes, CD8+ and CD4+ T-cells, and B-cells) and the neutrophile-lymphocyte ratio were estimated using methylation data as a measure for immunological analysis and a cell type-specific epigenome wide association study was conducted between medication-overuse headache-patients and controls, and longitudinally for reduction in headache days/month among medication-overuse headache-patients.

RESULTS

We found a higher neutrophile-lymphocyte ratio in medication-overuse headache-patients compared to controls, indicating a higher immunological response in medication-overuse headache-patients (false discovery rate (adjusted p-value)<0.001). Reduction in headache days/month (9.8; 95% CI 8.1-11.5) was associated with lower neutrophile-lymphocyte ratio (false discovery rate adjusted p-value = 0.041).Three genes (CORIN, CCKBR and CLDN9) were hypermethylated in specific cell types in medication-overuse headache-patients compared to controls. No methylation differences were associated with reduction in headache days in medication-overuse headache-patients after six months.

CONCLUSION

This pilot study was consistent with higher immunological response in medication-overuse headache-patients which decreased with a reduction in headache days in longitudinal analysis. medication-overuse headache-patients exhibited differential methylation in innate immune cells but did not exhibit longitudinal differences with alterations in headache days. Our study creates hypotheses for further biomarker searches.ClinicalTrials.gov Identifier: NCT02993289.

Distinct genome-wide DNA methylation and gene expression signatures in classical monocytes from African American patients with systemic sclerosis

BACKGROUND

Systemic sclerosis (SSc) is a multisystem autoimmune disorder that has an unclear etiology and disproportionately affects women and African Americans. Despite this, African Americans are dramatically underrepresented in SSc research. Additionally, monocytes show heightened activation in SSc and in African Americans relative to European Americans. In this study, we sought to investigate DNA methylation and gene expression patterns in classical monocytes in a health disparity population.

METHODS

Classical monocytes (CD14+ + CD16-) were FACS-isolated from 34 self-reported African American women. Samples from 12 SSc patients and 12 healthy controls were hybridized on MethylationEPIC BeadChip array, while RNA-seq was performed on 16 SSc patients and 18 healthy controls. Analyses were computed to identify differentially methylated CpGs (DMCs), differentially expressed genes (DEGs), and CpGs associated with changes in gene expression (eQTM analysis).

RESULTS

We observed modest DNA methylation and gene expression differences between cases and controls. The genes harboring the top DMCs, the top DEGs, as well as the top eQTM loci were enriched for metabolic processes. Genes involved in immune processes and pathways showed a weak upregulation in the transcriptomic analysis. While many genes were newly identified, several other have been previously reported as differentially methylated or expressed in different blood cells from patients with SSc, supporting for their potential dysregulation in SSc.

CONCLUSIONS

While contrasting with results found in other blood cell types in largely European-descent groups, the results of this study support that variation in DNA methylation and gene expression exists among different cell types and individuals of different genetic, clinical, social, and environmental backgrounds. This finding supports the importance of including diverse, well-characterized patients to understand the different roles of DNA methylation and gene expression variability in the dysregulation of classical monocytes in diverse populations, which might help explaining the health disparities.

Associations of peripheral blood DNA methylation and estimated monocyte proportion differences during infancy with toddler attachment style

Attachment is a motivational system promoting felt security to a caregiver resulting in a persistent internal working model of interpersonal behavior. Attachment styles are developed in early social environments and predict future health and development outcomes with potential biological signatures, such as epigenetic modifications like DNA methylation (DNAm). Thus, we hypothesized infant DNAm would associate with toddler attachment styles. An epigenome-wide association study (EWAS) of blood DNAm from 3-month-old infants was regressed onto children's attachment style from the Strange Situation Procedure at 22-months at multiple DNAm Cytosine-phosphate-Guanine (CpG) sites. The 26 identified CpGs associated with proinflammatory immune phenotypes and cognitive development. In post-hoc analyses, only maternal cognitive-growth fostering, encouraging intellectual exploration, contributed. For disorganized children, DNAm-derived cell-type proportions estimated higher monocytes -cells in immune responses hypothesized to increase with early adversity. Collectively, these findings suggested the potential biological embedding of both adverse and advantageous social environments as early as 3-months-old.

Dementia with Lewy bodies post-mortem brains reveal differentially methylated CpG sites with biomarker potential

Dementia with Lewy bodies (DLB) is a common form of dementia with known genetic and environmental interactions. However, the underlying epigenetic mechanisms which reflect these gene-environment interactions are poorly studied. Herein, we measure genome-wide DNA methylation profiles of post-mortem brain tissue (Broadmann area 7) from 15 pathologically confirmed DLB brains and compare them with 16 cognitively normal controls using Illumina MethylationEPIC arrays. We identify 17 significantly differentially methylated CpGs (DMCs) and 17 differentially methylated regions (DMRs) between the groups. The DMCs are mainly located at the CpG islands, promoter and first exon regions. Genes associated with the DMCs are linked to "Parkinson's disease" and "metabolic pathway", as well as the diseases of "severe intellectual disability" and "mood disorders". Overall, our study highlights previously unreported DMCs offering insights into DLB pathogenesis with the possibility that some of these could be used as biomarkers of DLB in the future.

Designing Studies for Epigenetic Biomarker Development in Autoimmune Rheumatic Diseases

In just a few years, the number of epigenetic studies in autoimmune rheumatic and inflammatory diseases has greatly increased. This is in part due to the need of identifying additional determinants to genetics to explain the pathogenesis and development of these disorders. In this regard, epigenetics provides potential mechanisms that determine gene function, are linked to environmental factors, and could explain a wide range of phenotypic variability among patients with these diseases. Despite the high interest and number of studies describing epigenetic alterations under these conditions and exploring their relationship to various clinical aspects, few of the proposed biomarkers have yet reached clinical practice. The potential of epigenetic markers is high, as these alterations link measurable features with a number of biological traits. In the present article, we present published studies in the field, discuss some frequent limitations in the existing research, and propose a number of considerations that should be taken into account by those starting new projects in the field, with an aim to generate biomarkers that could make it into the clinics.

DNA Methylation and Ischemic Stroke Risk: An Epigenome-Wide Association Study

BACKGROUND

Ischemic stroke (IS) risk heritability is partly explained by genetics. Other heritable factors, such as epigenetics, could explain an unknown proportion of the IS risk. The objective of this study is to evaluate DNA methylation association with IS using epigenome-wide association studies (EWAS).

METHODS

 We performed a two-stage EWAS comprising 1,156 subjects. Differentially methylated positions (DMPs) and differentially methylated regions (DMRs) were assessed using the Infinium 450K and EPIC BeadChip in the discovery cohort (252 IS and 43 controls). Significant DMPs were replicated in an independent cohort (618 IS and 243 controls). Stroke subtype associations were also evaluated. Differentially methylated cell-type (DMCT) was analyzed in the replicated CpG sites using EpiDISH. We additionally performed pathway enrichment analysis and causality analysis with Mendelian randomization for the replicated CpG sites.

RESULTS

 A total of 957 CpG sites were epigenome-wide-significant (p ≤ 10^-7) in the discovery cohort, being CpG sites in the top signals (logFC = 0.058, p = 2.35 × 10^-22; logFC = 0.035, p = 3.22 × 10^-22, respectively). ZFHX3 and MAP3K1 were among the most significant DMRs. In addition, 697 CpG sites were replicated considering Bonferroni-corrected p-values (p < 5.22 × 10^-5). All the replicated DMPs were associated with risk of cardioembolic, atherothrombotic, and undetermined stroke. The DMCT analysis demonstrated that the significant associations were driven by natural killer cells. The pathway enrichment analysis showed overrepresentation of genes belonging to certain pathways including oxidative stress. ZFHX3 and MAP3K1 methylation was causally associated with specific stroke-subtype risk.

CONCLUSION

 Specific DNA methylation pattern is causally associated with IS risk. These results could be useful for specifically predicting stroke occurrence and could potentially be evaluated as therapeutic targets.

Social Factors, Epigenomics and Lupus in African American Women (SELA) Study: protocol for an observational mechanistic study examining the interplay of multiple individual and social factors on lupus outcomes in a health disparity population

Introduction

Despite the disproportional impact of SLE on historically marginalised communities, the individual and sociocultural factors underlying these health disparities remain elusive. We report the design and methods for a study aimed at identifying epigenetic biomarkers associated with racism and resiliency that affect gene function and thereby influence SLE in a health disparity population.

Methods and analysis

The Social Factors, Epigenomics and Lupus in African American Women (SELA) Study is a cross-sectional, case–control study. A total of 600 self-reported African American women will be invited to participate. All participants will respond to questionnaires that capture detailed sociodemographic and medical history, validated measures of racial discrimination, social support, as well as disease activity and damage for cases. Participants who wish will receive their genetic ancestry estimates and be involved in research. Blood samples are required to provide peripheral blood mononuclear cell counts, DNA and RNA. The primary goals of SELA are to identify variation in DNA methylation (DNAm) associated with self-reported exposure to racial discrimination and social support, to evaluate whether social DNAm sites affect gene expression, to identify the synergistic effects of social factors on DNAm changes on SLE and to develop a social factors-DNAm predictive model for disease outcomes. This study is conducted in cooperation with the Sea Island Families Project Citizen Advisory Committee.

Discussion and dissemination

SELA will respond to the pressing need to clarify the interplay and regulatory mechanism by which various positive and negative social exposures influence SLE. Results will be published and shared with patients and the community. Knowledge of the biological impact of social exposures on SLE, as informed by the results of this study, can be leveraged by advocacy efforts to develop psychosocial interventions that prevent or mitigate risk exposures, and services or interventions that promote positive exposures. Implementation of such interventions is paramount to the closure of the health disparities gap.

Data Analysis of DNA Methylation Epigenome-Wide Association Studies (EWAS): A Guide to the Principles of Best Practice

Array-based EWAS have become an increasingly popular technique to identify population epigenetic effects, particularly in humans. With the arrival of nonhuman species arrays, such as the mouse, this is likely to become an even more widely used technology. This chapter provides the less experienced researcher a guide to the analysis of data from the most widely used platform, the Illumina Infinium Methylation assay. This includes an overview of quality filtering, data normalization, analysis options, and techniques to improve the interpretation of results.

pwrBRIDGE: a user-friendly web application for power and sample size estimation in batch-confounded microarray studies with dependent samples

Batch effect Reduction of mIcroarray data with Dependent samples usinG Empirical Bayes (BRIDGE) is a recently developed statistical method to address the issue of batch effect correction in batch-confounded microarray studies with dependent samples. The key component of the BRIDGE methodology is the use of samples run as technical replicates in two or more batches, "bridging samples", to inform batch effect correction/attenuation. While previously published results indicate a relationship between the number of bridging samples, M, and the statistical power of downstream statistical testing on the batch-corrected data, there is of yet no formal statistical framework or user-friendly software, for estimating M to achieve a specific statistical power for hypothesis tests conducted on the batch-corrected data. To fill this gap, we developed pwrBRIDGE, a simulation-based approach to estimate the bridging sample size, M, in batch-confounded longitudinal microarray studies. To illustrate the use of pwrBRIDGE, we consider a hypothetical, longitudinal batch-confounded study whose goal is to identify Alzheimer's disease (AD) progression-associated genes from amnestic mild cognitive impairment (aMCI) to AD in human blood after a 5-year follow-up. pwrBRIDGE helps researchers design and plan batch-confounded microarray studies with dependent samples to avoid over- or under-powered studies.

An investigation into DNA methylation patterns associated with risk preference in older individuals

Risk preference is a complex trait governed by psycho-social, environmental and genetic determinants. We aimed to examine how an individual’s risk preference associates with their epigenetic profile.

Risk preferences were ascertained by asking participants of the Northern Ireland COhort for the Longitudinal study of Ageing to make a series of choices between hypothetical income scenarios. From these, four risk preference categories were derived, ranging from risk-averse to risk-seeking. Illumina’s Infinium High-Density Methylation Assay was used to evaluate the status of 862,927 CpGs.

Risk preference and DNA methylation data were obtained for 1,656 individuals. The distribution of single-site DNA methylation levels between risk-averse and risk-seeking individuals was assessed whilst adjusting for age, sex and peripheral white cell counts. In this discovery cohort, 55 CpGs were identified with significantly different levels of methylation (p≤x10⁻⁵) between risk-averse and risk-seeking individuals when adjusting for the maximum number of covariates. No CpGs were significantly differentially methylated in any of the risk preference groups at an epigenome-wide association level (p<9x10⁻⁸) following covariate adjustment.

Protein-coding genes NWD1 and LRP1 were among the genes in which the top-ranked dmCpGs were located for all analyses conducted. Mutations in these genes have previously been linked to neurological conditions.

Epigenetic modifications have not previously been linked to risk-aversion using a population cohort, but may represent important biomarkers of accumulated, complex determinants of this trait. Several striking results from this study support further analysis of DNA methylation as an important link between measurable biomarkers and health behaviours.

Using Monozygotic Twins to Dissect Common Genes in Posttraumatic Stress Disorder and Migraine

Epigenetic mechanisms have been associated with genes involved in Posttraumatic stress disorder (PTSD). PTSD often co-occurs with other health conditions such as depression, cardiovascular disorder and respiratory illnesses. PTSD and migraine have previously been reported to be symptomatically positively correlated with each other, but little is known about the genes involved. The aim of this study was to understand the comorbidity between PTSD and migraine using a monozygotic twin disease discordant study design in six pairs of monozygotic twins discordant for PTSD and 15 pairs of monozygotic twins discordant for migraine. DNA from peripheral blood was run on Illumina EPIC arrays and analyzed. Multiple testing correction was performed using the Bonferroni method and 10% false discovery rate (FDR). We validated 11 candidate genes previously associated with PTSD including DOCK2, DICER1, and ADCYAP1. In the epigenome-wide scan, seven novel CpGs were significantly associated with PTSD within/near IL37, WNT3, ADNP2, HTT, SLFN11, and NQO2, with all CpGs except the IL37 CpG hypermethylated in PTSD. These results were significantly enriched for genes whose DNA methylation was previously associated with migraine (p-value = 0.036). At 10% FDR, 132 CpGs in 99 genes associated with PTSD were also associated with migraine in the migraine twin samples. Genes associated with PTSD were overrepresented in vascular smooth muscle, axon guidance and oxytocin signaling pathways, while genes associated with both PTSD and migraine were enriched for AMPK signaling and longevity regulating pathways. In conclusion, these results suggest that common genes and pathways are likely involved in PTSD and migraine, explaining at least in part the co-morbidity between the two disorders.

DNA methylation in human sperm: a systematic review

BACKGROUND

Studies in non-human mammals suggest that environmental factors can influence spermatozoal DNA methylation, and some research suggests that spermatozoal DNA methylation is also implicated in conditions such as subfertility and imprinting disorders in the offspring. Together with an increased availability of cost-effective methods of interrogating DNA methylation, this premise has led to an increasing number of studies investigating the DNA methylation landscape of human spermatozoa. However, how the human spermatozoal DNA methylome is influenced by environmental factors is still unclear, as is the role of human spermatozoal DNA methylation in subfertility and in influencing offspring health.

OBJECTIVE AND RATIONALE

The aim of this systematic review was to critically appraise the quality of the current body of literature on DNA methylation in human spermatozoa, summarize current knowledge and generate recommendations for future research.

SEARCH METHODS

A comprehensive literature search of the PubMed, Web of Science and Cochrane Library databases was conducted using the search terms 'semen' OR 'sperm' AND 'DNA methylation'. Publications from 1 January 2003 to 2 March 2020 that studied human sperm and were written in English were included. Studies that used sperm DNA methylation to develop methodologies or forensically identify semen were excluded, as were reviews, commentaries, meta-analyses or editorial texts. The Grading of Recommendations, Assessment, Development and Evaluations (GRADE) criteria were used to objectively evaluate quality of evidence in each included publication.

OUTCOMES

The search identified 446 records, of which 135 were included in the systematic review. These 135 studies were divided into three groups according to area of research; 56 studies investigated the influence of spermatozoal DNA methylation on male fertility and abnormal semen parameters, 20 studies investigated spermatozoal DNA methylation in pregnancy outcomes including offspring health and 59 studies assessed the influence of environmental factors on spermatozoal DNA methylation. Findings from studies that scored as 'high' and 'moderate' quality of evidence according to GRADE criteria were summarized. We found that male subfertility and abnormal semen parameters, in particular oligozoospermia, appear to be associated with abnormal spermatozoal DNA methylation of imprinted regions. However, no specific DNA methylation signature of either subfertility or abnormal semen parameters has been convincingly replicated in genome-scale, unbiased analyses. Furthermore, although findings require independent replication, current evidence suggests that the spermatozoal DNA methylome is influenced by cigarette smoking, advanced age and environmental pollutants. Importantly however, from a clinical point of view, there is no convincing evidence that changes in spermatozoal DNA methylation influence pregnancy outcomes or offspring health.

WIDER IMPLICATIONS

Although it appears that the human sperm DNA methylome can be influenced by certain environmental and physiological traits, no findings have been robustly replicated between studies. We have generated a set of recommendations that would enhance the reliability and robustness of findings of future analyses of the human sperm methylome. Such studies will likely require multicentre collaborations to reach appropriate sample sizes, and should incorporate phenotype data in more complex statistical models.

Genome-wide associations between alcohol consumption and blood DNA methylation: evidence from twin study

Aim: Alcohol intake alters DNA methylation profiles and methylation might mediate the association between alcohol and disease, but limited number of positive CpG sites repeatedly replicated. Materials & methods: In total, 57 monozygotic (MZ) twin pairs discordant for alcohol drinking from the Chinese National Twin Registry and 158 MZ and dizygotic twin pairs in the Swedish Adoption/Twin Study of Aging were evaluated. DNA methylation was detected using the Infinium HumanMethylation450 BeadChip. Results: Among candidate CpG sites, cg07326074 was significantly correlated with drinking after adjusting for covariates in MZ twins in both datasets but not in the entire sample or dizygotic twins. Conclusion: The hypermethylation of cg07326074, located in the tumor-promoting gene C16orf59, was associated with alcohol consumption.

Prenatal medication exposure and epigenetic outcomes: a systematic literature review and recommendations for prenatal pharmacoepigenetic studies

When used during pregnancy, analgesics and psychotropics pass the placenta to enter the foetal circulation and may induce epigenetic modifications. Where such modifications occur and whether they disrupt normal foetal developme nt, are currently unanswered questions. This field of prenatal pharmacoepigenetics has received increasing attention, with several studies reporting associations between in utero medication exposure and offspring epigenetic outcomes. Nevertheless, no recent systematic review of the literature is available. Therefore, the objectives of this review were to (i) provide an overview of the literature on the association of prenatal exposure to psychotropics a nd analgesics with epigenetic outcomes, and (ii) suggest recommendations for future studies within prenatal pharmacoepigenetics. We performed systematic literature searches in five databases. The eligible studies assessed human prenatal exposure to psychotropics or analgesics, with epigenetic analyses of offspring tissue as an outcome. We identified 18 eligible studies including 4,419 neonates exposed to either antidepressants, antiepileptic drugs, paracetamol, acetylsalicylic acid, or methadone. The epigenetic outcome in all studies was DNA methylation in cord blood, placental tissue or buccal cells. Although most studies found significant differences in DNA methylation upon medication exposure, almost no differences were persistent across studies for similar medications and sequencing methods. The reviewed studies were challenging to compare due to poor transparency in reporting, and heterogeneous methodology, design, genome coverage, and statistical modelling. We propose 10 recommendations for future prenatal pharmacoepigenetic studies considering both epidemiological and epigenetic perspectives. These recommendations may improve the quality, comparability, and clinical relevance of such studies. PROSPERO registration ID: CRD42020166675.

Individuals Diagnosed with Binge-Eating Disorder Have DNA Hypomethylated Sites in Genes of the Metabolic System: A Pilot Study

Binge-eating disorder, recently accepted as a diagnostic category, is differentiated from bulimia nervosa in that the former shows the presence of binge-eating episodes and the absence of compensatory behavior. Epigenetics is a conjunct of mechanisms (like DNA methylation) that regulate gene expression, which are dependent on environmental changes. Analysis of DNA methylation in eating disorders shows that it is reduced. The present study aimed to analyze the genome-wide DNA methylation differences between individuals diagnosed with BED and BN. A total of 46 individuals were analyzed using the Infinium Methylation EPIC array. We found 11 differentially methylated sites between BED- and BN-diagnosed individuals, with genome-wide significance. Most of the associations were found in genes related to metabolic processes (ST3GAL4, PRKAG2, and FRK), which are hypomethylated genes in BED. Cg04781532, located in the body of the PRKAG2 gene (protein kinase AMP-activated non-catalytic subunit gamma 2), was hypomethylated in individuals with BED. Agonists of PRKAG2, which is the subunit of AMPK (AMP-activated protein kinase), are proposed to treat obesity, BED, and BN. The present study contributes important insights into the effect that BED could have on PRKAG2 activation.

Integrative Genomic-Epigenomic Analysis of Clozapine-Treated Patients with Refractory Psychosis

Clozapine (CLZ) is the only antipsychotic drug that has been proven to be effective in patients with refractory psychosis, but it has also been proposed as an effective mood stabilizer; however, the complex mechanisms of action of CLZ are not yet fully known. To find predictors of CLZ-associated phenotypes (i.e., the metabolic ratio, dosage, and response), we explore the genomic and epigenomic characteristics of 44 patients with refractory psychosis who receive CLZ treatment based on the integration of polygenic risk score (PRS) analyses in simultaneous methylome profiles. Surprisingly, the PRS for bipolar disorder (BD-PRS) was associated with the CLZ metabolic ratio (pseudo-R2 = 0.2080, adjusted p-value = 0.0189). To better explain our findings in a biological context, we assess the protein-protein interactions between gene products with high impact variants in the top enriched pathways and those exhibiting differentially methylated sites. The GABAergic synapse pathway was found to be enriched in BD-PRS and was associated with the CLZ metabolic ratio. Such interplay supports the use of CLZ as a mood stabilizer and not just as an antipsychotic. Future studies with larger sample sizes should be pursued to confirm the findings of this study.

Differential methylation of imprinting genes and MHC locus in 22q11.2 deletion syndrome-related schizophrenia spectrum disorders

OBJECTIVES

22q11.2 deletion syndrome (DS) is the strongest known genetic risk for schizophrenia. Methylome screening was conducted to elucidate possible involvement of epigenetic alterations in the emergence of schizophrenia spectrum disorders (SZ-SD) in 22q11.2DS.

METHODS

Sixteen adult men with/without SZ-SD were recruited from a 22q11.2DS cohort and underwent genome-wide DNA methylation profile analysis. Differentially methylated probes (DMPs) and regions (DMRs) were analysed using the ChAMP software.

RESULTS

The DMPs (p-value <10^-6) and DMRs (p-valueArea <0.01) were enriched in two gene sets, 'imprinting genes' and 'chr6p21', a region overlapping the MHC locus. Most of the identified imprinting genes are involved in neurodevelopment and located in clusters under imprinting control region (ICR) regulation, including PEG10, SGCE (7q21.3), GNAS, GNAS-AS1 (20q13.32) and SNHG14, SNURF-SNRPN, SNORD115 (15q11.2). The differentially methylated genes from the MHC locus included immune HLA-genes and non-immune genes, RNF39, PPP1R18 and NOTCH4, implicated in neurodevelopment and synaptic plasticity. The most significant DMR is located in MHC locus and covered the transcription regulator ZFP57 that is required for control and maintenance of gene imprinting at multiple ICRs.

CONCLUSIONS

The differential methylation in imprinting genes and in chr6p21-22 indicate the neurodevelopmental nature of 22q11.2DS-related SZ and the major role of MHC locus in the risk to develop SZ.

An epigenome-wide association study of metabolic syndrome and its components

The role of metabolic syndrome (MetS) as a preceding metabolic state for type 2 diabetes and cardiovascular disease is widely recognised. To accumulate knowledge of the pathological mechanisms behind the condition at the methylation level, we conducted an epigenome-wide association study (EWAS) of MetS and its components, testing 1187 individuals of European ancestry for approximately 470 000 methylation sites throughout the genome. Methylation site cg19693031 in gene TXNIP —previously associated with type 2 diabetes, glucose and lipid metabolism, associated with fasting glucose level (P = 1.80 × 10⁻⁸). Cg06500161 in gene ABCG1 associated both with serum triglycerides (P = 5.36 × 10⁻⁹) and waist circumference (P = 5.21 × 10⁻⁹). The previously identified type 2 diabetes–associated locus cg08309687 in chromosome 21 associated with waist circumference for the first time (P = 2.24 × 10⁻⁷). Furthermore, a novel HDL association with cg17901584 in chromosome 1 was identified (P = 7.81 × 10⁻⁸). Our study supports previous genetic studies of MetS, finding that lipid metabolism plays a key role in pathology of the syndrome. We provide evidence regarding a close interplay with glucose metabolism. Finally, we suggest that in attempts to identify methylation loci linking separate MetS components, cg19693031 appears to represent a strong candidate.

Epigenetic profiling of Italian patients identified methylation sites associated with hereditary transthyretin amyloidosis

Hereditary transthyretin (TTR) amyloidosis (hATTR) is a rare life-threatening disorder caused by amyloidogenic coding mutations located in TTR gene. To understand the high phenotypic variability observed among carriers of TTR disease-causing mutations, we conducted an epigenome-wide association study (EWAS) assessing more than 700,000 methylation sites and testing epigenetic difference of TTR coding mutation carriers vs. non-carriers. We observed a significant methylation change at cg09097335 site located in Beta-secretase 2 (BACE2) gene (standardized regression coefficient = -0.60, p = 6.26 × 10-8). This gene is involved in a protein interaction network enriched for biological processes and molecular pathways related to amyloid-beta metabolism (Gene Ontology: 0050435, q = 0.007), amyloid fiber formation (Reactome HSA-977225, q = 0.008), and Alzheimer's disease (KEGG hsa05010, q = 2.2 × 10-4). Additionally, TTR and BACE2 share APP (amyloid-beta precursor protein) as a validated protein interactor. Within TTR gene region, we observed that Val30Met disrupts a methylation site, cg13139646, causing a drastic hypomethylation in carriers of this amyloidogenic mutation (standardized regression coefficient = -2.18, p = 3.34 × 10-11). Cg13139646 showed co-methylation with cg19203115 (Pearson's r2 = 0.32), which showed significant epigenetic differences between symptomatic and asymptomatic carriers of amyloidogenic mutations (standardized regression coefficient = -0.56, p = 8.6 × 10-4). In conclusion, we provide novel insights related to the molecular mechanisms involved in the complex heterogeneity of hATTR, highlighting the role of epigenetic regulation in this rare disorder.

DNA Methylation Associated With Diabetic Kidney Disease in Blood-Derived DNA

A subset of individuals with type 1 diabetes will develop diabetic kidney disease (DKD). DKD is heritable and large-scale genome-wide association studies have begun to identify genetic factors that influence DKD. Complementary to genetic factors, we know that a person’s epigenetic profile is also altered with DKD. This study reports analysis of DNA methylation, a major epigenetic feature, evaluating methylome-wide loci for association with DKD. Unique features (n = 485,577; 482,421 CpG probes) were evaluated in blood-derived DNA from carefully phenotyped White European individuals diagnosed with type 1 diabetes with (cases) or without (controls) DKD (n = 677 samples). Explicitly, 150 cases were compared to 100 controls using the 450K array, with subsequent analysis using data previously generated for a further 96 cases and 96 controls on the 27K array, and de novo methylation data generated for replication in 139 cases and 96 controls. Following stringent quality control, raw data were quantile normalized and beta values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls; resultant P-values for array-based data were adjusted for multiple testing. Genes with significantly increased (hypermethylated) and/or decreased (hypomethylated) levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways. Twenty-two loci demonstrated statistically significant fold changes associated with DKD and additional support for these associated loci was sought using independent samples derived from patients recruited with similar inclusion criteria. Markers associated with CCNL1 and ZNF187 genes are supported as differentially regulated loci (P < 10^–8), with evidence also presented for AFF3, which has been identified from a meta-analysis and subsequent replication of genome-wide association studies. Further supporting evidence for differential gene expression in CCNL1 and ZNF187 is presented from kidney biopsy and blood-derived RNA in people with and without kidney disease from NephroSeq. Evidence confirming that methylation sites influence the development of DKD may aid risk prediction tools and stimulate research to identify epigenomic therapies which might be clinically useful for this disease.

"Epigenome-wide methylation profile of chronic kidney disease-derived arterial DNA uncovers novel pathways in disease-associated cardiovascular pathology."

Chronic kidney disease (CKD) related cardiovascular disease (CVD) is characterized by vascular remodelling with well-established structural and functional changes in the vascular wall such as arterial stiffness, matrix deposition, and calcification. These phenotypic changes resemble pathology seen in ageing, and are likely to be mediated by sustained alterations in gene expression, which may be caused by epigenetic changes such as tissue-specific DNA methylation. We aimed to investigate tissue specific changes in DNA methylation that occur in CKD-related CVD. Genome-wide DNA methylation changes were examined in bisulphite converted genomic DNA isolated from the vascular media of CKD and healthy arteries. Methylation-specific PCR was used to validate the array data, and the association between DNA methylation and gene and protein expression was examined. The DNA methylation age was compared to the chronological age in both cases and controls. Three hundred and nineteen differentially methylated regions (DMR) were identified spread across the genome. Pathway analysis revealed that DMRs associated with genes were involved in embryonic and vascular development, and signalling pathways such as TGFβ and FGF. Expression of top differentially methylated gene HOXA5 showed a significant negative correlation with DNA methylation. Interestingly, DNA methylation age and chronological age were highly correlated, but there was no evidence of accelerated age-related DNA methylation in the arteries of CKD patients. In conclusion, we demonstrated that differential DNA methylation in the arterial tissue of CKD patients represents a potential mediator of arterial pathology and may be used to uncover novel pathways in the genesis of CKD-associated complications.

Genome-wide associations of human gut microbiome variation and implications for causal inference analyses

Recent population-based1-4 and clinical studies5 have identified a range of factors associated with human gut microbiome variation. Murine quantitative trait loci6, human twin studies7 and microbiome genome-wide association studies1,3,8-12 have provided evidence for genetic contributions to microbiome composition. Despite this, there is still poor overlap in genetic association across human studies. Using appropriate taxon-specific models, along with support from independent cohorts, we show an association between human host genotype and gut microbiome variation. We also suggest that interpretation of applied analyses using genetic associations is complicated by the probable overlap between genetic contributions and heritable components of host environment. Using faecal 16S ribosomal RNA gene sequences and host genotype data from the Flemish Gut Flora Project (n = 2,223) and two German cohorts (FoCus, n = 950; PopGen, n = 717), we identify genetic associations involving multiple microbial traits. Two of these associations achieved a study-level threshold of P = 1.57 × 10-10; an association between Ruminococcus and rs150018970 near RAPGEF1 on chromosome 9, and between Coprococcus and rs561177583 within LINC01787 on chromosome 1. Exploratory analyses were undertaken using 11 other genome-wide associations with strong evidence for association (P < 2.5 × 10-8) and a previously reported signal of association between rs4988235 (MCM6/LCT) and Bifidobacterium. Across these 14 single-nucleotide polymorphisms there was evidence of signal overlap with other genome-wide association studies, including those for age at menarche and cardiometabolic traits. Mendelian randomization analysis was able to estimate associations between microbial traits and disease (including Bifidobacterium and body composition); however, in the absence of clear microbiome-driven effects, caution is needed in interpretation. Overall, this work marks a growing catalogue of genetic associations that will provide insight into the contribution of host genotype to gut microbiome. Despite this, the uncertain origin of association signals will likely complicate future work looking to dissect function or use associations for causal inference analysis.

DNA methylation biomarkers of future health outcomes in children

Biomarkers which predict future health outcomes are key to the goals of precision health. Such biomarkers do not have to be involved in the causal pathway of a disease, and their performance is best assessed using statistical tests of clinical performance and evaluation of net health impact. DNA methylation is the most commonly studied epigenetic process and represents a potential biomarker of future health outcomes. We review 25 studies in non-oncological paediatric conditions where DNA methylation biomarkers of future health outcomes are assessed. Whilst a number of positive findings have been described, the body of evidence is severely limited by issues with outcome measures, tissue-specific samples, accounting for sample cell type heterogeneity, lack of appropriate statistical testing, small effect sizes, limited validation, and no assessment of net health impact. Future studies should concentrate on careful study design to overcome these issues, and integration of DNA methylation data with other 'omic', clinical, and environmental data to generate the most clinically useful biomarkers of paediatric disease.

DNA Methylation Associated with Mitochondrial Dysfunction in a South African Autism Spectrum Disorder Cohort

Autism spectrum disorder (ASD) is characterized by phenotypic heterogeneity and a complex genetic architecture which includes distinctive epigenetic patterns. We report differential DNA methylation patterns associated with ASD in South African children. An exploratory whole‐epigenome methylation screen using the Illumina 450 K MethylationArray identified differentially methylated CpG sites between ASD and controls that mapped to 898 genes (P ≤ 0.05) which were enriched for nine canonical pathways converging on mitochondrial metabolism and protein ubiquitination. Targeted Next Generation Bisulfite Sequencing of 27 genes confirmed differential methylation between ASD and control in our cohort. DNA pyrosequencing of two of these genes, the mitochondrial enzyme Propionyl‐CoA Carboxylase subunit Beta (PCCB) and Protocadherin Alpha 12 (PCDHA12), revealed a wide range of methylation levels (9–49% and 0–54%, respectively) in both ASD and controls. Three CpG loci were differentially methylated in PCCB (P ≤ 0.05), while PCDHA12, previously linked to ASD, had two significantly different CpG sites (P ≤ 0.001) between ASD and control. Differentially methylated CpGs were hypomethylated in ASD. Metabolomic analysis of urinary organic acids revealed that three metabolites, 3‐hydroxy‐3‐methylglutaric acid (P = 0.008), 3‐methyglutaconic acid (P = 0.018), and ethylmalonic acid (P = 0.043) were significantly elevated in individuals with ASD. These metabolites are directly linked to mitochondrial respiratory chain disorders, with a putative link to PCCB, consistent with impaired mitochondrial function. Our data support an association between DNA methylation and mitochondrial dysfunction in the etiology of ASD. Autism Res 2020, 13: 1079‐1093. © 2020 The Authors. Autism Research published by International Society for Autism Research published by Wiley Periodicals, Inc.

Lay Summary

Epigenetic changes are chemical modifications of DNA which can change gene function. DNA methylation, a type of epigenetic modification, is linked to autism. We examined DNA methylation in South African children with autism and identified mitochondrial genes associated with autism. Mitochondria are power‐suppliers in cells and mitochondrial genes are essential to metabolism and energy production, which are important for brain cells during development. Our findings suggest that some individuals with ASD also have mitochondrial dysfunction.

Intraspecific and interspecific investigations of skeletal DNA methylation and femur morphology in primates

OBJECTIVES

Epigenetic mechanisms influence the development and maintenance of complex phenotypes and may also contribute to the evolution of species-specific phenotypes. With respect to skeletal traits, little is known about the gene regulation underlying these hard tissues or how tissue-specific patterns are associated with bone morphology or vary among species. To begin exploring these topics, this study evaluates one epigenetic mechanism, DNA methylation, in skeletal tissues from five nonhuman primate species which display anatomical and locomotor differences representative of their phylogenetic groups.

MATERIALS AND METHODS

First, we test whether intraspecific variation in skeletal DNA methylation is associated with intraspecific variation in femur morphology. Second, we identify interspecific differences in DNA methylation and assess whether these lineage-specific patterns may have contributed to species-specific morphologies. Specifically, we use the Illumina Infinium MethylationEPIC BeadChip to identify DNA methylation patterns in femur trabecular bone from baboons (n = 28), macaques (n = 10), vervets (n = 10), chimpanzees (n = 4), and marmosets (n = 6).

RESULTS

Significant differentially methylated positions (DMPs) were associated with a subset of morphological variants, but these likely have small biological effects and may be confounded by other variables associated with morphological variation. Conversely, several species-specific DMPs were identified, and these are found in genes enriched for functions associated with complex skeletal traits.

DISCUSSION

Overall, these findings reveal that while intraspecific epigenetic variation is not readily associated with skeletal morphology differences, some interspecific epigenetic differences in skeletal tissues exist and may contribute to evolutionarily distinct phenotypes. This work forms a foundation for future explorations of gene regulation and skeletal trait evolution in primates.

Evaluating liquid biopsies for methylomic profiling of prostate cancer

Background: Liquid biopsies offer significant potential for informing on cancer progression and therapeutic resistance via minimally invasive serial monitoring of genetic alterations. Although the cancer epigenome is a central driving force in most neoplasia, the accuracy of monitoring the tumor methylome using liquid biopsies remains relatively unknown. Objectives: to investigate how well two types of liquid biopsy (urine and blood) capture the prostate cancer methylome, and may thus serve as a non-invasive surrogate for studying the tumor epigenome. Methods: A cohort of four metastatic treatment naïve prostate cancer (PCa) patients was selected. Matched biopsy cores (tumor and histologically matched-normal), post-DRE, pre-biopsy urine, and peripheral blood plasma were available for each subject. DNA methylation was profiled utilizing the Infinium® MethylationEPIC BeadChip (Illumina) and analysed using the RnBeads software. Significantly (FDR adjusted P < 0.05) differentially methylated probes (DMPs) between tumor and MN were identified and examined in the liquids (done at a grouped and individual subject level). Results: DNA methylation analysis of urine and blood in men with metastatic PCa showed highly correlated patterns between the different liquid types (ρ = 0.93, P < 0.0001), with large contributions from non-tumor sources. DNA methylation profiles of liquids were more similar between subjects, than intra-individual liquid-tumor correlations. Overall, both urine and plasma are viable surrogates for tumor tissue biopsies, capturing up to 39.40% and 64.14% of tumor-specific methylation alterations, respectively. Conclusion: We conclude that both urine and blood plasma are easily accessible and sensitive biofluids for the study of PCa epigenomic alterations.

Tobacco smoking induces changes in true DNA methylation, hydroxymethylation and gene expression in bronchoalveolar lavage cells

Background

While smoking is known to associate with development of multiple diseases, the underlying mechanisms are still poorly understood. Tobacco smoking can modify the chemical integrity of DNA leading to changes in transcriptional activity, partly through an altered epigenetic state. We aimed to investigate the impact of smoking on lung cells collected from bronchoalveolar lavage (BAL).

Methods

We profiled changes in DNA methylation (5mC) and its oxidised form hydroxymethylation (5hmC) using conventional bisulphite (BS) treatment and oxidative bisulphite treatment with Illumina Infinium MethylationEPIC BeadChip, and examined gene expression by RNA-seq in healthy smokers.

Findings

We identified 1667 total 5mC + 5hmC, 1756 5mC and 67 5hmC differentially methylated positions (DMPs) between smokers and non-smokers (FDR-adjusted P <.05, absolute Δβ >0.15). Both 5mC DMPs and to a lesser extent 5mC + 5hmC were predominantly hypomethylated. In contrast, almost all 5hmC DMPs were hypermethylated, supporting the hypothesis that smoking-associated oxidative stress can lead to DNA demethylation, via the established sequential oxidation of which 5hmC is the first step. While we confirmed differential methylation of previously reported smoking-associated 5mC + 5hmC CpGs using former generations of BeadChips in alveolar macrophages, the large majority of identified DMPs, 5mC + 5hmC (1639/1667), 5mC (1738/1756), and 5hmC (67/67), have not been previously reported. Most of these novel smoking-associating sites are specific to the EPIC BeadChip and, interestingly, many of them are associated to FANTOM5 enhancers. Transcriptional changes affecting 633 transcripts were consistent with DNA methylation profiles and converged to alteration of genes involved in migration, signalling and inflammatory response of immune cells.

Interpretation

Collectively, these findings suggest that tobacco smoke exposure epigenetically modifies BAL cells, possibly involving a continuous active demethylation and subsequent increased activity of inflammatory processes in the lungs.

Fund

The study was supported by the Swedish Research Council, the Swedish Heart-Lung Foundation, the Stockholm County Council (ALF), the King Gustav's and Queen Victoria's Freemasons' Foundation, Knut and Alice Wallenberg Foundation, Neuro Sweden, and the Swedish MS foundation.