Ewastools: Infinium Human Methylation BeadChip pipeline for population epigenetics integrated into Galaxy

Comparison of Infinium MethylationEPIC v2.0 to v1.0 for human population epigenetics: considerations for addressing EPIC version differences in DNA methylation-based tools

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.

Findings

We profiled DNA methylation on both EPIC versions using matched venous blood samples from individuals spanning early to late adulthood across three cohorts. On combining the DNA methylomes of the cohorts, we observed that samples primarily clustered by the EPIC version they were measured on. Within each cohort, when we calculated cell type proportions, epigenetic age acceleration (EAA), rate of aging estimates, and biomarker scores for the matched samples on each version, we noted significant differences between EPICv1 and EPICv2 in the majority of these estimates. These differences were not significant, however, when estimates were adjusted for EPIC version or when EAAs were calculated separately for each EPIC version.

Conclusions

Our findings indicate that EPIC version differences predominantly explain DNA methylation variation and influence estimates of DNA methylation-based tools, and therefore we recommend caution when combining cohorts run on different versions. We demonstrate the importance of calculating DNA methylation-based estimates separately for each EPIC version or accounting for EPIC version either as a covariate in statistical models or by using version correction algorithms.

Telehealth Parenting Program and Salivary Epigenetic Biomarkers in Preschool Children With Developmental Delay: NIMHD Social Epigenomics Program

Importance

Children with developmental delays are at a heightened risk of experiencing mental health challenges, and this risk is exacerbated among racially minoritized children who face disproportionate adversity. Understanding the impact of parenting interventions on biological markers associated with these risks is crucial for mitigating long-term health disparities.

Objective

To examine the effect of 20 weeks of an internet-based parent-child interaction training (iPCIT) program on biomarkers associated with aging and chronic inflammation among preschoolers with developmental delay at 12-month follow-up.

Design, Setting, and Participants

An observational secondary analysis of data from a randomized clinical trial conducted from March 17, 2016, to December 15, 2020, to assess changes in salivary DNA methylation (DNAm)-derived biomarkers following iPCIT intervention. Participants were recruited from 3 Part C early intervention sites in a large southeastern US city. Eligible participants included children recruited within 3 months of their third birthday who had a Child Behavior Checklist Externalizing Problems T score greater than 60 and provided saliva in at least 1 study wave. Data analysis was conducted May 2023 to April 2024.

Intervention

Participants received either iPCIT (a telehealth therapeutic intervention focused on enhancing the parent-child relationship and addressing behavioral challenges in young children) or referrals as usual.

Main Outcomes and Measures

DNAm at the 12-month follow-up was assessed using the Infinium HumanMethylationEPIC Bead Chip Assay to derive biomarkers DunedinPACE, C-reactive protein (CRP), and interleukin-6 (IL-6). Analyses were intent-to-treat and used path analysis.

Results

A total of 71 children (mean [SD] age, 36.27 [0.61] months 51 male [71.8%] and 20 female [28.2%]) were analyzed, of whom 34 received iPCIT and 37 received referrals as usual. The iPCIT group had a slower pace of aging (β = 0.26; 95% CI, 0.06 to 0.50; P = .03) and less DNAm-derived CRP (β = 0.27; 95% CI, 0.05 to 0.49; P = .01) relative to the control condition at the 12-month follow-up. These associations remained significant after accounting for baseline DNAm score, child demographics, and symptom severity, and were independent of predicted buccal epithelial cell proportion for both DunedinPACE and CRP. There was no association with DNAm-derived IL-6 (β = 0.14; 95% CI, -0.08 to 0.36; P = .21).

Conclusions and Relevance

In this study of a parenting intervention, iPCIT, the association of intervention with decreased molecular markers of inflammation and biological aging suggests their potential to modify aspects of the biological embedding of stress. Understanding the systemic biological impact of such interventions offers insights into addressing health disparities and promoting resilience among vulnerable populations.

Trial Registration

ClinicalTrials.gov Identifier: NCT03260816.

regionalpcs: improved discovery of DNA methylation associations with complex traits

We have developed the regional principal components (rPCs) method, a novel approach for summarizing gene-level methylation. rPCs address the challenge of deciphering complex epigenetic mechanisms in diseases like Alzheimer's disease (AD). In contrast to traditional averaging, rPCs leverage principal components analysis to capture complex methylation patterns across gene regions. Our method demonstrated a 54% improvement in sensitivity over averaging in simulations, offering a robust framework for identifying subtle epigenetic variations. Applying rPCs to the AD brain methylation data in ROSMAP, combined with cell type deconvolution, we uncovered 838 differentially methylated genes associated with neuritic plaque burden-significantly outperforming conventional methods. Integrating methylation quantitative trait loci (meQTL) with genome-wide association studies (GWAS) identified 17 genes with potential causal roles in AD, including MS4A4A and PICALM. Our approach is available in the Bioconductor package regionalpcs, opening avenues for research and facilitating a deeper understanding of the epigenetic landscape in complex diseases.

Serotonin Transporter (SLC6A4) and FK506-Binding Protein 5 (FKBP5) Genotype and Methylation Relationships with Response to Meditation in Veterans with PTSD

Meditation-based interventions are novel and effective non-pharmacologic treatments for veterans with PTSD. We examined relationships between treatment response, early life trauma exposure, DNA polymorphisms, and methylation in the serotonin transporter (SLC6A4) and FK506-binding protein 5 (FKBP5) genes. DNA samples and clinical outcomes were examined in 72 veterans with PTSD who received meditation-based therapy in two separate studies of mindfulness-based stress reduction (MBSR) and Transcendental Meditation (TM). The PTSD Checklist was administered to assess symptoms at baseline and after 9 weeks of meditation intervention. We examined the SLC6A4 promoter (5HTTLPR_L/S insertion/deletion + rs25531_A/G) polymorphisms according to previously defined gene expression groups, and the FKBP5 variant rs1360780 previously associated with PTSD disease risk. Methylation for CpG sites of SLC6A4 (28 sites) and FKBP5 (45 sites) genes was quantified in DNA samples collected before and after treatment. The 5HTTLPR LALA high expression genotype was associated with greater symptom improvement in participants exposed to early life trauma (p = 0.015). Separately, pre to post-treatment change of DNA methylation in a group of nine FKBP5 CpG sites was associated with greater symptom improvement (OR = 2.8, 95% CI 1.1-7.1, p = 0.027). These findings build on a wealth of existing knowledge regarding epigenetic and genetic relationships with PTSD disease risk to highlight the potential importance of SLC6A4 and FKBP5 for treatment mechanisms and as biomarkers of symptom improvement.

The association between neighborhood deprivation and DNA methylation in an autopsy cohort

Previous research has found that living in a disadvantaged neighborhood is associated with poor health outcomes. Living in disadvantaged neighborhoods may alter inflammation and immune response in the body, which could be reflected in epigenetic mechanisms such as DNA methylation (DNAm). We used robust linear regression models to conduct an epigenome-wide association study examining the association between neighborhood deprivation (Area Deprivation Index; ADI), and DNAm in brain tissue from 159 donors enrolled in the Emory Goizueta Alzheimer's Disease Research Center (Georgia, USA). We found one CpG site (cg26514961, gene PLXNC1) significantly associated with ADI after controlling for covariates and multiple testing (p-value=5.0e-8). Effect modification by APOE ε4 was statistically significant for the top ten CpG sites from the EWAS of ADI, indicating that the observed associations between ADI and DNAm were mainly driven by donors who carried at least one APOE ε4 allele. Four of the top ten CpG sites showed a significant concordance between brain tissue and tissues that are easily accessible in living individuals (blood, buccal cells, saliva), including DNAm in cg26514961 (PLXNC1). Our study identified one CpG site (cg26514961, PLXNC1 gene) that was significantly associated with neighborhood deprivation in brain tissue. PLXNC1 is related to immune response, which may be one biological pathway how neighborhood conditions affect health. The concordance between brain and other tissues for our top CpG sites could make them potential candidates for biomarkers in living individuals.

A novel approach toward optimal workflow selection for DNA methylation biomarker discovery

DNA methylation is a major epigenetic modification involved in many physiological processes. Normal methylation patterns are disrupted in many diseases and methylation-based biomarkers have shown promise in several contexts. Marker discovery typically involves the analysis of publicly available DNA methylation data from high-throughput assays. Numerous methods for identification of differentially methylated biomarkers have been developed, making the need for best practices guidelines and context-specific analyses workflows exceedingly high. To this end, here we propose TASA, a novel method for simulating methylation array data in various scenarios. We then comprehensively assess different data analysis workflows using real and simulated data and suggest optimal start-to-finish analysis workflows. Our study demonstrates that the choice of analysis pipeline for DNA methylation-based marker discovery is crucial and different across different contexts.

Associations between the New DNA-Methylation-Based Telomere Length Estimator, the Mediterranean Diet and Genetics in a Spanish Population at High Cardiovascular Risk

Biological aging is a relevant risk factor for chronic diseases, and several indicators for measuring this factor have been proposed, with telomere length (TL) among the most studied. Oxidative stress may regulate telomere shortening, which is implicated in the increased risk. Using a novel estimator for TL, we examined whether adherence to the Mediterranean diet (MedDiet), a highly antioxidant-rich dietary pattern, is associated with longer TL. We determined TL using DNA methylation algorithms (DNAmTL) in 414 subjects at high cardiovascular risk from Spain. Adherence to the MedDiet was assessed by a validated score, and genetic variants in candidate genes and at the genome-wide level were analyzed. We observed several significant associations (p < 0.05) between DNAmTL and candidate genes (TERT, TERF2, RTEL1, and DCAF4), contributing to the validity of DNAmTL as a biomarker in this population. Higher adherence to the MedDiet was associated with lower odds of having a shorter TL in the whole sample (OR = 0.93; 95% CI: 0.85-0.99; p = 0.049 after fully multivariate adjustment). Nevertheless, this association was stronger in women than in men. Likewise, in women, we observed a direct association between adherence to the MedDiet score and DNAmTL as a continuous variable (beta = 0.015; SE: 0.005; p = 0.003), indicating that a one-point increase in adherence was related to an average increase of 0.015 ± 0.005 kb in TL. Upon examination of specific dietary items within the global score, we found that fruits, fish, "sofrito", and whole grains exhibited the strongest associations in women. The novel score combining these items was significantly associated in the whole population. In the genome-wide association study (GWAS), we identified ten polymorphisms at the suggestive level of significance (p < 1 × 10-5) for DNAmTL (intergenics, in the IQSEC1, NCAPG2, and ABI3BP genes) and detected some gene-MedDiet modulations on DNAmTL. As this is the first study analyzing the DNAmTL estimator, genetics, and modulation by the MedDiet, more studies are needed to confirm these findings.

The application of epiphenotyping approaches to DNA methylation array studies of the human placenta

BACKGROUND

Genome-wide DNA methylation (DNAme) profiling of the placenta with Illumina Infinium Methylation bead arrays is often used to explore the connections between in utero exposures, placental pathology, and fetal development. However, many technical and biological factors can lead to signals of DNAme variation between samples and between cohorts, and understanding and accounting for these factors is essential to ensure meaningful and replicable data analysis. Recently, "epiphenotyping" approaches have been developed whereby DNAme data can be used to impute information about phenotypic variables such as gestational age, sex, cell composition, and ancestry. These epiphenotypes offer avenues to compare phenotypic data across cohorts, and to understand how phenotypic variables relate to DNAme variability. However, the relationships between placental epiphenotyping variables and other technical and biological variables, and their application to downstream epigenome analyses, have not been well studied.

RESULTS

Using DNAme data from 204 placentas across three cohorts, we applied the PlaNET R package to estimate epiphenotypes gestational age, ancestry, and cell composition in these samples. PlaNET ancestry estimates were highly correlated with independent polymorphic ancestry-informative markers, and epigenetic gestational age, on average, was estimated within 4 days of reported gestational age, underscoring the accuracy of these tools. Cell composition estimates varied both within and between cohorts, as well as over very long placental processing times. Interestingly, the ratio of cytotrophoblast to syncytiotrophoblast proportion decreased with increasing gestational age, and differed slightly by both maternal ethnicity (lower in white vs. non-white) and genetic ancestry (lower in higher probability European ancestry). The cohort of origin and cytotrophoblast proportion were the largest drivers of DNAme variation in this dataset, based on their associations with the first principal component.

CONCLUSIONS

This work confirms that cohort, array (technical) batch, cell type proportion, self-reported ethnicity, genetic ancestry, and biological sex are important variables to consider in any analyses of Illumina DNAme data. We further demonstrate the specific utility of epiphenotyping tools developed for use with placental DNAme data, and show that these variables (i) provide an independent check of clinically obtained data and (ii) provide a robust approach to compare variables across different datasets. Finally, we present a general framework for the processing and analysis of placental DNAme data, integrating the epiphenotype variables discussed here.

The application of epiphenotyping approaches to DNA methylation array studies of the human placenta

Background : Genome-wide DNA methylation (DNAme) profiling of the placenta with Illumina Infinium Methylation bead arrays is often used to explore the connections between in utero exposures, placental pathology, and fetal development. However, many technical and biological factors can lead to signals of DNAme variation between samples and between cohorts, and understanding and accounting for these factors is essential to ensure meaningful and replicable data analysis. Recently, "epiphenotyping" approaches have been developed whereby DNAme data can be used to impute information about phenotypic variables such as gestational age, sex, cell composition, and ancestry. These epiphenotypes offer avenues to compare phenotypic data across cohorts, and to understand how phenotypic variables relate to DNAme variability. However, the relationships between placental epiphenotyping variables and other technical and biological variables, and their application to downstream epigenome analyses, have not been well studied. Results : Using DNAme data from 204 placentas across three cohorts, we applied the PlaNET R package to estimate epiphenotypes gestational age, ancestry, and cell composition in these samples. PlaNET ancestry estimates were highly correlated with independent polymorphic ancestry informative markers, and epigenetic gestational age, on average, was estimated within 4 days of reported gestational age, underscoring the accuracy of these tools. Cell composition estimates varied both within and between cohorts, but reassuringly were robust to placental processing time. Interestingly, the ratio of cytotrophoblast to syncytiotrophoblast proportion decreased with increasing gestational age, and differed slightly by both maternal ethnicity (lower in white vs. non-white) and genetic ancestry (lower in higher probability European ancestry). The cohort of origin and cytotrophoblast proportion were the largest drivers of DNAme variation in this dataset, based on their associations with the first principal component. Conclusions : This work confirms that cohort, array (technical) batch, cell type proportion, self-reported ethnicity, genetic ancestry, and biological sex are important variables to consider in any analyses of Illumina DNAme data. Further, we demonstrate that estimating epiphenotype variables from the DNAme data itself, when possible, provides both an independent check of clinically-obtained data and can provide a robust approach to compare variables across different datasets.

DNA-Methylation Signatures of Tobacco Smoking in a High Cardiovascular Risk Population: Modulation by the Mediterranean Diet

Biomarkers based on DNA methylation are relevant in the field of environmental health for precision health. Although tobacco smoking is one of the factors with a strong and consistent impact on DNA methylation, there are very few studies analyzing its methylation signature in southern European populations and none examining its modulation by the Mediterranean diet at the epigenome-wide level. We examined blood methylation smoking signatures on the EPIC 850 K array in this population (n = 414 high cardiovascular risk subjects). Epigenome-wide methylation studies (EWASs) were performed analyzing differential methylation CpG sites by smoking status (never, former, and current smokers) and the modulation by adherence to a Mediterranean diet score was explored. Gene-set enrichment analysis was performed for biological and functional interpretation. The predictive value of the top differentially methylated CpGs was analyzed using receiver operative curves. We characterized the DNA methylation signature of smoking in this Mediterranean population by identifying 46 differentially methylated CpGs at the EWAS level in the whole population. The strongest association was observed at the cg21566642 (p = 2.2 × 10-32) in the 2q37.1 region. We also detected other CpGs that have been consistently reported in prior research and discovered some novel differentially methylated CpG sites in subgroup analyses. In addition, we found distinct methylation profiles based on the adherence to the Mediterranean diet. Particularly, we obtained a significant interaction between smoking and diet modulating the cg5575921 methylation in the AHRR gene. In conclusion, we have characterized biomarkers of the methylation signature of tobacco smoking in this population, and suggest that the Mediterranean diet can increase methylation of certain hypomethylated sites.

DNA quality evaluation of formalin-fixed paraffin-embedded heart tissue for DNA methylation array analysis

Archived formalin-fixed and paraffin-embedded (FFPE) heart tissue from autopsied individuals represents an important resource for investigating the DNA methylation of heart tissue of deceased individuals. The DNA quality of FFPE tissue from autopsies may be decreased, affecting the DNA methylation measurements. Therefore, inexpensive screening methods for estimating DNA quality are valuable. We investigated the correlation between the DNA quality of archived FFPE heart tissue examined with the Illumina Infinium HD FFPE QC assay (Infinium QC) and Thermo Fisher's Quantifiler Trio DNA Quantification kit (QuantifilerTrio), respectively, and the amount of usable DNA methylation data as measured by the probe detection rate (probe DR) obtained with the Illumina Infinium MethylationEPIC array. We observed a high correlation (r² = 0.75; p < 10^-11) between the QuantifilerTrio degradation index, DI, and the amount of usable DNA methylation data analysed with SeSAMe, whereas a much weaker correlation was observed between the Infinium QC and SeSAMe probe DR (r² = 0.17; p < 0.05). Based on the results, QuantifilerTrio DI seems to predict the proportion of usable DNA methylation data analysed with the Illumina Infinium MethylationEPIC array and SeSAMe by a linear model: SeSAMe probe DR = 0.80-log₁₀(DI) × 0.25.

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.

Methylation differences in Alzheimer's disease neuropathologic change in the aged human brain

Alzheimer's disease (AD) is the most common cause of dementia with advancing age as its strongest risk factor. AD neuropathologic change (ADNC) is known to be associated with numerous DNA methylation changes in the human brain, but the oldest old (> 90 years) have so far been underrepresented in epigenetic studies of ADNC. Our study participants were individuals aged over 90 years (n = 47) from The 90+ Study. We analyzed DNA methylation from bulk samples in eight precisely dissected regions of the human brain: middle frontal gyrus, cingulate gyrus, entorhinal cortex, dentate gyrus, CA1, substantia nigra, locus coeruleus and cerebellar cortex. We deconvolved our bulk data into cell-type-specific (CTS) signals using computational methods. CTS methylation differences were analyzed across different levels of ADNC. The highest amount of ADNC related methylation differences was found in the dentate gyrus, a region that has so far been underrepresented in large scale multi-omic studies. In neurons of the dentate gyrus, DNA methylation significantly differed with increased burden of amyloid beta (Aβ) plaques at 5897 promoter regions of protein-coding genes. Amongst these, higher Aβ plaque burden was associated with promoter hypomethylation of the Presenilin enhancer 2 (PEN-2) gene, one of the rate limiting genes in the formation of gamma-secretase, a multicomponent complex that is responsible in part for the endoproteolytic cleavage of amyloid precursor protein into Aβ peptides. In addition to novel ADNC related DNA methylation changes, we present the most detailed array-based methylation survey of the old aged human brain to date. Our open-sourced dataset can serve as a brain region reference panel for future studies and help advance research in aging and neurodegenerative diseases.

Prenatal exposure to phthalates and peripheral blood and buccal epithelial DNA methylation in infants: An epigenome-wide association study

BACKGROUND

Prenatal exposure to phthalates has been associated with adverse health and neurodevelopmental outcomes. DNA methylation (DNAm) alterations may be a mechanism underlying these effects, but prior investigations of prenatal exposure to phthalates and neonatal DNAm profiles are limited to placental tissue and umbilical cord blood.

OBJECTIVE

Conduct an epigenome-wide association study (EWAS) of the associations between prenatal exposure to phthalates and DNAm in two accessible infant tissues, venous buffy coat blood and buccal epithelial cells (BECs).

METHODS

Participants included 152 maternal-infant pairs from the Alberta Pregnancy Outcomes and Nutrition (APrON) study. Maternal second trimester urine samples were analyzed for nine phthalate metabolites. Blood (n = 74) or BECs (n = 78) were collected from 3-month-old infants and profiled for DNAm using the Infinium HumanMethylation450 (450K) BeadChip. Robust linear regressions were used to investigate the associations between high (HMWPs) and low molecular weight phthalates (LMWPs) and change in methylation levels at variable Cytosine-phosphate-Guanine (CpG) sites in infant tissues, as well as the sensitivity of associations to potential confounders.

RESULTS

One candidate CpG in gene RNF39 reported by a previous study examining prenatal exposure to phthalates and cord blood DNAm was replicated. The EWAS identified 12 high-confidence CpGs in blood and another 12 in BECs associated with HMWPs and/or LMWPs. Prenatal exposure to bisphenol A (BPA) associated with two of the CpGs associated with HMWPs in BECs.

DISCUSSION

Prenatal exposure to phthalates was associated with DNAm variation at CpGs annotated to genes associated with endocrine hormone activity (i.e., SLCO4A1, TPO), immune pathways and DNA damage (i.e., RASGEF1B, KAZN, HLA-A, MYO18A, DIP2C, C1or109), and neurodevelopment (i.e., AMPH, NOTCH3, DNAJC5). Future studies that characterize the stability of these associations in larger samples, multiple cohorts, across tissues, and investigate the potential associations between these biomarkers and relevant health and neurodevelopmental outcomes are needed.

Multi-Omics Profiling Approach to Asthma: An Evolving Paradigm

Asthma is a complex multifactorial and heterogeneous respiratory disease. Although genetics is a strong risk factor of asthma, external and internal exposures and their interactions with genetic factors also play important roles in the pathophysiology of asthma. Over the past decades, the application of high-throughput omics approaches has emerged and been applied to the field of asthma research for screening biomarkers such as genes, transcript, proteins, and metabolites in an unbiased fashion. Leveraging large-scale studies representative of diverse population-based omics data and integrating with clinical data has led to better profiling of asthma risk. Yet, to date, no omic-driven endotypes have been translated into clinical practice and management of asthma. In this article, we provide an overview of the current status of omics studies of asthma, namely, genomics, transcriptomics, epigenomics, proteomics, exposomics, and metabolomics. The current development of the multi-omics integrations of asthma is also briefly discussed. Biomarker discovery following multi-omics profiling could be challenging but useful for better disease phenotyping and endotyping that can translate into advances in asthma management and clinical care, ultimately leading to successful precision medicine approaches.

Ewastools: Infinium Human Methylation BeadChip pipeline for population epigenetics integrated into Galaxy

BACKGROUND

Infinium Human Methylation BeadChip is an array platform for complex evaluation of DNA methylation at an individual CpG locus in the human genome based on Illumina's bead technology and is one of the most common techniques used in epigenome-wide association studies. Finding associations between epigenetic variation and phenotype is a significant challenge in biomedical research. The newest version, HumanMethylationEPIC, quantifies the DNA methylation level of 850,000 CpG sites, while the previous versions, HumanMethylation450 and HumanMethylation27, measured >450,000 and 27,000 loci, respectively. Although a number of bioinformatics tools have been developed to analyse this assay, they require some programming skills and experience in order to be usable.

RESULTS

We have developed a pipeline for the Galaxy platform for those without experience aimed at DNA methylation analysis using the Infinium Human Methylation BeadChip. Our tool is integrated into Galaxy (http://galaxyproject.org), a web-based platform. This allows users to analyse data from the Infinium Human Methylation BeadChip in the easiest possible way.

CONCLUSIONS

The pipeline provides a group of integrated analytical methods wrapped into an easy-to-use interface. Our tool is available from the Galaxy ToolShed, GitHub repository, and also as a Docker image. The aim of this project is to make Infinium Human Methylation BeadChip analysis more flexible and accessible to everyone.