Comprehensive Evaluation of The Infinium Human MethylationEPIC v2 BeadChip

UHRF1 ubiquitin ligase activity supports the maintenance of low-density CpG methylation

The RING E3 ubiquitin ligase UHRF1 is an established cofactor for DNA methylation inheritance. The model posits that nucleosomal engagement through histone and DNA interactions directs UHRF1 ubiquitin ligase activity toward lysines on histone H3 tails, creating binding sites for DNMT1 through ubiquitin interacting motifs (UIM1 and UIM2). However, the extent to which DNMT1 relies on ubiquitin signaling through UHRF1 in support of DNA methylation maintenance remains unclear. Here, with integrative epigenomic and biochemical analyses, we reveal that DNA methylation maintenance at low-density cytosine-guanine dinucleotides (CpGs) is particularly vulnerable to disruption of UHRF1 ubiquitin ligase activity and DNMT1 ubiquitin reading activity through UIM1. Hypomethylation of low-density CpGs in this manner induces formation of partially methylated domains (PMDs), a methylation signature observed across human cancers. In contrast, UIM2 disruption completely abolishes the DNA methylation maintenance function of DNMT1 in a CpG density-independent manner. In the context of DNA methylation recovery following acute DNMT1 depletion, we further reveal a 'bookmarking' function for UHRF1 ubiquitin ligase activity in support of DNA re-methylation. Collectively, these studies show that DNMT1-dependent DNA methylation inheritance is a ubiquitin-regulated process that is partially reliant on UHRF1 and suggest a disrupted UHRF1-DNMT1 ubiquitin signaling axis contributes to PMD formation in cancers.

Role of epigenetic mechanisms in the pathogenesis of chronic respiratory diseases and response to inhaled exposures: From basic concepts to clinical applications

Epigenetic modifications are chemical groups in our DNA (and chromatin) that determine which genes are active and which are shut off. Importantly, they integrate environmental signals to direct cellular function. Upon chronic environmental exposures, the epigenetic signature of lung cells gets altered, triggering aberrant gene expression programs that can lead to the development of chronic lung diseases. In addition to driving disease, epigenetic marks can serve as attractive lung disease biomarkers, due to early onset, disease specificity, and stability, warranting the need for more epigenetic research in the lung field. Despite substantial progress in mapping epigenetic alterations (mostly DNA methylation) in chronic lung diseases, the molecular mechanisms leading to their establishment are largely unknown. This review is meant as a guide for clinicians and lung researchers interested in epigenetic regulation with a focus on DNA methylation. It provides a short introduction to the main epigenetic mechanisms (DNA methylation, histone modifications and non-coding RNA) and the machinery responsible for their establishment and removal. It presents examples of epigenetic dysregulation across a spectrum of chronic lung diseases and discusses the current state of epigenetic therapies. Finally, it introduces the concept of epigenetic editing, an exciting novel approach to dissecting the functional role of epigenetic modifications. The promise of this emerging technology for the functional study of epigenetic mechanisms in cells and its potential future use in the clinic is further discussed.

DNA hypomethylation promotes UHRF1-and SUV39H1/H2-dependent crosstalk between H3K18ub and H3K9me3 to reinforce heterochromatin states

Mono-ubiquitination of lysine 18 on histone H3 (H3K18ub), catalyzed by UHRF1, is a DNMT1 docking site that facilitates replication-coupled DNA methylation maintenance. Its functions beyond this are unknown. Here, we genomically map simultaneous increases in UHRF1-dependent H3K18ub and SUV39H1/H2-dependent H3K9me3 following DNMT1 inhibition. Mechanistically, transient accumulation of hemi-methylated DNA at CpG islands facilitates UHRF1 recruitment and E3 ligase activity toward H3K18. Notably, H3K18ub enhances SUV39H1/H2 methyltransferase activity and, in colon cancer cells, nucleates new H3K9me3 domains at CpG island promoters of DNA methylation-silenced tumor suppressor genes (TSGs). Disrupting UHRF1 enzyme activity prevents H3K9me3 accumulation while promoting PRC2-dependent H3K27me3 as a tertiary layer of gene repression in these regions. By contrast, disrupting H3K18ub-dependent SUV39H1/H2 activity enhances the transcriptional activating and antiproliferative effects of DNMT1 inhibition. Collectively, these findings reveal roles for UHRF1 and H3K18ub in regulating a hierarchy of repressive histone methylation signaling and rationalize a combination strategy for epigenetic cancer therapy.

Technical variability across the 450K, EPICv1, and EPICv2 DNA methylation arrays: lessons learned for clinical and longitudinal studies

DNA methylation (DNAm) is the most commonly measured epigenetic mechanism in human populations, with most studies using Illumina arrays to assess DNAm levels. In 2023, Illumina updated their DNAm arrays to the EPIC version 2 (EPICv2), building on prior iterations, namely the EPIC version 1 (EPICv1) and 450K arrays. Whether DNAm measurements are stable across these three generations of arrays has yet not been investigated, limiting the ability of researchers-especially those with longitudinal data-to compare and replicate results across arrays. Here, we present results from a study of 30 child participants (15 male; 15 female) from the Drakenstein Child Health Study, who had DNAm measured on all three of the latest arrays: 450K, EPICv1, and EPICv2. Using these data, we created an annotation of probe quality across arrays, which includes the intraclass correlations, interquartile ranges, correlations, and array bias (i.e., the extent to which DNAm levels were explained by array type) of all CpGs. We also present results from an analysis of sex differences, where we found that CpGs with lower replicability across arrays had higher array-based variance, suggesting this variance metric help guide replication efforts. We also showed that epigenetic age estimates across arrays were more stable when using the principal component versions of epigenetic clocks. Ultimately, this collection of results provides a framework for investigating the replicability and longitudinal stability of epigenetic changes across multiple versions of Illumina DNAm arrays.

MethylCallR : a comprehensive analysis framework for Illumina Methylation Beadchip

DNA methylation is a molecular process that mediates gene-environment interactions. Epigenome-wide association studies (EWAS) using the Illumina Human Methylation BeadChip are powerful tools for quantifying the relationship between DNA methylation and phenotypes. Recently, the Illumina Methylation EPICv2 BeadChip (EPICv2) was released, which includes new features, such as duplicated probes and changed probe names. Several published algorithms have been updated to address these features in EPICv2. However, appropriate EPICv2 preprocessing and integration with previous microarray versions remain complex. Therefore, MethylCallR, an open-source R package designed to provide standard procedures for performing EWAS using Illumina methylation microarrays including EPICv2, was developed. MethylCallR can be used to control duplicated probes in EPICv2, by using pre-set data implemented in MethylCallR or new customized data. MethylCallR includes a straightforward conversion function between different types of Illumina Human Methylation BeadChips. Using MethylCallR, potential outlier sample detection and statistical power estimation were conducted and used to select meaningful probes. Publicly available data was analyzed using MethylCallR and the findings were compared to that of a previous study.

Ultrasensitive Amplification-Free Quantification of a Methyl CpG-Rich Cancer Biomarker by Single-Molecule Kinetic Fingerprinting

The most well-studied epigenetic marker in humans is the 5-methyl modification of cytosine in DNA, which has great potential as a disease biomarker. Currently, quantification of DNA methylation relies heavily on bisulfite conversion followed by PCR amplification and NGS or microarray analysis. PCR is subject to potential bias in differential amplification of bisulfite-converted methylated versus unmethylated sequences. Here, we combine bisulfite conversion with single-molecule kinetic fingerprinting to develop an amplification-free assay for DNA methylation at the branched-chain amino acid transaminase 1 (BCAT1) promoter. Our assay selectively responds to methylated sequences with a limit of detection below 1 fM and a specificity of 99.9999%. Evaluating complex genomic DNA matrices, we reliably distinguish <5% DNA methylation at the BCAT1 promoter in whole blood DNA from completely unmethylated whole-genome amplified DNA. Taken together, these results demonstrate the feasibility and sensitivity of our amplification-free, single-molecule quantification approach to improve the early detection of methylated cancer DNA biomarkers.

Robust estimation of cancer and immune cell-type proportions from bulk tumor ATAC-Seq data

Assay for Transposase-Accessible Chromatin sequencing (ATAC-Seq) is a widely used technique to explore gene regulatory mechanisms. For most ATAC-Seq data from healthy and diseased tissues such as tumors, chromatin accessibility measurement represents a mixed signal from multiple cell types. In this work, we derive reliable chromatin accessibility marker peaks and reference profiles for most non-malignant cell types frequently observed in the microenvironment of human tumors. We then integrate these data into the EPIC deconvolution framework (Racle et al., 2017) to quantify cell-type heterogeneity in bulk ATAC-Seq data. Our EPIC-ATAC tool accurately predicts non-malignant and malignant cell fractions in tumor samples. When applied to a human breast cancer cohort, EPIC-ATAC accurately infers the immune contexture of the main breast cancer subtypes.

Applicability of epigenetic age models to next-generation methylation arrays

BACKGROUND

Epigenetic clocks are mathematical models used to estimate epigenetic age based on DNA methylation at specific CpG sites. As new methylation microarrays are developed and older models discontinued, existing epigenetic clocks might become obsolete. Here, we explored the effects of the changes introduced in the new EPICv2 DNA methylation array on existing epigenetic clocks.

METHODS

We tested the performance of four epigenetic clocks on the probeset of the EPICv2 array using a dataset of 10,835 samples. We developed a new epigenetic age prediction model compatible across the 450 k, EPICv1, and EPICv2 microarrays and validated it on 2095 samples. We estimated technical noise and intra-subject variation using two datasets with repeated sampling. We used data from (i) cancer survivors who had undergone different therapies, (ii) breast cancer patients and controls, and (iii) an exercise-based interventional study, to test the ability of our model to detect alterations in epigenetic age acceleration in response to theoretically antiaging interventions.

RESULTS

The results of the four epiclocks tested are significantly distorted by the EPICv2 probeset, causing an average difference of up to 25 years. Our new model produced highly accurate chronological age predictions, comparable to a state-of-the-art epiclock. The model reported the lowest epigenetic age acceleration in normal populations, as well as the lowest variation across technical replicates and repeated samples from the same subjects. Finally, our model reproduced previous results of increased epigenetic age acceleration in cancer patients and in survivors treated with radiation therapy, and no changes from exercise-based interventions.

CONCLUSION

Existing epigenetic clocks require updates for full EPICv2 compatibility. Our new model translates the capabilities of state-of-the-art epigenetic clocks to the EPICv2 platform and is cross-compatible with older microarrays. The characterization of epigenetic age prediction variation provides useful metrics to contextualize the relevance of epigenetic age alterations. The analysis of data from subjects influenced by radiation, cancer, and exercise-based interventions shows that despite being good predictors of chronological age, neither a pathological state like breast cancer, a hazardous environmental factor (radiation), nor exercise (a beneficial intervention) caused significant changes in the values of the "epigenetic age" determined by these first-generation models.

Retro-age: A unique epigenetic biomarker of aging captured by DNA methylation states of retroelements

Reactivation of retroelements in the human genome has been linked to aging. However, whether the epigenetic state of specific retroelements can predict chronological age remains unknown. We provide evidence that locus-specific retroelement DNA methylation can be used to create retroelement-based epigenetic clocks that accurately measure chronological age in the immune system, across human tissues, and pan-mammalian species. We also developed a highly accurate retroelement epigenetic clock compatible with EPICv.2.0 data that was constructed from CpGs that did not overlap with existing first- and second-generation epigenetic clocks, suggesting a unique signal for epigenetic clocks not previously captured. We found retroelement-based epigenetic clocks were reversed during transient epigenetic reprogramming, accelerated in people living with HIV-1, and responsive to antiretroviral therapy. Our findings highlight the utility of retroelement-based biomarkers of aging and support a renewed emphasis on the role of retroelements in geroscience.

Discrepancies in readouts between Infinium MethylationEPIC v2.0 and v1.0 reflected in DNA methylation-based tools: implications and considerations for human population epigenetic studies

Background

The recently launched DNA methylation profiling platform, Illumina MethylationEPIC BeadChip Infinium microarray v2.0 (EPICv2), is highly correlated with measurements obtained from its predecessor MethylationEPIC BeadChip Infinium microarray v1.0 (EPICv1). However, the concordance between the two versions in the context of DNA methylation-based tools, including cell type deconvolution algorithms, epigenetic clocks, and inflammation and lifestyle biomarkers has not yet been investigated. To address this, we profiled DNA methylation on both EPIC versions using matched venous blood samples from individuals spanning early to late adulthood across four cohorts.

Findings

Within each cohort, samples primarily clustered by the EPIC version they were measured on. High concordance between EPIC versions at the array level, but variable concordance at the individual probe level was noted. Significant differences between versions in estimates from DNA methylation-based tools were observed, irrespective of the normalization method, with some nuanced differences across cohorts and tools. Adjusting for EPIC version or calculating estimates separately for each version largely mitigated these version-specific discordances.

Conclusions

Our work illustrates the importance of accounting for EPIC version differences in research scenarios, especially in meta-analyses and longitudinal studies, when samples profiled across different versions are harmonized. Alongside DNA methylation-based tools, our observations also have implications in interpretation of epigenome-wide association studies (EWAS) findings, when results obtained from one version are compared to another, particularly for probes that are poorly concordant between versions.

Cross-Tissue Comparison of Epigenetic Aging Clocks in Humans

Epigenetic clocks are a common group of tools used to measure biological aging - the progressive deterioration of cells, tissues and organs. Epigenetic clocks have been trained almost exclusively using blood-based tissues but there is growing interest in estimating epigenetic age using less-invasive oral-based tissues (i.e., buccal or saliva) in both research and commercial settings. However, differentiated cell types across body tissues exhibit unique DNA methylation landscapes and age-related alterations to the DNA methylome. Applying epigenetic clocks derived from blood-based tissues to estimate epigenetic age of oral-based tissues may introduce biases. We tested the within-person comparability of common epigenetic clocks across five tissue types: buccal epithelial, saliva, dry blood spots, buffy coat (i.e., leukocytes), and peripheral blood mononuclear cells. We tested 284 distinct tissue samples from 83 individuals aged 9-70 years. Overall, there were significant within-person differences in epigenetic clock estimates from oral-based versus blood-based tissues, with average differences of almost 30 years observed in some age clocks. In addition, most epigenetic clock estimates of blood-based tissues exhibited low correlation with estimates from oral-based tissues despite controlling for cellular proportions and other technical factors. Our findings indicate that application of blood-derived epigenetic clocks in oral-based tissues may not yield comparable estimates of epigenetic age, highlighting the need for careful consideration of tissue type when estimating epigenetic age.

mLiftOver: harmonizing data across Infinium DNA methylation platforms

MOTIVATION

Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip.

RESULTS

We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells.

AVAILABILITY AND IMPLEMENTATION

mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license.

Gestational exposure to environmental chemicals and epigenetic alterations in the placenta and cord blood mononuclear cells

Background

Exposure to environmental chemicals such as phthalates, phenols, and polycyclic aromatic hydrocarbons (PAHs) during pregnancy can increase the risk of adverse newborn outcomes. We explored the associations between maternal exposure to select environmental chemicals and DNA methylation in cord blood mononuclear cells (CBMC) and placental tissue (maternal and fetal sides) to identify potential mechanisms underlying these associations.

Method

This study included 75 pregnant individuals who planned to give birth at the University of Cincinnati Hospital between 2014 and 2017. Maternal urine samples during the delivery visit were collected and analyzed for 37 biomarkers of phenols (12), phthalates (13), phthalate replacements (4), and PAHs (8). Cord blood and placenta tissue (maternal and fetal sides) were also collected to measure the DNA methylation intensities using the Infinium HumanMethylation450K BeadChip. We used linear regression, adjusting for potential confounders, to assess CpG-specific methylation changes in CBMC (n = 54) and placenta [fetal (n = 67) and maternal (n = 68) sides] associated with gestational chemical exposures (29 of 37 biomarkers measured in this study). To account for multiple testing, we used a false discovery rate q-values < 0.05 and presented results by limiting results with a genomic inflation factor of 1±0.5. Additionally, gene set enrichment analysis was conducted using the Kyoto Encyclopedia of Genes and Genomics pathways.

Results

Among the 29 chemical biomarkers assessed for differential methylation, maternal concentrations of PAH metabolites (1-hydroxynaphthalene, 2-hydroxyfluorene, 4-hydroxyphenanthrene, 1-hydroxypyrene), monocarboxyisononyl phthalate, mono-3-carboxypropyl phthalate, and bisphenol A were associated with altered methylation in placenta (maternal or fetal side). Among exposure biomarkers associated with epigenetic changes, 1-hydroxynaphthalene, and mono-3-carboxypropyl phthalate were consistently associated with differential CpG methylation in the placenta. Gene enrichment analysis indicated that maternal 1-hydroxynaphthalene was associated with lipid metabolism and cellular processes of the placenta. Additionally, mono-3-carboxypropyl phthalate was associated with organismal systems and genetic information processing of the placenta.

Conclusion

Among the 29 chemical biomarkers assessed during delivery, 1-hydroxynaphthalene and mono-3-carboxypropyl phthalate were associated with DNA methylation in the placenta.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43682-024-00027-7.

MSA: scalable DNA methylation screening BeadChip for high-throughput trait association studies

The Infinium DNA Methylation BeadChips have significantly contributed to population-scale epigenetics research by enabling epigenome-wide trait association discoveries. Here, we design, describe, and experimentally verify a new iteration of this technology, the Methylation Screening Array (MSA), to focus on human trait screening and discovery. This array utilizes extensive data from previous Infinium platform-based epigenome-wide association studies (EWAS). It incorporates knowledge from the latest single-cell and cell type-resolution whole genome methylome profiles. The MSA is engineered to achieve scalable screening of epigenetics-trait association in an ultra-high sample throughput. Our design encompassed diverse human trait associations, including those with genetic, cellular, environmental, and demographical variables and human diseases such as genetic, neurodegenerative, cardiovascular, infectious, and immune diseases. We comprehensively evaluated this array's reproducibility, accuracy, and capacity for cell-type deconvolution and supporting 5-hydroxymethylation profiling in diverse human tissues. Our first atlas data using this platform uncovered the complex chromatin and tissue contexts of DNA modification variations and genetic variants linked to human phenotypes.

Genome-Wide Methylation Analysis in Two Wild-Type Non-Small Cell Lung Cancer Subgroups with Negative and High PD-L1 Expression

We conducted a pilot study to analyze the differential methylation status of 20 primary acinar adenocarcinomas of the lungs. These adenocarcinomas had to be wild type in mutation analysis and had either high (TPS > 50%; n = 10) or negative (TPS < 1%; n = 10) PD-L1 status to be integrated into our study. To examine the methylation of 866,895 specific sites, we utilized the Illumina Infinium EPIC bead chip array. Both hypermethylation and hypomethylation play significant roles in tumor development, progression, and metastasis. They also impact the formation of the tumor microenvironment, which plays a decisive role in tumor differentiation, epigenetics, dissemination, and immune evasion. The gained methylation patterns were correlated with PD-L1 expression. Our analysis has identified distinct methylation patterns in lung adenocarcinomas with high and negative PD-L1 expression. After analyzing the correlation between the methylation results of genes and promoters with their pathobiology, we found that tumors with high expression of PD-L1 tend to exhibit oncogenic effects through hypermethylation. On the other hand, tumors with negative PD-L1 expression show loss of their suppressor functions through hypomethylation. The suppressor functions of hypermethylated genes and promoters are ineffective compared to simultaneously activated dominant oncogenic mechanisms. The tumor microenvironment supports tumor growth in both groups.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

Ultrasensitive amplification-free quantification of a methyl CpG-rich cancer biomarker by single-molecule kinetic fingerprinting

The most well-studied epigenetic marker in humans is the 5-methyl modification of cytosine in DNA, which has great potential as a disease biomarker in liquid biopsies of cell-free DNA. Currently, quantification of DNA methylation relies heavily on bisulfite conversion followed by PCR amplification and NGS or microarray analysis. PCR is subject to potential bias in differential amplification of bisulfite-converted methylated versus unmethylated sequences. Here, we combine bisulfite conversion with single-molecule kinetic fingerprinting to develop an amplification-free assay for DNA methylation at the branched-chain amino acid transaminase 1 (BCAT1) promoter. Our assay selectively responds to methylated sequences with a limit of detection below 1 fM and a specificity of 99.9999%. Evaluating complex genomic DNA matrices, we reliably distinguish 2-5% DNA methylation at the BCAT1 promoter in whole blood DNA from completely unmethylated whole-genome amplified DNA. Taken together, these results demonstrate the feasibility and sensitivity of our amplification-free, single-molecule quantification approach to improve the early detection of methylated cancer DNA biomarkers.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling

BACKGROUND

The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples.

RESULTS

We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding.

CONCLUSIONS

Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip.

Creation and Validation of the First Infinium DNA Methylation Array for the Human Imprintome

Background

Differentially methylated imprint control regions (ICRs) regulate the monoallelic expression of imprinted genes. Their epigenetic dysregulation by environmental exposures throughout life results in the formation of common chronic diseases. Unfortunately, existing Infinium methylation arrays lack the ability to profile these regions adequately. Whole genome bisulfite sequencing (WGBS) is the unique method able to profile these regions, but it is very expensive and it requires not only a high coverage but it is also computationally intensive to assess those regions.

Findings

To address this deficiency, we developed a custom methylation array containing 22,819 probes. Among them, 9,757 probes map to 1,088 out of the 1,488 candidate ICRs recently described. To assess the performance of the array, we created matched samples processed with the Human Imprintome array and WGBS, which is the current standard method for assessing the methylation of the Human Imprintome. We compared the methylation levels from the shared CpG sites and obtained a mean R 2 = 0.569. We also created matched samples processed with the Human Imprintome array and the Infinium Methylation EPIC v2 array and obtained a mean R 2 = 0.796. Furthermore, replication experiments demonstrated high reliability (ICC: 0.799-0.945).

Conclusions

Our custom array will be useful for replicable and accurate assessment, mechanistic insight, and targeted investigation of ICRs. This tool should accelerate the discovery of ICRs associated with a wide range of diseases and exposures, and advance our understanding of genomic imprinting and its relevance in development and disease formation throughout the life course.