An evaluation of processing methods for HumanMethylation450 BeadChip data

Methods in DNA methylation array dataset analysis: A review

Understanding the intricate relationships between gene expression levels and epigenetic modifications in a genome is crucial to comprehending the pathogenic mechanisms of many diseases. With the advancement of DNA Methylome Profiling techniques, the emphasis on identifying Differentially Methylated Regions (DMRs/DMGs) has become crucial for biomarker discovery, offering new insights into the etiology of illnesses. This review surveys the current state of computational tools/algorithms for the analysis of microarray-based DNA methylation profiling datasets, focusing on key concepts underlying the diagnostic/prognostic CpG site extraction. It addresses methodological frameworks, algorithms, and pipelines employed by various authors, serving as a roadmap to address challenges and understand changing trends in the methodologies for analyzing array-based DNA methylation profiling datasets derived from diseased genomes. Additionally, it highlights the importance of integrating gene expression and methylation datasets for accurate biomarker identification, explores prognostic prediction models, and discusses molecular subtyping for disease classification. The review also emphasizes the contributions of machine learning, neural networks, and data mining to enhance diagnostic workflow development, thereby improving accuracy, precision, and robustness.

Predicting type 2 diabetes via machine learning integration of multiple omics from human pancreatic islets

Type 2 diabetes (T2D) is the fastest growing non-infectious disease worldwide. Impaired insulin secretion from pancreatic beta-cells is a hallmark of T2D, but the mechanisms behind this defect are insufficiently characterized. Integrating multiple layers of biomedical information, such as different Omics, may allow more accurate understanding of complex diseases such as T2D. Our aim was to explore and use Machine Learning to integrate multiple sources of biological/molecular information (multiOmics), in our case RNA-sequening, DNA methylation, SNP and phenotypic data from islet donors with T2D and non-diabetic controls. We exploited Machine Learning to perform multiOmics integration of DNA methylation, expression, SNPs, and phenotypes from pancreatic islets of 110 individuals, with ~ 30% being T2D cases. DNA methylation was analyzed using Infinium MethylationEPIC array, expression was analyzed using RNA-sequencing, and SNPs were analyzed using HumanOmniExpress arrays. Supervised linear multiOmics integration via DIABLO based on Partial Least Squares (PLS) achieved an accuracy of 91 ± 15% of T2D prediction with an area under the curve of 0.96 ± 0.08 on the test dataset after cross-validation. Biomarkers identified by this multiOmics integration, including SACS and TXNIP DNA methylation, OPRD1 and RHOT1 expression and a SNP annotated to ANO1, provide novel insights into the interplay between different biological mechanisms contributing to T2D. This Machine Learning approach of multiOmics cross-sectional data from human pancreatic islets achieved a promising accuracy of T2D prediction, which may potentially find broad applications in clinical diagnostics. In addition, it delivered novel candidate biomarkers for T2D and links between them across the different Omics.

Toward a more holistic approach to the study of exposures and child outcomes

Aim: The current work was designed to demonstrate the application of the exposome framework in examining associations between exposures and children's long-term neurodevelopmental and behavioral outcomes. Methods: Longitudinal data were collected from birth through age 6 from 402 preterm infants. Three statistical methods were utilized to demonstrate the exposome framework: exposome-wide association study, cumulative exposure and machine learning models, with and without epigenetic data. Results: Each statistical approach answered a distinct research question regarding the impact of exposures on longitudinal child outcomes. Findings highlight associations between exposures, epigenetics and executive function. Conclusion: Findings demonstrate how an exposome-based approach can be utilized to understand relationships between internal (e.g., DNA methylation) and external (e.g., prenatal risk) exposures and long-term developmental outcomes in preterm children.

Epigenome-wide association study identifies neonatal DNA methylation associated with two-year attention problems in children born very preterm

Prior research has identified epigenetic predictors of attention problems in school-aged children but has not yet investigated these in young children, or children at elevated risk of attention problems due to preterm birth. The current study evaluated epigenome-wide associations between neonatal DNA methylation and attention problems at age 2 years in children born very preterm. Participants included 441 children from the Neonatal Neurobehavior and Outcomes in Very Preterm Infants (NOVI) Study, a multi-site study of infants born < 30 weeks gestational age. DNA methylation was measured from buccal swabs collected at NICU discharge using the Illumina MethylationEPIC Bead Array. Attention problems were assessed at 2 years of adjusted age using the attention problems subscale of the Child Behavior Checklist (CBCL). After adjustment for multiple testing, DNA methylation at 33 CpG sites was associated with child attention problems. Differentially methylated CpG sites were located in genes previously linked to physical and mental health, including several genes associated with ADHD in prior epigenome-wide and genome-wide association studies. Several CpG sites were located in genes previously linked to exposure to prenatal risk factors in the NOVI sample. Neonatal epigenetics measured at NICU discharge could be useful in identifying preterm children at risk for long-term attention problems and related psychiatric disorders, who could benefit from early prevention and intervention efforts.

Epigenetic signature of human immune aging in the GESTALT study

Age-associated DNA methylation in blood cells convey information on health status. However, the mechanisms that drive these changes in circulating cells and their relationships to gene regulation are unknown. We identified age-associated DNA methylation sites in six purified blood-borne immune cell types (naive B, naive CD4+ and CD8+ T cells, granulocytes, monocytes, and NK cells) collected from healthy individuals interspersed over a wide age range. Of the thousands of age-associated sites, only 350 sites were differentially methylated in the same direction in all cell types and validated in an independent longitudinal cohort. Genes close to age-associated hypomethylated sites were enriched for collagen biosynthesis and complement cascade pathways, while genes close to hypermethylated sites mapped to neuronal pathways. In silico analyses showed that in most cell types, the age-associated hypo- and hypermethylated sites were enriched for ARNT (HIF1β) and REST transcription factor (TF) motifs, respectively, which are both master regulators of hypoxia response. To conclude, despite spatial heterogeneity, there is a commonality in the putative regulatory role with respect to TF motifs and histone modifications at and around these sites. These features suggest that DNA methylation changes in healthy aging may be adaptive responses to fluctuations of oxygen availability.

DNA methylation-based classifier differentiates intrahepatic pancreato-biliary tumours

BACKGROUND

Differentiating intrahepatic cholangiocarcinomas (iCCA) from hepatic metastases of pancreatic ductal adenocarcinoma (PAAD) is challenging. Both tumours have similar morphological and immunohistochemical pattern and share multiple driver mutations. We hypothesised that DNA methylation-based machine-learning algorithms may help perform this task.

METHODS

We assembled genome-wide DNA methylation data for iCCA (n = 259), PAAD (n = 431), and normal bile duct (n = 70) from publicly available sources. We split this cohort into a reference (n = 399) and a validation set (n = 361). Using the reference cohort, we trained three machine learning models to differentiate between these entities. Furthermore, we validated the classifiers on the technical validation set and used an internal cohort (n = 72) to test our classifier.

FINDINGS

On the validation cohort, the neural network, support vector machine, and the random forest classifiers reached accuracies of 97.68%, 95.62%, and 96.5%, respectively. Filtering by anomaly detection and thresholds improved the accuracy to 99.07% (37 samples excluded by filtering), 96.22% (17 samples excluded), and 100% (44 samples excluded) for the neural network, support vector machine and random forest, respectively. Because of best balance between accuracy and number of predictable cases we tested the neural network with applied filters on the in-house cohort, obtaining an accuracy of 95.45%.

INTERPRETATION

We developed a classifier that can differentiate between iCCAs, intrahepatic metastases of a PAAD, and normal bile duct tissue with high accuracy. This tool can be used for improving the diagnosis of pancreato-biliary cancers of the liver.

FUNDING

This work was supported by Berlin Institute of Health (JCS Program), DKTK Berlin (Young Investigator Grant 2022), German Research Foundation (493697503 and 314905040 - SFB/TRR 209 Liver Cancer B01), and German Cancer Aid (70113922).

Who's afraid of the X? Incorporating the X and Y chromosomes into the analysis of DNA methylation array data

BACKGROUND

Many human disease phenotypes manifest differently by sex, making the development of methods for incorporating X and Y-chromosome data into analyses vital. Unfortunately, X and Y chromosome data are frequently excluded from large-scale analyses of the human genome and epigenome due to analytical complexity associated with sex chromosome dosage differences between XX and XY individuals, and the impact of X-chromosome inactivation (XCI) on the epigenome. As such, little attention has been given to considering the methods by which sex chromosome data may be included in analyses of DNA methylation (DNAme) array data.

RESULTS

With Illumina Infinium HumanMethylation450 DNAme array data from 634 placental samples, we investigated the effects of probe filtering, normalization, and batch correction on DNAme data from the X and Y chromosomes. Processing steps were evaluated in both mixed-sex and sex-stratified subsets of the analysis cohort to identify whether including both sexes impacted processing results. We found that identification of probes that have a high detection p-value, or that are non-variable, should be performed in sex-stratified data subsets to avoid over- and under-estimation of the quantity of probes eligible for removal, respectively. All normalization techniques investigated returned X and Y DNAme data that were highly correlated with the raw data from the same samples. We found no difference in batch correction results after application to mixed-sex or sex-stratified cohorts. Additionally, we identify two analytical methods suitable for XY chromosome data, the choice between which should be guided by the research question of interest, and we performed a proof-of-concept analysis studying differential DNAme on the X and Y chromosome in the context of placental acute chorioamnionitis. Finally, we provide an annotation of probe types that may be desirable to filter in X and Y chromosome analyses, including probes in repetitive elements, the X-transposed region, and cancer-testis gene promoters.

CONCLUSION

While there may be no single "best" approach for analyzing DNAme array data from the X and Y chromosome, analysts must consider key factors during processing and analysis of sex chromosome data to accommodate the underlying biology of these chromosomes, and the technical limitations of DNA methylation arrays.

Comparative analysis of the DNA methylation landscape in CD4, CD8, and B memory lineages

BACKGROUND

There is considerable evidence that epigenetic mechanisms and DNA methylation are critical drivers of immune cell lineage differentiation and activation. However, there has been limited coordinated investigation of common epigenetic pathways among cell lineages. Further, it remains unclear if long-lived memory cell subtypes differentiate distinctly by cell lineages.

RESULTS

We used the Illumina EPIC array to investigate the consistency of DNA methylation in B cell, CD4 T, and CD8 T naïve and memory cells states. In the process of naïve to memory activation across the three lineages, we identify considerable shared epigenetic regulation at the DNA level for immune memory generation. Further, in central to effector memory differentiation, our analyses revealed specific CpG dinucleotides and genes in CD4 T and CD8 T cells with DNA methylation changes. Finally, we identified unique DNA methylation patterns in terminally differentiated effector memory (TEMRA) CD8 T cells compared to other CD8 T memory cell subtypes.

CONCLUSIONS

Our data suggest that epigenetic alterations are widespread and essential in generating human lymphocyte memory. Unique profiles are involved in methylation changes that accompany memory genesis in the three subtypes of lymphocytes.

CpG Site-Based Signature Predicts Survival of Colorectal Cancer

BACKGROUND

A critical unmet medical need in clinical management of colorectal cancer (CRC) pivots around lack of noninvasive and or minimally invasive techniques for early diagnosis and prognostic prediction of clinical outcomes. Because DNA methylation can capture the regulatory landscape of tumors and can be measured in body fluids, it provides unparalleled opportunities for the discovery of early diagnostic and prognostics markers predictive of clinical outcomes. Here we investigated use of DNA methylation for the discovery of potential clinically actionable diagnostic and prognostic markers for predicting survival in CRC.

METHODS

We analyzed DNA methylation patterns between tumor and control samples to discover signatures of CpG sites and genes associated with CRC and predictive of survival. We conducted functional analysis to identify molecular networks and signaling pathways driving clinical outcomes.

RESULTS

We discovered a signature of aberrantly methylated genes associated with CRC and a signature of thirteen (13) CpG sites predictive of survival. We discovered molecular networks and signaling pathways enriched for CpG sites likely to drive clinical outcomes.

CONCLUSIONS

The investigation revealed that CpG sites can predict survival in CRC and that DNA methylation can capture the regulatory state of tumors through aberrantly methylated molecular networks and signaling pathways.

Traffic-related air pollution and genome-wide DNA methylation: A randomized, crossover trial

BACKGROUND

Traffic-related air pollution (TRAP) has been associated with changes in gene-specific DNA methylation. However, few studies have investigated impact of TRAP exposure on genome-wide DNA methylation in circulating blood of human.

OBJECTIVE

To explore the association between TRAP exposure and genome-wide DNA methylation.

METHODS

We conducted a randomized, crossover exposure trial among 35 healthy adults in Shanghai, China. All subjects were randomly allocated to a traffic-free park or a main road for consecutive 4 h, respectively. Blood genome-wide DNA methylation after each exposure session was measured by the Infinium Methylation EPIC BeadChip (850K). The differentially methylated CpGs loci associated with TRAP exposure were identified using linear mixed-effect model.

RESULTS

The average concentrations of traffic-related air pollutants including black carbon, ultrafine particles, carbon dioxide, and nitrogen dioxide were 2-3 times higher in the road compared to those in the park. Methylation levels of 68 CpG loci were significantly changed (false discovery rate < 0.05) following TRAP exposure, among which 49 were hypermethylated and 19 were hypomethylated. The annotated genes based on the differential CpGs loci were related to pathways in cardiovascular signaling, cytokine signaling, immune response, nervous system signaling, and metabolism.

CONCLUSIONS

We found that TRAP exposure was associated with DNA methylation in dozens of genes concerning cardiometabolic health. This trial for the first-time profiled genome-wide methylation changes induced by TRAP exposure using the 850K assay, providing epigenetic insights in understanding the cardiometabolic effects of TRAP exposure.

Updates to data versions and analytic methods influence the reproducibility of results from epigenome-wide association studies

Biomedical research has grown increasingly cooperative through the sharing of consortia-level epigenetic data. Since consortia preprocess data prior to distribution, new processing pipelines can lead to different versions of the same dataset. Similarly, analytic frameworks evolve to incorporate cutting-edge methods and best practices. However, it remains unknown how different data and analytic versions alter the results of epigenome-wide analyses, which could influence the replicability of epigenetic associations. Thus, we assessed the impact of these changes using data from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. We analysed DNA methylation from two data versions, processed using separate preprocessing and analytic pipelines, examining associations between seven childhood adversities or prenatal smoking exposure and DNA methylation at age 7. We performed two sets of analyses: (1) epigenome-wide association studies (EWAS); (2) Structured Life Course Modelling Approach (SLCMA), a two-stage method that models time-dependent effects. SLCMA results were also compared across two analytic versions. Data version changes impacted both EWAS and SLCMA analyses, yielding different associations at conventional p-value thresholds. However, the magnitude and direction of associations was generally consistent between data versions, regardless of p-values. Differences were especially apparent in analyses of childhood adversity, while smoking associations were more consistent using significance thresholds. SLCMA analytic versions similarly altered top associations, but time-dependent effects remained concordant. Alterations to data and analytic versions influenced the results of epigenome-wide analyses. Our findings highlight that magnitude and direction are better measures for replication and stability than p-value thresholds.

DNA methylation as a mediator of genetic and environmental influences on Parkinson's disease susceptibility: Impacts of alpha-Synuclein, physical activity, and pesticide exposure on the epigenome

Parkinson's disease (PD) is a neurodegenerative disorder with a complex etiology and increasing prevalence worldwide. As PD is influenced by a combination of genetic and environment/lifestyle factors in approximately 90% of cases, there is increasing interest in identification of the interindividual mechanisms underlying the development of PD as well as actionable lifestyle factors that can influence risk. This narrative review presents an outline of the genetic and environmental factors contributing to PD risk and explores the possible roles of cytosine methylation and hydroxymethylation in the etiology and/or as early-stage biomarkers of PD, with an emphasis on epigenome-wide association studies (EWAS) of PD conducted over the past decade. Specifically, we focused on variants in the SNCA gene, exposure to pesticides, and physical activity as key contributors to PD risk. Current research indicates that these factors individually impact the epigenome, particularly at the level of CpG methylation. There is also emerging evidence for interaction effects between genetic and environmental contributions to PD risk, possibly acting across multiple omics layers. We speculated that this may be one reason for the poor replicability of the results of EWAS for PD reported to date. Our goal is to provide direction for future epigenetics studies of PD to build upon existing foundations and leverage large datasets, new technologies, and relevant statistical approaches to further elucidate the etiology of this disease.

Evaluation of different computational methods for DNA methylation-based biological age

In recent years there has been a widespread interest in researching biomarkers of aging that could predict physiological vulnerability better than chronological age. Aging, in fact, is one of the most relevant risk factors for a wide range of maladies, and molecular surrogates of this phenotype could enable better patients stratification. Among the most promising of such biomarkers is DNA methylation-based biological age. Given the potential and variety of computational implementations (epigenetic clocks), we here present a systematic review of such clocks. Furthermore, we provide a large-scale performance comparison across different tissues and diseases in terms of age prediction accuracy and age acceleration, a measure of deviance from physiology. Our analysis offers both a state-of-the-art overview of the computational techniques developed so far and a heterogeneous picture of performances, which can be helpful in orienting future research.

Accelerated epigenetic aging in newborns with Down syndrome

Accelerated aging is a hallmark of Down syndrome (DS), with adults experiencing early-onset Alzheimer's disease and premature aging of the skin, hair, and immune and endocrine systems. Accelerated epigenetic aging has been found in the blood and brain tissue of adults with DS but when premature aging in DS begins remains unknown. We investigated whether accelerated aging in DS is already detectable in blood at birth. We assessed the association between age acceleration and DS using five epigenetic clocks in 346 newborns with DS and 567 newborns without DS using Illumina MethylationEPIC DNA methylation array data. We compared two epigenetic aging clocks (DNAmSkinBloodClock and pan-tissue DNAmAge) and three epigenetic gestational age clocks (Haftorn, Knight, and Bohlin) between DS and non-DS newborns using linear regression adjusting for observed age, sex, batch, deconvoluted blood cell proportions, and genetic ancestry. Targeted sequencing of GATA1 was performed in a subset of 184 newborns with DS to identify somatic mutations associated with transient abnormal myelopoiesis. DS was significantly associated with increased DNAmSkinBloodClock (effect estimate = 0.2442, p < 0.0001), with an epigenetic age acceleration of 244 days in newborns with DS after adjusting for potential confounding factors (95% confidence interval: 196-292 days). We also found evidence of epigenetic age acceleration associated with somatic GATA1 mutations among newborns with DS (p = 0.015). DS was not associated with epigenetic gestational age acceleration. We demonstrate that accelerated epigenetic aging in the blood of DS patients begins prenatally, with implications for the pathophysiology of immunosenescence and other aging-related traits in DS.

Effect of excess weight and insulin resistance on DNA methylation in prepubertal children

Epigenetic mechanisms, such as DNA methylation, regulate gene expression and play a role in the development of insulin resistance. This study evaluates how the BMI z-score (BMIz) and the homeostatic model assessment of insulin resistance (HOMA-IR), alone or in combination, relate to clinical outcomes and DNA methylation patterns in prepubertal children. DNA methylation in peripheral blood mononuclear cells (PBMCs) and clinical outcomes were measured in a cohort of 41 prepubertal children. Children with higher HOMA-IR had higher blood pressure and plasma lactate levels while children with higher BMIz had higher triglycerides levels. Moreover, the DNA methylation analysis demonstrated that a 1 unit increase in the BMIz was associated with a 0.41 (95% CI: 0.29, 0.53) increase in methylation of a CpG near the PPP6R2 gene. This gene is important in the regulation of NF-kB expression. However, there was no strong evidence that the BMIz and the HOMA-IR were synergistically related to any clinical or DNA methylation outcomes. In summary, the results suggest that obesity and insulin resistance may impact metabolic health both independently in prepubertal children. In addition, obesity also has an impact on the DNA methylation of the PPP6R2 gene. This may be a novel underlying starting point for the systemic inflammation associated with obesity and insulin resistance, in this population.

Batch-effect detection, correction and characterisation in Illumina HumanMethylation450 and MethylationEPIC BeadChip array data

Background

Genomic technologies can be subject to significant batch-effects which are known to reduce experimental power and to potentially create false positive results. The Illumina Infinium Methylation BeadChip is a popular technology choice for epigenome-wide association studies (EWAS), but presently, little is known about the nature of batch-effects on these designs. Given the subtlety of biological phenotypes in many EWAS, control for batch-effects should be a consideration.

Results

Using the batch-effect removal approaches in the ComBat and Harman software, we examined two in-house datasets and compared results with three large publicly available datasets, (1214 HumanMethylation450 and 1094 MethylationEPIC BeadChips in total), and find that despite various forms of preprocessing, some batch-effects persist. This residual batch-effect is associated with the day of processing, the individual glass slide and the position of the array on the slide. Consistently across all datasets, 4649 probes required high amounts of correction. To understand the impact of this set to EWAS studies, we explored the literature and found three instances where persistently batch-effect prone probes have been reported in abstracts as key sites of differential methylation. As well as batch-effect susceptible probes, we also discover a set of probes which are erroneously corrected. We provide batch-effect workflows for Infinium Methylation data and provide reference matrices of batch-effect prone and erroneously corrected features across the five datasets spanning regionally diverse populations and three commonly collected biosamples (blood, buccal and saliva).

Conclusions

Batch-effects are ever present, even in high-quality data, and a strategy to deal with them should be part of experimental design, particularly for EWAS. Batch-effect removal tools are useful to reduce technical variance in Infinium Methylation data, but they need to be applied with care and make use of post hoc diagnostic measures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01277-9.

Characterization of the Illumina EPIC Array for Optimal Applications in Epigenetic Research Targeting Diverse Human Populations

The Illumina EPIC array is widely used for high-throughput profiling of DNA cytosine modifications in human samples, covering more than 850,000 modification sites across various genomic features. The application of this platform is expected to provide novel insights into the epigenetic contribution to human complex traits and diseases. Considering the diverse inter-population genetic and epigenetic variation, it will benefit the research community with a comprehensive characterization of this platform for its applicability to major global populations. Specifically, we mapped 866,836 CpG probes from the EPIC array to the human genome reference. We detected 91,034 CpG probes that did not align reliably to the human genome reference. In addition, 21,256 CpG probes were found to ambiguously map to multiple loci in the human genome, and 448 probes showing inaccurate genomic information from the original Illumina annotations. We further characterized those uniquely mapped CpG probes in terms of whether they contained common genetic variants, i.e., single nucleotide polymorphisms (SNPs), in major global populations, by utilizing the 1000 Genomes Project data. A list of optimal CpG probes on the EPIC array was generated for major global populations, with the aim of providing a resource to facilitate future studies of diverse human populations. In conclusion, our analysis indicated that studies of diverse human populations using the EPIC array would be benefited by taking into account of the technical features of this platform.

DNA methylation signatures reveal that distinct combinations of transcription factors specify human immune cell epigenetic identity

Epigenetic reprogramming underlies specification of immune cell lineages, but patterns that uniquely define immune cell types and the mechanisms by which they are established remain unclear. Here, we identified lineage-specific DNA methylation signatures of six immune cell types from human peripheral blood and determined their relationship to other epigenetic and transcriptomic patterns. Sites of lineage-specific hypomethylation were associated with distinct combinations of transcription factors in each cell type. By contrast, sites of lineage-specific hypermethylation were restricted mostly to adaptive immune cells. PU.1 binding sites were associated with lineage-specific hypo- and hypermethylation in different cell types, suggesting that it regulates DNA methylation in a context-dependent manner. These observations indicate that innate and adaptive immune lineages are specified by distinct epigenetic mechanisms via combinatorial and context-dependent use of key transcription factors. The cell-specific epigenomics and transcriptional patterns identified serve as a foundation for future studies on immune dysregulation in diseases and aging.

NEOage clocks - epigenetic clocks to estimate post-menstrual and postnatal age in preterm infants

Epigenetic clocks based on DNA methylation (DNAm) can accurately predict chronological age and are thought to capture biological aging. A variety of epigenetic clocks have been developed for different tissue types and age ranges, but none have focused on postnatal age prediction for preterm infants. Epigenetic estimators of biological age might be especially informative in epidemiologic studies of neonates since DNAm is highly dynamic during the neonatal period and this is a key developmental window. Additionally, markers of biological aging could be particularly important for those born preterm since they are at heightened risk of developmental impairments. We aimed to fill this gap by developing epigenetic clocks for neonatal aging in preterm infants. As part of the Neonatal Neurobehavior and Outcomes in Very Preterm Infants (NOVI) study, buccal cells were collected at NICU discharge to profile DNAm levels in 542 very preterm infants. We applied elastic net regression to identify four epigenetic clocks (NEOage Clocks) predictive of post-menstrual and postnatal age, compatible with the Illumina EPIC and 450K arrays. We observed high correlations between predicted and reported ages (0.93 - 0.94) with root mean squared errors (1.28 - 1.63 weeks). Epigenetic estimators of neonatal aging in preterm infants can be useful tools to evaluate biological maturity and associations with neonatal and long-term morbidities.

Prenatal risk factors and neonatal DNA methylation in very preterm infants

Background

Prenatal risk factors are related to poor health and developmental outcomes for infants, potentially via epigenetic mechanisms. We tested associations between person-centered prenatal risk profiles, cumulative prenatal risk models, and epigenome-wide DNA methylation (DNAm) in very preterm neonates.

Methods

We studied 542 infants from a multi-center study of infants born < 30 weeks postmenstrual age. We assessed 24 prenatal risk factors via maternal report and medical record review. Latent class analysis was used to define prenatal risk profiles. DNAm was quantified from neonatal buccal cells using the Illumina MethylationEPIC Beadarray.

Results

We identified three latent profiles of women: a group with few risk factors (61%) and groups with elevated physical (26%) and psychological (13%) risk factors. Neonates born to women in higher risk subgroups had differential DNAm at 2 CpG sites. Higher cumulative prenatal risk was associated with methylation at 15 CpG sites, 12 of which were located in genes previously linked to physical and mental health and neurodevelopment.

Conclusion

We observed associations between prenatal risk factors and DNAm in very preterm infants using both person-centered and cumulative risk approaches. Epigenetics offers a potential biological indicator of prenatal risk exposure.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01164-9.

Redox Imbalance and Methylation Disturbances in Early Childhood Obesity

Obesity is increasing worldwide in prepubertal children, reducing the age of onset of associated comorbidities, including type 2 diabetes. Sulfur-containing amino acids, methionine, cysteine, and their derivatives play important roles in the transmethylation and transsulfuration pathways. Dysregulation of these pathways leads to alterations in the cellular methylation patterns and an imbalanced redox state. Therefore, we tested the hypothesis that one-carbon metabolism is already dysregulated in prepubertal children with obesity. Peripheral blood was collected from 64 children, and the plasma metabolites from transmethylation and transsulfuration pathways were quantified by HPLC. The cohort was stratified by BMI z-scores and HOMA-IR indices into healthy lean (HL), healthy obese (HO), and unhealthy obese (UHO). Fasting insulin levels were higher in the HO group compared to the HL, while the UHO had the highest. All groups presented normal fasting glycemia. Furthermore, high-density lipoprotein (HDL) was lower while triglycerides and lactate levels were higher in the UHO compared to HO subjects. S-adenosylhomocysteine (SAH) and total homocysteine levels were increased in the HO group compared to HL. Additionally, glutathione metabolism was also altered. Free cystine and oxidized glutathione (GSSG) were increased in the HO as compared to HL subjects. Importantly, the adipocyte secretory function was already compromised at this young age. Elevated circulating leptin and decreased adiponectin levels were observed in the UHO as compared to the HO subjects. Some of these alterations were concomitant with alterations in the DNA methylation patterns in the obese group, independent of the impaired insulin levels. In conclusion, our study informs on novel and important metabolic alterations in the transmethylation and the transsulfuration pathways in the early stages of obesity. Moreover, the altered secretory function of the adipocyte very early in life may be relevant in identifying early metabolic markers of disease that may inform on the increased risk for specific future comorbidities in this population.

Association of mammographic density with blood DNA methylation

Background: Altered DNA methylation may be an intermediate phenotype between breast cancer risk factors and disease. Mammographic density is a strong risk factor for breast cancer. However, no studies to date have identified an epigenetic signature of mammographic density. We performed an epigenome-wide association study of mammographic density.

Methods: White blood cell DNA methylation was measured for 385 postmenopausal women using the Illumina Infinium MethylationEPIC BeadChip array. Differential methylation was assessed using genome-wide, probe-level, and regional analyses. We implemented a resampling-based approach to improve the stability of our findings.

Results: On average, women with elevated mammographic density exhibited DNA hypermethylation within CpG islands and gene promoters compared to women with lower mammographic density. We identified 250 CpG sites for which DNA methylation was significantly associated with mammographic density. The top sites were located within genes associated with cancer, including HDLBP, TGFB2, CCT4, and PAX8, and were more likely to be located in regulatory regions of the genome. We also identified differential DNA methylation in 37 regions, including within the promoters of PAX8 and PF4, a gene involved in the regulation of angiogenesis. Overall, our results paint a picture of epigenetic dysregulation associated with mammographic density.

Conclusion: Mammographic density is associated with differential DNA methylation throughout the genome, including within genes associated with cancer. Our results suggest the potential involvement of several genes in the biological mechanisms behind differences in breast density between women. Further studies are warranted to explore these potential mechanisms and potential links to breast cancer risk.

Genetic impacts on DNA methylation: research findings and future perspectives

Multiple recent studies highlight that genetic variants can have strong impacts on a significant proportion of the human DNA methylome. Methylation quantitative trait loci, or meQTLs, allow for the exploration of biological mechanisms that underlie complex human phenotypes, with potential insights for human disease onset and progression. In this review, we summarize recent milestones in characterizing the human genetic basis of DNA methylation variation over the last decade, including heritability findings and genome-wide identification of meQTLs. We also discuss challenges in this field and future areas of research geared to generate insights into molecular processes underlying human complex traits.

Elucidation of the Genomic-Epigenomic Interaction Landscape of Aggressive Prostate Cancer

Background

Majority of prostate cancer (PCa) deaths are attributed to localized high-grade aggressive tumours which progress rapidly to metastatic disease. A critical unmet need in clinical management of PCa is discovery and characterization of the molecular drivers of aggressive tumours. The development and progression of aggressive PCa involve genetic and epigenetic alterations occurring in the germline, somatic (tumour), and epigenomes. To date, interactions between genes containing germline, somatic, and epigenetic mutations in aggressive PCa have not been characterized. The objective of this investigation was to elucidate the genomic-epigenomic interaction landscape in aggressive PCa to identify potential drivers aggressive PCa and the pathways they control. We hypothesized that aggressive PCa originates from a complex interplay between genomic (both germline and somatic mutations) and epigenomic alterations. We further hypothesized that these complex arrays of interacting genomic and epigenomic factors affect gene expression, molecular networks, and signaling pathways which in turn drive aggressive PCa.

Methods

We addressed these hypotheses by performing integrative data analysis combining information on germline mutations from genome-wide association studies with somatic and epigenetic mutations from The Cancer Genome Atlas using gene expression as the intermediate phenotype.

Results

The investigation revealed signatures of genes containing germline, somatic, and epigenetic mutations associated with aggressive PCa. Aberrant DNA methylation had effect on gene expression. In addition, the investigation revealed molecular networks and signalling pathways enriched for germline, somatic, and epigenetic mutations including the STAT3, PTEN, PCa, ATM, AR, and P53 signalling pathways implicated in aggressive PCa.

Conclusions

The study demonstrated that integrative analysis combining diverse omics data is a powerful approach for the discovery of potential clinically actionable biomarkers, therapeutic targets, and elucidation of oncogenic interactions between genomic and epigenomic alterations in aggressive PCa.

Bioinformatic Estimation of DNA Methylation and Hydroxymethylation Proportions

Simultaneous measurement of 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) at the single-nucleotide level can be obtained by combining data from DNA processing methods including traditional bisulfite (BS), oxidative bisulfite (oxBS), or Tet-assisted (TAB) bisulfite conversion. Array-based technologies have been widely used in this task, due to their time and cost efficiency. For methylation studies using BS data, many protocols and related packages have been suggested in the literature to deal with limitations and confounders that arise from array data. In this chapter, we illustrate how the reader can make small adjustments to these protocols to obtain estimates of methylation and hydroxymethylation proportions.

Before and After: Comparison of Legacy and Harmonized TCGA Genomic Data Commons' Data

We present a systematic analysis of the effects of synchronizing a large-scale, deeply characterized, multi-omic dataset to the current human reference genome, using updated software, pipelines, and annotations. For each of 5 molecular data platforms in The Cancer Genome Atlas (TCGA)-mRNA and miRNA expression, single nucleotide variants, DNA methylation and copy number alterations-comprehensive sample, gene, and probe-level studies were performed, towards quantifying the degree of similarity between the 'legacy' GRCh37 (hg19) TCGA data and its GRCh38 (hg38) version as 'harmonized' by the Genomic Data Commons. We offer gene lists to elucidate differences that remained after controlling for confounders, and strategies to mitigate their impact on biological interpretation. Our results demonstrate that the hg19 and hg38 TCGA datasets are very highly concordant, promote informed use of either legacy or harmonized omics data, and provide a rubric that encourages similar comparisons as new data emerge and reference data evolve.

Unraveling the Genomic-Epigenomic Interaction Landscape in Triple Negative and Non-Triple Negative Breast Cancer

Background: The recent surge of next generation sequencing of breast cancer genomes has enabled development of comprehensive catalogues of somatic mutations and expanded the molecular classification of subtypes of breast cancer. However, somatic mutations and gene expression data have not been leveraged and integrated with epigenomic data to unravel the genomic-epigenomic interaction landscape of triple negative breast cancer (TNBC) and non-triple negative breast cancer (non-TNBC). Methods: We performed integrative data analysis combining somatic mutation, epigenomic and gene expression data from The Cancer Genome Atlas (TCGA) to unravel the possible oncogenic interactions between genomic and epigenomic variation in TNBC and non-TNBC. We hypothesized that within breast cancers, there are differences in somatic mutation, DNA methylation and gene expression signatures between TNBC and non-TNBC. We further hypothesized that genomic and epigenomic alterations affect gene regulatory networks and signaling pathways driving the two types of breast cancer. Results: The investigation revealed somatic mutated, epigenomic and gene expression signatures unique to TNBC and non-TNBC and signatures distinguishing the two types of breast cancer. In addition, the investigation revealed molecular networks and signaling pathways enriched for somatic mutations and epigenomic changes unique to each type of breast cancer. The most significant pathways for TNBC were: retinal biosynthesis, BAG2, LXR/RXR, EIF2 and P2Y purigenic receptor signaling pathways. The most significant pathways for non-TNBC were: UVB-induced MAPK, PCP, Apelin endothelial, Endoplasmatic reticulum stress and mechanisms of viral exit from host signaling Pathways. Conclusion: The investigation revealed integrated genomic, epigenomic and gene expression signatures and signing pathways unique to TNBC and non-TNBC, and a gene signature distinguishing the two types of breast cancer. The study demonstrates that integrative analysis of multi-omics data is a powerful approach for unravelling the genomic-epigenomic interaction landscape in TNBC and non-TNBC.

Evaluation of pre-processing on the meta-analysis of DNA methylation data from the Illumina HumanMethylation450 BeadChip platform

Introduction

Meta-analysis is a powerful means for leveraging the hundreds of experiments being run worldwide into more statistically powerful analyses. This is also true for the analysis of omic data, including genome-wide DNA methylation. In particular, thousands of DNA methylation profiles generated using the Illumina 450k are stored in the publicly accessible Gene Expression Omnibus (GEO) repository. Often, however, the intensity values produced by the BeadChip (raw data) are not deposited, therefore only pre-processed values -obtained after computational manipulation- are available. Pre-processing is possibly different among studies and may then affect meta-analysis by introducing non-biological sources of variability.

Material and methods

To systematically investigate the effect of pre-processing on meta-analysis, we analysed four different collections of DNA methylation samples (datasets), each composed of two subsets, for which raw data from controls (i.e. healthy subjects) and cases (i.e. patients) are available. We pre-processed the data from each dataset with nine among the most common pipelines found in literature. Moreover, we evaluated the performance of regRCPqn, a modification of the RCP algorithm that aims to improve data consistency. For each combination of pre-processing (9 × 9), we first evaluated the between-sample variability among control subjects and, then, we identified genomic positions that are differentially methylated between cases and controls (differential analysis).

Results and conclusion

The pre-processing of DNA methylation data affects both the between-sample variability and the loci identified as differentially methylated, and the effects of pre-processing are strongly dataset-dependent. By contrast, application of our renormalization algorithm regRCPqn: (i) reduces variability and (ii) increases agreement between meta-analysed datasets, both critical components of data harmonization.

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.

Longitudinal DNA methylation changes at MET may alter HGF/c-MET signalling in adolescents at risk for depression

Unrecognized depression during adolescence can result in adult suicidal behaviour. The aim of this study was to identify, replicate and characterize DNA methylation (DNAm) shifts in depression aetiology, using a longitudinal, multi-tissue (blood and brain) and multi-layered (genetics, epigenetics, transcriptomics) approach. We measured genome-wide blood DNAm data at baseline and one-year follow-up, and imputed genetic variants, in 59 healthy adolescents comprising the discovery cohort. Depression and suicidal symptoms were determined using the Development and Well-Being Assessment (DAWBA) depression band, Montgomery-Åsberg Depression Rating Scale-Self (MADRS-S) and SUicide Assessment Scale (SUAS). DNAm levels at follow-up were regressed against depression scores, adjusting for sex, age and the DNAm residuals at baseline. Higher methylation levels of 5% and 13% at cg24627299 within the MET gene were associated with higher depression scores (p_raw<1e-4) and susceptibility for suicidal symptoms (p_adj.<0.005). The nearby rs39748 was discovered to be a methylation and expression quantitative trait locus in blood cells. mRNA levels of hepatocyte growth factor (HGF) expression, known to strongly interact with MET, were inversely associated with methylation levels at cg24627299, in an independent cohort of 1180 CD14+ samples. In an open-access dataset of brain tissue, lower methylation at cg24627299 was found in 45 adults diagnosed with major depressive disorder compared with matched controls (p_adj.<0.05). Furthermore, lower MET expression was identified in the hippocampus of depressed individuals compared with controls in a fourth, independent cohort. Our findings reveal methylation changes at MET in the pathology of depression, possibly involved in downregulation of HGF/c-MET signalling the hippocampal region.

Deconvolution of the Genomic and Epigenomic Interaction Landscape of Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is the most aggressive form of breast cancer. Emerging evidenced suggests that both genetics and epigenetic factors play a role in the pathogenesis of TNBC. However, oncogenic interactions and cooperation between genomic and epigenomic variation have not been characterized. The objective of this study was to deconvolute the genomic and epigenomic interaction landscape in TNBC using an integrative genomics approach, which integrates information on germline, somatic, epigenomic and gene expression variation. We hypothesized that TNBC originates from a complex interplay between genomic (both germline and somatic variation) and epigenomic variation. We further hypothesized that these complex arrays of interacting genomic and epigenomic factors affect entire molecular networks and signaling pathways which, in turn, drive TNBC. We addressed these hypotheses using germline variation from genome-wide association studies and somatic, epigenomic and gene expression variation from The Cancer Genome Atlas (TCGA). The investigation revealed signatures of functionally related genes containing germline, somatic and epigenetic variations. DNA methylation had an effect on gene expression. Network and pathway analysis revealed molecule networks and signaling pathways enriched for germline, somatic and epigenomic variation, among them: Role of BRCA1 in DNA Damage Response, Hereditary Breast Cancer Signaling, Molecular Mechanisms of Cancer, Estrogen-Dependent Breast Cancer, p53, MYC Mediated Apoptosis, and PTEN Signaling pathways. The investigation revealed that integrative genomics is a powerful approach for deconvoluting the genomic-epigenomic interaction landscape in TNBC. Further studies are needed to understand the biological mechanisms underlying oncogenic interactions between genomic and epigenomic factors in TNBC.

Methylomic profiles reveal sex-specific differences in leukocyte composition associated with post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) is a debilitating mental disorder precipitated by trauma exposure. However, only some persons exposed to trauma develop PTSD. There are sex differences in risk; twice as many women as men develop a lifetime diagnosis of PTSD. Methylomic profiles derived from peripheral blood are well-suited for investigating PTSD because DNA methylation (DNAm) encodes individual response to trauma and may play a key role in the immune dysregulation characteristic of PTSD pathophysiology. In the current study, we leveraged recent methodological advances to investigate sex-specific differences in DNAm-based leukocyte composition that are associated with lifetime PTSD. We estimated leukocyte composition on a combined methylation array dataset (483 participants, ∼450 k CpG sites) consisting of two civilian cohorts, the Detroit Neighborhood Health Study and Grady Trauma Project. Sex-stratified Mann-Whitney U test and two-way ANCOVA revealed that lifetime PTSD was associated with significantly higher monocyte proportions in males, but not in females (Holm-adjusted p-val < 0.05). No difference in monocyte proportions was observed between current and remitted PTSD cases in males, suggesting that this sex-specific difference may reflect a long-standing trait of lifetime history of PTSD, rather than current state of PTSD. Associations with lifetime PTSD or PTSD status were not observed in any other leukocyte subtype and our finding in monocytes was confirmed using cell estimates based on a different deconvolution algorithm, suggesting that our sex-specific findings are robust across cell estimation approaches. Overall, our main finding of elevated monocyte proportions in males, but not in females with lifetime history of PTSD provides evidence for a sex-specific difference in peripheral blood leukocyte composition that is detectable in methylomic profiles and that may reflect long-standing changes associated with PTSD diagnosis.

SeSAMe: reducing artifactual detection of DNA methylation by Infinium BeadChips in genomic deletions

We report a new class of artifacts in DNA methylation measurements from Illumina HumanMethylation450 and MethylationEPIC arrays. These artifacts reflect failed hybridization to target DNA, often due to germline or somatic deletions and manifest as incorrectly reported intermediate methylation. The artifacts often survive existing preprocessing pipelines, masquerade as epigenetic alterations and can confound discoveries in epigenome-wide association studies and studies of methylation-quantitative trait loci. We implement a solution, P-value with out-of-band (OOB) array hybridization (pOOBAH), in the R package SeSAMe. Our method effectively masks deleted and hyperpolymorphic regions, reducing or eliminating spurious reports of epigenetic silencing at oft-deleted tumor suppressor genes such as CDKN2A and RB1 in cases with somatic deletions. Furthermore, our method substantially decreases technical variation whilst retaining biological variation, both within and across HM450 and EPIC platform measurements. SeSAMe provides a light-weight, modular DNA methylation data analysis suite, with a performant implementation suitable for efficient analysis of thousands of samples.

Different Methylation of CpG-SNPs in Behcet's Disease

Purpose

We recently performed an Epigenome-Wide Association Studies (EWAS) study in Behcet's disease (BD) and identified various cytosine–phosphate–guanine (CpG) loci that were aberrantly methylated. In the current study, we wanted to investigate whether these sites contained genetic polymorphisms and whether the frequency of these polymorphisms was altered in BD.

Methods

A two-stage study was performed. The first stage involved 358 BD patients and 704 healthy controls to investigate genetic variants of 10 CpG-SNPs (rs10454134, rs176249, rs3808620, rs10176517, rs11247118, rs78016579, rs9461624, rs10492166, rs34929465, and rs6507921) using an iPLEX Gold genotyping assay and a Sequenom MassARRAY. In the second stage, an additional 172 independent BD patients and 330 healthy individuals are to confirm trends found in the first stage.

Results

A higher frequency of both the rs10454134 AG genotypes (p = 0.008, OR = 1.413, 95% CI = 1.094-1.826) and a lower GG genotype frequency (p = 0.003, OR = 0.630, 95% CI = 0.465-0.854) were found in BD patients compared to the controls in the first stage. However, after correcting for multiple comparisons, all associations identified in the first stage lost statistical significance. The frequencies of the other CpG-SNPs investigated were not different between BD patients and controls. The second stage was designed using an additional cohort to confirm the association with CpG-SNP, rs10454134. The data failed to confirm the association between this CpG-SNP and BD.

Conclusions

This study did not show an association between BD and CpG-SNPs in gene sites that were earlier shown to be aberrantly methylated.

Genome-wide methylation profiling identified novel differentially hypermethylated biomarker MPPED2 in colorectal cancer

BACKGROUND

Epigenetic alternation is a common contributing factor to neoplastic transformation. Although previous studies have reported a cluster of aberrant promoter methylation changes associated with silencing of tumor suppressor genes, little is known concerning their sequential DNA methylation changes during the carcinogenetic process. The aim of the present study was to address a genome-wide search for identifying potentially important methylated changes and investigate the onset and pattern of methylation changes during the progression of colorectal neoplasia.

METHODS

A three-phase design was employed in this study. In the screening phase, DNA methylation profile of 12 pairs of colorectal cancer (CRC) and adjacent normal tissues was analyzed by using the Illumina MethylationEPIC BeadChip. Significant CpG sites were selected based on a cross-validation analysis from The Cancer Genome Atlas (TCGA) database. Methylation levels of candidate CpGs were assessed using pyrosequencing in the training dataset (tumor lesions and adjacent normal tissues from 46 CRCs) and the validation dataset (tumor lesions and paired normal tissues from 13 hyperplastic polyps, 129 adenomas, and 256 CRCs). A linear mixed-effects model was used to examine the incremental changes of DNA methylation during the progression of colorectal neoplasia.

RESULTS

The comparisons between normal and tumor samples in the screening phase revealed an extensive CRC-specific methylomic pattern with 174,006 (21%) methylated CpG sites, of which 22,232 (13%) were hyermethylated and 151,774 (87%) were hypomethylated. Hypermethylation mostly occurred in CpG islands with an overlap of gene promoters, while hypomethylation tended to be mapped far away from functional regions. Further cross validation analysis from TCGA dataset confirmed 265 hypermethylated promoters coupling with downregulated gene expression. Among which, hypermethylated changes in MEEPD2 promoter was successfully replicated in both training and validation phase. Significant hypermethylation appeared since precursor lesions with an extensive modification in CRCs. The linear mixed-effects modeling analysis found that a cumulative pattern of MPPED2 methylation changes from normal mucosa to hyperplastic polyp to adenoma, and to carcinoma (P < 0.001).

CONCLUSIONS

Our findings indicate that epigenetic alterations of MPPED2 promoter region appear sequentially during the colorectal neoplastic progression. It might be able to serve as a promising biomarker for early diagnosis and stage surveillance of colorectal tumorigenesis.

Meffil: efficient normalization and analysis of very large DNA methylation datasets

Motivation

DNA methylation datasets are growing ever larger both in sample size and genome coverage. Novel computational solutions are required to efficiently handle these data.

Results

We have developed meffil, an R package designed for efficient quality control, normalization and epigenome-wide association studies of large samples of Illumina Methylation BeadChip microarrays. A complete re-implementation of functional normalization minimizes computational memory without increasing running time. Incorporating fixed and random effects within functional normalization, and automated estimation of functional normalization parameters reduces technical variation in DNA methylation levels, thus reducing false positive rates and improving power. Support for normalization of datasets distributed across physically different locations without needing to share biologically-based individual-level data means that meffil can be used to reduce heterogeneity in meta-analyses of epigenome-wide association studies.

Availability and implementation

https://github.com/perishky/meffil/.

Supplementary information

Supplementary data are available at Bioinformatics online.

LogLoss-BERAF: An ensemble-based machine learning model for constructing highly accurate diagnostic sets of methylation sites accounting for heterogeneity in prostate cancer

Although modern methods of whole genome DNA methylation analysis have a wide range of applications, they are not suitable for clinical diagnostics due to their high cost and complexity and due to the large amount of sample DNA required for the analysis. Therefore, it is crucial to be able to identify a relatively small number of methylation sites that provide high precision and sensitivity for the diagnosis of pathological states. We propose an algorithm for constructing limited subsamples from high-dimensional data to form diagnostic panels. We have developed a tool that utilizes different methods of selection to find an optimal, minimum necessary combination of factors using cross-entropy loss metrics (LogLoss) to identify a subset of methylation sites. We show that the algorithm can work effectively with different genome methylation patterns using ensemble-based machine learning methods. Algorithm efficiency, precision and robustness were evaluated using five genome-wide DNA methylation datasets (totaling 626 samples), and each dataset was classified into tumor and non-tumor samples. The algorithm produced an AUC of 0.97 (95% CI: 0.94-0.99, 9 sites) for prostate adenocarcinoma and an AUC of 1.0 (from 2 to 6 sites) for urothelial bladder carcinoma, two types of kidney carcinoma and colorectal carcinoma. For prostate adenocarcinoma we showed that identified differential variability methylation patterns distinguish cluster of samples with higher recurrence rate (hazard ratio for recurrence = 0.48, 95% CI: 0.05-0.92; log-rank test, p-value < 0.03). We also identified several clusters of correlated interchangeable methylation sites that can be used for the elaboration of biological interpretation of the resulting models and for further selection of the sites most suitable for designing diagnostic panels. LogLoss-BERAF is implemented as a standalone python code and open-source code is freely available from https://github.com/bioinformatics-IBCH/logloss-beraf along with the models described in this article.

Methylation-derived Neutrophil-to-Lymphocyte Ratio and Lung Cancer Risk in Heavy Smokers

The neutrophil-to-lymphocyte ratio (NLR) is a biomarker that indicates systemic inflammation and can be estimated using array-based DNA methylation data as methylation-derived NLR (mdNLR). We assessed the relationship between prediagnosis mdNLR and lung cancer risk in a nested case-control study in the β-Carotene and Retinol Efficacy Trial (CARET) of individuals at high risk for lung cancer due to heavy smoking or substantial occupational asbestos exposure. We matched 319 incident lung cancer cases to controls based on age at blood draw, smoking, sex, race, asbestos, enrollment year, and time at risk. We computed mdNLR using the ratio of predicted granulocyte and lymphocyte proportions derived from DNA methylation signatures in whole blood collected prior to diagnosis (median 4.4 years in cases). Mean mdNLR was higher in cases than controls (2.06 vs. 1.86, P = 0.03). Conditional logistic regression models adjusted for potential confounders revealed a 21% increased risk of lung cancer per unit increase in mdNLR [OR 1.21; 95% confidence interval (CI) 1.01-1.45]. A 30% increased risk of non-small cell lung cancer (NSCLC) was observed for each unit increase in mdNLR (n = 240 pairs; OR 1.30, 95% CI, 1.03-1.63), and there was no statistically significant association between mdNLR and small-cell lung cancer risk. The mdNLR-NSCLC association was most pronounced in those with asbestos exposure (n = 42 male pairs; OR 3.39; 95% CI, 1.32-8.67). A better understanding of the role of mdNLR in lung cancer etiology may improve prevention and detection of lung cancer. Cancer Prev Res; 11(11); 727-34. ©2018 AACR.

Comparison of pre-processing methods for Infinium HumanMethylation450 BeadChip array

Motivation

Microarrays are widely used to quantify DNA methylation because they are economical, require only small quantities of input DNA and focus on well-characterized regions of the genome. However, pre-processing of methylation microarray data is challenging because of confounding factors that include background fluorescence, dye bias and the impact of germline polymorphisms. Therefore, we present valuable insights and a framework for those seeking the most optimal pre-processing method through a data-driven approach.

Results

Here, we show that Dasen is the optimal pre-processing methodology for the Infinium HumanMethylation450 BeadChip array in prostate cancer, a frequently employed platform for tumour methylome profiling in both the TCGA and ICGC consortia. We evaluated the impact of 11 pre-processing methods on batch effects, replicate variabilities, sensitivities and sample-to-sample correlations across 809 independent prostate cancer samples, including 150 reported for the first time in this study. Overall, Dasen is the most effective for removing artefacts and detecting biological differences associated with tumour aggressivity. Relative to the raw dataset, it shows a reduction in replicate variances of 67% and 76% for β- and M-values, respectively. Our study provides a unique pre-processing benchmark for the community with an emphasis on biological implications.

Availability and implementation

All software used in this study are publicly available as detailed in the article.

Contact

paul.boutros@oicr.on.ca.

Supplementary information

Supplementary data are available at Bioinformatics online.

Investigating the Epigenetic Discrimination of Identical Twins Using Buccal Swabs, Saliva, and Cigarette Butts in the Forensic Setting

Monozygotic (MZ) twins are typically indistinguishable via forensic DNA profiling. Recently, we demonstrated that epigenetic differentiation of MZ twins is feasible; however, proportions of twin differentially methylated CpG sites (tDMSs) identified in reference-type blood DNA were not replicated in trace-type blood DNA. Here we investigated buccal swabs as typical forensic reference material, and saliva and cigarette butts as commonly encountered forensic trace materials. As an analog to a forensic case, we analyzed one MZ twin pair. Epigenome-wide microarray analysis in reference-type buccal DNA revealed 25 candidate tDMSs with >0.5 twin-to-twin differences. MethyLight quantitative PCR (qPCR) of 22 selected tDMSs in trace-type DNA revealed in saliva DNA that six tDMSs (27.3%) had >0.1 twin-to-twin differences, seven (31.8%) had smaller (<0.1) but robustly detected differences, whereas for nine (40.9%) the differences were in the opposite direction relative to the microarray data; for cigarette butt DNA, results were 50%, 22.7%, and 27.3%, respectively. The discrepancies between reference-type and trace-type DNA outcomes can be explained by cell composition differences, method-to-method variation, and other technical reasons including bisulfite conversion inefficiency. Our study highlights the importance of the DNA source and that careful characterization of biological and technical effects is needed before epigenetic MZ twin differentiation is applicable in forensic casework.

Association between fetal exposure to phthalate endocrine disruptor and genome-wide DNA methylation at birth

BACKGROUND

Phthalic acid esters are ubiquitous and antiandrogenic, and may cause systemic effects in humans, particularly with in utero exposure. Epigenetic modification, such as DNA methylation, has been hypothesized to be an important mechanism that mediates certain biological processes and pathogenic effects of in utero phthalate exposure.

OBJECTIVE

The aim of this study was to examine the association between genome-wide DNA methylation at birth and prenatal exposure to phthalate.

METHODS

We studied 64 infant-mother pairs included in TMICS (Taiwan Maternal and Infant Cohort Study), a long-term follow-up birth cohort from the general population. DNA methylation levels at more than 450,000 CpG sites were measured in cord blood samples using Illumina Infinium HumanMethylation450 BeadChips. The concentrations of three metabolites of di-(2-ethylhexyl) phthalate (DEHP) were measured using liquid chromatography tandem-mass spectrometry (LC-MS/MS) in urine samples collected from the pregnant women during 28-36 weeks gestation.

RESULTS

We identified 25 CpG sites whose methylation levels in cord blood were significantly correlated with prenatal DEHP exposure using a false discovery rate (FDR) of 5% (q-value < 0.05). Via gene-set enrichment analysis (GSEA), we also found that there was significant enrichment of genes involved in the androgen response, estrogen response, and spermatogenesis within those genes showing DNA methylation changes in response to exposure. Specifically, PA2G4, HMGCR, and XRCC6 genes were involved in genes in response to androgen.

CONCLUSIONS

Phthalate exposure in utero may cause significant alterations in the DNA methylation in cord blood. These changes in DNA methylation might serve as biomarkers of maternal exposure to phthalate in infancy and potential candidates for studying mechanisms via which phthalate may impact on health in later life. Future investigations are warranted.

Agreement in DNA methylation levels from the Illumina 450K array across batches, tissues, and time

Epigenome-wide association studies (EWAS) have focused primarily on DNA methylation as a chemically stable and functional epigenetic modification. However, the stability and accuracy of the measurement of methylation in different tissues and extraction types is still being actively studied, and the longitudinal stability of DNA methylation in commonly studied peripheral tissues is of great interest. Here, we used data from two studies, three tissue types, and multiple time points to assess the stability of DNA methylation measured with the Illumina Infinium HumanMethylation450 BeadChip array. Redundancy analysis enabled visual assessment of agreement of replicate samples overall and showed good agreement after removing effects of tissue type, age, and sex. At the probe level, analysis of variance contrasts separating technical and biological replicates clearly showed better agreement between technical replicates versus longitudinal samples, and suggested increased stability for buccal cells versus blood or blood spots. Intraclass correlations (ICCs) demonstrated that inter-individual variability is of similar magnitude to within-sample variability at many probes; however, as inter-individual variability increased, so did ICC. Furthermore, we were able to demonstrate decreasing agreement in methylation levels with time, despite a maximal sampling interval of only 576 days. Finally, at 6 popular candidate genes, there was a large range of stability across probes. Our findings highlight important sources of technical and biological variation in DNA methylation across different tissues over time. These data will help to inform longitudinal sampling strategies of future EWAS.

Epigenetic research in multiple sclerosis: progress, challenges, and opportunities

Multiple sclerosis (MS) is a chronic inflammatory and demyelinating disease of the central nervous system. MS likely results from a complex interplay between predisposing causal gene variants (the strongest influence coming from HLA class II locus) and environmental risk factors such as smoking, infectious mononucleosis, and lack of sun exposure/vitamin D. However, little is known about the mechanisms underlying MS development and progression. Moreover, the clinical heterogeneity and variable response to treatment represent additional challenges to a comprehensive understanding and efficient treatment of disease. Epigenetic processes, such as DNA methylation and histone posttranslational modifications, integrate influences from the genes and the environment to regulate gene expression accordingly. Studying epigenetic modifications, which are stable and reversible, may provide an alternative approach to better understand and manage disease. We here aim to review findings from epigenetic studies in MS and further discuss the challenges and clinical opportunities arising from epigenetic research, many of which apply to other diseases with similar complex etiology. A growing body of evidence supports a role of epigenetic processes in the mechanisms underlying immune pathogenesis and nervous system dysfunction in MS. However, disparities between studies shed light on the need to consider possible confounders and methodological limitations for a better interpretation of the data. Nevertheless, translational use of epigenetics might offer new opportunities in epigenetic-based diagnostics and therapeutic tools for a personalized care of MS patients.