A new approach to epigenome-wide discovery of non-invasive methylation biomarkers for colorectal cancer screening in circulating cell-free DNA using pooled samples

Clinical Epigenetics, Apr 2018

María Gallardo-Gómez, Sebastian Moran, María Páez de la Cadena, Vicenta Soledad Martínez-Zorzano, Francisco Javier Rodríguez-Berrocal, Mar Rodríguez-Girondo, Manel Esteller, Joaquín Cubiella, Luis Bujanda, Antoni Castells, et al.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1186%2Fs13148-018-0487-y.pdf

A new approach to epigenome-wide discovery of non-invasive methylation biomarkers for colorectal cancer screening in circulating cell-free DNA using pooled samples

Gallardo-Gómez et al. Clinical Epigenetics A new approach to epigenome-wide discovery of non-invasive methylation biomarkers for colorectal cancer screening in circulating cell-free DNA using pooled samples María Gallardo-Gómez 0 3 Sebastian Moran 2 María Páez de la Cadena 0 3 Vicenta Soledad Martínez-Zorzano 0 3 Francisco Javier Rodríguez-Berrocal 0 3 Mar Rodríguez-Girondo 1 7 Manel Esteller 2 Joaquín Cubiella 6 Luis Bujanda 5 Antoni Castells 4 Francesc Balaguer 4 Rodrigo Jover 8 Loretta De Chiara 0 3 0 Department of Biochemistry, Genetics and Immunology, Centro Singular de Investigación de Galicia (CINBIO), University of Vigo , Campus As Lagoas-Marcosende s/n, 36310 Vigo , Spain 1 Department of Medical Statistics and Bioinformatics, Leiden University Medical Centre , Leiden , The Netherlands 2 Cancer Epigenetics and Biology Program (PEBC), Bellvitge Biomedical Research Institute (IDIBELL) , Barcelona , Spain 3 Department of Biochemistry, Genetics and Immunology, Centro Singular de Investigación de Galicia (CINBIO), University of Vigo , Campus As Lagoas-Marcosende s/n, 36310 Vigo , Spain 4 Gastroenterology Department, Hospital Clínic, IDIBAPS, CIBERehd, University of Barcelona , Barcelona , Spain 5 Department of Gastroenterology, Instituto Biodonostia, Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd), Universidad del País Vasco (UPV/EHU) , San Sebastián , Spain 6 Department of Gastroenterology, Complexo Hospitalario Universitario de Ourense, Instituto de Investigación Biomédica Galicia Sur, Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd) , Ourense , Spain 7 SiDOR Research Group and Centro de Investigaciones Biomédicas (CINBIO), Faculty of Economics and Business Administration, University of Vigo , Vigo , Spain 8 Department of Gastroenterology, Hospital General Universitario de Alicante , Alicante , Spain Background: Colorectal cancer is the fourth cause of cancer-related deaths worldwide, though detection at early stages associates with good prognosis. Thus, there is a clear demand for novel non-invasive tests for the early detection of colorectal cancer and premalignant advanced adenomas, to be used in population-wide screening programs. Aberrant DNA methylation detected in liquid biopsies, such as serum circulating cell-free DNA (cfDNA), is a promising source of non-invasive biomarkers. This study aimed to assess the feasibility of using cfDNA pooled samples to identify potential serum methylation biomarkers for the detection of advanced colorectal neoplasia (colorectal cancer or advanced adenomas) using microarray-based technology. Results: cfDNA was extracted from serum samples from 20 individuals with no colorectal findings, 20 patients with advanced adenomas, and 20 patients with colorectal cancer (stages I and II). Two pooled samples were prepared for each pathological group using equal amounts of cfDNA from 10 individuals, sex-, age-, and recruitment hospitalmatched. We measured the methylation levels of 866,836 CpG positions across the genome using the MethylationEPIC array. Pooled serum cfDNA methylation data meets the quality requirements. The proportion of detected CpG in all pools (> 99% with detection p value < 0.01) exceeded Illumina Infinium methylation data quality metrics of the number of sites detected. The differential methylation analysis revealed 1384 CpG sites (5% false discovery rate) with at least 10% difference in the methylation level between no colorectal findings controls and advanced neoplasia, the majority of which were hypomethylated. Unsupervised clustering showed that cfDNA methylation patterns can distinguish advanced neoplasia from healthy controls, as well as separate tumor tissue from healthy mucosa in an independent dataset. We also observed that advanced adenomas and stage I/II colorectal cancer methylation profiles, grouped as advanced neoplasia, are largely homogenous and clustered close together. (Continued on next page) - (Continued from previous page) Conclusions: This preliminary study shows the viability of microarray-based methylation biomarker discovery using pooled serum cfDNA samples as an alternative approach to tissue specimens. Our strategy sets an open door for deciphering new non-invasive biomarkers not only for colorectal cancer detection, but also for other types of cancers. Background Colorectal cancer (CRC) is the fourth leading cause of cancer-related deaths worldwide, accounting for over 1.4 million new cases in 2012 [ 1, 2 ]. While diagnosis at early stages associates with good prognosis and reduced mortality rates, the detection and removal of premalignant advanced adenomas (AA) results in the reduction of CRC incidence [3]. Since neoplastic transformation can last decades, there is a broad time window for implementing screening strategies for the detection of advanced neoplasia (AN: CRC or AA) [ 3, 4 ]. Approaches for CRC screening can be divided into two groups. Invasive procedures like colonoscopy allow the examination of the entire colon and the removal of lesions (polypectomy); however, limitations of this strategy include considerably low participation rates and high cost [ 5 ]. On the other hand, non-invasive methods like fecal immunological test (FIT) have the advantage of increased acceptance and adequate specificity, though sensitivity for colorectal tumors, especially of proximal location, and AA is moderate to low [ 6, 7 ]. Blood-based markers are capable of improving CRC screening adherence, and a large number of candidates have been reported for CRC diagnosis, reviewed in [8]. Currently, the most promising is the SEPT9 methylation assay, though its performance for the detection of early-stage tumors and AA needs to be improved [ 9 ]. Therefore, there is an imperative need of finding new non-invasive biomarkers for CRC screening. Nowadays, it is well-established that not only genetic alterations but also epigenetic modifications are involved in CRC development and progression [ 10 ]. The abnormal methylation occurring during colorectal neoplasia is characterized by promoter hypermethylation and transcriptional silencing of tumor suppressor or DNA repair genes [ 11, 12 ], coexisting with a global loss of methylation that leads to chromosomal and microsatellite instability and oncogene activation [ 11, 13 ]. Both promoter hypermethylation and global hypomethylation are hallmarks of early stages of colorectal carcinogenesis [ 10, 14 ]. Several methodologies are suitable for genomewide methylation biomarker discovery, including whole and reduced genome bisulfite sequencing and array-based genotyping technology [ 15, 16 ]. These epigenome-wide measurements allow a more successful identification of methylation alterations related to complex diseases compared to target studies. The main drawback is the large sample size needed, which increases project costs. DNA sample pooling strategies represent an affordable approach for biomarker discovery, resulting in reduced costs and increased amount of input DNA when small amounts are available. Additionally, it has been reported that pooled samples provide similar results to individual samples in both genome-wide [ 17, 18 ] and epigenome-wide [ 19 ] association studies. During the last years, it has been demonstrated that circulating cell-free DNA (cfDNA) present in liquid biopsies reflects methylation changes originated in tumor cells [ 20–22 ]. Given the stability of DNA methylation in body fluids [ 23–25 ], the discovery of cfDNA methylation markers using serum samples seems a very attractive alternative to direct the search of non-invasive biomarkers. Taking advantage of this fact, we hypothesized that an array-based epigenome-wide analysis using serum cfDNA as input could be a novel and affordable approach for the discovery of a methylation marker panel with greater diagnostic value, compared to other indirect strategies using tumor tissue and mucosa as input DNA. Therefore, in the present study, we aim to assess the feasibility of hybridizing pooled serum cfDNA to the MethylationEPIC array to detect differentially methylated patterns between patients with advanced neoplasia and individuals with no colorectal findings. Methods Study population and serum samples Individuals were recruited from the following Spanish Hospitals: Hospital Donostia (San Sebastián), Complexo Hospitalario Universitario de Ourense (Ourense), Hospital Clínic de Barcelona (Barcelona), and Hospital General Universitario de Alicante (Alicante). Patients’ characteristics are described elsewhere [ 26 ]. We carried out a stratified random sampling using colorectal finding and gender as stratifying variables. Moreover, age was restricted to 50–75 years and strata were matched by recruitment hospital and age. We selected from this multicenter cohort 20 individuals with no colorectal findings (NCF), 20 individuals with AA (adenomas ≥ 10 mm, with villous component or high-grade dysplasia), and 20 CRC cases (7 stage I and 13 stage II, according to the AJCC staging system [ 27 ]). Individuals were classified according to the most advanced lesion after colonoscopy. Lesions were considered “proximal” when located only proximal to the splenic flexure of the colon and “distal” when lesions were found only in the distal colon or in both distal and proximal colons. Advanced neoplasia (AN) was defined as AA or CRC. Blood samples were obtained the same day of the colonoscopy, immediately prior to the procedure. Blood samples were coagulated and subsequently centrifuged according to the manufacturer’s instruction for serum collection. Serum samples were stored at − 20 °C until used. DNA extraction and sample pooling We extracted cfDNA from 0.5–1.5 mL of serum using a phenol-chloroform protocol as described by Clemens et al. [ 28 ], with minor modifications, and resuspended in 20 μL sterile water. DNA concentration was determined for each individual sample using the Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific, MA, USA), a fluorimetric assay specific for double-stranded DNA that gives an accurate measurement of DNA concentration. All cfDNA samples were stored at − 20 °C. Two independent pooled samples were constructed for each pathological group (NCF, AA, and CRC) using equal amounts of cfDNA from 10 individuals per pool. The factors considered to match between pools were gender, age, and recruitment hospital. Table 1 shows epidemiologic and clinical data of each individual included in pool A and B (NCF), pool C and D (AA), and pool E and F (CRC). Since the preparation of pooled samples is a critical step that requires high accuracy, cfDNA from each individual included in a pool was thawed, tempered, and re-quantified using the Qubit assay. As reported by previous DNA pooling protocols [ 18, 29 ], in order to avoid inaccuracies derived from pipetting small volumes, we decided to dilute by a factor of two samples with more than 10 ng/μL of DNA. Diluted DNA was measured again. Once the actual concentration of all the individual samples of a pool was available, we determined the sample containing the limiting ng of cfDNA (based on measured concentration and volume). Based on this limiting nanogram, we calculated for each of the nine remaining samples the volume containing the same nanogram of cfDNA as the limiting sample. Finally, the pool was constructed by incorporating into the tube the corresponding volume of each of the 10 individual samples of the pool. The cfDNA mixture was allowed to stand for 1 h, and then the DNA concentration was quantified with the Qubit Assay to ensure that the final DNA concentration of the pool was as expected according to the theoretical calculation: ðlimiting ngÞ n ðtotal volume of the poolÞ where n is the number of individuals included in each pool (10). Pools were considered valid for the Infinium Methylation Assay protocol when the difference between expected and measured concentration (Qubit) was less than 5%. A graphical description of the pooling protocol is presented [see Additional file 1]. This protocol was followed for each of the pools included in the study. The six pooled cfDNA samples were stored at − 20 °C and were submitted to the Cancer Epigenetics and Biology Program (PEBC) facilities at the Bellvitge Biomedical Research Institute for processing. Epigenome-wide methylation measurements DNA methylation was analyzed with the Infinium MethylationEPIC BeadChip microarray (EPIC; Illumina Inc., CA, USA), that quantitatively measures the methylation levels of more than 850,000 CpG sites across the genome [ 30 ], located in promoter regions and gene bodies, and also in intergenic enhancer regions identified by the ENCODE [ 31, 32 ] and FANTOM5 [33] projects. Pooled samples were bisulfite treated in the same batch, and MethylationEPIC arrays were hybridized according to manufacturer’s instructions. Data preprocessing and differential methylation analysis Data quality control was assessed with the GenomeStudio V2011.1, based on the internal control probes present on the array. The preprocessing, normalization, and correction steps were conducted using the R environment (versions 3. 3.3 and 3.4.0) with Bioconductor packages. The pipeline was a sequence of R functions adapted from the minfi [ 34 ] and ChAMP [ 35 ] Bioconductor packages. Our dataset was normalized using the Functional Normalization implemented in the minfi package. This algorithm does not rely on any biological assumption and therefore is suitable for cases where global changes in the methylation levels are expected, such as in cancer-normal comparisons [ 36 ]. Detection p values were computed with the minfi package, and mean detection p values were examined across all samples in order to identify any failed sample. Probes with a detection p value > 0.01 in at least one sample were discarded. We filtered out probes containing a single nucleotide polymorphism (SNP) at the CpG interrogation site and at the single nucleotide extension for any minor allele frequency (MAF), and probes containing a SNP at the probe body for a MAF >5 %, Genderb Age Lesion descriptionc Lesion locationd F F F M M M M M Age 61 60 55 67 63 61 59 57 T2N0, MD T3N0, MD T3N0, NA T2N0M0, WD T2N0, MD T3N0, MD T3N0, WD T2N0, MD 68 61 54 54 67 67 65 62 54 72 68 65 63 54 71 66 65 61 58 70 67 65 61 59 71 71 64 64 54 70 67 65 59 59 72 66 62 60 51 72 65 10 mm, T, LGD 25 mm, V, LGD 10 mm, TV, LGD 12 mm, T, LGD 10 mm, T, LGD 10 mm, T, LGD 10 mm, TV, LGD 20 mm, T, LGD 15 mm, T, LGD 3 mm, TV, LGD 30 mm, T, LGD 20 mm, TV, LGD 30 mm, TV, LGD 10 mm, TV, LGD 10 mm, TV, LGD 8 mm, TV, LGD 10 mm, V, LGD 20 mm, V, LGD 20 mm, T, LGD 5 mm, TV, LGD T3N0, WD T3N0, MD T3N0, MD T4N0, MD T3N0, WD T2N0, MD T2N0, NA T3N0, MD T1N0M0, WD T3N0, MD T3N0, WD T3N0, WD F F F F F M M M M M Proximal Distal Distal Distal Distal Distal Distal Distal Proximal Proximal Proximal Distal Proximal Distal Distal Distal Proximal Distal Proximal Distal Distal Distal Distal Distal Distal Distal Distal Distal Distal Proximal Distal Distal Age 63 62 59 54 54 71 68 66 60 53 because differential methylation levels can be confounded with actual polymorphisms in the DNA sequence. According to the list provided by Pidsley et al. [ 37 ], cross-reactive probes were removed. Probes targeting X and Y chromosomes were also discarded. In accordance with Du et al. [ 38 ], methylation levels were expressed as beta and M values. Beta values were used for visualization and intuitive interpretation of the results, and M values were used for the differential methylation analysis. Prior to differential methylation analysis, data was checked for batch effects across all array runs using the combat method implemented in the ChAMP package. Differentially methylated positions (DMP) between NCF and AN (AA or CRC) were detected with the dmpFinder function from the minfi package, which uses an F-test for categorical phenotype comparisons at a probe level. p values for each probe were corrected for multiple testing using the Benjamini-Hochberg procedure, with a false discovery rate (FDR) of 5% to determine significant DMPs. In silico evaluation of differential methylation We applied unsupervised clustering approaches to evaluate the differentially methylated patterns between AN and NCF pools in an independent dataset. The publicly available dataset GSE48684 that includes the methylation data of 64 colorectal tumor biopsies (adenocarcinomas) and 41 healthy mucosa biopsies, measured with the Infinium HumanMethylation450 BeadChip array (450K) [ 14 ], was used as a test cohort. This independent evaluation was limited to the probes shared by 450K and EPIC arrays due to the absence of colorectal tumor and mucosa EPIC public datasets. final cfDNA concentration of the six pooled samples ranged from 135 to 250 ng. Results and discussion DNA pooling methodology To our knowledge, this is the first pooling-based study that analyzes the methylation patterns in cfDNA, aiming to assess the feasibility of liquid biopsy methylation biomarker discovery using microarray technology in a more affordable manner compared to individual samples. DNA sample pooling has been reported as an efficient tool for genome-wide and high-throughput association studies [ 17, 18 ]. More recently, its potential utility was highlighted in microarray-based epigenome-wide association studies (EWAS), as Gallego-Fabrega et al. reported high correlation of the methylation levels between pools and individual DNA samples using the Infinium HumanMethylation450 BeadChip [19]. Taking into account the limitation that only mean methylation levels can be obtained from pooled samples, pooling strategy is an accurate and affordable alternative that can significantly reduce costs in large EWAS. DNA pooling is also an efficient alternative when small amounts of DNA are available and when working with precious samples. For sample pooling, accurate construction is critical, and each DNA sample must be equally represented in the pool. To guarantee the most precise pool construction, we first tested two different pooling strategies: diluting all samples to a common concentration and then mixing equal volumes in a tube as in previous works [ 19, 29 ] or directly adding the same nanogram of DNA (calculated corresponding volume from each sample) into the tube. Once test pools were constructed and DNA concentration measured, we checked for discrepancies between the actual and the expected concentration. Variations inferior to 5% of the expected pooled DNA concentration were found when using the second protocol; therefore, sample pooling was performed as described [see Additional file 1]. We prepared two pooled samples for each pathological group (two pools of individuals with NCF, two pools of AA patients, and two pools of CRC patients stages I and II). We included 10 individuals per pool to ensure an acceptable amount of DNA input for the microarray analyses and also to reduce population stratification and the presence of unobserved confounding variables. The categories considered to match between pools were gender, age (median 63.5, range 51–72 years), and recruitment hospital. The age range was selected based on the USPSTF guideline recommendation for CRC screening, targeting individuals from 50 to 75 years [ 39 ]. No statistically significant difference was found in the mean age between pools (ANOVA, p < 0.05). The Quality control of methylation data The methylation levels of 866,836 CpG positions across the genome in the six pooled samples were measured using the MethylationEPIC BeadChip. The quality control based on the internal control probes present on the array, which include bisulfite conversion efficiency, hybridization, extension, and staining, among others, indicates that pooled serum cfDNA methylation data meets the quality requirements. The QC report is presented [see Additional file 2] and shows that the signals observed are much higher than the background signal, coinciding with what is expected for high-quality DNA. In relation to CpG detection, all the pools showed more than 99% of CpG detected correctly (only 2811 probes presented a detection p value > 0.01 at least in one sample, and were discarded). The number of probes detected in each pool was 866,497; 866,021; 865,865; 865,501; 865,778; and 866,463 for pools A, B, C, D, E, and F, respectively. These results are indicative of a uniform amplification and hybridization in all the pooled samples. The proportion of CpG detection observed in our samples exceeded Illumina Infinium methylation data quality metrics of the number of sites detected (> 96% for genomic DNA and > 90% for FFPE samples). A probable explanation could be the pooling design, as measuring methylation of 10 individuals in the same assay would increase the representation of each CpG in the input DNA. Therefore, no samples were discarded due to QC issues. The distribution of methylation levels in pooled samples presented the expected bimodal distribution for both beta and M values, with the two peaks indicating fully methylated and unmethylated states characteristic of DNA methylation data (Fig. 1a). Then, we evaluated the distribution of beta values by type I and II probes separately. As observed in Fig. 1b, all the pools showed the distribution of type II probes shifted in relation to type I, as previously reported in the 450K and EPIC data [ 37, 40 ]. Differential methylation Once the technical quality of the six pooled sample data is verified, we performed a differential methylation analysis. In order to detect differentially methylated positions (DMPs), we compared the two NCF pools (colonoscopically confirmed controls) with the two AA pools together with the two CRC pools (considered as AN). The differential methylation analysis was performed on the 703,653 probes left after the filtering step (see the “Methods” section). We first assessed the global methylation in the NCF and AN groups, and in the AA and CRC groups. As shown in Fig. 2a, a lower content of global methylation was observed in AN, AA, and CRC compared to NCF. In addition, we found there is no difference in terms of global methylation between AA and CRC cfDNA pooled samples. Since the purpose of screening programs include the detection of early stage CRC together with the identification and removal of premalignant AA [ 4 ], we grouped AA and CRC in the single group AN for the analyses. We found a total of 5808 significant DMPs between the NCF and AN groups, identified with a FDR of 5% (Fig. 2b). Of these, 1384 presented at least 10% difference in the methylation level between NCF and AN (|Δbeta| > 0.1): 135 (9.75%) were found hypermethylated in AN, while 1249 (90.25%) appeared hypomethylated (Fig. 2c). The distribution of the DMPs identified according to their location relative to CpG islands (CGI) and promoter regions is represented in Fig. 2d. The majority of the differentially hypomethylated CpGs in AN are located in opensea (78.08%), outside gene promoters, within regions with no enrichment in CpG content followed by CGI-shore (10.23%), CGI-shelf (7.50%), and CGI (4.19%). DNA hypomethylation was the first aberrant methylation alteration described in several human cancers (reviewed in [ 41 ]). This global loss of genome-wide methylation was also described long time ago in both CRC and colorectal adenomas [ 42 ], indicating that global hypomethylation is characteristic of early stages of colorectal carcinogenesis [ 10, 14, 43 ]. Large hypomethylated blocks were also identified by Timp and colleagues (2014) using the 450K array. Among the six different tumor types analyzed, colon cancer tissue showed the highest proportions of hypomethylation in opensea, CGI-shelf, CGI-shore, and CGI [ 44 ]. In our work, using pooled serum samples, we found that more than 90% of the DMPs appeared hypomethylated in AN, agreeing with these previous reports, and suggesting that perhaps efforts should be centered on hypomethylated candidates to accomplish a greater discrimination capacity. Unsupervised clustering performed on DNA methylation values for the top 1384 DMPs identified is presented in Fig. 3a, b. These results highlight the differences between AN and NCF pooled samples and suggest that differential cfDNA methylation profiles obtained with pooled samples can discriminate AN from NCF controls. We further evaluated the DMPs identified in our pooled serum cfDNA samples with dataset GSE48684 consisting of tumor and mucosa tissue samples [ 14 ] as a test cohort, restricting the analysis to the 518 DMPs between AN and NCF with |Δbeta| > 0.1 targeted by probes shared by 450K and EPIC arrays. The unsupervised clustering shown in Fig. 4, performed on tumor and mucosa samples from GSE48684 based on our DMPs, reveals that the differential methylation patterns found between AN and NCF cfDNA can also separate tumor tissue from healthy mucosa samples. It is worth to mention that 24 healthy mucosa samples from GSE48684 were normal colon concurrent with CRC and were obtained from the normal-appearing resection margin of the colorectal tumor biopsy [ 14 ], as represented in Fig. 4b. This can be related to the fact that a subcluster of mucosa samples partially overlaps with CRC samples. Though this in silico verification is limited, a considerable degree of concordance can be deduced. It should also be mentioned that discrepancy in the frequencies of methylation alterations have been reported in tumor and cfDNA, showing the latter considerably lower frequencies [ 45 ]. Hence, arraybased strategies that rely on tissue samples for cfDNA methylation marker discovery have the inconvenience of resulting in decreased sensitivity of the selected candidate markers once tested in serum or plasma, limiting their utility as non-invasive tests [ 46, 47 ]. An alternative approach for biomarker discovery was accomplished by Heiss et al. that used whole blood. However, these authors indicate that the methylation signature identified in leukocyte DNA may not be specific for CRC, reflecting immune responses [48]. The exploratory nature of this study, with a reduced number of samples, limits further analyses, but offers a new affordable strategy for biomarker discovery, providing an alternative approach to tissue biopsy, reducing costs in microarray-based EWAS. This work should be followed by new studies that include a greater number of pooled serum cfDNA samples and a greater range of colorectal pathologies, allowing a more robust comparison between methylation profiles. Furthermore, differential methylation profiles must be validated in independent serum cfDNA individual samples, using quantitative realtime techniques, with the aim of finding a serum methylation panel for CRC diagnosis and screening. Conclusion As far as we are concerned, this proof-of-principle study is the first to evaluate pooled serum cfDNA profiling on an epigenome-wide scale for CRC biomarker discovery using the MethylationEPIC array. Our data, although preliminary, revealed that the whole epigenome is represented in pooled serum cfDNA samples and that differentially methylated cfDNA profiles can discriminate NCF controls from AN cases (AA or CRC). These results suggest that a pooling strategy using cfDNA may be a valuable source of novel non-invasive methylation biomarkers for CRC early detection and screening. Also, our approach can be translated to the search of biomarkers for other types of tumors, as an affordable alternative approach to tissue biopsy. Additional files Additional file 1: Graphical description of the protocol for DNA sample pooling. aExpected DNA concentration of the pool was calculated as follows: ðtotal ðvloimluimtinegonfgtÞhen poolÞ where n is the number of individuals included in each pool (10). (PDF 311 kb) Additional file 2: GenomeStudio software quality control report based on the internal control probes present on the array. (PDF 97 kb) Abbreviations AA: Advanced adenoma; AN: Advanced neoplasia; cfDNA: Circulating cell-free DNA; CRC: Colorectal cancer; DMP: Differentially methylated position; FDR: False discovery rate; NCF: No colorectal findings Acknowledgements We would like to thank Leticia Barcia for her support in daily laboratory work and also Dr. Jezabel Varadé for her advice and tips about the results presentation. Funding This work received funding from Plan Nacional I +D +I 2015-2018 (Acción Estratégica en Salud) Instituto de Salud Carlos III (Spain)-FEDER (PI15/02007), “Fundación Científica de la Asociación Española contra el Cáncer” (GCB13131592CAST), and support from Centro Singular de Investigación de Galicia (Consellería de Cultura, Educación e Ordenación Universitaria) (ED431G/02, Xunta de Galicia and FEDER-European Union). María Gallardo-Gómez is supported by a predoctoral fellowship from Ministerio de Educación, Cultura y Deporte (Spanish Government) (FPU15/02350). Availability of data and materials The EPIC data from all the pooled samples generated and analyzed during this study has been deposited in the NCBI Gene Expression Omnibus (GEO) (www.ncbi.nlm.nih.gov/geo) and are accessible through GEO Series accession number GSE110185. Authors’ contributions LD, VSMZ, and SM conceived and designed the study. LD, MP, FJRB, and ME supervised the study. JC, LB, AC, FB, and RJ clinical advise for the study design, collection, and management of clinical data. MGG, MRG, and MP contributed to the experimental design. MGG, LD, and SM contributed to the sample preparation and data acquisition. LD, MGG, VSMZ, and SM performed the analysis and interpretation of data. MGG, LD, FJRB, and MP prepared the manuscript. All authors critically reviewed and approved the final manuscript. Ethics approval and consent to participate All individuals provided written informed consent, and the study followed the ethical and clinical practices of the Spanish Government and the Helsinki Declaration, and was approved by the Galician Ethical Committee for Clinical Research. Competing interests The authors declare that they have no competing interests. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 1. Ferlay J , Soerjomataram I , Dikshit R , Eser S , Mathers C , Rebelo M , et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012 . Int J Cancer . 2015 ; 136 : E359 - 86 . 2. Arnold M , Sierra MS , Laversanne M , Soerjomataram I , Jemal A , Bray F . Global patterns and trends in colorectal cancer incidence and mortality . Gut . 2017 ; 66 : 683 - 91 . https://doi.org/10.1136/gutjnl-2015-310912. 3. Ng SC , Lau JYW , Chan FKL , Suen BY , Leung WK , Tse YK , et al. Increased risk of advanced neoplasms among asymptomatic siblings of patients with colorectal cancer . Gastroenterology . 2013 ; 144 : 544 - 50 . 4. Brenner H , Stock C , Hoffmeister M. Colorectal cancer screening: the time to act is now . BMC Med . 2015 ; 13 : 262 . https://doi.org/10.1186/ s12916-015-0498-x. 5. Salas Trejo D , Portillo Villares I , Espinàs Piñol JA , Ibáñez Cabanell J , Vanaclocha Espí M , Pérez Riquelme F , et al. Implementation of colorectal cancer screening in Spain . Eur J Cancer Prev . 2017 ; 26 : 17 - 26 . https://doi.org/10.1097/CEJ.0000000000000232. 6. Castro I , Cubiella J , Rivera C , González-Mao C , Vega P , Soto S , et al. Fecal immunochemical test accuracy in familial risk colorectal cancer screening . Int J Cancer . 2014 ; 134 : 367 - 75 . https://doi.org/10.1002/ijc.28353. 7. Kim NH , Yang H-J , Park S-K , Park JH , Park DI , Sohn CI , et al. Does low threshold value use improve proximal neoplasia detection by fecal immunochemical test? Dig Dis Sci . 2016 ; 61 : 2685 - 93 . https://doi.org/10.1007/s10620-016-4169-3. 8. Shah R , Jones E , Vidart V , Kuppen PJK , Conti JA , Francis NK . Biomarkers for early detection of colorectal cancer and polyps: systematic review . Cancer Epidemiol Biomark Prev . 2014 ; 23 : 1712 - 28 . https://doi.org/10.1158/ 1055 - 9965 .EPI- 14 -0412. 9. Song L , Li Y. Progress on the clinical application of the SEPT9 gene methylation assay in the past 5 years . Biomark Med . 2017 ; 11 : 415 - 8 . https://doi.org/10.2217/bmm-2017-0091. 10. Lao VV , Grady WM . Epigenetics and colorectal cancer . Nat Rev Gastroenterol Hepatol . 2012 ; 8 : 686 - 700 . 11. Bariol C , Suter C , Cheong K , Ku SL , Meagher A , Hawkins N , et al. The relationship between hypomethylation and CpG island methylation in colorectal neoplasia . Am J Pathol . 2003 ; 162 : 1361 - 71 . 12. Øster B , Thorsen K , Lamy P , Wojdacz TK , Hansen LL , Birkenkamp-Demtröder K , et al. Identification and validation of highly frequent CpG island hypermethylation in colorectal adenomas and carcinomas . Int J Cancer . 2011 ; 129 : 2855 - 66 . 13. Beggs AD , Jones A , El-Bahwary M , Abulafi M , Hodgson SV , Tomlinson IPM . Whole-genome methylation analysis of benign and malignant colorectal tumours . J Pathol . 2013 ; 229 : 697 - 704 . 14. Luo Y , Wong CJ , Kaz AM , Dzieciatkowski S , Carter KT , Morris SM , et al. Differences in DNA methylation signatures reveal multiple pathways of progression from adenoma to colorectal cancer . Gastroenterology . 2014 ; 147 : 418 - 29 . https://doi.org/10.1053/j.gastro. 2014 . 04 .039. 15. Laird PW . Principles and challenges of genome-wide DNA methylation analysis . Nat Rev Genet . 2010 ; 11 : 191 - 203 . https://doi.org/10.1038/nrg2732. 16. Fan S , Chi W. Methods for genome-wide DNA methylation analysis in human cancer . Brief Funct Genomics . 2016 ; 15 :elw010. https://doi.org/10.1093/bfgp/elw010. 17. Norton N , Williams NM , O'Donovan MC , Owen MJ . DNA pooling as a tool for large-scale association studies in complex traits . Ann Med . 2004 ; 36 : 146 - 52 . 18. Pearson JV , Huentelman MJ , Halperin RF , Tembe WD , Melquist S , Homer N , et al. Identification of the genetic basis for complex disorders by use of pooling-based genomewide single-nucleotide-polymorphism association studies . Am J Hum Genet . 2007 ; 80 : 126 - 39 . https://doi.org/10.1086/510686. 19. Gallego-Fabrega C , Carrera C , Muiño E , Montaner J , Krupinski J , Fernandez-Cadenas I. DNA methylation levels are highly correlated between pooled samples and averaged values when analysed using the Infinium HumanMethylation450 BeadChip array . Clin Epigenetics . 2015 ; 7 : 78 . https://doi.org/10.1186/s13148-015-0097-x. 20. Schwarzenbach H , Hoon DSB , Pantel K. Cell-free nucleic acids as biomarkers in cancer patients . Nat Rev Cancer . 2011 ; 11 : 426 - 37 . https://doi.org/10.1038/nrc3066. 21. Krishnamurthy N , Spencer E , Torkamani A , Nicholson L . Liquid biopsies for cancer: coming to a patient near you . J Clin Med . 2017 ; 6:3 . https://doi.org/ 10.3390/jcm6010003. 22. Zhai R , Zhao Y , Su L , Cassidy L , Liu G , Christiani DC . Genome-wide DNA methylation profiling of cell-free serum DNA in esophageal adenocarcinoma and Barrett esophagus . Neoplasia . 2012 ; 14 : 29 - 33 . https:// doi.org/10.1593/neo.111626. 23. Bosch LJW , Mongera S , sive Droste JST , Oort FA , van Turenhout ST , Penning MT , et al. Analytical sensitivity and stability of DNA methylation testing in stool samples for colorectal cancer detection . Cell Oncol . 2012 ; 35 : 309 - 15 . https://doi.org/10.1007/s13402-012-0092-6. 24. Forat S , Huettel B , Reinhardt R , Fimmers R , Haidl G , Denschlag D , et al. Methylation markers for the identification of body fluids and tissues from forensic trace evidence . PLoS One . 2016 ; 11 :e0147973. https://doi.org/10.1371/journal.pone. 0147973 . 25. Galanopoulos M , Tsoukalas N , Papanikolaou IS , Tolia M , Gazouli M , Mantzaris GJ . Abnormal DNA methylation as a cell-free circulating DNA biomarker for colorectal cancer detection: a review of literature . World J Gastrointest Oncol . 2017 ; 9 : 142 - 52 . https://doi.org/10.4251/wjgo.v9. i4 . 142 . 26. Quintero E , Castells A , Bujanda L , Cubiella J , Salas D , Lanas Á , et al. Colonoscopy versus fecal immunochemical testing in colorectal-cancer screening . N Engl J Med . 2012 ; 366 : 697 - 706 . https://doi.org/10.1056/NEJMoa1108895. 27. Edge SB , Compton CC . The American Joint Committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM . Ann Surg Oncol. 2010 ; 17 : 1471 - 4 . https://doi.org/10.1245/s10434-010-0985-4. 28. Clemens H , Markus S , Martin M , Roland G , Richard G. A modified phenolchloroform extraction method for isolating circulating cell free DNA of tumor patients . J Nucleic Acids Investig . 2013 ; 4 : 1 . 29. Sham P , Bader JS , Craig I , O'Donovan M , Owen M. DNA pooling: a tool for large-scale association studies . Nat Rev Genet . 2002 ; 3 : 862 - 71 . 30. Moran S , Arribas C , Esteller M. Validation of a DNA methylation microarray for 850,000 CpG sites of the human genome enriched in enhancer sequences . Epigenomics . 2016 ; 8 : 389 - 99 . https://doi.org/10.2217/epi.15.114. 31. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome . Nature . 2012 ; 489 : 57 - 74 . https://doi.org/10.1038/nature11247. 32. Siggens L , Ekwall K . Epigenetics, chromatin and genome organization: recent advances from the ENCODE project . J Intern Med . 2014 ; 276 : 201 - 14 . 33. Lizio M , Harshbarger J , Shimoji H , Severin J , Kasukawa T , Sahin S , et al. Gateways to the FANTOM5 promoter level mammalian expression atlas . Genome Biol . 2015 ; 16 : 22 . https://doi.org/10.1186/s13059-014-0560-6. 34. Aryee MJ , Jaffe AE , Corrada-Bravo H , Ladd-Acosta C , Feinberg AP , Hansen KD , et al. Minfi: a flexible and comprehensive bioconductor package for the analysis of Infinium DNA methylation microarrays . Bioinformatics . 2014 ; 30 : 1363 - 9 . 35. Morris TJ , Butcher LM , Feber A , Teschendorff AE , Chakravarthy AR , Wojdacz TK , et al. ChAMP: 450k chip analysis methylation pipeline . Bioinformatics . 2014 ; 30 : 428 - 30 . 36. Fortin J-P , Labbe A , Lemire M , Zanke BW , Hudson TJ , Fertig EJ , et al. Functional normalization of 450k methylation array data improves replication in large cancer studies . Genome Biol . 2014 ; 15 : 503 . https://doi.org/10.1186/s13059-014-0503-2. 37. Pidsley R , Zotenko E , Peters TJ , Lawrence MG , Risbridger GP , Molloy P , et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling . Genome Biol . 2016 ; 17 : 208 . https://doi.org/10.1186/s13059-016-1066-1. 38. Du P , Zhang X , Huang C-C , Jafari N , Kibbe WA , Hou L , et al. Comparison of beta-value and M-value methods for quantifying methylation levels by microarray analysis . BMC Bioinformatics . 2010 ; 11 : 587 . https://doi.org/10. 1186/ 1471 -2105-11-587. 39. Bibbins-Domingo K , Grossman DC , Curry SJ , Davidson KW , Epling JW , García FAR , et al. Screening for colorectal cancer . JAMA . 2016 ; 315 : 2564 - 75 . https://doi.org/10.1001/jama. 2016 . 5989 . 40. Dedeurwaerder S , Defrance M , Calonne E , Denis H , Sotiriou C , Fuks F. Evaluation of the Infinium methylation 450K technology . Epigenomics . 2011 ; 3 : 771 - 84 . 41. Ehrlich M. DNA hypomethylation in cancer cells . Epigenomics . 2009 ; 1 : 239 - 59 . https://doi.org/10.2217/epi.09.33. 42. Goelz SE , Vogelstein B , Hamilton SR , Feinberg AP . Hypomethylation of DNA from benign and malignant human colon neoplasms . Science . 1985 ; 228 : 187 - 90 . https://doi.org/10.1126/SCIENCE.2579435. 43. Okugawa Y , Grady WM , Goel A . Epigenetic alterations in colorectal cancer: emerging biomarkers . Gastroenterology . 2015 ; 149 : 1204 - 1225 . e12 . https://doi.org/10.1053/j.gastro. 2015 . 07 .011. 44. Timp W , Bravo HC , McDonald OG , Goggins M , Umbricht C , Zeiger M , et al. Large hypomethylated blocks as a universal defining epigenetic alteration in human solid tumors . Genome Med . 2014 ;6 https://doi.org/10.1186/s13073-014-0061-y. 45. Jung K , Fleischhacker M , Rabien A . Cell-free DNA in the blood as a solid tumor biomarker-a critical appraisal of the literature . Clin Chim Acta . 2010 ; 411 : 1611 - 24 . https://doi.org/10.1016/j.cca. 2010 . 07 .032. 46. Roperch JP , Incitti R , Forbin S , Bard F , Mansour H , Mesli F , et al. Aberrant methylation of NPY, PENK, and WIF1 as a promising marker for blood-based diagnosis of colorectal cancer . BMC Cancer . 2013 ; 13 : 566 . 47. Barault L , Amatu A , Siravegna G , Ponzetti A , Moran S , Cassingena A , et al. Discovery of methylated circulating DNA biomarkers for comprehensive non-invasive monitoring of treatment response in metastatic colorectal cancer . Gut . 2017 ; 0 : 1 - 11 . https://doi.org/10.1136/gutjnl-2016-313372. 48. Heiss JA , Brenner H . Epigenome-wide discovery and evaluation of leukocyte DNA methylation markers for the detection of colorectal cancer in a screening setting . Clin Epigenetics . 2017 ; 9 : 24 . https://doi.org/10.1186/s13148-017-0322-x.


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1186%2Fs13148-018-0487-y.pdf

María Gallardo-Gómez, Sebastian Moran, María Páez de la Cadena, Vicenta Soledad Martínez-Zorzano, Francisco Javier Rodríguez-Berrocal, Mar Rodríguez-Girondo, Manel Esteller, Joaquín Cubiella, Luis Bujanda, Antoni Castells, Francesc Balaguer, Rodrigo Jover, Loretta De Chiara. A new approach to epigenome-wide discovery of non-invasive methylation biomarkers for colorectal cancer screening in circulating cell-free DNA using pooled samples, Clinical Epigenetics, 2018, 53, DOI: 10.1186/s13148-018-0487-y