Skip to main content

Evaluation of a panel of tumor-specific differentially-methylated DNA regions in IRF4, IKZF1 and BCAT1 for blood-based detection of colorectal cancer

Abstract

Background

Differentially-methylated regions (DMRs) are characteristic of colorectal cancer (CRC) and some occur more frequently than common mutations. This study aimed to evaluate the clinical utility of assaying circulating cell-free DNA for methylation in BCAT1, IKZF1 and IRF4 for detection of CRC.

Methods

A multiplexed real-time PCR assay targeting DMRs in each of the three genes was developed. Assay accuracy was explored in plasma specimens banked from observational cross-sectional trials or from volunteers scheduled for colonoscopy or prior to CRC surgery.

Results

1620 specimens were suitable for study inclusion including 184 and 616 cases with CRC and adenomas, respectively, and 820 cases without neoplasia (overall median age, 63.0 years; 56% males). Combining the PCR signals for all targeted DMRs returned the best sensitivity for CRC (136/184, 73.9%, 95% CI 67.1–79.7), advanced adenomas (53/337, 15.7%, 95% CI 12.0–20.1) and high-grade dysplastic (HGD) adenomas (9/35, 25.7%, 95% CI 14.0–42.3) with a 90.1%, specificity for neoplasia (739/820, 95% CI 87.9–92.0, p < 0.01). Detection of methylation in all three genes were more likely in CRC cases than those without it (OR 28.5, 95% CI 7.3–121.2, p < 0.0001). Of the 81 positive cases without neoplasia, 62 (76.5%) were positive by a single PCR replicate only and predominantly due to detection of methylated BCAT1 (53.2%). Single replicate positivity was significantly higher than that in CRC (26/136, 19.1%, p < 0.0001), and single BCAT1 replicate positivity was more likely in cases without neoplasia than in CRC (OR 17.7, 95% CI 6.6–43.3, p < 0.0001). When a positive result was limited to those with ≥ 1 PCR replicate positive for either IKZF1 or IRF4, or at least two replicates positive for BCAT1, the multi-panel test maintained a high sensitivity for CRC (131/184, 71.2%, 95% CI 64.3–77.3) and HGD adenomas (8/35, 22.9%, 95% CI 11.8–39.3, p = 0.029) but improved specificity significantly (772/820, 94.1%, 95% CI 92.3–95.6, p < 0.0001 vs. any PCR replicate positive).

Conclusion

The multi-panel methylation assay differentiates cases with CRC from those without it and does so with high specificity when criteria for BCAT1 detection are applied. The marker panel is flexible and studies in those at average risk for CRC are now warranted to determine which panel configuration best suits screening goals.

Trial registration: ACTRN12611000318987. Registered 25 March 2011, https://www.anzctr.org.au/ ACTRN12611000318987.

Introduction

Screening programs for colorectal cancer (CRC) are near universal in developed countries, but suboptimal participation rates are commonly reported, especially for stool-based screening tests [1]. A blood test might overcome the behavioral barriers observed with stool-based screening tests [2].

Circulating tumor DNA (ctDNA) is a promising biomarker of cancer, including CRC [3, 4], but reliable detection of ctDNA is subject to the frequency of the targeted tumor-specific sequence(s), whether or not the tumor DNA enters into the circulation, and to its fragmented state in circulation.

Successful implementation of a ctDNA-based test in a screening setting is technically challenging as ctDNA comprises as little as 0.01% of the total amount of circulating cell-free DNA (ccfDNA) [3]. Hence, there is a growing interest in using multiple markers for ctDNA detection to increase sensitivity for early stage neoplasia which is the ideal target for screening.

Virtually all colorectal tumors have thousands of abnormally methylated DNA regions that seem likely to be independent of various molecular sub-types of CRC [4, 5]. Many of these epigenetic changes occur more frequently and earlier in tumorigenesis than most mutations [6]. Therefore, assaying ccfDNA for CRC-specific hypermethylation in multiple regions might be a better universal identifier of CRC-derived DNA in circulation than mutations. Further, methylation-based ctDNA tests are simple in construction and not confounded by the need to cover multiple and often large regions for possible mutations.

The process of defining and validating panels of biomarkers is complex and links biomarker discovery and assay configuration with clinical validation and adaptation of a platform suitable for a clinical setting [7]. Following marker discovery and initial validation, biomarker tests will often benefit from methodological refinement to optimise assay sensitivity and specificity. We have previously reported on multiple differentially-methylated regions (DMRs) which are methylated with high frequency in CRC and the identification of DMRs residing in four genes (BCAT1, IKZF1, IRF4 and GRASP) which exhibited low to no methylation in ccfDNA isolated from healthy subjects [8, 9]. An initial epigenetic ctDNA test targeting single DMRs in BCAT1 and IKZF1 resulted in a 62–64% sensitivity for CRC with a 92–94% specificity [10, 11]. To explore whether increasing the number of DMRs improves sensitivity for early stage CRC, we redesigned the methylation-specific ctDNA qPCR assay to include an additional DMR target in IRF4 and changed the IKZF1 PCR assay component to detect the targeted DMR on both strands [11, 12].

This study evaluated performance of the multi-panel assay for presence of targeted methylated regions in biobanked specimens obtained from a cohort where clinical phenotype was colonoscopy-confirmed by exploring the ways in which the targeted regions methylated with high frequency in colorectal neoplastic tissues might be utilized for blood-based detection of CRC.

Materials and methods

Study overview

This study was performed retrospectively using stored plasma samples that were collected from cases scheduled for colonoscopy (indicated by a wide range of clinical indications including screening and symptoms) or prior to colonic surgery. Circulating cell-free DNA (ccfDNA) was isolated from plasma, bisulphite converted and the resulting DNA was assayed in triplicate using a real-time multiplexed PCR for detection of DNA methylation in BCAT1, IKZF1 and IRF4. The qPCR assay also targeted a region in ACTB as a quality control for bisulphite converted DNA. True- and false-positive results, using findings at colonoscopy as the diagnostic standard, were determined for various combinations of the three genes.

The biobanked samples

Specimens were from plasma biobanks established from observational, predominantly prospective, cross-sectional trials undertaken at Flinders Medical Centre, Bedford Park, South Australia; Amsterdam University Medical Centers at Amsterdam Medical Center and Flevo Hospital, Almere, The Netherlands; and Hvidovre Hospital, Hvidovre, Denmark.

The Danish plasma biobank received specimens from volunteers participating in the Danish National CRC screening program. Diagnostic information was available for those undergoing colonoscopy following a fecal immunochemical test (FIT) positive result. The Dutch and Australian specimens were collected either prior to colonoscopy (for standard clinical indications including positive FIT, symptoms, surveillance due to family/personal history of neoplasia or for inflammatory bowel disease) or from cases shown at colonoscopy to have CRC and who had not yet received treatment.

An additional biobank of plasma specimens sourced through Proteogenex (CA, USA) was also included. These specimens were collected from volunteers 3–10 days after diagnostic colonoscopy (polyps were not removed) undertaken in Moscow, Russia.

The clinical trials were approved by the Institutional Review Boards of the respective sites and written informed consent was obtained from all cases for sample collection, storage and testing for research purposes. No study-wide control of colonoscopy or pathology procedures was undertaken as these specimens were collected for biomarker evaluation studies aimed to assess biomarker performance relative to outcomes determined in usual clinical practice. All venous blood was collected in either K2- or K3 EDTA tubes and processed to plasma using a 2-spin centrifugation approach (1500–3000 g for 10 min, 4–25 °C, lowest deceleration setting).

Clinical classification and specimen selection

Clinical phenotype was determined using clinicopathological findings by experts at each site, with main outcomes categorized as CRC, adenoma or no neoplasia. CRC was further subcategorized into stages according to the AJCC 7th Edition [13]. Advanced adenoma was defined as adenoma with any of the following characteristics: (a) ≥ 10 mm in size, (b) villous histology (> 20% villous component), (c) high-grade dysplasia (HGD) and/or (d) the presence of ≥ 3 tubular adenomas (< 10 mm and with low-grade dysplasia (LGD). The presence of multiple adenomas was included with the advanced adenoma classification as previous studies have shown that these lesions are associated with an increased risk for future advanced neoplasia [14, 15]. Non-advanced adenoma refers to those not meeting the characteristics of an advanced adenoma. Adenomas were also classified separately into dysplasia status. Stage 0 CRC, where there was severe cellular atypia or marked architectural distortion but no evidence of invasion, was included in advanced adenoma (as HGD). The most advanced neoplasm was used as the principal diagnosis when multiple colorectal pathologies were present.

Specimens were selected for inclusion in the analysis on the basis of sufficient plasma available for testing (3.9–4.5 mL) provided that complete clinical and demographic data were available. Cases excluded were those with known or suspected cancer of another organ at the time of collection, familial adenomatous polyposis or hereditary non-polyposis CRC syndrome (Lynch syndrome) or incomplete diagnostic information.

Detection of methylated DNA in plasma

All frozen plasma samples were couriered to Clinical Genomics Technologies for storage and subsequent testing (Sydney, NSW, Australia). ccfDNA was isolated from plasma using the QS DSP Circulating DNA Kit (Qiagen) on a QIASymphony SP instrument as per manufacturer’s instruction (Qiagen) and bisulphite converted using the EpiTect Fast 96 DNA Bisulfite Conversion kit (Qiagen) on a QIACube HT instrument as previously described [11]. The resulting purified bisulphite-converted DNA (~ 45μL) was assayed as triplicates of 12μL in a total PCR volume of 30μL on a Light Cycler 480 II (Roche Diagnostics, IN, USA) as previously described [11] (see Additional file 1: Fig. S1 for further details). Cycle threshold (Ct) values were calculated using the second derivative maximum algorithm provided with the LC480 software. The ACTB assay component was used for estimation of ccfDNA yield as well as a quality control parameter. Samples with mean ACTB Ct values ≥ 36.6 were not accepted for analysis unless positive for one or more of the methylation targets. Plasma specimens were processing by staff who were blinded to the associated clinical and demographic data.

Statistical methods

Detection rates for each DMR were determined for each clinical phenotype and assay positivity was based on using various combinations of the DMRs. Sensitivity for a colorectal neoplastic condition was estimated from the positivity rate in the presence of that clinical phenotype. Specificity for neoplasia was estimated from the positivity rate in the absence of neoplasia. Medians and interquartile ranges (IQR) were determined where appropriate. Clinical sub-populations were compared using two correlated proportion methodologies for discrete (Z score two-population proportion test and Chi Square test, or Fisher’s exact test when sample size was small) and continuous data (Wilcoxon rank-sum test). One-way ANOVA was used to compare three or more independent populations. McNemar’s test was used for concordance analyses. Multivariable logistic regression was used to explore the impact of gender, age, other comorbidities and yield on assay result. All statistical tests were two-sided and a p value of < 0.05 determined statistical significance. All analyses were performed using GraphPad Prism version 8.2.0 as well as the online tool (http://graphpad.com/scientific-software).

Results

Study population

In total, 1620 specimens were suitable for assay and obtained from patients that met the inclusion criteria. Patient characteristics for each of the four collection sites are shown in Table 1.

Table 1 Demographic characteristics and clinical findings by collection site

The specimens included 184 cases with CRC [median 67.2 years (35–88), 52.7% males], 616 with adenomas [median 63.5 years (33.7–85.4), 62.8% males] and 820 cases without neoplasia [median 60.0 years (18.0–85.0), 51.0% males]. Compared to the other three collection sites, the specimens from Proteogenex had a lower content of males and were obtained from patients 5–6 years younger.

Detection of methylated DNA by clinical phenotype

Of the 1620 plasma specimens, 296 had at least one PCR replicate (18.3%, 95% CI 16.4–20.2) positive for methylation in any one of the three genes. The detection rate was highest in those with CRC (136/184, 73.9% sensitivity), followed by a 15.7% detection rate for advanced adenomas (53/337), which was also significantly higher than that of non-advanced adenomas, 9.3% (26/279, Z score t test p = 0.018) and those without neoplasia [81/820, 9.9% (specificity 90.1%), p = 0.005], Fig. 1.

Fig. 1
figure1

Detection of methylated DNA by clinical phenotype. A specimen was deemed positive (Pos) if at least one PCR replicate was positive for DNA methylation in at least one of the 3 genes. Black closed circles, calculated mean detection rates (%); horizontal bars, 95% CI; severity of dysplasia was available for 584 of 616 adenomas, HGD high grade dysplasia, LGD low grade dysplasia; non-neoplastic pathologies included benign polyps (hyperplastic, unspecified, inflammatory, other polyps), inflammatory bowel disease, diverticular disease, angiodysplasia, hemorrhoids; Odds ratios (95% CI) are relative to cases without neoplasia, **p < 0.0001; *p < 0.05. Two of the 184 CRC cases were unstaged [2/2, 100% (15.8–100)] and were omitted from the figure

When the multi-target assay was positive for methylation in any one of the targeted genes, the odds ratio for presence of CRC was high compared to cases without neoplasia [odds ratio (OR) for CRC 25.9, 95% confidence interval (95% CI): 17.3–38.4, p < 0.0001], Fig. 1. The odds of advanced adenoma was also higher (OR 1.7, 95% CI 1.2–2.5, p = 0.006), whereas the odds ratio for non-advanced adenoma was not different from those without neoplasia (OR 0.9, 95% CI 0.6–1.5, p = 0.907).

The individual markers of methylation in BCAT1, IKZF1 and IRF4 also discriminated well between those with cancer or advanced adenomas, and those with no neoplasia, Table 2.

Table 2 Positivity rate for methylated BCAT1, IKZF1 and IRF4 in each clinical phenotype

Positivity rate by stage of neoplasia

The detection rate of methylated DNA increased significantly from non-advanced adenoma to stage IV CRC (Fig. 1; one-way ANOVA, post-test for trend, p < 0.0001). The increase by stage was also observed for each of the targeted genes, Table 2. The detection rate (sensitivity) for CRC increased significantly with increasing depth of invasion as assessed by T-stage (Fig. 2, one-way ANOVA p < 0.0001). The detection rate in Stage I T2N0M0 cancers was more than two fold higher than that for Stage I T1N0M0 (p = 0.029).

Fig. 2
figure2

Positivity by progression of cellular atypia (dysplasia in adenomas) and degree of invasion (T stage in cancer), in cases without disseminated disease. LGD low grade dysplasia, HGD high grade dysplasia. Y axis: average detection rates, % (sensitivity). Closed black circles—at least one PCR replicate positive for methylation in any of the three genes; blue, red and black symbols—detection rates of the individual methylation markers, IRF4, IKZF1 and BCAT1, respectively. Additional file 1: Table S1 provides details for count of positives, % detected and 95% CI. *Only 53 of the 57 Stage II cases had full TNM information

When classifying adenomas by grade of dysplasia, the detection rate was significantly higher for those with HGD (9/35, 25.7%, 95% CI 14.0–42.3) than for those with LGD (65/549, 11.8% (9.4–14.8), p = 0.017), Fig. 2 and Additional file 1: Table S1. Adenomas with HGD and T1N0M0 cancers returned similar detection rates, 25.7% and 25.0%, respectively. The odds of adenoma with HGD being present given a positive multi-target assay was high compared to no neoplasia (OR 3.2, 95% CI 1.4–6.9, p = 0.007). The odds of an adenoma with LGD was not different from those with no neoplasia (OR 1.2, 95% CI 0.9–1.7, p = 0.284). These positivity patterns were mirrored in each of the targeted three genes, Fig. 2 and Additional file 1: Table S1.

Concordance between methylated genes

Of the 296 specimens positive for at least one of the biomarkers, 69 (23.3%) were positive for methylation in all three targeted genes (i.e. concordance was positive), Fig. 3. The proportion of cases positive for methylation in all three genes was highest for CRC (57/136 positives, 41.9%) and lowest for cases without neoplasia (2/81 positives, 2.5%). When the multi-target assay was positive for DNA methylation in all 3 genes, the odds ratio for presence of cancer compared to those with no neoplasia was 28.5 (95% CI 7.3–121.2, p < 0.0001).

Fig. 3
figure3

Concordance between positive methylated BCAT1, IKZF1 and IRF4 results in plasma. The number of positive specimens is shown for each phenotype, together with the numbers returning a positive for each BCAT1, IKZF1 and IRF4 combination.1There were two unstaged cancer cases, both positive for methylation, which are not included in the CRC Stage Venn diagrams. Severity of dysplasia was available for 584 of 616 adenomas, HGD high grade dysplasia, LGD low grade dysplasia; non-neoplastic pathologies included benign polyps (hyperplastic, unspecified, inflammatory, other polyps), inflammatory bowel disease, diverticular disease, angiodysplasia, hemorrhoids

The proportion of CRC cases positive for methylation in any of the targeted genes rose with increasing stage: Stage I—3/16 (18.8%), Stage II—17/50 (34.0%), Stage III—17/40 (42.5%), Stage IV—19/28 (67.9%), ANOVA post trend test, p = 0.005. Similarly, the three genes were positive in 4/9 (44.4%) adenomas with HGD compared to 5/65 (7.7%; p = 0.01, Fisher’s exact test) with LGD.

Test performance refinement based on different methylation target combinations

Each of the targeted methylated regions contributed to the detection of cancer (true-positives), and different combinations of these targets are shown in Table 3. The best specificity, 97.8%, was achieved using IRF4 as the sole marker of methylated DNA in circulation with a 50% sensitivity. The best sensitivity (73.9%) was achieved using all three genes with a specificity of 90.1%. A similar trend was observed in advanced adenomas and adenomas with HGD, with the best sensitivity for these being achieved using all three genes, Additional file 1: Table S2.

Table 3 Test accuracy for detection of colorectal cancer based on all possible gene combinations

Of the 81 false-positive cases without neoplasia, 67 (82.7%) were positive for methylation in a single gene only, and methylation in BCAT1 was the most frequent cause of the three genes (34/67, 50.7%).

PCR replicate positivity by clinical status

In view of the effect of BCAT1 detection on specificity and as a total of nine PCR replicate results were generated (three for each gene target) per assayed specimen, the relationship between PCR replicate count and clinical status was examined. A comparison of PCR replicate positivity between cases with CRC and those without neoplasia showed that the number of positive PCR replicates was much higher in cases with CRC, Fig. 4. Of the 136 CRC cases with a DNA methylation signal (i.e. at least 1 PCR replicate positive for methylation in either BCAT1, IKZF1 and/or IRF4), 35 (25.7%) cases were methylation positive in all 9 PCR replicates and 26 (19.1%) were methylation positive by just a single PCR replicate. In contrast, none of the 81 positive cases without neoplasia returned 9 positive PCR replicates. Of the 67 cases without neoplasia that were positive for methylation in a single gene only, 62 (92.5%) were positive by just a single PCR replicate. Single PCR replicate positivity was significantly higher than that in CRC (62/81, 76.5% vs. 26/136, 19.1%, Z score p < 0.0001).

Fig. 4
figure4

PCR replicate positivity rates in positive specimens. a Frequency distribution of the number of PCR replicates positive for methylation in patients without neoplasia (n = 81) and with cancer (n = 136). Each assayed specimen generates a total of 9 PCR replicates, 3 for each gene. b In cases (cancer or no neoplasia) where only one replicate was positive, the proportion for each gene responsible for the positive result. FP false positive, TP true positive, TN true negative, FN false negative

Evaluation of a replicate rule for panel positivity

More than half of the single PCR replicate positive cases without neoplasia were due to detection of methylation in BCAT1 (33/62, 53.2%, Fig. 4b). As the odds ratio for a single BCAT1 positive PCR replicate representing a false-positive (i.e. a case without neoplasia) compared to a true-positive (i.e. a case with CRC) was 17.7 (95% CI6.6–43.3), p < 0.0001), a ‘BCAT1 replicate rule’ for assay positivity requiring at least two positive BCAT1 replicates was evaluated. For the ‘BCAT1 replicate rule’, a specimen was deemed positive when at least one PCR replicate was positive for methylation in either IKZF1 or IRF4, or at least two PCR replicates were positive for methylation in BCAT1. With this rule, the false-positive rate for cancer was significantly reduced from 9.9% (81/820) to 5.9% (48/820, 95% CI4.4–7.8, p < 0.0001), while the sensitivity was not significantly affected; 73.9% (136/184, 95% CI 67.1–79.7) versus 71.2% (131/184, 95% CI 64.1–77.6; p = 0.562). To explore the dynamic nature of assay response to number of replicates, a ROC curve was created that yielded an area under the curve (AUC) of 0.827 (0.786–0.867) and the LRP was 12.2 (9.1–16.3). The odds ratio for presence of CRC compared to those without neoplasia was 25.9 (95% CI 17.3–38.4) for the multi-target assay (any gene positive) compared to 39.8 (95% CI 25.7–61.5) if using the ‘BCAT1 replicate rule’. Figure 5 compares ROC analysis of possible combinations of genes and the effect of applying the ‘BCAT1 replicate’ rule for detection of CRC.

Fig. 5
figure5

ROC analysis for potential combinations of different genes for cancer detection. Black dots—one methylation target; grey squares—two methylation targets; black triangle—all three methylation targets. Red triangle: resulting performance specification when applying the ‘BCAT1 replicate rule’ to the 3-gene combination (see text for details)

When applying the replicate rule, the true-positive rates (sensitivity) for detection of advanced adenomas compared to a three gene panel without the rule was significantly lower (53/337, 15.7%, 95% CI 12.2–20.0 vs. 37/337, 11.0%, 95% CI 8.1–14.8, p = 0.048); as was the case for adenomas with HGD (9/35, 25.7%, 95% CI 14.0–42.3 vs. 8/35, 22.9%, 95% CI 11.8–39.3, p = 0.029). Additional file 1: Table S2.

Confounding variables associated with assay positivity in cases without neoplasia

A variety of non-malignant pathologies, comorbidities and demographic factors, such as age and gender, are known to be associated the detection of methylated DNA in blood [16,17,18]. Logistic regression was used to explore the impact of clinical variables on assay positivity in the subgroup without neoplasia, Table 4. Yield of ccfDNA was a significant independent predictor of multi-target assay positivity in cases without neoplasia (OR 1.89, 95% CI 1.30–2.75, p = 0.0008). Male gender and aging showed a positive association with assay positivity but failed to reach significance (p = 0.07–0.10).

Table 4 Predictors of test positivity (any marker detected) in cases without neoplasia

As the yield of ccfDNA was the only independent predictor of assay positivity in cases without neoplasia, we looked at known sources of ccfDNA variance [19, 20]. There was a significant relationship between yield of ccfDNA and increasing age (p = 0.0001), inflammatory disorders (p = 0.008) and male gender (p = 0.015).

Including all 1620 specimens, ccfDNA levels varied greatly (median 3.8 ng/mL, IQR 2.5–5.5) but were only significantly elevated in cases with stage III (4.7 ng/mL, IQR 2.6–9.1) or IV (5.6 ng/mL, IQR 3.8–9.9) cancer compared to cases with no evidence of disease (3.6 ng/mL, IQR 2.4–5.2, p < 0.0001).

Discussion

Following initial discovery of potential biomarkers for the detection of CRC, the next step is to determine the efficient and effective biomarker combinations that are best able to discriminate between those with and without cancer. This evaluation is an important step before moving to translational research in a screening context [21]. We have previously reported on genes, including BCAT1 and IKZF1, that contain regions methylated with high frequency in neoplastic colonic lesions and detection of such DMRs in DNA in circulation for indication of CRC [10, 11]. In an endeavor to achieve better sensitivity, we developed a multi-target qPCR assay detecting an additional DMR in IRF4 as well as detecting an additional DMR in IKZF1. The resulting multi-target qPCR assay was evaluated in a population comprising the spectrum of pathologies typically encountered when screening for CRC and which was more diverse than the cohort used in our initial discovery [10].

We found that each of the targeted DMRs residing in BCAT1, IKZF1 and IRF4 contributed to detection of CRC and advanced adenomas (especially those with HGD) with a significantly higher odds ratio for either disease state being present. The multi-target qPCR assay discriminated between CRC or advanced adenomas compared to those without neoplasia and was more sensitive than previously observed when just 2 DMRs were included.

Positive concordance between the three genes was highest in those with cancer and also higher for those with HGD adenomas compared to cases without neoplasia. Overall, this concordance reflects an increasing rate of aberrant methylation as neoplasia progresses.

As expected, increasing the number of targets resulted in a lower specificity than previously observed [10, 11]. Detection of circulating DNA methylated in only BCAT1 was the most prevalent event leading to a false-positive assay result.

The assay positivity rate increased as one progressed from adenoma through to late stage cancer indicating that presence of hypermethylated tumor DNA in plasma was dependent on stage of cancer, especially T stage, but also on degree of dysplasia (the initial morphological change characterizing neoplasia) in adenomas. Adenomas with HGD were detected at a similar rate to T1N0M0 CRC, which was a significantly higher rate than those with non-advanced adenomas or adenomas with LGD. The detection rate of advanced adenomas (defined by commonly accepted clinical criteria) was significantly higher than the detection rate of those with non-advanced adenomas but significantly lower those with HGD. This positivity pattern is not unexpected since the clinical definition for advanced adenoma include states that relate to risk of developing metachronous neoplasia at a later stage and are not restricted to morphological features characterising progressive cellular atypia at the time of testing for the methylation biomarkers. These findings for HGD adenomas are novel and point to potential for using this class of ctDNA biomarkers to target clinically relevant adenomas as well as cancer.

The relationship between marker detection and cancer stage agrees with previous observations, which indicates that as neoplastic lesions progress and especially as they invade, the number of ctDNA molecules increases due to egress into the circulation [22, 23]. Each of the targeted ctDNA methylation markers had a 39–52% sensitivity for early-stage CRC (Table 2) which aligns with performance reported for other single-target somatic- or epigenetic-based ctDNA tests [3, 24]. Enhancing the analytical sensitivity of the assay by targeting just 44 methylated CpG sites in 4 DMRs residing in 3 genes (Figure S1) resulted in a 67.4% sensitivity for early stage CRC which is comparable to the 60.4% reported sensitivity for early stage CRC using a ctDNA blood test targeting more than 28,000 DMRs [25].

This study demonstrates that, with the right biomarkers, it is possible to design a relatively simple ctDNA test that can detect CRC with good sensitivity without the need for large dimensional data sets such as those involving next generation sequencing. Importantly, our strategy avoids the inflexibility of complex marker algorithms and shows that by using different marker combinations, an end-user has flexibility to adjust the assay output to suit the clinical goals (especially important in screening), such as maximising detection (sensitivity) while controlling cost-effectiveness (specificity) or feasibility. For instance, if a high specificity is required but a 50% sensitivity for CRC is acceptable, as is already the case in some screening programs around the world that use a conservative fecal haemoglobin concentration cut-off for FIT [26], this is achievable by detection of just methylated IRF4 DNA in blood alone, although this approach would require further validation in a true screening context. But other jurisdictions consider a higher sensitivity for CRC, as observed with the multitarget assay here, as being desirable. The risk of adding biomarkers is that specificity falls as one strives for higher sensitivity. By considering the specificity of each individual biomarker, however, we have been able to improve specificity of the panel from 90.1% to 94.1% with little compromise in sensitivity for CRC (73.9% vs. 71.2%). This improvement was achieved by defining a specimen positive for methylation if there was at least one PCR replicate positive for methylation in IKZF1 or IRF4 (irrespective of BCAT1 positivity), or at least two PCR replicates positive for methylation in BCAT1 (the “replicate rule”). Using this rule, the odds ratio for presence of CRC was 39.8 when positive for the rule compared to an OR of 25.9 for the full biomarker panel when not applying the replicate rule. This improvement was because BCAT1 was the major contributor to deterioration in specificity. A high replicate count was characteristic of patients with CRC (a quarter of cancer cases were positive in all nine replicates) while in the majority of cases without neoplasia, only one replicate was positive, usually BCAT1. Based on the cohort tested herein, the odds ratio of a single BCAT1 positive PCR replicate being a false-positive compared to a true positive was 17.7.

The three targeted genes examined in this study may play a role in the tumorigenesis of CRC. Differential methylation in the BCAT1 promoters alters the ratio of generated BCAT1 protein isoforms [27], and aberrant expression of BCAT1 has been associated with CRC [28]. The promoters of IKZF1 and IRF4 are silenced when hypermethylated [29, 30]. The transcriptional factors, IKZF1 and IRF4 are important transcriptional regulators of notch and c-myc [29, 31,32,33], and c-MYC have been linked to BCAT1 expression [34]. As well-described for c-MYC, both IKZF1 and IRF4 are involved in the development of cancer, including CRC [29, 35, 36].

This study has strengths and limitations. As the estimates of test accuracy are derived from biobanked specimens from four sites where patients were undergoing diagnostic assessment by colonoscopy for a wide range of clinical applications, the actual accuracy estimates might not be the same in an unbiased typical screening population. Thus, it will be important to proceed with application of the panel described here in a prospective population screening study, ideally compared with another proven screening test such as the fecal immunochemical test. The diverse nature of the population, however, has advantages given the broad spread of colorectal pathologies. We were able to comprehensively assess the impact of non-neoplastic conditions and other possible confounding variables associated with detection of these biomarkers that might affect specificity in a large number of subjects. The concentration of ccfDNA was the only significant factor affecting the multi-target assay response in the subgroup of cases without neoplasia. Male gender and those aged over 80 years showed a weak positive association with assay result but failed to reach significance.

As observed by others, the levels of ccfDNA varied greatly in cases without neoplasia [18, 37,38,39], and there was a significant relationship between higher ccfDNA concentrations and increasing age.

When screening for CRC, it is ideal to use a test capable of detecting advanced adenomas given that their removal reduces incidence of CRC. Thus, endoscopic screening for CRC, or use of high-sensitivity FIT (with sensitivity for advanced adenomas in the order of 40%), is adopted in some jurisdictions. By targeting multiple regions methylated with high frequency and at the earliest onset of colorectal neoplasia, we have achieved a sensitivity of 25% for adenomas exhibiting HGD which raises the prospect of further success with future modifications of the panel.

Similar observations of increased adenoma detection have been observed for the well-studied methylation biomarker, SEPT9 when studied in a small cohort of 76 cases in conjunction with another methylation marker, ALX4. The additional of ALX4 increased sensitivity for advanced precancerous colorectal lesions from 12 to 45% (6/49 vs. 22/49) [24, 40] compared to SEPT9 alone. Nevertheless, the specificity in that study also fell significantly from 95.5% (1/22) to 82% (4/22) with the additional biomarker demonstrating that the higher positivity was not specific for neoplasia. The present study did not suffer such a large increase in false positives.

Our findings show that detection of methylated BCAT1, IKZF1 and IRF4 in circulating ccfDNA differentiates cases with CRC from those without neoplasia, and that the specificity of the multi-target assay can be substantially improved with no significant effect on sensitivity by applying a PCR replicate rule to BCAT1. This panel of markers should now be prospectively evaluated in a typical screening population to clarify accuracy of different marker configurations against that of FIT. As there is considerable global variation in what constitutes acceptable sensitivity and/or specificity [1], the flexibility which is provided by using different configurations will allow health care provides to choose a performance that suits the goals of specific screening programs.

Availability of data and materials

The dataset analyzed during this study are available from the corresponding author on reasonable requests.

References

  1. 1.

    Young GP, Rabeneck L, Winawer SJ. The global paradigm shift in screening for colorectal cancer. Gastroenterology. 2019;156:843–51.

    Article  Google Scholar 

  2. 2.

    Young GP, Symonds EL, Allison SJ, Cole SR, Fraser CG, Halloran SP, Kuipers EJ, Seaman HE. Advances in fecal occult blood tests: the FIT revolution. Dig Dis Sci. 2014;60:1–14.

    Google Scholar 

  3. 3.

    Bettegowda C, Sausen M, Leary RJ, Kinde I, Wang Y, Agrawal N, et al. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci Transl Med. 2014;6:224.

    Article  Google Scholar 

  4. 4.

    Lao VV, Grady WM. Epigenetics and colorectal cancer. Nat Rev Gastroenterol Hepatol . 2011;8:686–700.

    CAS  Article  Google Scholar 

  5. 5.

    Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA, Kinzler KW. Cancer genome landscapes. Science. 2013;339:1546–58.

    CAS  Article  Google Scholar 

  6. 6.

    Weisenberger DJ. Characterizing DNA methylation alterations from the cancer genome atlas. J Clin Invest. 2014;124:17–23.

    CAS  Article  Google Scholar 

  7. 7.

    Pepe MS, Etzioni R, Feng Z, Potter JD, Thompson ML, Thornquist M, et al. Phases of biomarker development for early detection of cancer. JNCI. 2001;93:1054–61.

    CAS  Article  Google Scholar 

  8. 8.

    Mitchell SM, Ross JP, Drew HR, Ho T, Brown GS, Saunders NFR, et al. A panel of genes methylated with high frequency in colorectal cancer. BMC Cancer. 2014;14:54.

    Article  Google Scholar 

  9. 9.

    Mitchell SM, Ho T, Brown G, Baker R, Thomas M, McEvoy A, et al. Evaluation of methylation biomarkers for detection of circulating tumor DNA and application to colorectal cancer. Genes. 2016;7:125.

    Article  Google Scholar 

  10. 10.

    Pedersen SK, Symonds EL, Baker RT, Murray DH, McEvoy A, van Doorn SC, et al. Evaluation of an assay for methylated BCAT1 and IKZF1 in plasma for detection of colorectal neoplasia. BMC Cancer. 2015;15:654.

    Article  Google Scholar 

  11. 11.

    Symonds EL, Pedersen SK, Baker RT, Murray DH, Gaur S, Cole SR, et al. A blood test for methylated BCAT1 and IKZF1 vs. a fecal immunochemical test for detection of colorectal neoplasia. Clin Transl Gastroenterol. 2016;7:137.

    Article  Google Scholar 

  12. 12.

    Redshaw N, Huggett JF, Taylor MS, Foy CA, Devonshire AS. Quantification of epigenetic biomarkers: an evaluation of established and emerging methods for DNA methylation analysis. BMC Genomics. 2014;15:1174.

    Article  Google Scholar 

  13. 13.

    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.

    Article  Google Scholar 

  14. 14.

    Chung SJ, Kim YS, Yang SY, Song JH, Kim D, Park MJ, et al. Five-year risk for advanced colorectal neoplasia after initial colonoscopy according to the baseline risk stratification: a prospective study in 2452 asymptomatic Koreans. Gut. 2011;60:1537–43.

    Article  Google Scholar 

  15. 15.

    Martinez ME, Baron JA, Lieberman DA, Schatzkin A, Lanza E, Winawer SJ, et al. A pooled analysis of advanced colorectal neoplasia diagnoses after colonoscopic polypectomy. Gastroenterology. 2009;136:832–41.

    Article  Google Scholar 

  16. 16.

    Tsang SY, Ahmad T, Mat FWK, Zhao C, Xiao S, Xia C, Xue H. Variation of global DNA methylation levels with age and in autistic children. Hum Genomics. 2016;10:31–6.

    Article  Google Scholar 

  17. 17.

    Zhang FF, Cardarelli R, Carroll J, Fulda KG, Kaur M, Gonzalez K, et al. Significant differences in global genomic DNA methylation by gender and race/ethnicity in peripheral blood. Epigenetics. 2011;6:623–9.

    CAS  Article  Google Scholar 

  18. 18.

    Teo YV, Capri M, Morsiani C, Pizza G, Faria AMC, Franceschi C, Neretti N. Cell-free DNA as a biomarker of aging. Aging Cell. 2019;18:e12890.

    Article  Google Scholar 

  19. 19.

    Kustanovich A, Schwartz R, Peretz T, Grinshpun A. Life and death of circulating cell-free DNA. Cancer Biol Therapy. 2019;20:1057–67.

    CAS  Article  Google Scholar 

  20. 20.

    Volik S, Alcaide M, Morin RD, Collins C. Cell-free DNA (ccfDNA): Clinical Significance and utility in cancer shaped by emerging technologies. Mol Cancer Res. 2016;14:898–908.

    CAS  Article  Google Scholar 

  21. 21.

    Young GP, Senore C, Mandel JS, Allison JE, Atkin WS, Benamouzig R, et al. Recommendations for a step-wise comparative approach to the evaluation of new screening tests for colorectal cancer. Cancer. 2016. https://doi.org/10.1002/cncr.29865.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Diehl F, Li M, Dressman D, He Y, Shen D, Szabo S, Diaz LA, et al. Detection and quantification of mutations in the plasma of patients with colorectal tumors. PNAS. 2005;102:16368–73.

    CAS  Article  Google Scholar 

  23. 23.

    Fiala C, Diamandis EP. Utility of circulating tumor DNA in cancer diagnostics with emphasis on early detection. BMC Med. 2018;16:166.

    CAS  Article  Google Scholar 

  24. 24.

    Church TR, Wandell M, Lofton-Day C, Mongin SJ, Burger M, Payne SR, et al. Prospective evaluation of methylated SEPT9 in plasma for detection of asymptomatic colorectal cancer. Gut. 2014;63:317–25.

    CAS  Article  Google Scholar 

  25. 25.

    Liu MC, Oxnard GR, Klein EA, Swanton C, Seiden MC, et al. Sensitive and specific multi-cancer detection and localization using methylation signatures in cell-free DNA. Ann Oncol. 2020;31:745–59.

    Article  Google Scholar 

  26. 26.

    Clark G, Strachan JA, Carey FA, Godfrey T, Irvine A, McPherson A, et al. Transition to quantitative faecal immunochemical testing from guaiac faecal occult blood testing in a fully rolled-out population-based national bowel screening programme. Gut. 2020. https://doi.org/10.1136/gutjnl-2019-320297.

    Article  PubMed  Google Scholar 

  27. 27.

    Tonjes M, Barbus S, Park YJ, Wang W, Schlotter M, Lindroth AM, et al. BCAT1 promotes cell proliferation through amino acid catabolism in gliomas carrying wild-type IDH1. Nature. 2013;19:901–8.

    Google Scholar 

  28. 28.

    Yoshikawa R, Yanagi H, Shen CS, Fujiwara Y, Noda M, Yagyu T, et al. ECA39 is a novel distant metastasis-related biomarker in colorectal cancer. WJG. 2006;12:5884–9.

    CAS  Article  Google Scholar 

  29. 29.

    Javierre BM, Rodriguez-Ubreva J, Al-Shahrour F, Corominas M, Graña O, Ciudad L, et al. Long-range epigenetic silencing associates with deregulation of Ikaros Targets in colorectal cancer cells. Mol Cancer Res. 2011;9:1139–51.

    CAS  Article  Google Scholar 

  30. 30.

    Ortmann CA, Burchert A, Hölzle K, Nitsche A, Wittig B, Neubauer A, et al. Down-regulation of interferon regulatory factor 4 gene expression in leukemic cells due to hypermethylation of CpG motifs in the promoter region. NAR. 2005;33:6895–905.

    CAS  Article  Google Scholar 

  31. 31.

    Patil A, Manzano M, Gottwein E. CK1α and IRF4 are essential and independent effectors of immunomodulatory drugs in primary effusion lymphoma. Blood. 2018;132:577–86.

    CAS  Article  Google Scholar 

  32. 32.

    Pathak S, Ma S, Trinh L, Eudy J, Wagner KU, Joshi SS, et al. IRF4 is a suppressor of c-Myc induced B cell leukemia. PLoS ONE. 2011;6:e22628.

    CAS  Article  Google Scholar 

  33. 33.

    Palomero T, Lim WK, Odom DT, Sulis ML, Real PJ, Margolin A, et al. NOTCH1 directly regulates c-MYC and activates a feed-forward-loop transcriptional network promoting leukemic cell growth. PNAS. 2006;103:18261–6.

    CAS  Article  Google Scholar 

  34. 34.

    Ben-Yosef T, Eden A, Benvenisty N. Characterization of murine BCAT genes: Bcat1, a c-Myc target, and its homolog, Bcat2. Mamm Genome. 1998;9:595–7.

    CAS  Article  Google Scholar 

  35. 35.

    Slattery ML, Lundgreen A, Bondurant KL, Wolff RK. Interferon-signaling pathway: associations with colon and rectal cancer risk and subsequent survival. Carcinogenesis. 2011;32:1660–7.

    CAS  Article  Google Scholar 

  36. 36.

    Sharrard RM, Royds JA, Rogers S, Shorthouse AJ. Patterns of methylation of the c-myc gene in human colorectal cancer progression. Br J Cancer. 1992;65:667–72.

    CAS  Article  Google Scholar 

  37. 37.

    Duvvuri B, Lood C. Cell-Free DNA as a Biomarker in Autoimmune Rheumatic Diseases. Front Immunol. 2019;10:502.

    CAS  Article  Google Scholar 

  38. 38.

    Frank MO. Circulating cell-free DNA differentiates severity of inflammation. Biol Res Nurs. 2016;18:477–88.

    CAS  Article  Google Scholar 

  39. 39.

    van der Vaart M, Pretorius PJ. Characterization of circulating DNA in healthy human plasma. Clin Chim Acta. 2008;395:186.

    Article  Google Scholar 

  40. 40.

    Tänzer M, Balluff B, Distler J, Hale K, Leodolter A, Röcken C, et al. Performance of epigenetic markers SEPT9 and ALX4 in plasma for detection of colorectal precancerous lesions. PLoS ONE. 2010;5:e9061.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported in part by Clinical Genomics Pty Ltd and by grants received by the National Health & Medical Research Council (APP1006242 and APP1017083).

Author information

Affiliations

Authors

Contributions

GPY and SKP conducted the study. ES, HJN, LF, ED, MV, IJC, LF, MV, ES and RMH were responsible for the data collection and auditing of cases. NB, BY, MC and GT were responsible for assay testing and reporting. LCL, SKP and GPY were responsible for study design. GPY and SKP drafted the manuscript. All authors read and approved the final manuscript to be published.

Corresponding author

Correspondence to Susanne K. Pedersen.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Southern Adelaide Clinical Human Research Ethics Committee (April 4, 2005; ethics number 134.045), the Medical Ethical Board of Academic Medical Centre Amsterdam (July 12, 2011), and the Danish Regional Ethics Committee (H-4-2013-050) and Data Protection Agency (2007-58-0015/HVH-2013-022).

Competing interests

GPY is a paid consultant to Clinical Genomics. SKP, BY, MC, GT, LCL, NB are paid employees of Clinical Genomics. The other authors declare that they have no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1

. The multi-panel real-time PCR asssay. Table S1. Positivity by progression of cellular atypia (dysplasia in adenomas) and degree of invasion (T stage in cancer) in cases without disseminated disease. Table S2. Test accuracy for detection of advanced adenomas and adenomas with high grade dysplasia based on all possible combinations of genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Young, G.P., Symonds, E.L., Nielsen, H.J. et al. Evaluation of a panel of tumor-specific differentially-methylated DNA regions in IRF4, IKZF1 and BCAT1 for blood-based detection of colorectal cancer. Clin Epigenet 13, 14 (2021). https://doi.org/10.1186/s13148-020-00999-y

Download citation

Keywords

  • Methylated circulating tumor DNA
  • BCAT1
  • IKZF1
  • IRF4
  • Colorectal cancer screening
\