Sex differences in DNA methylation assessed by 450 K BeadChip in newborns

A systematic review on the contribution of DNA methylation to hearing loss

BACKGROUND

DNA methylation may have a regulatory role in monogenic sensorineural hearing loss and complex, polygenic phenotypic forms of hearing loss, including age-related hearing impairment or Meniere disease. The purpose of this systematic review is to critically assess the evidence supporting a functional role of DNA methylation in phenotypes associated with hearing loss.

RESULTS

The search strategy yielded a total of 661 articles. After quality assessment, 25 records were selected (12 human DNA methylation studies, 5 experimental animal studies and 8 studies reporting mutations in the DNMT1 gene). Although some methylation studies reported significant differences in CpG methylation in diverse gene promoters associated with complex hearing loss phenotypes (ARHI, otosclerosis, MD), only one study included a replication cohort that supported a regulatory role for CpG methylation in the genes TCF25 and POLE in ARHI. Conversely, several studies have independently confirmed pathogenic mutations within exon 21 of the DNMT1 gene, which encodes the DNA (cytosine-5)-methyltransferase 1 enzyme. This methylation enzyme is strongly associated with a rare disease defined by autosomal dominant cerebellar ataxia, deafness and narcolepsy (ADCA-DN). Of note, rare variants in DNMT1 and DNMT3A genes have also been reported in noise-induced hearing loss.

CONCLUSIONS

Evidence supporting a functional role for DNA methylation in hearing loss is limited to few genes in complex disorders such as ARHI. Mutations in the DNMT1 gene are associated with ADCA-DN, suggesting the CpG methylation in hearing loss genes deserves further attention in hearing research.

Prenatal exposure to environmental phenols and phthalates and altered patterns of DNA methylation in childhood

INTRODUCTION

Epigenetic marks are key biomarkers linking the prenatal environment to health and development. However, DNA methylation associations and persistence of marks for prenatal exposure to multiple Endocrine Disrupting Chemicals (EDCs) in human populations have not been examined in great detail.

METHODS

We measured Bisphenol-A (BPA), triclosan, benzophenone-3 (BP3), methyl-paraben, propyl-paraben, and butyl-paraben, as well as 11 phthalate metabolites, in two pregnancy urine samples, at approximately 13 and 26 weeks of gestation in participants of the Center for the Health Assessment of Mothers and Children of Salinas (CHAMACOS) study (N = 309). DNA methylation of cord blood at birth and child peripheral blood at ages 9 and 14 years was measured with 450K and EPIC arrays. Robust linear regression was used to identify differentially methylated probes (DMPs), and comb-p was used to identify differentially methylated regions (DMRs) in association with pregnancy-averaged EDC concentrations. Quantile g-computation was used to assess associations of the whole phenol/phthalate mixture with DMPs and DMRs.

RESULTS

Prenatal BPA exposure was associated with 1 CpG among males and Parabens were associated with 10 CpGs among females at Bonferroni-level significance in cord blood. Other suggestive DMPs (unadjusted p-value < 1 × 10-6) and several DMRs associated with the individual phenols and whole mixture were also identified. A total of 10 CpG sites at least suggestively associated with BPA, Triclosan, BP3, Parabens, and the whole mixture in cord blood were found to persist into adolescence in peripheral blood.

CONCLUSIONS

We found sex-specific associations between prenatal phenol exposure and DNA methylation, particularly with BPA in males and Parabens in females. Additionally, we found several DMPs that maintained significant associations with prenatal EDC exposures at age 9 and age 14 years.

Sex-related DNA methylation is associated with inflammation and gene expression in the lungs of healthy individuals

Lung cancer exhibits sex-biased molecular characteristics and epidemiological trends, suggesting a need for sex-specific approaches to understanding its etiology and treatment. DNA methylation alterations play critical roles in lung carcinogenesis and may serve as valuable biomarkers for precision medicine strategies. We employed the Infinium MethylationEPIC array to identify autosomal sex-related differentially methylated CpG sites (DM-CpGs) in lung epithelium of healthy individuals (32 females and 37 males) while controlling for age, BMI, and tobacco use. We correlated DM-CpGs with gene expression in lung epithelium and immune responses in bronchoalveolar lavage. We validated these DM-CpGs in lung tumors and adjacent normal tissue from The Cancer Genome Atlas (TCGA). Among 522 identified DM-CpGs, 61% were hypermethylated in females, predominantly located in promoter regions. These DM genes were implicated in cell-to-cell signaling, cellular function, transport, and lipid metabolism. Correlation analysis revealed sex-specific patterns between DM-CpGs and gene expression. Additionally, several DM-CpGs were correlated significantly with cytokines (IL-1β, IL-4, IL-12p70, and IFN-γ), macrophage, and lymphocyte counts. Also, some DM-CpGs were observed in TCGA lung adenocarcinoma, squamous cell carcinoma, and adjacent normal tissues. Our findings highlight sex-specific DNA methylation patterns in healthy lung epithelium and their associations with lung gene expression and lung immune biomarkers. These findings underscore the potential role of lung sex-related CpGs as epigenetic predispositions influencing sex disparities in lung cancer risk and outcomes, warranting further investigation for personalized lung cancer management strategies.

Maternal sensitivity and child internalizing and externalizing behavior: a mediating role for glucocorticoid receptor gene (NR3C1) methylation?

The early caregiving environment can have lasting effects on child mental health. Animal models suggest that glucocorticoid receptor gene (NR3C1) DNA methylation plays a mediating role in linking more responsive caregiving to improved behavioral outcomes by its impact on the stress regulatory system. In this longitudinal study, we examined whether children's NR3C1 methylation levels mediate an effect of maternal sensitivity in infancy on levels of child internalizing and externalizing behavior in a community sample. Maternal sensitivity of 145 mothers was rated at infant age 5 weeks, 12 months, and 30 months by observing mother-infant interactions. Buccal DNA methylation was assessed in the same children at age 6 years and maternal-reported internalizing and externalizing behavior was assessed at age 6 and 10 years. Higher sensitivity at age 5 weeks significantly predicted lower DNA methylation levels at two NR3C1 CpG loci, although methylation levels at these loci did not mediate an effect of maternal sensitivity on levels of child internalizing and externalizing behavior. Overall, the study provides evidence that maternal sensitivity in early infancy is associated with DNA methylation levels at loci involved in stress regulation, but the significance of this finding for child mental health remains unclear.

Sex differences in DNA methylation across gestation: a large scale, cross-cohort, multi-tissue analysis

Biological sex is a key variable influencing many physiological systems. Disease prevalence as well as treatment success can be modified by sex. Differences emerge already early in life and include pregnancy complications and adverse birth outcomes. The placenta is a critical organ for fetal development and shows sex-based differences in the expression of hormones and cytokines. Epigenetic regulation, such as DNA methylation (DNAm), may underlie the previously reported placental sexual dimorphism. We associated placental DNAm with fetal sex in three cohorts. Individual cohort results were meta-analyzed with random-effects modelling. CpG-sites differentially methylated with sex were further investigated regarding pathway enrichment, overlap with methylation quantitative trait loci (meQTLs), and hits from phenome-wide association studies (PheWAS). We evaluated the consistency of findings across tissues (CVS, i.e. chorionic villus sampling from early placenta, and cord blood) as well as with gene expression. We identified 10,320 epigenome-wide significant sex-differentially methylated probes (DMPs) spread throughout the epigenome of the placenta at birth. Most DMPs presented with lower DNAm levels in females. DMPs mapped to genes upregulated in brain, were enriched for neurodevelopmental pathways and significantly overlapped with meQTLs and PheWAS hits. Effect sizes were moderately correlated between CVS and placenta at birth, but only weakly correlated between birth placenta and cord blood. Sex differential gene expression in birth placenta was less pronounced and implicated genetic regions only marginally overlapped with those associated with differential DNAm. Our study provides an integrative perspective on sex-differential DNAm in perinatal tissues underscoring the possible link between placenta and brain.

Donor whole blood DNA methylation is not a strong predictor of acute graft versus host disease in unrelated donor allogeneic haematopoietic cell transplantation

Allogeneic hematopoietic cell transplantation (HCT) is used to treat many blood-based disorders and malignancies, however it can also result in serious adverse events, such as the development of acute graft-versus-host disease (aGVHD). This study aimed to develop a donor-specific epigenetic classifier to reduce incidence of aGVHD by improving donor selection. Genome-wide DNA methylation was assessed in a discovery cohort of 288 HCT donors selected based on recipient aGVHD outcome; this cohort consisted of 144 cases with aGVHD grades III-IV and 144 controls with no aGVHD. We applied a machine learning algorithm to identify CpG sites predictive of aGVHD. Receiver operating characteristic (ROC) curve analysis of these sites resulted in a classifier with an encouraging area under the ROC curve (AUC) of 0.91. To test this classifier, we used an independent validation cohort (n = 288) selected using the same criteria as the discovery cohort. Attempts to validate the classifier failed with the AUC falling to 0.51. These results indicate that donor DNA methylation may not be a suitable predictor of aGVHD in an HCT setting involving unrelated donors, despite the initial promising results in the discovery cohort. Our work highlights the importance of independent validation of machine learning classifiers, particularly when developing classifiers intended for clinical use.

Sex effects on DNA methylation affect discovery in epigenome-wide association study of schizophrenia

Sex differences in the epidemiology and clinical characteristics of schizophrenia are well-known; however, the molecular mechanisms underlying these differences remain unclear. Further, the potential advantages of sex-stratified meta-analyses of epigenome-wide association studies (EWAS) of schizophrenia have not been investigated. Here, we performed sex-stratified EWAS meta-analyses to investigate whether sex stratification improves discovery, and to identify differentially methylated regions (DMRs) in schizophrenia. Peripheral blood-derived DNA methylation data from 1519 cases of schizophrenia (male n = 989, female n = 530) and 1723 controls (male n = 997, female n = 726) from three publicly available datasets, and the TOP cohort were meta-analyzed to compare sex-specific, sex-stratified, and sex-adjusted EWAS. The predictive power of each model was assessed by polymethylation score (PMS). The number of schizophrenia-associated differentially methylated positions identified was higher for the sex-stratified model than for the sex-adjusted one. We identified 20 schizophrenia-associated DMRs in the sex-stratified analysis. PMS from sex-stratified analysis outperformed that from sex-adjusted analysis in predicting schizophrenia. Notably, PMSs from the sex-stratified and female-only analyses, but not those from sex-adjusted or the male-only analyses, significantly predicted schizophrenia in males. The findings suggest that sex-stratified EWAS meta-analyses improve the identification of schizophrenia-associated epigenetic changes and highlight an interaction between sex and schizophrenia status on DNA methylation. Sex-specific DNA methylation may have potential implications for precision psychiatry and the development of stratified treatments for schizophrenia.

Sex-specific DNA methylation in saliva from the multi-ethnic Future of Families and Child Wellbeing Study

The prevalence and severity of many diseases differs by sex, potentially due to sex-specific patterns in DNA methylation. Autosomal sex-specific differences in DNA methylation have been observed in cord blood and placental tissue but are not well studied in saliva or in diverse populations. We sought to characterize sex-specific DNA methylation on autosomal chromosomes in saliva samples from children in the Future of Families and Child Wellbeing Study, a multi-ethnic prospective birth cohort containing an oversampling of Black, Hispanic and low-income families. DNA methylation from saliva samples was analysed on 796 children (50.6% male) at both ages 9 and 15 with DNA methylation measured using the Illumina HumanMethylation 450k array. An epigenome-wide association analysis of the age 9 samples identified 8,430 sex-differentiated autosomal DNA methylation sites (P < 2.4 × 10-7), of which 76.2% had higher DNA methylation in female children. The strongest sex-difference was in the cg26921482 probe, in the AMDHD2 gene, with 30.6% higher DNA methylation in female compared to male children (P < 1 × 10-300). Treating the age 15 samples as an internal replication set, we observed highly consistent results between the ages 9 and 15 measurements, indicating stable and replicable sex-differentiation. Further, we directly compared our results to previously published DNA methylation sex differences in both cord blood and saliva and again found strong consistency. Our findings support widespread and robust sex-differential DNA methylation across age, human tissues, and populations. These findings help inform our understanding of potential biological processes contributing to sex differences in human physiology and disease.

The methylome of buccal epithelial cells is influenced by age, sex, and physiological properties

Epigenetic modifications, particularly DNA methylation, have emerged as regulators of gene expression and are implicated in various biological processes and disease states. Understanding the factors influencing the epigenome is essential for unraveling its complexity. In this study, we aimed to identify how the methylome of buccal epithelial cells, a noninvasive and easily accessible tissue, is associated with demographic and health-related variables commonly used in clinical settings, such as age, sex, blood immune composition, hemoglobin levels, and others. We developed a model to assess the association of multiple factors with the human methylome and identify the genomic loci significantly impacted by each trait. We demonstrated that DNA methylation variation is accurately modeled by several factors. We confirmed the well-known impact of age and sex and unveiled novel clinical factors associated with DNA methylation, such as blood neutrophils, hemoglobin, red blood cell distribution width, high-density lipoprotein cholesterol, and urea. Genomic regions significantly associated with these traits were enriched in relevant transcription factors, drugs, and diseases. Among our findings, we showed that neutrophil-impacted loci were involved in neutrophil functionality and maturation. Similarly, hemoglobin-influenced sites were associated with several diseases, including aplastic anemia, and the genomic loci affected by urea were related to congenital anomalies of the kidney and urinary tract. Our findings contribute to a better understanding of the human methylome plasticity and provide insights into novel factors shaping DNA methylation patterns, highlighting their potential clinical implications as biomarkers and the importance of considering these physiological traits in future medical epigenomic investigations.NEW & NOTEWORTHY We have developed a quantitative model to assess how the human methylome is associated with several factors and to identify the genomic loci significantly impacted by each trait. We reported novel health-related factors driving DNA methylation patterns and new site-specific regulations that further elucidate methylome dynamics. Our study contributes to a better understanding of the plasticity of the human methylome and unveils novel physiological traits with a potential role in future medical epigenomic investigations.

A higher dysregulation burden of brain DNA methylation in female patients implicated in the sex bias of Schizophrenia

Sex differences are pervasive in schizophrenia (SCZ), but the extent and magnitude of DNA methylation (DNAm) changes underlying these differences remain uncharacterized. In this study, sex-stratified differential DNAm analysis was performed in postmortem brain samples from 117 SCZ and 137 controls, partitioned into discovery and replication datasets. Three differentially methylated positions (DMPs) were identified (adj.p < 0.05) in females and 29 DMPs in males without overlap between them. Over 81% of these sex-stratified DMPs were directionally consistent between sexes but with different effect sizes. Females experienced larger magnitude of DNAm changes and more DMPs (based on data of equal sample size) than males, contributing to a higher dysregulation burden of DNAm in females SCZ. Additionally, despite similar proportions of female-related DMPs (fDMPs, 8%) being under genetic control compared with males (10%), significant enrichment of DMP-related single nucleotide polymorphisms (SNPs) in signals of genome-wide association studies was identified only in fDMPs. One DMP in each sex connected the SNPs and gene expression of CALHM1 in females and CCDC149 in males. PPI subnetworks revealed that both female- and male-related differential DNAm interacted with synapse-related dysregulation. Immune-related pathways were unique for females and neuron-related pathways were associated with males. This study reveals remarkable quantitative differences in DNAm-related sexual dimorphism in SCZ and that females have a higher dysregulation burden of SCZ-associated DNAm than males.

Maternal adverse childhood experiences (ACEs) and DNA methylation of newborns in cord blood

BACKGROUND

Adverse childhood experiences (ACEs) increase the risk of poor health outcomes later in life. Psychosocial stressors may also have intergenerational health effects by which parental ACEs are associated with mental and physical health of children. Epigenetic programming may be one mechanism linking parental ACEs to child health. This study aimed to investigate epigenome-wide associations of maternal preconception ACEs with DNA methylation patterns of children. In the Center for the Health Assessment of Mothers and Children of Salinas study, cord blood DNA methylation was measured using the Illumina HumanMethylation450 BeadChip. Preconception ACEs, which occurred during the mothers' childhoods, were collected using a standard ACE questionnaire including 10 ACE indicators. Maternal ACE exposures were defined in this study as (1) the total number of ACEs; (2) the total number of ACEs categorized as 0, 1-3, and > 4; and (3) individual ACEs. Associations of ACE exposures with differential methylated positions, regions, and CpG modules determined using weighted gene co-expression network analysis were evaluated adjusting for covariates.

RESULTS

Data on maternal ACEs and cord blood DNA methylation were available for 196 mother/newborn pairs. One differential methylated position was associated with maternal experience of emotional abuse (cg05486260/FAM135B gene; q value < 0.05). Five differential methylated regions were significantly associated with the total number of ACEs, and 36 unique differential methylated regions were associated with individual ACEs (Šidák p value < 0.05). Fifteen CpG modules were significantly correlated with the total number of ACEs or individual ACEs, of which 8 remained significant in fully adjusted models (p value < 0.05). Significant modules were enriched for pathways related to neurological and immune development and function.

CONCLUSIONS

Maternal ACEs prior to conception were associated with cord blood DNA methylation of offspring at birth. Although there was limited overlap between differential methylated regions and CpGs in modules associated with ACE exposures, statistically significant regions and networks were related to genes involved in neurological and immune function. Findings may provide insights to pathways linking psychosocial stressors to health. Further research is needed to understand the relationship between changes in DNA methylation and child health.

Downstream Assays for Variant Resolution: Epigenetics, RNA Sequnncing, and Metabolomics

As the availability of advanced molecular testing like whole exome and genome sequencing expands, it comes with the added complication of interpreting inconclusive results, including determining the relevance of variants of uncertain significance or failing to find a variant in an otherwise suspected specific genetic disorder. This complication necessitates the use of alternative testing methods to gather more information in support of, or against, a particular genetic diagnosis. Therefore, new genome-wide approaches, including DNA epigenetic testing, RNA sequencing, and metabolomics, are increasingly being used to increase the diagnostic yield when used in conjunction with more conventional genetic tests.

A permutation-based approach using a rank-based statistic to identify sex differences in epigenetics

Epigenetic sex differences and their resulting implications for human health have been studied for about a decade. The objective of this paper is to use permutation-based inference and a new ranked-based test statistic to identify sex-based epigenetic differences in the human DNA methylome. In particular, we examine whether we could identify separations between the female and male distributions of DNA methylation across hundred of thousands CpG sites in two independent cohorts, the Swedish Adoption Twin study and the Lamarck study. Based on Fisherian p-values, we set a threshold for methylation differences "worth further scrutiny". At this threshold, we were able to confirm previously-found CpG sites that stratify with respect to sex. These CpG sites with sex differences in DNA methylation should be further investigated for their possible contribution to various physiological and pathological functions in the human body. We followed-up our statistical analyses with a literature review in order to inform the proposed disease implications for the loci we uncovered.

Sexually dimorphic methylation patterns characterize the placenta and blood from extremely preterm newborns

BACKGROUND

Health outcomes among children born prematurely are known to be sexually dimorphic, with male infants often more affected, yet the mechanism behind this observation is not clear. CpG methylation levels in the placenta and blood also differ by sex and are associated with adverse health outcomes. We contrasted CpG methylation levels in the placenta and neonatal blood (n = 358) from the Extremely Low Gestational Age Newborn (ELGAN) cohort based on the EPIC array, which assays over 850,000 CpG sites across the epigenome. Sex-specific epigenome-wide association analyses were conducted for the placenta and neonatal blood samples independently, and the results were compared to determine tissue-specific differences between the methylation patterns in males and females. All models were adjusted for cell type heterogeneity. Enrichment pathway analysis was performed to identify the biological functions of genes related to the sexually dimorphic CpG sites.

RESULTS

Approximately 11,500 CpG sites were differentially methylated in relation to sex. Of these, 5949 were placenta-specific and 5361 were blood-specific, with only 233 CpG sites overlapping in both tissues. For placenta-specific CpG sites, 90% were hypermethylated in males. For blood-specific CpG sites, 95% were hypermethylated in females. In the placenta, keratinocyte differentiation biological pathways were enriched among the differentially methylated genes. No enrichment pathways were observed for blood.

CONCLUSIONS

Distinct methylation patterns were observed between male and female children born extremely premature, and keratinocyte differentiation pathways were enriched in the placenta. These findings provide new insights into the epigenetic mechanisms underlying sexually dimorphic health outcomes among extremely premature infants.

Associations between prenatal phthalate exposure and childhood epigenetic age acceleration

BACKGROUND

Phthalates, a group of pervasive endocrine-disrupting chemicals found in plastics and personal care products, have been associated with a wide range of developmental and health outcomes. However, their impact on biomarkers of aging has not been characterized. We tested associations between prenatal exposure to 11 phthalate metabolites on epigenetic aging in children at birth, 7, 9, and 14 years of age. We hypothesized that prenatal phthalate exposure will be associated with epigenetic age acceleration measures at birth and in early childhood, with patterns dependent on sex and timing of DNAm measurement.

METHODS

Among 385 mother-child pairs from the CHAMACOS cohort, we measured DNAm at birth, 7, 9, and 14 years of age, and utilized adjusted linear regression to assess the association between prenatal phthalate exposure and Bohlin's Gestational Age Acceleration (GAA) at birth and Intrinsic Epigenetic Age Acceleration (IEAA) throughout childhood. Additionally, quantile g-computation was utilized to assess the effect of the phthalate mixture on GAA at birth and IEAA throughout childhood.

RESULTS

We found a negative association between prenatal di (2-ethylhexyl) phthalate (DEHP) exposure and IEAA among males at age 7 (-0.62 years; 95% CI:-1.06 to -0.18), and a marginal negative association between the whole phthalate mixture and GAA among males at birth (-1.54 days, 95% CI: -2.79 to -0.28), while most other associations were nonsignificant.

CONCLUSIONS

Our results suggest that prenatal exposure to certain phthalates is associated with epigenetic aging in children. Additionally, our findings suggest that the influence of prenatal exposures on epigenetic age may only manifest during specific periods of child development, and studies relying on DNAm measurements solely from cord blood or single time points may overlook potential relationships.

Epigenome-Wide Associations of Placental DNA Methylation and Behavioral and Emotional Difficulties in Children at 3 Years of Age

The placenta is a key organ for fetal and brain development. Its epigenome can be regarded as a biochemical record of the prenatal environment and a potential mechanism of its association with the future health of the fetus. We investigated associations between placental DNA methylation levels and child behavioral and emotional difficulties, assessed at 3 years of age using the Strengths and Difficulties Questionnaire (SDQ) in 441 mother-child dyads from the EDEN cohort. Hypothesis-driven and exploratory analyses (on differentially methylated probes (EWAS) and regions (DMR)) were adjusted for confounders, technical factors, and cell composition estimates, corrected for multiple comparisons, and stratified by child sex. Hypothesis-driven analyses showed an association of cg26703534 (AHRR) with emotional symptoms, and exploratory analyses identified two probes, cg09126090 (intergenic region) and cg10305789 (PPP1R16B), as negatively associated with peer relationship problems, as well as 33 DMRs, mostly positively associated with at least one of the SDQ subscales. Among girls, most associations were seen with emotional difficulties, whereas in boys, DMRs were as much associated with emotional than behavioral difficulties. This study provides the first evidence of associations between placental DNA methylation and child behavioral and emotional difficulties. Our results suggest sex-specific associations and might provide new insights into the mechanisms of neurodevelopment.

Genetic strategies for sex-biased persistence of gut microbes across human life

Although compositional variation in the gut microbiome during human development has been extensively investigated, strain-resolved dynamic changes remain to be fully uncovered. In the current study, shotgun metagenomic sequencing data of 12,415 fecal microbiomes from healthy individuals are employed for strain-level tracking of gut microbiota members to elucidate its evolving biodiversity across the human life span. This detailed longitudinal meta-analysis reveals host sex-related persistence of strains belonging to common, maternally-inherited species, such as Bifidobacterium bifidum and Bifidobacterium longum subsp. longum. Comparative genome analyses, coupled with experiments including intimate interaction between microbes and human intestinal cells, show that specific bacterial glycosyl hydrolases related to host-glycan metabolism may contribute to more efficient colonization in females compared to males. These findings point to an intriguing ancient sex-specific host-microbe coevolution driving the selective persistence in women of key microbial taxa that may be vertically passed on to the next generation.

DNA methylation markers for kidney function and progression of diabetic kidney disease

Epigenetic markers are potential biomarkers for diabetes and related complications. Using a prospective cohort from the Hong Kong Diabetes Register, we perform two independent epigenome-wide association studies to identify methylation markers associated with baseline estimated glomerular filtration rate (eGFR) and subsequent decline in kidney function (eGFR slope), respectively, in 1,271 type 2 diabetes subjects. Here we show 40 (30 previously unidentified) and eight (all previously unidentified) CpG sites individually reach epigenome-wide significance for baseline eGFR and eGFR slope, respectively. We also develop a multisite analysis method, which selects 64 and 37 CpG sites for baseline eGFR and eGFR slope, respectively. These models are validated in an independent cohort of Native Americans with type 2 diabetes. Our identified CpG sites are near genes enriched for functional roles in kidney diseases, and some show association with renal damage. This study highlights the potential of methylation markers in risk stratification of kidney disease among type 2 diabetes individuals.

The Epigenetics of Animal Personality

Animal personality, consistent individual differences in behaviour, is an important concept for understanding how individuals vary in how they cope with environmental challenges. In order to understand the evolutionary significance of animal personality, it is crucial to understand the underlying regulatory mechanisms. Epigenetic marks such as DNA methylation are hypothesised to play a major role in explaining variation in phenotypic changes in response to environmental alterations. Several characteristics of DNA methylation also align well with the concept of animal personality. In this review paper, we summarise the current literature on the role that molecular epigenetic mechanisms may have in explaining personality variation. We elaborate on the potential for epigenetic mechanisms to explain behavioural variation, behavioural development and temporal consistency in behaviour. We then suggest future routes for this emerging field and point to potential pitfalls that may be encountered. We conclude that a more inclusive approach is needed for studying the epigenetics of animal personality and that epigenetic mechanisms cannot be studied without considering the genetic background.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Sex differences in the intergenerational link between maternal and neonatal whole blood DNA methylation: a genome-wide analysis in 2 birth cohorts

BACKGROUND

The mother-child inheritance of DNA methylation (DNAm) variations could contribute to the inheritance of disease susceptibility across generations. However, no study has investigated patterns of mother-child associations in DNAm at the genome-wide scale. It remains unknown whether there are sex differences in mother-child DNAm associations.

RESULTS

Using genome-wide DNAm profiling data (721,331 DNAm sites, including 704,552 on autosomes and 16,779 on the X chromosome) of 396 mother-newborn pairs (54.5% male) from the Boston Birth Cohort, we found significant sex differences in mother-newborn correlations in genome-wide DNAm patterns (Spearman's rho = 0.91-0.98; p = 4.0 × 10-8), with female newborns having stronger correlations. Sex differences in correlations were attenuated but remained significant after excluding X-chromosomal DNAm sites (Spearman's rho = 0.91-0.98; p = 0.035). Moreover, 89,267 DNAm sites (12.4% of all analyzed, including 88,051 [12.5% of analyzed] autosomal and 1,216 [7.2% of analyzed] X-chromosomal sites) showed significant mother-newborn associations in methylation levels, and the top autosomal DNAm sites had high heritability than the genome-wide background (e.g., the top 100 autosomal DNAm sites had a medium h2 of 0.92). Additionally, significant interactions between newborn sex and methylation levels were observed for 11 X-chromosomal and 4 autosomal DNAm sites that were mapped to genes that have been associated with sex-specific disease/traits or early development (e.g., EFHC2, NXY, ADCYAP1R1, and BMP4). Finally, 18,769 DNAm sites (14,482 [77.2%] on the X chromosome) showed mother-newborn differences in methylation levels that were significantly associated with newborn sex, and the top autosomal DNAm sites had relatively small heritability (e.g., the top 100 autosomal DNAm sites had a medium h2 of 0.23). These DNAm sites were mapped to 2,532 autosomal genes and 978 X-chromosomal genes with significant enrichment in pathways involved in neurodegenerative and psychological diseases, development, neurophysiological process, immune response, and sex-specific cancers. Replication analysis in the Isle of Wight birth cohort yielded consistent results.

CONCLUSION

In two independent birth cohorts, we demonstrated strong mother-newborn correlations in whole blood DNAm on both autosomes and ChrX, and such correlations vary substantially by sex. Future studies are needed to examine to what extent our findings contribute to developmental origins of pediatric and adult diseases with well-observed sex differences.

Data processing choices can affect findings in differential methylation analyses: an investigation using data from the LIMIT RCT

Objective

A wide array of methods exist for processing and analysing DNA methylation data. We aimed to perform a systematic comparison of the behaviour of these methods, using cord blood DNAm from the LIMIT RCT, in relation to detecting hypothesised effects of interest (intervention and pre-pregnancy maternal BMI) as well as effects known to be spurious, and known to be present.

Methods

DNAm data, from 645 cord blood samples analysed using Illumina 450K BeadChip arrays, were normalised using three different methods (with probe filtering undertaken pre- or post- normalisation). Batch effects were handled with a supervised algorithm, an unsupervised algorithm, or adjustment in the analysis model. Analysis was undertaken with and without adjustment for estimated cell type proportions. The effects estimated included intervention and BMI (effects of interest in the original study), infant sex and randomly assigned groups. Data processing and analysis methods were compared in relation to number and identity of differentially methylated probes, rankings of probes by p value and log-fold-change, and distributions of p values and log-fold-change estimates.

Results

There were differences corresponding to each of the processing and analysis choices. Importantly, some combinations of data processing choices resulted in a substantial number of spurious 'significant' findings. We recommend greater emphasis on replication and greater use of sensitivity analyses.

Beadchip technology to detect DNA methylation in mouse faithfully recapitulates whole-genome bisulfite sequencing

Aim: To facilitate wide-scale implementation of Illumina Mouse Methylation BeadChip (MMB) technology, array-based measurement of cytosine methylation was compared with the gold-standard assessment of DNA methylation by whole-genome bisulfite sequencing (WGBS). Methods: DNA methylation across two mouse strains (C57B6 and C3H) and both sexes was assessed using the MMB and compared with previously existing deep-coverage WGBS of mice of the same strain and sex. Results & conclusion: The findings demonstrated that 93.3-99.2% of sites had similar measurements of methylation across technologies and that differentially methylated cytosines and regions identified by each technology overlap and enrich for similar biological functions, suggesting that the MMB faithfully recapitulates the findings of WGBS.

A Higher Dysregulation Burden of Brain DNA Methylation in Female Patients Implicated in the Sex Bias of Schizophrenia

Sex differences are pervasive in schizophrenia (SCZ), but the extent and magnitude of DNA methylation (DNAm) changes underlying these differences remain uncharacterized. In this study, sex-stratified differential DNAm analysis was performed in postmortem brain samples from 117 SCZ and 137 controls, partitioned into discovery and replication datasets. Three differentially methylated positions (DMPs) were identified (adj. p  < 0.05) in females and 29 DMPs in males without overlap between them. Over 81% of these sex-stratified DMPs were directionally consistent between sexes but with different effect sizes. Down-sampling analysis revealed more DMPs in females than in males when the sample sizes matched. Females had higher DNAm levels in healthy individuals and larger magnitude of DNAm changes in patients than males. Despite similar proportions of female-related DMPs (fDMPs, 8%) being under genetic control compared with males (10%), significant enrichment of DMP-related SNPs in signals of genome-wide association studies was identified only in fDMPs. One DMP in each sex connected the SNPs and gene expression of CALHM1 in females and CCDC149 in males. PPI subnetworks revealed that both female- and male-related differential DNAm interacted with synapse-related dysregulation. Immune-related pathways were unique for females and neuron-related pathways were associated with males. This study reveals remarkable quantitative differences in DNAm-related sexual dimorphism in SCZ and that females have a higher dysregulation burden of SCZ-associated DNAm than males.

Association between child maltreatment and depressive symptoms in emerging adulthood: The mediating and moderating roles of DNA methylation

Prospective studies suggest that child maltreatment substantially increases the risk for depression in adulthood. However, the mechanisms underlying this association require further elucidation. In recent years, DNA methylation has emerged as a potential mechanism by which maltreatment experiences (a) could partly explain the emergence or aggravation of depressive symptoms (i.e., mediation) and/or (b) could increase (or decrease) the risk for depressive symptoms (i.e., moderation). The present study tested whether the methylation levels of nine candidate genes mediated and/or moderated the association between maltreatment experiences in childhood and depressive symptoms in emerging adulthood. The sample comprised 156 men aged between 18 and 35 years. Maltreatment experiences and depressive symptoms were assessed retrospectively using self-reported questionnaires. Methylation levels of nine candidate genes (COMT, FKBP5, IL6, IL10, MAOA, NR3C1, OXTR, SLC6A3 and SLC6A4), previously reported to be sensitive to early-life stress, were quantified from saliva samples. Maltreatment experiences in childhood were significantly associated with depressive symptoms in emerging adulthood. Both maltreatment experiences and depressive symptoms were associated with the methylation levels of two genomic sites, which cumulatively, but not individually, explained 16% of the association between maltreatment experiences in childhood and depressive symptoms in emerging adulthood. Moreover, maltreatment experiences in childhood interacted with the methylation levels of fourteen genomic sites, which cumulatively, but not individually, modulated the level of depressive symptoms in young male adults who were maltreated as children. However, none of these effects survived multiple testing correction. These findings bring attention to the cumulative effects of DNA methylation measured in several candidate genes on the risk of reporting depressive symptoms following maltreatment experiences in childhood. Nonetheless, future studies need to clarify the robustness of these putative cumulative effects in larger samples and longitudinal cohorts.

Sex-specific DNA methylation and associations with in utero tobacco smoke exposure at nuclear-encoded mitochondrial genes

Sex-linked differences in mitochondrial ATP production, enzyme activities, and reactive oxygen species generation have been reported in multiple tissue and cell types. While the effects of reproductive hormones underlie many of these differences, regulation of sexually dimorphic mitochondrial function has not been fully characterized. We hypothesized that sex-specific DNA methylation contributes to sex-specific expression of nuclear genes that influence mitochondrial function. Herein, we analysed DNA methylation data specifically focused on nuclear-encoded mitochondrial genes in 191 males and 190 females. We found 596 differentially methylated sites (DMSs) (FDR p < 0.05), corresponding to 324 genes, with at least a 1% difference in methylation between sexes. To investigate the potential functional significance, we utilized gene expression microarray data. Of the 324 genes containing DMSs, 17 showed differences in gene expression by sex. Particularly striking was that ATP5G2, encoding subunit C of ATP synthase, contains seven DMSs and exhibits a sex difference in expression (p = 0.04). Finally, we also found that alterations in DNA methylation associated with in utero tobacco smoke exposure were sex-specific in these nuclear-encoded mitochondrial genes. Interestingly, the level of sex differences in DNA methylation at nuclear-encoded mitochondrial genes and the level of methylation changes associated with smoke exposure were less prominent than that of other genes. This suggests more conservative regulation of DNA methylation at these nuclear-encoded mitochondrial genes as compared to others. Overall, our findings suggest that sex-specific DNA methylation may help establish sex differences in expression and function and that sex-specific alterations in DNA methylation in response to exposures could contribute to sex-variable toxicological responses.

Childhood Trauma and Epigenetics: State of the Science and Future

PURPOSE OF REVIEW

There is a great deal of interest regarding the biological embedding of childhood trauma and social exposures through epigenetic mechanisms, including DNA methylation (DNAm), but a comprehensive understanding has been hindered by issues of limited reproducibility between studies. This review presents a summary of the literature on childhood trauma and DNAm, highlights issues in the field, and proposes some potential solutions.

RECENT FINDINGS

Investigations of the associations between DNAm and childhood trauma are commonly performed using candidate gene approaches, specifically involving genes related to neurological and stress pathways. Childhood trauma is defined in a wide range of ways in several societal contexts. However, although variations in DNAm are frequently found in stress-related genes, unsupervised epigenome-wide association studies (EWAS) have shown limited reproducibility both between studies and in relating these changes to exposures. The reproducibility of childhood trauma DNAm studies, and the field of social epigenetics in general, may be improved by increasing sample sizes, standardizing variables, making use of effect size thresholds, collecting longitudinal and intervention samples, appropriately accounting for known confounding factors, and applying causal analysis wherever possible, such as "two-step epigenetic Mendelian randomization."

An examination of mediation by DNA methylation on birthweight differences induced by assisted reproductive technologies

BACKGROUND

Children born after assisted reproductive technologies (ART) differ in birthweight from those naturally conceived. It has been hypothesized that this might be explained by epigenetic mechanisms. We examined whether cord blood DNA methylation mediated the birthweight difference between 890 newborns conceived by ART (764 by fresh embryo transfer and 126 frozen thawed embryo transfer) and 983 naturally conceived newborns from the Norwegian Mother, Father, and Child Cohort Study (MoBa). DNA methylation was measured by the Illumina Infinium MethylationEPIC array. We conducted mediation analyses to assess whether differentially methylated CpGs mediated the differences in birthweight observed between: (1) fresh embryo transfer and natural conception and (2) frozen and fresh embryo transfer.

RESULTS

We observed a difference in birthweight between fresh embryo transfer and naturally conceived offspring of - 120 g. 44% (95% confidence interval [CI] 26% to 81%) of this difference in birthweight between fresh embryo transfer and naturally conceived offspring was explained by differences in methylation levels at four CpGs near LOXL1, CDH20, and DRC1. DNA methylation differences at two CpGs near PTGS1 and RASGRP4 jointly mediated 22% (95% CI 8.1% to 50.3%) of the birthweight differences between fresh and frozen embryo transfer.

CONCLUSION

Our findings suggest that DNA methylation is an important mechanism in explaining birthweight differences according to the mode of conception. Further research should examine how gene regulation at these loci influences fetal growth.

The immune factors driving DNA methylation variation in human blood

Epigenetic changes are required for normal development, yet the nature and respective contribution of factors that drive epigenetic variation in humans remain to be fully characterized. Here, we assessed how the blood DNA methylome of 884 adults is affected by DNA sequence variation, age, sex and 139 factors relating to life habits and immunity. Furthermore, we investigated whether these effects are mediated or not by changes in cellular composition, measured by deep immunophenotyping. We show that DNA methylation differs substantially between naïve and memory T cells, supporting the need for adjustment on these cell-types. By doing so, we find that latent cytomegalovirus infection drives DNA methylation variation and provide further support that the increased dispersion of DNA methylation with aging is due to epigenetic drift. Finally, our results indicate that cellular composition and DNA sequence variation are the strongest predictors of DNA methylation, highlighting critical factors for medical epigenomics studies.

Epigenome-wide association study of global cortical volumes in generation Scotland: Scottish family health study

A complex interplay of genetic and environmental risk factors influence global brain structural alterations associated with brain health and disease. Epigenome-wide association studies (EWAS) of global brain imaging phenotypes have the potential to reveal the mechanisms of brain health and disease and can lead to better predictive analytics through the development of risk scores.We perform an EWAS of global brain volumes in Generation Scotland using peripherally measured whole blood DNA methylation (DNAm) from two assessments, (i) at baseline recruitment, ~6 years prior to MRI assessment (N = 672) and (ii) concurrent with MRI assessment (N=565). Four CpGs at baseline were associated with global cerebral white matter, total grey matter, and whole-brain volume (Bonferroni p≤7.41×10-8, βrange = -1.46x10-6 to 9.59 × 10-7). These CpGs were annotated to genes implicated in brain-related traits, including psychiatric disorders, development, and ageing. We did not find significant associations in the meta-analysis of the EWAS of the two sets concurrent with imaging at the corrected level.These findings reveal global brain structural changes associated with DNAm measured ~6 years previously, indicating a potential role of early DNAm modifications in brain structure. Although concurrent DNAm was not associated with global brain structure, the nominally significant findings identified here present a rationale for future investigation of associations between DNA methylation and structural brain phenotypes in larger population-based samples.

DMRscaler: a scale-aware method to identify regions of differential DNA methylation spanning basepair to multi-megabase features

BACKGROUND

Pathogenic mutations in genes that control chromatin function have been implicated in rare genetic syndromes. These chromatin modifiers exhibit extraordinary diversity in the scale of the epigenetic changes they affect, from single basepair modifications by DNMT1 to whole genome structural changes by PRM1/2. Patterns of DNA methylation are related to a diverse set of epigenetic features across this full range of epigenetic scale, making DNA methylation valuable for mapping regions of general epigenetic dysregulation. However, existing methods are unable to accurately identify regions of differential methylation across this full range of epigenetic scale directly from DNA methylation data.

RESULTS

To address this, we developed DMRscaler, a novel method that uses an iterative windowing procedure to capture regions of differential DNA methylation (DMRs) ranging in size from single basepairs to whole chromosomes. We benchmarked DMRscaler against several DMR callers in simulated and natural data comparing XX and XY peripheral blood samples. DMRscaler was the only method that accurately called DMRs ranging in size from 100 bp to 1 Mb (pearson's r = 0.94) and up to 152 Mb on the X-chromosome. We then analyzed methylation data from rare-disease cohorts that harbor chromatin modifier gene mutations in NSD1, EZH2, and KAT6A where DMRscaler identified novel DMRs spanning gene clusters involved in development.

CONCLUSION

Taken together, our results show DMRscaler is uniquely able to capture the size of DMR features across the full range of epigenetic scale and identify novel, co-regulated regions that drive epigenetic dysregulation in human disease.

Toxicoepigenetics and Environmental Health: Challenges and Opportunities

The rapidly growing field of toxicoepigenetics seeks to understand how toxicant exposures interact with the epigenome to influence disease risk. Toxicoepigenetics is a promising field of environmental health research, as integrating epigenetics into the field of toxicology will enable a more thorough evaluation of toxicant-induced disease mechanisms as well as the elucidation of the role of the epigenome as a biomarker of exposure and disease and possible mediator of exposure effects. Likewise, toxicoepigenetics will enhance our knowledge of how environmental exposures, lifestyle factors, and diet interact to influence health. Ultimately, an understanding of how the environment impacts the epigenome to cause disease may inform risk assessment, permit noninvasive biomonitoring, and provide potential opportunities for therapeutic intervention. However, the translation of research from this exciting field into benefits for human and animal health presents several challenges and opportunities. Here, we describe four significant areas in which we see opportunity to transform the field and improve human health by reducing the disease burden caused by environmental exposures. These include (1) research into the mechanistic role for epigenetic change in environment-induced disease, (2) understanding key factors influencing vulnerability to the adverse effects of environmental exposures, (3) identifying appropriate biomarkers of environmental exposures and their associated diseases, and (4) determining whether the adverse effects of environment on the epigenome and human health are reversible through pharmacologic, dietary, or behavioral interventions. We then highlight several initiatives currently underway to address these challenges.

Interrelationships and determinants of aging biomarkers in cord blood

BACKGROUND

Increasing evidence supports the concept of prenatal programming as an early factor in the aging process. DNA methylation age (DNAm age), global genome-wide DNA methylation (global methylation), telomere length (TL), and mitochondrial DNA content (mtDNA content) have independently been shown to be markers of aging, but their interrelationship and determinants at birth remain uncertain.

METHODS

We assessed the inter-correlation between the aging biomarkers DNAm age, global methylation, TL and mtDNA content using Pearson's correlation in 190 cord blood samples of the ENVIRONAGE birth cohort. TL and mtDNA content was measured via qPCR, while the DNA methylome was determined using the human 450K methylation Illumina microarray. Subsequently, DNAm age was calculated according to Horvath's epigenetic clock, and mean global, promoter, gene-body, and intergenic DNA methylation were determined. Path analysis, a form of structural equation modeling, was performed to disentangle the complex causal relationships among the aging biomarkers and their potential determinants.

RESULTS

DNAm age was inversely correlated with global methylation (r = -0.64, p < 0.001) and mtDNA content (r = - 0.16, p = 0.027). Cord blood TL was correlated with mtDNA content (r = 0.26, p < 0.001) but not with global methylation or DNAm age. Path analysis showed the strongest effect for global methylation on DNAm age with a decrease of 0.64 standard deviations (SD) in DNAm age for each SD (0.01%) increase in global methylation (p < 0.001). Among the applied covariates, newborn sex and season of delivery were the strongest determinants of aging biomarkers.

CONCLUSIONS

We provide insight into molecular aging signatures at the start of life, including their interrelations and determinants, showing that cord blood DNAm age is inversely associated with global methylation and mtDNA content but not with newborn telomere length. Our findings demonstrate that cord blood TL and DNAm age relate to different pathways/mechanisms of biological aging and can be influenced by environmental factors already at the start of life. These findings are relevant for understanding fetal programming and for the early prevention of noncommunicable diseases.

Prenatal and birth associations of epigenetic gestational age acceleration in the Center for the Health Assessment of Mothers and Children of Salinas (CHAMACOS) cohort

Gestational age (GA) is an important determinant of child health and disease risk. Two epigenetic GA clocks have been developed using DNA methylation (DNAm) patterns in cord blood. We investigate the accuracy of GA clocks and determinants of epigenetic GA acceleration (GAA), a biomarker of biological ageing. We hypothesize that prenatal and birth characteristics are associated with altered GAA, thereby disrupting foetal biological ageing. We examined 372 mother-child pairs from the Center for the Health Assessment of Mothers and Children of Salinas study of primarily Latino farmworkers in California. Chronological GA was robustly correlated with epigenetic GA (DNAm GA) estimated by the Knight (r = 0.48, p < 2.2x10^-16) and Bohlin clocks (r = 0.67, p < 2.2x10^-16) using the Illumina 450K array in cord blood samples collected at birth. GA clock performance was robust, though slightly lower, using DNAm profiles from the Illumina EPIC array in a smaller subsample (Knight: r = 0.39, p < 3.5x10^-5; Bohlin: r = 0.60, p < 7.7x10^-12). After adjusting for confounders, high maternal serum triglyceride levels (Bohlin: β = -0.01 days per mg/dL, p = 0.03), high maternal serum lipid levels (Bohlin: β = -4.31x10^-3 days per mg/dL, p = 0.04), preterm delivery (Bohlin: β = -4.03 days, p = 9.64x10^-4), greater maternal parity (Knight: β = -4.07 days, p = 0.01; Bohlin: β = -2.43 days, p = 0.01), and male infant sex (Knight: β = -3.15 days, p = 3.10x10^-3) were associated with decreased GAA.Prenatal and birth characteristics affect GAA in newborns. Understanding factors that accelerate or delay biological ageing at birth may identify early-life targets for disease prevention and improve ageing across the life-course. Future research should test the impact of altered GAA on the long-term burden of age-related diseases.

Tissue- and ethnicity-independent hypervariable DNA methylation states show evidence of establishment in the early human embryo

We analysed DNA methylation data from 30 datasets comprising 3474 individuals, 19 tissues and 8 ethnicities at CpGs covered by the Illumina450K array. We identified 4143 hypervariable CpGs ('hvCpGs') with methylation in the top 5% most variable sites across multiple tissues and ethnicities. hvCpG methylation was influenced but not determined by genetic variation, and was not linked to probe reliability, epigenetic drift, age, sex or cell heterogeneity effects. hvCpG methylation tended to covary across tissues derived from different germ-layers and hvCpGs were enriched for proximity to ERV1 and ERVK retrovirus elements. hvCpGs were also enriched for loci previously associated with periconceptional environment, parent-of-origin-specific methylation, and distinctive methylation signatures in monozygotic twins. Together, these properties position hvCpGs as strong candidates for studying how stochastic and/or environmentally influenced DNA methylation states which are established in the early embryo and maintained stably thereafter can influence life-long health and disease.

InterpolatedXY: a two-step strategy to normalize DNA methylation microarray data avoiding sex bias

MOTIVATION

Data normalization is an essential step to reduce technical variation within and between arrays. Due to the different karyotypes and the effects of X chromosome inactivation, females and males exhibit distinct methylation patterns on sex chromosomes; thus, it poses a significant challenge to normalize sex chromosome data without introducing bias. Currently, existing methods do not provide unbiased solutions to normalize sex chromosome data, usually, they just process autosomal and sex chromosomes indiscriminately.

RESULTS

Here, we demonstrate that ignoring this sex difference will lead to introducing artificial sex bias, especially for thousands of autosomal CpGs. We present a novel two-step strategy (interpolatedXY) to address this issue, which is applicable to all quantile-based normalization methods. By this new strategy, the autosomal CpGs are first normalized independently by conventional methods, such as funnorm or dasen; then the corrected methylation values of sex chromosome-linked CpGs are estimated as the weighted average of their nearest neighbors on autosomes. The proposed two-step strategy can also be applied to other non-quantile-based normalization methods, as well as other array-based data types. Moreover, we propose a useful concept: the sex explained fraction of variance, to quantitatively measure the normalization effect.

AVAILABILITY AND IMPLEMENTATION

The proposed methods are available by calling the function 'adjustedDasen' or 'adjustedFunnorm' in the latest wateRmelon package (https://github.com/schalkwyk/wateRmelon), with methods compatible with all the major workflows, including minfi.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Illuminating the Mechanisms Underlying Sex Differences in Cardiovascular Disease

Sex is a key risk factor for many types of cardiovascular disease. It is imperative to understand the mechanisms underlying sex differences to devise optimal preventive and therapeutic approaches for all individuals. Both biological sex (determined by sex chromosomes and gonadal hormones) and gender (social and cultural behaviors associated with femininity or masculinity) influence differences between men and women in disease susceptibility and pathology. Here, we focus on the application of experimental mouse models that elucidate the influence of 2 components of biological sex-sex chromosome complement (XX or XY) and gonad type (ovaries or testes). These models have revealed that in addition to well-known effects of gonadal hormones, sex chromosome complement influences cardiovascular risk factors, such as plasma cholesterol levels and adiposity, as well as the development of atherosclerosis and pulmonary hypertension. One mechanism by which sex chromosome dosage influences cardiometabolic traits is through sex-biased expression of X chromosome genes that escape X inactivation. These include chromatin-modifying enzymes that regulate gene expression throughout the genome. The identification of factors that determine sex-biased gene expression and cardiometabolic traits will expand our mechanistic understanding of cardiovascular disease processes and provide insight into sex differences that remain throughout the lifespan as gonadal hormone levels alter with age.

Large-scale placenta DNA methylation integrated analysis reveals fetal sex-specific differentially methylated CpG sites and regions

Although male–female differences in placental structure and function have been observed, little is understood about their molecular underpinnings. Here, we present a mega-analysis of 14 publicly available placenta DNA methylation (DNAm) microarray datasets to identify individual CpGs and regions associated with fetal sex. In the discovery dataset of placentas from full term pregnancies (N = 532 samples), 5212 CpGs met genome-wide significance (p < 1E−8) and were enriched in pathways such as keratinization (FDR p-value = 7.37E−14), chemokine activity (FDR p-value = 1.56E−2), and eosinophil migration (FDR p-value = 1.83E−2). Nine differentially methylated regions were identified (fwerArea < 0.1) including a region in the promoter of ZNF300 that showed consistent differential DNAm in samples from earlier timepoints in pregnancy and appeared to be driven predominately by effects in the trophoblast cell type. We describe the largest study of fetal sex differences in placenta DNAm performed to date, revealing genes and pathways characterizing sex-specific placenta function and health outcomes later in life.

Challenges in Analyzing Functional Epigenetic Data in Perspective of Adolescent Psychiatric Health

The formative period of adolescence plays a crucial role in the development of skills and abilities for adulthood. Adolescents who are affected by mental health conditions are at risk of suicide and social and academic impairments. Gene–environment complementary contributions to the molecular mechanisms involved in psychiatric disorders have emphasized the need to analyze epigenetic marks such as DNA methylation (DNAm) and non-coding RNAs. However, the large and diverse bioinformatic and statistical methods, referring to the confounders of the statistical models, application of multiple-testing adjustment methods, questions regarding the correlation of DNAm across tissues, and sex-dependent differences in results, have raised challenges regarding the interpretation of the results. Based on the example of generalized anxiety disorder (GAD) and depressive disorder (MDD), we shed light on the current knowledge and usage of methodological tools in analyzing epigenetics. Statistical robustness is an essential prerequisite for a better understanding and interpretation of epigenetic modifications and helps to find novel targets for personalized therapeutics in psychiatric diseases.

Maternal prenatal depressive symptoms and toddler behavior: an umbilical cord blood epigenome-wide association study

Children of mothers with prenatal depressive symptoms (PND) have a higher risk of behavioral problems; fetal programming through DNA methylation is a possible underlying mechanism. This study investigated DNA methylation in cord blood to identify possible "at birth" signatures that may indicate susceptibility to behavioral problems at 18 months of age. Cord blood was collected from 256 children of mothers who had self-reported on symptoms of depression during pregnancy and the behavior of their child at 18 months of age. Whole genome DNA methylation was assessed using Illumina MethylationEPIC assay. The mother and child pairs were categorized into four groups, based on both self-reported depressive symptoms, PND or Healthy control (HC), and scores from the Child Behavior checklist (high or low for internalizing, externalizing, and total scores). Adjustments were made for batch effects, cell-type, and clinical covariates. Differentially methylated sites were identified using Kruskal-Wallis test, and Benjamini-Hochberg adjusted p values < 0.05 were considered significant. The analysis was also stratified by sex of the child. Among boys, we observed higher and correlated DNA methylation of one CpG-site in the promoter region of TPP1 in the HC group, with high externalizing scores compared to HC with low externalizing scores. Boys in the PND group showed lower DNA methylation in NUDT15 among those with high, compared to low, internalizing scores; the DNA methylation levels of CpGs in this gene were positively correlated with the CBCL scores. Hence, the differentially methylated CpG sites could be of interest for resilience, regardless of maternal mental health during pregnancy. The findings are in a relatively healthy study cohort, thus limiting the possibility of detecting strong effects associated with behavioral difficulties. This is the first investigation of cord blood DNA methylation signs of fetal programming of PND on child behavior at 18 months of age and thus calls for independent replications.

Differential Methylation Analysis of Suicidal Ideation Severity in Schizophrenia with the Illumina MethylationEPIC Array

There is a multitude of factors that makes difficult to identify those at risk for suicide, especially among schizophrenia patients. Suicide cannot be explained by genetics alone, therefore epigenetic mechanisms including DNA methylation are thought to play a role. DNA methylation could be a valuable tool in helping predict those at-risk individuals. This cross-sectional study comprised 112 subjects diagnosed with schizophrenia spectrum disorders, and were grouped according to the current suicidal ideation severity. DNA methylation across the genome was measured with the Infinium® MethylationEPIC BeadChip. We utilized the dmpFinder and bumphunter functions within the Bioconductor minfi package to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs), respectively. Following quality control, we removed one sample from the analysis and reported the most significant DMPs and DMRs associated with suicidal ideation severity. All positions and regions identified in this analysis were only found to have suggestive levels of significance at the genome-wide level. The present study was one of the first to investigate genome-wide methylation and suicidal ideation severity. While there were many strengths of our study, including investigating both differentially methylated positions and regions, further larger-scale studies are necessary to replicate, support, and validate our findings presented here.

Epigenetic Signatures Discriminate Patients With Primary Sclerosing Cholangitis and Ulcerative Colitis From Patients With Ulcerative Colitis

Background

Primary sclerosing cholangitis (PSC) is a chronic inflammatory liver disease affecting the intra- and extrahepatic bile ducts, and is strongly associated with ulcerative colitis (UC). In this study, we explored the peripheral blood DNA methylome and its immune cell composition in patients with PSC-UC, UC, and healthy controls (HC) with the aim to develop a predictive assay in distinguishing patients with PSC-UC from those with UC alone.

Methods

The peripheral blood DNA methylome of male patients with PSC and concomitant UC, UC and HCs was profiled using the Illumina HumanMethylation Infinium EPIC BeadChip (850K) array. Differentially methylated CpG position (DMP) and region (DMR) analyses were performed alongside gradient boosting classification analyses to discern PSC-UC from UC patients. As observed differences in the DNA methylome could be the result of differences in cellular populations, we additionally employed mass cytometry (CyTOF) to characterize the immune cell compositions.

Results

Genome wide methylation analysis did not reveal large differences between PSC-UC and UC patients nor HCs. Nonetheless, using gradient boosting we were capable of discerning PSC-UC from UC with an area under the receiver operator curve (AUROC) of 0.80. Four CpG sites annotated to the NINJ2 gene were found to strongly contribute to the predictive performance. While CyTOF analyses corroborated the largely similar blood cell composition among patients with PSC-UC, UC and HC, a higher abundance of myeloid cells was observed in UC compared to PSC-UC patients.

Conclusion

DNA methylation enables discerning PSC-UC from UC patients, with a potential for biomarker development.

A robust mean and variance test with application to high-dimensional phenotypes

Most studies of continuous health-related outcomes examine differences in mean levels (location) of the outcome by exposure. However, identifying effects on the variability (scale) of an outcome, and combining tests of mean and variability (location-and-scale), could provide additional insights into biological mechanisms. A joint test could improve power for studies of high-dimensional phenotypes, such as epigenome-wide association studies of DNA methylation at CpG sites. One possible cause of heterogeneity of variance is a variable interacting with exposure in its effect on outcome, so a joint test of mean and variability could help in the identification of effect modifiers. Here, we review a scale test, based on the Brown-Forsythe test, for analysing variability of a continuous outcome with respect to both categorical and continuous exposures, and develop a novel joint location-and-scale score (JLSsc) test. These tests were compared to alternatives in simulations and used to test associations of mean and variability of DNA methylation with gender and gestational age using data from the Accessible Resource for Integrated Epigenomics Studies (ARIES). In simulations, the Brown-Forsythe and JLSsc tests retained correct type I error rates when the outcome was not normally distributed in contrast to the other approaches tested which all had inflated type I error rates. These tests also identified > 7500 CpG sites for which either mean or variability in cord blood methylation differed according to gender or gestational age. The Brown-Forsythe test and JLSsc are robust tests that can be used to detect associations not solely driven by a mean effect.

Permutation-based significance analysis reduces the type 1 error rate in bisulfite sequencing data analysis of human umbilical cord blood samples

DNA methylation patterns are largely established in-utero and might mediate the impacts of in-utero conditions on later health outcomes. Associations between perinatal DNA methylation marks and pregnancy-related variables, such as maternal age and gestational weight gain, have been earlier studied with methylation microarrays, which typically cover less than 2% of human CpG sites. To detect such associations outside these regions, we chose the bisulphite sequencing approach. We collected and curated clinical data on 200 newborn infants; whose umbilical cord blood samples were analysed with the reduced representation bisulphite sequencing (RRBS) method. A generalized linear mixed-effects model was fit for each high coverage CpG site, followed by spatial and multiple testing adjustment of P values to identify differentially methylated cytosines (DMCs) and regions (DMRs) associated with clinical variables, such as maternal age, mode of delivery, and birth weight. Type 1 error rate was then evaluated with a permutation analysis. We discovered a strong inflation of spatially adjusted P values through the permutation analysis, which we then applied for empirical type 1 error control. The inflation of P values was caused by a common method for spatial adjustment and DMR detection, implemented in tools comb-p and RADMeth. Based on empirically estimated significance thresholds, very little differential methylation was associated with any of the studied clinical variables, other than sex. With this analysis workflow, the sex-associated differentially methylated regions were highly reproducible across studies, technologies, and statistical models.

Modeling dependency structures in 450k DNA methylation data

MOTIVATION

DNA methylation has been shown to be spatially dependent across chromosomes. Previous studies have focused on the influence of genomic context on the dependency structure, while not considering differences in dependency structure between individuals.

RESULTS

We modeled spatial dependency with a flexible framework to quantify the dependency structure, focusing on inter-individual differences by exploring the association between dependency parameters and technical and biological variables. The model was applied to a subset of the Finnish Twin Cohort study (N = 1611 individuals). The estimates of the dependency parameters varied considerably across individuals, but were generally consistent across chromosomes within individuals. The variation in dependency parameters was associated with bisulfite conversion plate, zygosity, sex and age. The age differences presumably reflect accumulated environmental exposures and/or accumulated small methylation differences caused by stochastic mitotic events, establishing recognizable, individual patterns more strongly seen in older individuals.

AVAILABILITY AND IMPLEMENTATION

The twin dataset used in the current study are located in the Biobank of the National Institute for Health and Welfare, Finland. All the biobanked data are publicly available for use by qualified researchers following a standardized application procedure (https://thl.fi/en/web/thl-biobank/for-researchers). A R-script for fitting the dependency structure to publicly available DNA methylation data with the software used in this article is provided in supplementary data.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Variation in Genotype and DNA Methylation Patterns Based on Alcohol Use and CVD in the Korean Genome and Epidemiology Study (KoGES)

Alcohol consumption can increase the risk of chronic diseases, such as myocardial infarction, coronary artery disease, hyperlipidemia, and hypertension. We aimed to assess the association between genotype, DNA methylation patterns, alcohol consumption, and chronic diseases in Korean population. We analyzed 8840 subjects for genotypes and 446 for DNA methylation among the 9351 subjects from the Korean Genome and Epidemiology Study (KoGES). We further divided both groups into two sub-groups according to the presence/absence of chronic diseases. We selected genes whose methylation varied significantly with alcohol consumption, and visualized genotype and DNA methylation patterns specific to each group. Genome-wide association study (GWAS) revealed single nucleotide polymorphisms (SNPs) rs2074356 and rs11066280 in HECT domain E3 ubiquitin protein ligase 4 (HECTD4) to be significantly associated with alcohol consumption in both the presence. The rs12229654 genotype also displayed significantly different patterns with alcohol consumption. Furthermore, we retrieved differentially methylated regions (DMRs) from four groups based on sex and chronic diseases and compared them by drinking status. In genotype analysis, cardiovascular diseases (CVDs) showed a higher proportion in drinker than in non-drinker, but not in DMR analysis. Additionally, we analyzed the enriched Gene Ontology terms and Kyoto Gene and Genome Encyclopedia (KEGG) pathways and visualized the network, heatmap, and upset plot. We show that the pattern of DNA methylation associated with CVD is strongly influenced by alcoholism. Overall, this study identified genetic and epigenetic variants influenced by alcohol consumption and chronic diseases.

Estimating and accounting for unobserved covariates in high-dimensional correlated data

Many high dimensional and high-throughput biological datasets have complex sample correlation structures, which include longitudinal and multiple tissue data, as well as data with multiple treatment conditions or related individuals. These data, as well as nearly all high-throughput 'omic' data, are influenced by technical and biological factors unknown to the researcher, which, if unaccounted for, can severely obfuscate estimation of and inference on the effects of interest. We therefore developed CBCV and CorrConf: provably accurate and computationally efficient methods to choose the number of and estimate latent confounding factors present in high dimensional data with correlated or nonexchangeable residuals. We demonstrate each method's superior performance compared to other state of the art methods by analyzing simulated multi-tissue gene expression data and identifying sex-associated DNA methylation sites in a real, longitudinal twin study.

Meta-analysis of epigenome-wide association studies in newborns and children show widespread sex differences in blood DNA methylation

BACKGROUND

Among children, sex-specific differences in disease prevalence, age of onset, and susceptibility have been observed in health conditions including asthma, immune response, metabolic health, some pediatric and adult cancers, and psychiatric disorders. Epigenetic modifications such as DNA methylation may play a role in the sexual differences observed in diseases and other physiological traits.

METHODS

We performed a meta-analysis of the association of sex and cord blood DNA methylation at over 450,000 CpG sites in 8438 newborns from 17 cohorts participating in the Pregnancy And Childhood Epigenetics (PACE) Consortium. We also examined associations of child sex with DNA methylation in older children ages 5.5-10 years from 8 cohorts (n = 4268).

RESULTS

In newborn blood, sex was associated at Bonferroni level significance with differences in DNA methylation at 46,979 autosomal CpG sites (p < 1.3 × 10-7) after adjusting for white blood cell proportions and batch. Most of those sites had lower methylation levels in males than in females. Of the differentially methylated CpG sites identified in newborn blood, 68% (31,727) met look-up level significance (p < 1.1 × 10-6) in older children and had methylation differences in the same direction.

CONCLUSIONS

This is a large-scale meta-analysis examining sex differences in DNA methylation in newborns and older children. Expanding upon previous studies, we replicated previous findings and identified additional autosomal sites with sex-specific differences in DNA methylation. Differentially methylated sites were enriched in genes involved in cancer, psychiatric disorders, and cardiovascular phenotypes.