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

* Amy P. Webster
* Simone Ecker
* Ismail Moghul
* Pawan Dhami
* Sarah Marzi
* Dirk S. Paul
* Michelle Kuxhausen
* Stephanie J. Lee
* Stephen R. Spellman
* Tao Wang
* Andrew Feber
* Vardhman Rakyan
* Karl S. Peggs
* Stephan Beck

## Abstract

Allogeneic hematopoietic cell transplantation (HCT) is used to treat many blood-based disorders and malignancies. While this is an effective treatment, it can result in serious adverse events, such as the development of acute graft-versus-host disease (aGVHD). This study aimed to develop a donor-specific epigenetic classifier that could be used in donor selection in HCT to reduce the incidence of aGVHD.

The discovery cohort of the study consisted of 288 donors from a population receiving HLA-A, -B, -C and -DRB1 matched unrelated donor HCT with T cell replete peripheral blood stem cell grafts for treatment of acute leukaemia or myelodysplastic syndromes after myeloablative conditioning. Donors were selected based on recipient aGVHD outcome; this cohort consisted of 144 cases with aGVHD grades III-IV and 144 controls with no aGVHD that survived at least 100 days post-HCT matched for sex, age, disease and GVHD prophylaxis.

Genome-wide DNA methylation was assessed using the Infinium Methylation EPIC BeadChip (Illumina), measuring CpG methylation at >850,000 sites across the genome. Following quality control, pre-processing and exploratory analyses, we applied a machine learning algorithm (Random Forest) to identify CpG sites predictive of aGVHD. Receiver operating characteristic (ROC) curve analysis of these sites resulted in a classifier with an encouraging area under the ROC curve (AUC) of 0.91.

To test this classifier, we used an independent validation cohort (n=288) selected using the same criteria as the discovery cohort. Different attempts to validate the classifier using the independent validation cohort failed with the AUC falling to 0.51. These results indicate that donor DNA methylation may not be a suitable predictor of aGVHD in an HCT setting involving unrelated donors, despite the initial promising results in the discovery cohort.

Our work highlights the importance of independent validation of machine learning classifiers, particularly when developing classifiers intended for clinical use.

## Introduction

In the past six decades, allogeneic hematopoietic cell transplantation (HCT) has become a cornerstone of treatment for haematological malignancies and is still often considered the only curative option1. Despite advances in the precision of HLA matching in unrelated donor selection and supportive care leading to ongoing improvements in HCT outcomes, severe graft versus host disease (GVHD) regularly occurs, increasing the risk of morbidity and mortality2. Acute GVHD occurs when the donor immune cells attack healthy tissue in the graft recipient, causing a range of inflammatory lesions which primarily affect the skin and digestive organs. Acute GVHD (aGVHD) typically occurs within 100 days of transplant. While the incidence has decreased in the last decade due to better HLA matching of donors, aGVHD still affects ∼30-50% of allogeneic HCT recipients3, making the prevention of aGVHD an important area of research.

DNA methylation is a stable modification of the DNA which can influence gene expression without altering the underlying genetic sequence. DNA methylation has an emerging role in precision medicine due to the environmental and developmental exposures it can capture. Several factors associated with the development of aGVHD are also known to influence the epigenome, including age4,5, sex6 and viral infections7. Despite the relative infancy of the field, DNA methylation classifiers predictive of clinical outcome are now being used in the clinic, notably in oncology to guide treatment of brain tumours8,9. The development of machine learning algorithms and increasing size of datasets has also allowed improvement in the development of such classifiers for early diagnosis and determining subtypes of disease10.

In 2015, we published a pilot study investigating DNA methylation as a potential classifier of aGVHD in HCT of HLA matched sibling pairs11. In that study, we assessed DNA methylation in a cohort of 85 HCT donors selected based on recipient outcome, identifying 31 DNA methylation markers associated with aGVHD severity in graft recipients. In internal cross-validation these markers showed strong predictive performance (AUC=0.98) indicating the potential utility of DNA methylation in improving donor selection in sibling HCT. The purpose of the current study was to investigate if DNA methylation is also predictive of outcome in HLA matched unrelated donor-recipient pairs, which constitute a much greater proportion of HCTs. To do this, we assessed genome-wide DNA methylation of 576 individuals recruited from the Center for International Blood and Marrow Transplant Research (CIBMTR). The scale and quality of annotation of the CIBMTR donor collection allowed us to use stringent selection criteria to minimise confounding and increase our power to detect methylation differences which were predictive of the development of aGVHD following HCT.

## Methods

### Study population

The discovery study cohort consisted of 288 HLA-A, -B, -C and -DRB1 matched, unrelated donor transplants reported to the CIBMTR that had pre-transplant donor peripheral blood samples available through the CIBMTR Research Repository. Patients received a transplant between 2002 to 2017 for acute lymphoblastic leukemia (ALL), acute myelogenous leukemia (AML) and myelodysplastic syndromes (MDS) using T-cell replete peripheral blood stem cell grafts, myeloablative conditioning and tacrolimus with methotrexate or mycophenylate mofetil based GVHD prophylaxis. The population was selected as a case-control cohort with 144 cases that developed aGVHD III-IV and controls with no aGVHD. Cases and controls were matched for sex, age, disease and GVHD prophylaxis. Donors were all self-reported as Caucasian.

The validation cohort (n=288) was selected using the same criteria. Using a previously described method12, power calculations for the discovery study using the EPIC array for genome-wide methylation measurement were performed with genome-wide significance set at 1×10−6. Sample groups of 140 donors matched to recipients with grade III-IV aGvHD, and 140 donors matched to recipients with no aGVHD, would give us 88% power to detect a methylation difference of 10% between the groups, and 100% power to detect methylation differences of 25%. Several additional samples for each group were profiled to ensure adequate power even if samples were removed during quality control.

### Samples

Genomic DNA was extracted from whole blood samples obtained from CIBMTR using the QIAamp DNA Blood Mini Kit (Qiagen) at the UCL Pathology Department (discovery study) and the UCL Genomics facility (validation study). The quality and concentration of DNA was assessed using NanoDrop and Qubit (Thermo Fisher).

### Genome-wide DNA methylation profiling

For each sample, 500ng high-quality DNA was bisulphite converted using the EZ DNA methylation kit (Zymo Research), using alternative incubation conditions recommended for Illumina methylation arrays. Methylation was subsequently analysed using the Infinium MethylationEPIC array (Illumina) measuring CpG methylation at >850,000 sites across the genome. Array preparation was performed at the UCL Genomics facility using standard operating procedures. Discovery and validation cohorts were processed independent at different timepoints, but within each cohort batches were minimised by distributing comparison groups evenly across BeadChips and position on BeadChip.

### Analysis overview

All analyses were performed in R version 3.6. Samples remaining following quality control (n=282 for discovery cohort and 288 for validation cohorts) were normalised using SWAN, then problematic probes were removed including those with a detection P value >0.01, probes with a beadcount <3 in more than 5% of samples13, non-cg probes, probes containing any common SNPs in dbSNP14 and probes mapping to the X or Y chromosomes. Singular variable decomposition (SVD)15 and principal components analysis (PCA) were used to assess batch effects in the data, which were subsequently adjusted for using Combat16. Cell composition was estimated and adjusted for using the Houseman method17 as implemented in ChAMP18,19, estimating cell proportions using the Reinius reference dataset20. Differentially methylated positions (DMPs) were assessed using a linear model in Limma21.

Machine learning analysis was performed using the random forest method22 as implemented in the RandomForest package. Instead of using all CpG sites as input for the RandomForest analysis, a subset of 10,000 CpG sites were selected through feature selection.

A supervised approach was used, where DMPs were identified in the discovery cohort using a linear model and the top 10,000 ranked probes were used as input for the random forest analysis. An alternative unsupervised approach was also carried out where the top 10,000 probes with the largest overall beta variance across all samples in the discovery cohort were used as input for the random forest analysis.

In both cases, the classifiers were then tested on matched probe sets from the validation cohort, and sensitivity and specificity of the classifiers were calculated.

### Data availability

The participants involved in the study had been recruited under different consents which require different levels of data access. According to consent given, the corresponding data are being made available in a three-tiered data access approach:

1.  Processed data (beta matrix) for all individuals (n=570) are available from the open access ‘Gene Expression Omnibus’ under accession number GSE196696. To reduce the chance of reidentification, all non-cg probes, including SNP targeting rs probes have been removed. The data are provided in both raw (unnormalized) and SWAN normalised formats.

2.  Raw data (IDAT files) are available for individuals with appropriate consent (n=403 in total) from the controlled access ‘European Genome-Phenome Archive’ under accession number EGAS00001006033.

3.  Raw data (IDAT files) and associated phenotype information are available for all individuals included in this study (n=570) directly from CIBMTR. Data are available under controlled access release upon reasonable request and execution of a data use agreement. Requests should be submitted to CIBMTR at info-request{at}mcw.edu and include the study reference IB17-04.

## Results

### Study Design

Unrelated donor-recipient pairs undergoing HCT were selected from the CIBMTR Research Repository, based on the aGVHD outcomes in recipients (Figure 1). Blood-based DNA methylation from donors was assessed using the Illumina EPIC arrays. Methylation differences were assessed, and random forest analysis was used to test for the presence of a classifier of aGVHD outcome.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/10/2022.03.08.22272071/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2022/03/10/2022.03.08.22272071/F1)

Figure 1: Study Design.
Unrelated donor-recipient pairs were selected based on the outcome of recipients following HCT. DNA methylation levels were assessed in donors associated with no (Grade 0) or severe (Grades 3-4) aGVHD in recipients. Donor-recipient pairs were HLA matched, and comparison groups were matched for sex, age, disease and GVHD prophylaxis. Feature selection reduced the number of probes in the discovery dataset to 10,000 for input to random forest analyses, and this classifier was subsequently tested in the validation cohort following pre-processing of data and refinement to the same set of probes.

### Study population

Unrelated donor-recipient pairs were selected by CIBMTR using stringent criteria as described in methods, resulting in 282 individuals in the discovery cohort following initial data quality control, and 288 individuals in the validation cohort. The resulting cohorts were well matched for characteristics that can influence DNA methylation profile, including age and sex (as shown in Tables 1 and 2).

View this table:
[Table 1:](http://medrxiv.org/content/early/2022/03/10/2022.03.08.22272071/T1)

Table 1: Discovery and validation cohort characteristics.
Characteristics of adult patients undergoing first allogeneic PB HCT for acute leukemia or MDS from an 8/8 HLA-matched unrelated donor between 2000-2016 with available donor blood samples, as reported to the CIBMTR. Restricted to Caucasian donors, myeloablative preparative regimens, **no ATG/Campath** and patients surviving >100 days with no aGVHD or those that developed grades III-IV aGVHD at any time post-HCT. Donors were matched between comparison groups based on sex and age by decade.

The discovery cohort was well matched for disease, with no significant difference in proportion of AML, ALL and MDS between comparison groups (p=0.339). Median recipient ages for the no/severe aGVHD groups were 45 (range 19-76) and 47 (range 18-72), respectively. There was no significant difference for recipient sex (p=0.716) or ethnicity (p=0.113) across comparison groups. Donors were well matched across comparison groups for sex (p=0.585), however, there was a difference in median age (p=0.003), though this was not apparent when individuals were stratified into age brackets (p=0.090). There were no significant differences across comparison groups for donor/recipient ABO type, blood type, Rh factor, CMV status or sex match.

The validation cohort had a significant difference in proportions of these diseases across comparison groups (p=0.02). The median recipient ages for the no/severe aGVHD groups in the validation cohort was 49 (range 20-75) and 50 (range 19-71). There was no significant difference in the recipient age distribution across comparison groups (p=0.998). There was no difference in recipient sex across the comparison groups (41% female recipients, p=1.0). Donors were well matched across comparison groups for sex (p=0.063) and median age (p=0.076). There were no significant differences in ethnicity, donor/recipient ABO type, blood type, Rh factor, CMV status or sex match across groups. There were differences in conditioning regimen across comparison groups (p<0.001).

### Data exploration and pre-processing

Following sample removal, quality control plots showed that the 282 individuals remaining in the discovery dataset and 288 individuals remaining in the validation dataset had very high quality methylation profiles (Supplementary Figure 1). Following probe filtering, 661,114 probes remained in the discovery dataset. Singular Value Decomposition (SVD) and principal components analysis (PCA) indicated that estimated ‘cell composition’, ‘Slide/BeadChip’ and ‘Array’ batch effects were having the largest impact on the data (Supplementary Figures 2-3), which were subsequently adjusted for using ChAMP cell composition correction and ComBat adjustment respectively. Cell type proportions were estimated for each group using the DNA methylation profiles and were found to be well balanced in each cohort with no significant difference between groups (Supplementary Table 1).

### Differential methylation analysis

No CpG sites passed a false discovery rate adjusted p-value significance threshold of 0.05 during DMP analysis when comparing the ‘no aGVHD’ group to the ‘severe aGVHD’ group. As the main batch and confounding effects of slide, array and cell composition had been previously adjusted in the dataset, no additional covariates were included during linear regression.

### Classifier generation from discovery data

Random forest analysis was performed on two sets of probes; the unsupervised analysis using the top-ranked 10,000 most variable probes, which all had a beta variance of >33% across all samples. The supervised analysis used the top 10,000 probes resulting from the linear model DMP analysis, though none passed statistical significance these were considered sites with putative methylation differences. Random forest analysis was run with 500 trees, with 100 variables tested at each split for both analysis approaches.

The high variability classifier showed very poor performance, with an out-of-bag (OOB) estimate of error rate of 45.39% and area under the curve (AUC) of 0.516 during internal cross-validation of the discovery dataset (Supplementary Figure 4). The differential methylation dataset produced an initially promising classifier with an OOB estimate of error rate of 14.89% and an AUC of 0.913 (Figure 2).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/10/2022.03.08.22272071/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2022/03/10/2022.03.08.22272071/F2)

Figure 2: ROC curve of classifier performance of the unsupervised Random Forest Classifier.
The figure shows the performance of the differential methylation (supervised approach) classifier which used the top 10,000 most differentially methylated CpG sites as input, during internal cross validation on the training dataset (blue line). The performance of the differential methylation classifier on the independent validation cohort is indicated by the orange line, which had an AUC of 0.508, a sensitivity of 90.97% with a very poor specificity of 6.25%. While initially this differential methylation-based classifier appeared encouraging with the discovery cohort, the classifier did not perform well during validation analyses.

During validation analysis, the matched CpG sites used as input to the original random forest training analysis were extracted from the validation dataset as all probes present in training analyses are required as input for validation. Validation analyses indicated that the differential methylation classifier had a sensitivity of 90.97% but a specificity of just 6.25%, and an AUC of 0.508. This is driven by an over-prediction of the ‘severe aGVHD’ group in the independent validation cohort, resulting in many false positive predictions. The unsupervised differential variability classifier also had an extremely poor performance in the validation cohort, with a sensitivity of just 50%, a specificity of 51.39% and an AUC of 0.523. As such, neither of these approaches yielded a useful classifier.

## Discussion

Recently developed predictors of aGVHD using clinical variables have had modest success with an AUC of ∼ 0.623, however this indicated that biological markers of gene expression, such as epigenetic markers, could provide additional insight to improve prediction of aGVHD. This was also supported by the recent finding that hypermethylation of the TP53 gene in HCT recipients was found to correlate with relapse of myelodysplastic syndromes following transplantation, indicating recipient-based DNA methylation could be predictive of outcomes during HCT24. As DNA methylation levels reflect both the underlying genetic sequence and factors known to be associated with aGVHD development (such as donor age, sex and cytomegalovirus serostatus), we hypothesised they would be a strong candidate for classifier identification. Our initial study focused on sibling donor-recipient pairs, in which a DNA methylation classifier of aGVHD development was identified in the blood of donors11. In the current study, we tested if DNA methylation as measured by EPIC arrays is also predictive of aGVHD in unrelated donor-recipient pairs and found that it is not.

There are several potential technical and biological reasons that a robust classifier of aGVHD was not identified in this study. Firstly, while the study performed was shown to have power to detect larger methylation differences of >10%, the relatively small sample size of the discovery cohort (n=280) and validation cohort (n=288) may have limited our ability to detect more subtle methylation differences. In the future, larger scale studies may provide increased power to detect such differences.

Secondly, the tissue we investigated was peripheral blood of donors which was intended to act as a surrogate tissue reflecting outcome. DNA methylation profiles are known to be highly cell type specific25, and while blood based DNA methylation may reflect certain exposures and factors associated with aGVHD development, it is possible that a specific cellular subtype which is not present in the whole blood of donors is responsible for the development of aGVHD and as such would not be reflected in the methylation profile.

Another possibility is that the specific cell type which is causing aGVHD could be present in whole blood, but in small proportions, making the signal significantly diluted by other more prominent cell types. Indeed, in the current analysis, cell composition was the biggest driver of variation in the data, and though this was balanced overall between the comparison groups and adjusted for in the data analysis, it could have been a confounding factor in the study, or subtle methylation effects could have been lost during adjustment. In the future, methylation analysis of individual cell types isolated from stem cell grafts may provide more insight into DNA methylation differences driving the development of aGVHD. While this approach would provide a more refined methylation measurement, it would be a significantly less practical approach for a clinical test, limiting the utility for optimising donor selection as usually these cells would only be collected once a donor is committed.

A classifier of aGVHD development was identified in our previously published work, which investigated donor DNA methylation from sibling HCT. A potential reason a similar biomarker was not identified in this cohort is that it could have been specific to sibling transplants, which generally have a lower incidence of aGVHD which may be driven more by extrinsic factors which influence DNA methylation, while aGVHD following HCT from an unrelated donor may be driven more by genetic factors. There may also be an issue of ‘epigenetic compatibility’, with donors and recipients varying in epigenetic profile inciting the initiation of aGVHD in certain individuals, without this being driven by a specifically differentially methylated gene or pathway. This would explain why a classifier was not identified in the current study, as the epigenetic marks conferring risk of aGVHD would be different for each individual. In the future, studies investigating the DNA methylation of both donors and recipients during HCT could provide more insight into this possibility.

When considering the clinical context of the development of aGVHD, it is likely the end result of a complicated clinical setting with multiple donor and recipient factors affecting the outcome. If the epigenetic pattern was highly predictive, it might infer that the occurrence of severe aGVHD is pre-ordained just by donor factors, which seems biologically unlikely.

On a technical level, this study has also demonstrated the importance of careful development and testing of analysis pipelines for methylation studies, in particular when applying complex machine learning methods to datasets. Our initial findings indicated a robust classifier might be present within the dataset, a finding which was amplified when data was pre-processed as a single batch with subsequent splitting of the dataset and internal cross validation. While our validation dataset was of exceptionally high quality and donors included were matched to a very high degree with the discovery cohort, the classifier was not validated even with extensive optimisation and testing of alternate pipeline settings. This demonstrates that even with the identification of a promising and robust classifier in a well-designed study, independent validation is critical26, and such validation datasets need to be generated completely independently with unique individuals and pre-processed separately to the training/discovery cohort. This also better mimics the experimental realities of clinical classifier use, making any findings that do stand up to the validation process more robust and clinically useful.

## Conclusions

In this study, we performed the definitive investigation of donor-derived blood-based DNA methylation as a classifier of aGVHD outcome in HCT and found that donor DNA methylation as assessed by methylation arrays is not a strong candidate for prediction of aGVHD. It is possible that other methylation signals exist which might improve our understanding of the development of aGVHD in these cohorts, which we plan to investigate in the future. We have also highlighted the importance of study design and well-designed independent validation of methylation differences especially when applying machine learning approaches.

## Supporting information

Supplementary Figure 1 [[supplements/272071_file03.docx]](pending:yes)

## Data Availability

The participants involved in the study had been recruited under different consents which require different levels of data access. According to consent given, the corresponding data are being made available in a three-tiered data access approach: 1. Processed data (beta matrix) for all individuals (n=570) are available from the open access 'Gene Expression Omnibus' under accession number GSE196696. To reduce the chance of reidentification, all non-cg probes, including SNP targeting rs probes have been removed. The data are provided in both raw (unnormalized) and SWAN normalised formats. 2. Raw data (IDAT files) are available for individuals with appropriate consent (n=403 in total) from the controlled access 'European Genome-Phenome Archive' under accession number EGAxxxxx. 3. Raw data (IDAT files) and associated phenotype information are available for all individuals included in this study (n=570) directly from CIBMTR. Data are available under controlled access release upon reasonable request and execution of a data use agreement. Requests should be submitted to CIBMTR at info-request@mcw.edu and include the study reference IB17-04.

[https://www.ncbi.nlm.nih.gov/geo/](https://www.ncbi.nlm.nih.gov/geo/) 

[https://ega-archive.org](https://ega-archive.org) 

## Conflict of Interest

The authors declare no conflicts of interest.

## Acknowledgements

The authors would like to thank UCL Genomics and UCL Pathology for their support with DNA extraction and array preparation.

This project was funded by the National Institute for Health Research (NIHR) Blood & Transplant Research Unit (BTRU) (NIHR-BTRU-2014-10074). The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care. The CIBMTR is supported primarily by Public Health Service U24CA076518 from the National Cancer Institute (NCI), the National Heart, Lung and Blood Institute (NHLBI) and the National Institute of Allergy and Infectious Diseases (NIAID); HHSH250201700006C from the Health Resources and Services Administration (HRSA); and N00014-20-1-2705 and N00014-20-1-2832 from the Office of Naval Research; Support is also provided by Be the Match Foundation, the Medical College of Wisconsin, and the National Marrow Donor Program.

Finally, the authors would like to sincerely thank all donors, patients and their families who contributed to CIBMTR cohort.

*   Received March 8, 2022.
*   Revision received March 8, 2022.
*   Accepted March 10, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/)

## References

1.  Duarte, R. F. et al. Indications for haematopoietic stem cell transplantation for haematological diseases, solid tumours and immune disorders: current practice in Europe, 2019. Bone Marrow Transplant 54, 1525–1552, doi:10.1038/s41409-019-0516-2 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41409-019-0516-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30953028&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

2.  McDonald-Hyman, C., Turka, L. A. & Blazar, B. R. Advances and challenges in immunotherapy for solid organ and hematopoietic stem cell transplantation. Sci Transl Med 7, 280rv282, doi:10.1126/scitranslmed.aaa6853 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/scitranslmed.aaa6853&link_type=DOI) 

3.  Al-Kadhimi, Z. et al. High incidence of severe acute graft-versus-host disease with tacrolimus and mycophenolate mofetil in a large cohort of related and unrelated allogeneic transplantation patients. Biol Blood Marrow Transplant 20, 979–985, doi:10.1016/j.bbmt.2014.03.016 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbmt.2014.03.016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24709007&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000337850500011&link_type=ISI) 

4.  Horvath, S. DNA methylation age of human tissues and cell types. Genome Biol 14, R115, doi:10.1186/gb-2013-14-10-r115 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2013-14-10-r115&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24138928&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

5.  Hannum, G. et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell 49, 359–367, doi:10.1016/j.molcel.2012.10.016 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2012.10.016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23177740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314379400018&link_type=ISI) 

6.  Yousefi, P. et al. Sex differences in DNA methylation assessed by 450 K BeadChip in newborns. BMC Genomics 16, 911, doi:10.1186/s12864-015-2034-y (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-015-2034-y&link_type=DOI) 

7.  Birdwell, C. E. et al. Genome-wide DNA methylation as an epigenetic consequence of Epstein-Barr virus infection of immortalized keratinocytes. J Virol 88, 11442–11458, doi:10.1128/JVI.00972-14 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjExOiI4OC8xOS8xMTQ0MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAzLzEwLzIwMjIuMDMuMDguMjIyNzIwNzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

8.  Capper, D. et al. DNA methylation-based classification of central nervous system tumours. Nature 555, 469–474, doi:10.1038/nature26000 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature26000&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29539639&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

9.  Koelsche, C. et al. Sarcoma classification by DNA methylation profiling. Nat Commun 12, 498, doi:10.1038/s41467-020-20603-4 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-20603-4&link_type=DOI) 

10. Maros, M. E. et al. Machine learning workflows to estimate class probabilities for precision cancer diagnostics on DNA methylation microarray data. Nat Protoc 15, 479–512, doi:10.1038/s41596-019-0251-6 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41596-019-0251-6&link_type=DOI) 

11. Paul, D. S. et al. A donor-specific epigenetic classifier for acute graft-versus-host disease severity in hematopoietic stem cell transplantation. Genome Med 7, 128, doi:10.1186/s13073-015-0246-z (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13073-015-0246-z&link_type=DOI) 

12. Tsai, P. C. & Bell, J. T. Power and sample size estimation for epigenome-wide association scans to detect differential DNA methylation. Int J Epidemiol 44, 1429–1441, doi:10.1093/ije/dyv041 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyv041&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25972603&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

13. Pidsley, R. et al. A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics 14, 293, doi:10.1186/1471-2164-14-293 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2164-14-293&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23631413&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

14. Zhou, W., Laird, P. W. & Shen, H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res 45, e22, doi:10.1093/nar/gkw967 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw967&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

15. Teschendorff, A. E. et al. An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS One 4, e8274, doi:10.1371/journal.pone.0008274 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0008274&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20019873&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

16. Johnson, W. E., Li, C. & Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8, 118–127, doi:10.1093/biostatistics/kxj037 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biostatistics/kxj037&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16632515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242715400008&link_type=ISI) 

17. Houseman, E. A. et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 13, 86, doi: 10.1186/1471-2105-13-86 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-13-86&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22568884&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

18. Morris, T. J. et al. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics 30, 428–430, doi:10.1093/bioinformatics/btt684 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt684&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24336642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

19. Tian, Y. et al. ChAMP: updated methylation analysis pipeline for Illumina BeadChips. Bioinformatics 33, 3982–3984, doi:10.1093/bioinformatics/btx513 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btx513&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

20. Reinius, L. E. et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One 7, e41361, doi:10.1371/journal.pone.0041361 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0041361&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22848472&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

21. Smyth, G. K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3, Article3, doi:10.2202/1544-6115.1027 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2202/1544-6115.1027&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16646809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 

22. Breiman, L. Random Forests. Machine Learning 45, 5–32, doi:10.1023/A:1010933404324 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1010933404324&link_type=DOI) 

23. Lee, C. et al. Prediction of absolute risk of acute graft-versus-host disease following hematopoietic cell transplantation. PLoS One 13, e0190610, doi:10.1371/journal.pone.0190610 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0190610&link_type=DOI) 

24. Wang, W. et al. Impact of Epigenomic Hypermethylation at TP53 on Allogeneic Hematopoietic Cell Transplantation Outcomes for Myelodysplastic Syndromes. Transplant Cell Ther, doi:10.1016/j.jtct.2021.04.027 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtct.2021.04.027&link_type=DOI) 

25. Ji, H. et al. Comprehensive methylome map of lineage commitment from haematopoietic progenitors. Nature 467, 338–342, doi: 10.1038/nature09367 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09367&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20720541&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281824900041&link_type=ISI) 

26. Ransohoff, D. F. Rules of evidence for cancer molecular-marker discovery and validation. Nat Rev Cancer 4, 309–314, doi: 10.1038/nrc1322 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrc1322&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15057290&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F10%2F2022.03.08.22272071.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220558900015&link_type=ISI)