Skip to main content

Chromatin conformation changes in peripheral blood can detect prostate cancer and stratify disease risk groups

Abstract

Background

Current diagnostic blood tests for prostate cancer (PCa) are unreliable for the early stage disease, resulting in numerous unnecessary prostate biopsies in men with benign disease and false reassurance of negative biopsies in men with PCa. Predicting the risk of PCa is pivotal for making an informed decision on treatment options as the 5-year survival rate in the low-risk group is more than 95% and most men would benefit from surveillance rather than active treatment. Three-dimensional genome architecture and chromosome structures undergo early changes during tumourigenesis both in tumour and in circulating cells and can serve as a disease biomarker.

Methods

In this prospective study we screened whole blood of newly diagnosed, treatment naïve PCa patients (n = 140) and cancer-free controls (n = 96) for the presence of 14,241 chromosomal loops in the loci of 425 genes.

Results

We have detected specific chromosome conformation changes in the loci of ETS1, MAP3K14, SLC22A3 and CASP2 genes in peripheral blood from PCa patients yielding PCa detection with 80% sensitivity and 80% specificity. Further analysis between PCa risk groups yielded prognostic validation sets consisting of HSD3B2, VEGFC, APAF1, BMP6, ERG, MSR1, MUC1, ACAT1 and DAPK1 genes that achieved 80% sensitivity and 93% specificity stratifying high-risk category 3 vs low risk category 1 and 84% sensitivity and 89% specificity stratifying high risk category 3 vs intermediate risk category 2 disease.

Conclusions

Our results demonstrate specific chromosome conformations in the blood of PCa patients that allow PCa diagnosis and risk stratification with high sensitivity and specificity.

Background

In the Western world prostate cancer (PCa) is now the most commonly diagnosed non-cutaneous cancer in men and is the second leading cause of cancer-related death [1, 2]. Many men as young as 30 show evidence of histological PCa [3, 4], most of which is microscopic and some may not manifest clinically. For the diagnosis and prognosis, prostate specific antigen (PSA), magnetic resonance imaging (MRI), an invasive needle biopsy, Gleason score and disease stage are used [5, 6].

The only available blood test for PCa in widespread clinical use involves measuring circulating levels of PSA (at 4 ng/ml has 21% sensitivity and 93% specificity, PCa prevention trial) [7], however, the prostate size, benign prostatic hyperplasia and prostatitis may also increase PSA levels. PSA cut off levels are age-specific. American urological association thresholds (ng/ml) are as follows: age 40–49 = 2.5; age 50–59 = 3.5; 60–69 = 4.5; 70–79 = 6.5. In the UK, in men aged 50–69 years, a reading of 3.0 ng/ml or higher leads to urgent referral for further investigation. However, men with low PSA can have cancer. In a large US study, 2950 men (age range, 62 to 91 years) with PSA < 4.0 ng/ml underwent biopsy. PCa was diagnosed in 449 (15.2%); 67 of these 449 cancers (14.9%) had a Gleason score of 7 or higher [8]. Furthermore, younger patients with Gleason 8–10 cancer and PSA levels of < 4.0 ng/ml have more aggressive disease than those with PSA levels of 4 to 9.9 ng/ml [9]. In early PCa, PSA testing is not specific enough to differentiate between early-stage invasive cancers and latent, non-lethal tumours that might otherwise have remained asymptomatic during a man’s lifetime. In advanced PCa, PSA kinetics are used as a clinical surrogate endpoint for outcome. However, while they do give a general prognosis they lack specificity for the individual [10]. A number of more specific blood tests are emerging for PCa detection including 4 K blood test (AUC 0.8) and PHI blood test (90% sensitivity, 17% specificity) [11]. Both PHI and 4 K tests have a proposed prognostic value [12], however, its utility is limited [13] and it is currently not widely used in clinical practice. PSA levels, clinical disease stage and Gleason score are used to establish the severity of PCa and stratify patients to risk groups [13]. To date, there is no prognostic blood test available that allows differentiation between low- and high-risk PCa.

Genome-wide association studies (GWAS) have identified more than 100 PCa susceptibility loci, explaining 30% of the familial risk for this disease [14]. However, most of PCas are not familial and are linked to unidentified multiple somatic mutations. These include mutations in p53 (up to 64% of tumours), p21 (up to 55%), p73 and MMAC1/PTEN tumour suppressor genes [15], as well as alterations in 25 PCa susceptibility loci [16], but these mutations do not explain all the observed effects on gene regulation. Indeed, one interesting observation that arose from GWAS, was that most of the loci within the genome that confer risk to diseases including cancer are located outside of known protein-coding regions [17].

Epigenetic modifications play an important role in PCa progression [18, 19]. Aberrant DNA methylation (hypo- and hypermethylation) is the best-characterized cancer-related epigenetic alteration in PCa [20, 21]. Histone modifications also contribute to PCa progression [22, 23]. We have previously found that in human cells epigenetic mechanisms involving dynamic and multi-layered chromosomal loop interactions are powerful regulators of gene expression [24]. Chromosome conformation capture (3C) technologies allow these signatures to be recorded and have gained considerable attention for disease diagnosis [25,26,27,28]. Prostate tumours have been shown to undergo long-range epigenetic alterations in 3-dimension chromosome conformations and distinct epigenetic signatures were found in circulating DNA from PCa patients [29]. Another study identified a strong association of single nucleotide polymorphism rs11672691 at 19q13 with clinical features indicative of aggressive PCa. Mechanistically, this locus resides in an enhancer element and alters the binding site of HOXA2, a novel oncogenic transcription factor [30].

We have developed a novel epigenetic assay, as a next generation of the 3C technique [31]. EpiSwitch™ technology employs 3C technique and algorithmic-based analysis to identify a panel of epigenetic differences capable of distinguishing between cancer-free controls and diseased tissue samples. It detects epigenetic regulatory signature changes in the higher order structures of human chromosomes at the loci implicated in the onset and progression of the disease. It offers highly effective means of screening; early detection; monitoring and prognostic analysis of major diseases associated with aberrant gene expression [24]. One of the main advantages of using 3C-based chromatin interactions as biomarkers is that DNA cross-linking is relatively stable, and following proximity ligation, give rise to a stable DNA product (Fig. 1) [32]. Using EpiSwitch™ technology, we have shown the presence of melanoma-specific chromatin conformations in peripheral blood mononuclear cells (PBMCs) and primary tumours of melanoma patients [33, 34]. Fractionation studies showed that the detected signature comes from lymphocytes and not circulating tumour cells [33, 34]. The concentration of circulating tumour cells in circulation is too low to allow detection of chromosomal conformations, while circulating free DNA does not retain 3D conformations [33, 34]. In this study, we used the EpiSwitch™ assay to screen for, define and evaluate specific chromosome conformations in the blood of PCa patients and to identify loci with potential to act as diagnostic and prognostic markers.

Fig. 1
figure 1

Chromosomal loops structure and conformation assay. a Chromosomal loop contains the enhancer region that increases target gene promoter activity. b During cancer progression there is increased looping in some tumour-related genes leading to abnormal gene expression. c Chromosome conformation assay: DNA is crosslinked, digested, ligated and new sequences (in places where loops were) are predicted using relevance machine vector. Loops presence is then detected using DNA microarray. Resulting markers are analysed using multivariate analysis yielding specific epigenetic signatures for selected patient cohorts

Methods

Patient population

A total of 140 PCa patients and 96 cancer-free controls were recruited, in two cohorts. Cohort 1: men with (n = 105) or without (n = 77) PCa diagnosis attending urology clinic at Imperial College Healthcare NHS Trust were prospectively recruited from October 2010 through September 2013. Cohort 2: Patients' samples (19 cancer-free controls and 35 PCa) obtained from Janssen (USA). Upon recruitment, a single blood sample (5 ml) was collected from PCa patients using the current practice for needle and blood collection methods into the BD Vacutainer® plastic EDTA tubes. Blood samples were passively frozen and stored at − 80 °C until processed. Prostate tumour samples were obtained from previously recruited patients (n = 5) that subsequently underwent radical prostatectomy at Imperial College Healthcare NHS Trust. Patient clinical characteristics are shown in Additional file 1: Table S1. The primary endpoint of this study was to detect changes in chromosomal conformations in PBMCs from PCa patients in comparison to cancer-free controls. Therefore, all treatment naïve PCa patients were eligible for this study irrespective of grade, stage and PSA levels. Patients that had previous chemotherapy or patients with other cancers were excluded from this study. PCa diagnosis was established as per clinical routine and patients were assigned to appropriate treatment. For prognostic study (secondary endpoint), patients were stratified according to the relevant National Comprehensive Cancer Network (NCCN) risk groups (Additional file 1: Table S2). No follow up study was conducted. All patients who participated in the trial have given informed consent. The study was approved by the local ethics committee (Imperial College NHS Trust REF 07/MRE09/54). All procedures and protocols were performed in accordance with the relevant guidelines and regulations. Based on the preliminary findings in melanoma [33, 34], an a priori power analysis was performed using the pwr.t.test function in the R package pwd.

Stepwise diagnostic and prognostic biomarker discovery process using EpiSwitch™ technology and data analysis

EpiSwitch™ technology platform (Oxford Biodynamics, Oxford, UK) pairs high resolution 3C results with regression analysis and a machine learning algorithm to develop disease classifications [33,34,35,36]. To select epigenetic biomarkers that can diagnose cancers, samples from patients suffering from cancer, in comparison to cancer-free control samples were screened for statistically significant differences in conditional and stable profiles of genome architecture (Fig. 1). All the samples used for the nested polymerase chain reaction (PCR) diagnostic and prognostic biomarker discovery were processed as per the manufacturer’s instructions using the proprietary EpiSwitch™ reagents and standard protocols (the standard protocol for nested PCR is Work Instruction 52 (WI52) and for real-time quantitative PCR (qPCR) is WI82 [36]). The samples for the development of the qPCR process were processed using the EpiSwitch™ reagents and protocols as described previously [36]. Briefly, 50 µl of whole blood was suspended in 450 µl PBS, formaldehyde containing solution was added and the mixture incubated at room temperature for 15 min to fix the samples. Glycine solution was added to bind any unreacted formaldehyde and quench the reaction. Cell membrane lysis was performed, and the intact nuclei purified using a density cushion centrifugation. TaqI restriction enzyme digestion (Fermentas FastDigest™ TaqI kit) was performed and proximity ligation with T4 DNA ligase (Takara T4 Ligase kit) was used to capture the inter/intra chromatin associations as chromosome conformation signatures. Proteinase K (Sigma) was added to the libraries to digest the constraining protein content of the samples and allow for greater accessibility in the PCR. The Proteinase K was denatured, and cross-links partially reversed at 90 °C. The protocols were repeated to produce the required volume of library for nested PCR investigation. The same process was applied to the matched tissue samples with the addition of an initial incubation in 0.125% (v/v) Collagenase (Sigma) for 37 °C with gentle agitation for 30 min until the tissue was suspended. The freed cells were pelleted by centrifugation at 800g for 5 min and the supernatant removed before being resuspended in 500 µl PBS and subjected to the same protocol as described previously [37]. In short, the assay was performed on a whole blood sample by first fixing chromatin with formaldehyde to capture intrachromatin associations, as described before [37]. The fixed chromatin was then digested into fragments with TaqI restriction enzyme, and the DNA strands were joined favouring cross-linked fragments. The cross-links were reversed and PCR performed using the primers previously established by the EpiSwitch™ software (Fig. 1).

The version 2 EpiSwitch™ reagents and protocols were used to generate the libraries for the qPCR translation stage. These follow the protocol as described previously [37] with an additional detergent-based step between the density cushion centrifugation and TaqI restriction to increase enzyme accessibility and restriction rates. In addition, the Proteinase K treatment was replaced with the use of the QIAmp® DNA FFPE Tissue Kit (Qiagen), starting at the step of resuspension of sample in ATL buffer.

EpiSwitch™ was used on blood samples in a three-step process to identify, evaluate, and validate statistically significant differences in chromosomal conformations between PCa patients and cancer-free controls (Fig. 2). For the first step, sequences from 425 manually curated PCa-related genes (obtained from the public databases (www.ensembl.org)) were used as templates for this computational probabilistic identification of regulatory signals involved in chromatin interaction (Additional file 1: Table S3). A customized CGH Agilent microarray (8 × 60 k) platform was designed to test technical and biological repeats for 14,241 potential chromosome conformations across 425 genetic loci. Eight PCa and eight cancer-free control samples were competitively hybridized to the array, and differential presence or absence of each chromosome conformation was defined by LIMMA linear modelling with empirical Bayes moderation of the standard errors, subsequent abundance filtering and cluster analysis. This initially revealed 60 chromosomal interactions with the ability to best discriminate PCa patients from cancer-free controls (Fig. 2). For the second evaluation stage, the 60 chromosomal interactions were translated into EpiSwitch™ PCR based-detection assays. Nested PCR primers were designed for the 60 chromosomal interactions and these tested using a pooled PCa and a pooled cancer-free control, to ascertain if they produced the expected PCR product. Of these 60 primer combinations, 7 failed to produce PCR products and were subsequently dropped from further investigation. The 53 successful primer combinations were screened on PCa and cancer-free controls (n = 6 in each group). Feature reduction of the 53 PCR was performed using an univariant statistical approach (Fishers Exact test), this analysis identified 15 nested PCR markers with the most statistical discernibility between PCa and cancer-free controls.

Fig. 2
figure 2

Three-step approach to identify, evaluate, and validate diagnostic and prognostic biomarkers for prostate cancer (PCa)

The known sample cohort used to identify the 15 markers was expanded to include 49 samples (PCa n = 24, cancer-free controls n = 25) to provide better statistical power for marker selection. With the expanded cohort univariate (Fishers Exact test) and multivariant (GLMNET, penalized logistic regression, alpha 1 (lasso) and lambda 0.01) statistical approaches were used to further identify markers with the greatest statistical discernibility. GLMNET was used in combination with permutation, where the cohort was split 500 times randomly into a training (0.667) and test (0.333). A GLMNET model is developed for each data split (500) and the markers used in the model are given a vote, the frequency of use of the marker is tallied, markers used in all data splits would get be voted 500 times. This resampling methodology counters false positives in terms of marker selection. In combination with the univariate test for independence (Fishers Exact test), the final model identified was made up of a five-marker signature (Table 1). Principal component analysis was also used to determine group variance and to identify potential outlier samples. The 49 known samples used for marker reduction were used to create a logistic model classifier with five-fold cross validation. The five markers were further statistically assessed using 500 randomised classifiers, 250 logistic model and 250 naïve bayes classifiers with data resampling in each iteration.

Table 1 Five-marker signature used for the diagnosis of prostate cancer (PCa)

For the last step, to further validate the chromosome conformation signature used to inform PCa diagnosis, the training sample cohort size was expanded to 95 PCa and 96 cancer-free control samples. The five-marker set was tested on a blinded, independent (n = 20) cohort of blood samples (Additional file 1: Table S4). The results were analysed using Bayesian Logistic modelling, p-value null hypothesis (Pr(N|z|) analysis, Fisher-Exact P test and Glmnet.

Data analysis and presentation were performed in accordance with CONSORT recommendations. All measurements were performed in a blinded manner. STARD criteria have been used to validate the analytical procedures. A similar stepwise approach was followed for the identification of prognostic markers. A subset of prognostic markers was translated from the nested PCR format to qPCR assays with hydrolysable probes, to the industry standards of MIQE-compliant qPCR. The alternative qPCR format assays were validated on a subset of the prognostic sample cohort.

Nested polymerase chain reaction (PCR)

Sequence specific oligonucleotides were designed around the chosen sites for screening potential markers by nested PCR using Primer3. All PCR amplified samples were visualized by electrophoresis in the LabChip GX, using the LabChip DNA 1 K Version2 kit (Perkin Elmer, Beaconsfield, UK) and internal DNA marker was loaded on the DNA chip according to the manufacturer’s protocol using fluorescent dyes. Fluorescence was detected by laser and electropherogram read-outs translated into a simulated band on gel picture using the instrument software. The threshold we set for a band to be deemed positive was 30 fluorescence units and above.

Real-time quantitative PCR (qPCR) assay development

Both the inner and outer primer combinations of primers for markers identified in the nested PCR investigation were subjected to a temperature gradient PCR to identify the optimum annealing temperature using EpiSwitch™ libraries as the template inputs. The PCR products were detected using SYBR green amplification and melt curves. All PCR amplified samples were visualised by electrophoresis. The PCR products were gel purified, sequenced and mapped to the human reference genome to validate the 3C products. The sequence verified PCR products were used to create standards for absolute quantification starting at a 106 copies standard. Hydrolysis probes were designed across the restriction/ligation junction of the 3C products to produce specific detection of the 3C products of interest. The hydrolysis probe detection was again subjected to a temperature gradient PCR to identify the optimum annealing/extension temperature for the specific detection of the PCR products. The generated 106 copy standard was used as the positive control alongside standard EpiSwitch template and no template controls.

Detection of similar epigenetic markers in blood and matching primary prostate tumours

Primary tumour samples were obtained from biopsies of selected patients from Imperial College Healthcare NHS Trust cohort (n = 5). The pulverized tissue samples were incubated in 0.125% collagenase at 37 °C with gentle agitation for 30 min. The resuspended cells (250 µl) were then centrifuged at 800g for 5 min at room temperature in a fixed arm centrifuge, supernatant removed, and the pellets resuspended in phosphate-buffered saline. Primary tumours and matching blood samples were analysed for the presence of the six-marker set for categories 3 vs 1 and 3 vs 2 at a fixed range of assay sensitivity (dilution factor 1:2). When matching PCR bands of the correct size were detected, a score of 1 was assigned, detection of no band was assigned a score of 0 (Table 2).

Table 2 Detection of similar epigenetic markers in blood and in matching primary prostate tumours at a fixed range of assay sensitivity

Results

Identification of the diagnostic markers

We have applied a stepwise diagnostic biomarker discovery process using EpiSwitch™ technology as described in methods. A customized CGH Agilent microarray (8 × 60 k) platform was designed to test technical and biological repeats for 14,241 potential chromosome conformations across 425 genetic loci (Additional file 1: Table S3) in eight PCa and eight cancer-free control samples (Figs. 1, 2). The presence or absence of each locus was defined by LIMMA linear modelling, subsequent binary filtering and cluster analysis. In the second evaluation stage, nested PCR was used for the 53 selected biomarkers further reducing them to 15 markers and finally to a five-marker signature (Fig. 2). This distinct chromosome conformation disease signature for PCa comprised of chromosomal interactions in five genomic loci: ETS proto-oncogene 1, transcription factor (ETS1), mitogen-activated protein kinase kinase kinase 14 (MAP3K14), solute carrier family 22 member 3 (SLC22A3) and caspase 2 (CASP2) (Table 1). The presence of an altered loop in each of the genes individually (except SLC22A3) was not statistically significant in predicting cancer. Only a combination of five of them provided required accuracy in the logistic model (please see the Stepwise Data Analysis section in Methods). The genomic locations of specific chromosomal loops in ETS1, MAP3K14, SLC22A3 and CASP2 genes in the chromosome conformation signature (Table 1) were mapped on their relative chromosomes in Fig. 3. The two genomic sites that corresponded to the junction of each chromosome conformation signature locus for ETS1, MAP3K14, SLC22A3 and CASP2 genes were mapped on chromosome 11 from 128,260,682 to 128,537,926; chromosome 17 from 43,303,603 to 43,432,282; chromosome 6 from 160,744,223 to 160,944,757 and chromosome 7 from 142,935,233 to 143,008,163, respectively, with the EpiSwitch™ sites marked with colour (Fig. 3a). Circos plots of ETS1, MAP3K14, SLC22A3 and CASP2 chromosome conformation signature markers showing the chromosomal loops are depicted in Fig. 3b.

Fig. 3
figure 3

Graphical representation of the genomic co-ordinates of the ETS1, MAP3K14, SLC22A3 and CASP2 chromosome conformation signature markers associated with PCa. a The Ensembl browser view of the ETS1, MAP3K14, SLC22A3 and CASP2 genes on chromosomes 11, 17, 6 and 7 with the EpiSwitch™ sites marked with green, pink, blue and red symbols, respectively. b Circos plots of ETS1 (green) MAP3K14 (pink), SLC22A3 (blue) and CASP2 (red) chromosome conformation signature markers showing the chromosomal loop

Principal component analysis for the five-markers was used to determine group variance and to identify potential outlier samples. This analysis was applied to 78 samples containing two groups. First group, 49 known samples (24 PCa and 25 cancer-free controls) combined with a second group of 29 samples including, 24 PCa samples and 5 cancer-free control samples demonstrating no outliers (Fig. 4). The final training set was built using 95 PCa and 96 cancer-free control samples and then tested on an independent blinded validation cohort of 20 samples (10 cancer-free controls and 10 PCa). The sensitivity and specificity for PCa detection using chromosomal interactions in five genomic loci were 80% (CI 44.39% to 97.48%) and 80% (CI 44.39% to 97.48%), respectively (Additional file 1: Table S4).

Fig. 4
figure 4

Principal component analysis for the five-markers applied to 78 samples containing two groups. First group, 49 known samples (24 PCa and 25 cancer-free controls (Cntrl)) combined with a second group of 29 samples including, 24 PCa samples and 5 healthy Cntrl samples. PCa samples in green and cancer-free Cntrl samples in red

Identification of the prognostic markers

To select epigenetic biomarkers that can stratify PCa, the samples from PCa patients categorised into risk group categories 1–3 (low (n = 42), intermediate (n = 18) and high (n = 25), respectively, Additional file 1: Table S2) were screened for statistically significant differences in conditional and stable profiles of genome architecture. EpiSwitch™ was used on blood samples in a three-step process to identify, evaluate, and validate statistically significant differences in chromosomal conformations between PCa patients at different stages of the disease (Fig. 2). For the first step, the array used covered 425 genetic loci, with testing probes for the total of 14,241 potential chromosomal conformations. Patients with high-risk PCa category 3 were compared to low-risk category 1 or intermediate-risk category 2. In total, 181 potential stratification marker leads for PCR evaluation were identified using enrichment statistics (Additional file 1: Table S5). The top 70 top markers were then taken to the next stage of PCR detection for further evaluation and stratification of high-risk category 3, vs low-risk category 1 patient samples. The 70 markers were reduced to 16 and finally a six-marker set for high category 3 vs low category 1 was established (Additional file 1: Table S6). The best markers were identified using Chi-square and then built into a classifier on a testing set of category 1 (n = 21) and category 3 (n = 19). An independent cohort of category 1 (n = 21) and category 3 (n = 6) which were not used for any marker reduction were then used for first round of blind validation. Similarly, a six-marker set was evaluated for high-risk category 3 vs intermediate-risk category 2 on a testing set of category 3 and category 2 including, 17 and 12 samples, respectively. An independent cohort of category 2 and category 3 (n = 6 in each group) which were not used for any marker reduction were then used for first round of blind validation.

For the last step, to further validate the chromosome conformation signature used to inform PCa prognosis, the six-marker set for high-risk category 3 vs low-risk category 1 was tested on a larger, more representative cohort. The original blind cohort was expanded to 67 samples, including 40 samples used in marker reduction (Additional file 1: Table S7). Similarly, the six-marker set for high-risk category 3 vs intermediate-risk category 2 was tested on a on a larger, more representative cohort. The original blind cohort was expanded to 43 samples (Additional file 1: Table S8).

A six-marker set for category 3 vs category 1 was established. This set contained bone morphogenetic protein 6 (BMP6), ETS transcription factor ERG (ERG), macrophage scavenger receptor 1 (MSR1), mucin 1 (MUC1), acetyl-CoA acetyltransferase 1 (ACAT1) and death-associated protein kinase 1 (DAPK1) genes (Additional file 1: Table S6). Six-biomarkers were identified for high-risk category 3 vs intermediate-risk category 2, including hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 2 (HSD3B2), vascular endothelial growth factor C (VEGFC), apoptotic peptidase activating factor 1 (APAF1), MUC1, ACAT1 and DAPK1. Notably, the last three-biomarkers (MUC1, ACAT1 and DAPK1) were common between categories 1 vs 3 and 3 vs 2 (Additional file 1: Table S6). Stratification of high-risk category 3 vs low-risk category 1 PCa using chromosomal interactions in six genomic loci showed sensitivity of 80% (CI 59.30% to 93.17%) and specificity of 93% (CI 80.52% to 98.50%) in the blind cohort of 67 samples (Additional file 1: Table S7). Similarly, the six-marker set for high-risk category 3 vs intermediate-risk category 2 was tested on a on a larger, more representative cohort of 43 samples demonstrating sensitivity of 84% (CI 63.92% to 95.46%), and specificity of 89% (CI 65.29% to 98.62%) (Additional file 1: Table S8).

Detection of similar epigenetic markers in blood and matching primary prostate tumours

Using five matching peripheral blood and primary tumour samples, we have compared the epigenetic markers identified in peripheral circulation (Additional file 1: Table S6) to the tumour tissue. Our results showed that a number of deregulation markers detected in the blood as part of stratifying signatures for category 1 vs 3 and category 2 vs 3 could be detected in the tumour tissue (Table 2). This demonstrates that the chromosome interactions that can be detected systemically could be detected under same conditions in the primary site of tumourigenesis.

Real-time quantitative PCR (qPCR) translation

The 16 statistically significant markers identified in the prognostic marker screen were translated to real-time qPCR assays. Due to the nature of 3C templates the sequence available for suitable primer and probe design for qPCR is constrained. This led to 6 of the markers failing to be translated to the qPCR assay. The assays were further assessed over the nested PCR assays by screening a 40-sample subset of the prognostic sample cohort. The cohort consisted of 10 cancer-free controls, 10 category 1, 10 category 2 and 10 category 3 samples. A new classifier was produced using the qPCR data and tested with resampling for the ability to distinguish between the cancer-free controls and PCa based on 6 markers, demonstrating average sensitivity of an 91% and specificity of 89.3% based on 500 randomisations for cross validation and a 66/34% train/test split.

Discussion

Timely diagnosis of PCa is crucial to reducing mortality. The European randomised study of screening for PCa has shown significant reduction in PCa mortality in men who underwent routine PSA screening [38, 39]. This notion was, however, not supported by the results of prostate, lung, colorectal, and ovarian (PLCO) cancer screening trial in the USA [40]. Total population screening leads to overdiagnosis of clinically insignificant disease and new less invasive tests capable of discriminating low- from high-risk disease are urgently required.

Our epigenetic analysis approach provides a potentially powerful means to address this need. The binary nature of the test (the chromosomal loop is either present or not) and the enormous combinatorial power (> 1010 combinations are possible with ~ 50,000 loops screened) may allow creating signatures that accurately fit clinically well-defined criteria. In PCa that would be discerning low-risk vs high-risk disease or identifying small but aggressive tumours and determining most appropriate therapeutic options. In addition, epigenetic changes are known to manifest early in tumourigenesis, making them useful for both diagnosis and prognosis [41].

In this study, we identified and validated chromosome conformations as distinctive biomarkers for a non-invasive blood-based epigenetic signature for PCa. Our data demonstrate the presence of stable chromatin loops in the loci of ETS1, MAP3K14, SLC22A3 and CASP2 genes present only in PCa patients (Table 1). Validation of these markers in an independent set of 20 blinded samples showed 80% sensitivity and 80% specificity (Additional file 1: Table S4), which is remarkable for a PCa blood test. Interestingly, the expression of some of these genes has already been linked to cancer pathophysiology. ETS1 is a member of ETS transcription factor family. ETS1‐overexpressing prostate tumours are associated with increased cell migration, invasion and induction of epithelial‐to‐mesenchymal transition [42]. MAP3K14 (also known as nuclear factor-kappa-beta (NF-kβ)-inducing kinase (NIK)) is a member of MAP3K group (or MEKK). Physiologically, MAP3K14/NIK can activate noncanonical NF-kβ signalling and induce canonical NF-kβ signalling, particularly when MAP3K14/NIK is overexpressed [43]. A novel role for MAP3K14/NIK in regulating mitochondrial dynamics to promote tumour cell invasion has been described [44]. SLC22A3 (also known as organic cation transporter 3) is a member of SLC group of membrane transport proteins. SLC22A3 expression is associated with PCa progression [45]. CASP2 is a member of caspase activation and recruitment domains group. Physiologically, CASP2 can act as an endogenous repressor of autophagy [46]. Two of the identified genes (SLC22A3 and CASP2) were previously shown to be inversely correlated with cancer progression [45, 47, 48]. Importantly, the presence of the chromatin loop can have indeterminate effect on gene expression.

To screen for PCa prognostic markers we performed the EpiSwitch™ custom array to analyse competitive hybridization of DNA from peripheral blood from patients with low-risk PCa (category 1) and high risk PCa (category 3). Six-marker set was identified for high-risk category 3 vs low-risk category 1, including BMP6, ERG, MSR1, MUC1, ACAT1 and DAPK1. Six-biomarkers were identified for high-risk category 3 vs intermediate-risk category 2, including HSD3B2, VEGFC, APAF1, MUC1, ACAT1 and DAPK1. Three of these biomarkers (MUC1, ACAT1 and DAPK1) were shared between these sets. Our data show high concordance between chromosomal conformations in the primary tumour and in the blood of matched PCa patients at stages 1 and 3 (Table 2). The prognostic significance and diagnostic value of some of these genes have previously been suggested [49,50,51,52]. BMP6 plays an important role in PCa bone metastasis [50]. In addition to ETS1, ERG is another member of the ETS family of transcription factors. Overwhelming evidence, reviewed in [49], suggested that ERG is implicated in several processes relevant to PCa progression including metastasis, epithelial–mesenchymal transition, epigenetic reprogramming, and inflammation [49]. MSR1 may confer a moderate risk for PCa [53]. MUC1 is a membrane-bound glycoprotein that belongs to the mucin family. MUC1 high expression in advanced PCa is associated with adverse clinicopathological tumour features and poor outcomes [51]. ACAT1 expression is elevated in high-grade and advanced PCa and acts as an indicator of reduced biochemical recurrence-free survival [52]. DAPK1 could function either as a tumour suppressor or as an oncogenic molecule in different cellular context [54]. HSD3B2 plays a crucial role in steroid hormone biosynthesis and it is up-regulated in a relevant fraction of PCa that are characterized by an adverse tumour phenotype, increased androgen receptor signalling and early biochemical recurrence. [55]. VEGFC is a member of VEGF family and its increased expression is associated with lymph node metastasis in PCa specimens [56]. In a comprehensive biochemical approach, APAF1 has been described as the core of the apoptosome [57].

Despite the identification of these loci, the mechanism of cancer-related epigenetic changes in PBMCs remains unidentified. The interaction, however, could be detected systemically and under same conditions in the primary site of tumourigenesis (Table 2). It is therefore assumed that the acquired changes must be directed by an external factor; presumably generated by the tumour cells. It is known that a significant proportion of chromosomal conformations are controlled by non-coding RNAs, which regulate the tumour-specific conformations [58, 59]. Tumour cells have been shown to secrete non-coding RNAs that are endocytosed by neighbouring or circulating cells and may change their chromosomal conformations [59, 60], and are possible regulators in this case.

RNA detection as a biomarker remains highly challenging (low stability, background drift, continuous variable for statistical stratification analysis). Chromosome conformation signatures offer well recognized stable binary advantages for the biomarker targeting use [25], specifically when tested in the nuclei [31], since the circulating DNA present in plasma does not retain 3D conformational topological structures present in the intact cellular nuclei. It is important to mention, that looking at one genetic locus does not equate to looking at one marker, as there may be multiple chromosome conformations present, representing parallel pathways of epigenetic regulation over the locus of interest.

Other technologies for plasma-based cancer detection such as using plasma cell-free DNA (cfDNA) methylomes [61] were recently introduced. The validity of this assay was tested to identify patients with renal cell carcinoma using urine cfDNA with area under the ROC curve (AUROC) of 0.86 [62]. It is worth noting that cfDNA is capturing post-apoptotic and necrotic passive distribution of free DNA, with significant variations, while EpiSwitch is measuring 3D genomic profiling in intact cells, capturing systemic surrogate readouts, which, for epigenetic modalities, have been shown to contain synchronized modulation at specific genetic loci concordant with primary sites of deregulation [63]. Systemic surrogate signatures at selective loci at the level of 3D genomics are sustained through exosome signalling and are not restricted to oncology [36].

One of the key challenges in the present clinical practice of PCa diagnosis is the time it takes to make a definitive diagnosis. So far, there is no single, definitive test for PCa. High levels of PSA will set the patient on a long journey of uncertainty where he will undergo a MRI scan followed by biopsy, if needed. Although a biopsy is more reliable than a PSA test, it is an invasive procedure at which targeting the cancer remains a significant issue. The five-set biomarker panel described here is based on a relatively inexpensive and well-established molecular biology technique (PCR). The samples are based on biofluid, which is simple to collect and provides clinicians with rapidly available clinical readouts within few hours. This in turn, offers a substantial time and cost savings and aids an informative diagnostic decision with the promise to help fill the gap in the current protocols for assertive diagnosis of PCa.

Predicting the risk of PCa is pivotal for making an informed decision on treatment options. Five-year survival rate in the low risk group is more than 95% and most men would benefit from less invasive therapy. Currently, PCa risk stratification is based on combined assessment of circulating PSA, tumour grade (from biopsy) and tumour stage (from imaging findings). The ability to derive similar information using a simple blood test would allow significant reduction in costs and would speed up the diagnostic process. Of particular importance in PCa treatment is identifying the few tumours that initially present as low-risk, but then progress to high-risk. This subset would therefore benefit from a quicker and more-radical intervention.

This study has several limitations including smaller number of patients in diagnostic validation cohort, unavailability of other clinical indices like PHI and 4 K (which are not part of the standard of care in the UK) and no follow-up (due to double anonymity). It is worth noting that the aim of this study was to demonstrate the potential capability of our general method and the new paradigm, rather than introduce a new clinical test for PCa that would require larger specialised studies in different settings. The exact mechanism of concordance between leukocyte 3D chromosome conformations and those in primary tumours remains unclear. We hypothesise that it is mediated by the horizontal transfer of epigenetic information, likely through non-coding RNAs, however, dedicated studies are required to investigate this further.

Conclusions

In conclusion, we have identified subsets of chromosomal conformations in patients’ PBMCs that are strongly indicative of PCa presence and prognosis. These signatures have significant potential for the development of quick diagnostic and prognostic blood tests for PCa and significantly exceed the specificity of the currently used PSA test.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

3C:

Chromosome conformation capture

ACAT1:

Acetyl-CoA acetyltransferase 1

APAF1:

Apoptotic peptidase activating factor 1

BMP6:

Morphogenetic protein 6

CASP2:

Caspase 2

DAPK1:

Death-associated protein kinase 1

ERG:

ETS transcription factor ERG

ETS1:

ETS proto-oncogene 1, transcription factor

GWAS:

Genome-wide association studies

HSD3B2:

Hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 2

MAP3K14:

Mitogen-activated protein kinase kinase kinase 14

MRI:

Magnetic resonance imaging

MSR1:

Macrophage scavenger receptor 1

MUC1:

Mucin 1

NF-kβ:

Nuclear factor-kappa-beta

NIK:

Nuclear factor-kappa-beta-inducing kinase

NCCN:

National Comprehensive Cancer Network

PBMCs:

Peripheral blood mononuclear cells

PCa:

Prostate cancer

PCR:

Polymerase chain reactions

PSA:

Prostate specific antigen

qPCR:

Quantitative PCR

SLC22A3:

Solute carrier family 22 member 3

VEGFC:

Vascular endothelial growth factor C

References

  1. ACS. American Cancer Society: http://www.cancer.org/Research/CancerFactsFigures/index; 2020. http://www.cancer.org/cancer/prostatecancer/overviewguide/prostate-cancer-overview-key-statistics.

  2. CRUK. Cancer Research UK; 2020. http://www.cancerresearchuk.org/home/.

  3. Sakr WA, Haas GP, Cassin BF, Pontes JE, Crissman JD. The frequency of carcinoma and intraepithelial neoplasia of the prostate in young male patients. J Urol. 1993;150(2 Pt 1):379–85.

    Article  CAS  PubMed  Google Scholar 

  4. Sanchez-Chapado M, Olmedilla G, Cabeza M, Donat E, Ruiz A. Prevalence of prostate cancer and prostatic intraepithelial neoplasia in Caucasian Mediterranean males: an autopsy study. Prostate. 2003;54(3):238–47.

    Article  PubMed  Google Scholar 

  5. Humphrey PA, Moch H, Cubilla AL, Ulbright TM, Reuter VE. The 2016 WHO classification of tumours of the urinary system and male genital organs-Part B: prostate and bladder tumours. Eur Urol. 2016;70(1):106–19.

    Article  PubMed  Google Scholar 

  6. EAU. Prostate cancer guidelines. 2020.

  7. Ankerst DP, Thompson IM. Sensitivity and specificity of prostate-specific antigen for prostate cancer detection with high rates of biopsy verification. Archivio italiano di urologia, andrologia. 2006;78(4):125–9.

    PubMed  Google Scholar 

  8. Thompson IM, Pauler DK, Goodman PJ, Tangen CM, Lucia MS, Parnes HL, et al. Prevalence of prostate cancer among men with a prostate-specific antigen level < or =4.0 ng per milliliter. N Engl J Med. 2004;350(22):2239–46.

    Article  CAS  PubMed  Google Scholar 

  9. Falchook AD, Martin NE, Basak R, Smith AB, Milowsky MI, Chen RC. Stage at presentation and survival outcomes of patients with Gleason 8–10 prostate cancer and low prostate-specific antigen. Urologic Oncol. 2016;34(3):119 e19-126.

    Article  Google Scholar 

  10. Nelius T, Filleur S. PSA surge/flare-up in patients with castration-refractory prostate cancer during the initial phase of chemotherapy. Prostate. 2009;69(16):1802–7 (Epub 2009/08/14. eng).

    Article  CAS  PubMed  Google Scholar 

  11. Punnen S, Pavan N, Parekh DJ. Finding the wolf in sheep’s clothing: the 4Kscore is a novel blood test that can accurately identify the risk of aggressive prostate cancer. Rev Urol. 2015;17(1):3–13.

    PubMed  PubMed Central  Google Scholar 

  12. Filella X, Fernandez-Galan E, Fernandez Bonifacio R, Foj L. Emerging biomarkers in the diagnosis of prostate cancer. Pharmacogenom Personal Med. 2018;11:83–94.

    Article  CAS  Google Scholar 

  13. NICE. Prostate Cancer: Diagnosis and management Evidence reviews for following-up people at risk of prostate cancer 2019.

  14. Eeles RA, Olama AA, Benlloch S, Saunders EJ, Leongamornlert DA, Tymrakiewicz M, et al. Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array. Nat Genet. 2013;45(4):385–91.

    Article  CAS  PubMed  Google Scholar 

  15. Reynolds MA. Molecular alterations in prostate cancer. Cancer Lett. 2008;271(1):13–24.

    Article  CAS  PubMed  Google Scholar 

  16. Amin Al, Olama A, Benlloch S, Antoniou AC, Giles GG, Severi G, Neal DE, et al. Risk analysis of prostate cancer in PRACTICAL, a multinational consortium, using 25 known prostate cancer susceptibility loci. Cancer Epidemiol Biomarkers Prev. 2015; 24(7): 1121–9.

  17. Frazer KA, Murray SS, Schork NJ, Topol EJ. Human genetic variation and its contribution to complex traits. Nat Rev Genet. 2009;10(4):241–51.

    Article  CAS  PubMed  Google Scholar 

  18. Jerónimo C, Bastian PJ, Bjartell A, Carbone GM, Catto JW, Clark SJ, et al. Epigenetics in prostate cancer: biologic and clinical relevance. European Urol. 2011;60(4):753–66 (Epub 2011/07/02. eng).

    Article  CAS  Google Scholar 

  19. Sugiura M, Sato H, Kanesaka M, Imamura Y, Sakamoto S, Ichikawa T, et al. Epigenetic modifications in prostate cancer. Int J Urol. 2020;14:100915.

    Google Scholar 

  20. Schulz WA, Elo JP, Florl AR, Pennanen S, Santourlidis S, Engers R, et al. Genomewide DNA hypomethylation is associated with alterations on chromosome 8 in prostate carcinoma. Genes Chromosomes Cancer. 2002;35(1):58–65 (Epub 2002/08/31. eng).

    Article  CAS  PubMed  Google Scholar 

  21. Yegnasubramanian S, Kowalski J, Gonzalgo ML, Zahurak M, Piantadosi S, Walsh PC, et al. Hypermethylation of CpG islands in primary and metastatic human prostate cancer. Cancer Res. 2004;64(6):1975–86 (Epub 2004/03/18. eng).

    Article  CAS  PubMed  Google Scholar 

  22. Chen H, Tu SW, Hsieh JT. Down-regulation of human DAB2IP gene expression mediated by polycomb Ezh2 complex and histone deacetylase in prostate cancer. J Biol Chem. 2005;280(23):22437–44 (Epub 2005/04/09. eng).

    Article  CAS  PubMed  Google Scholar 

  23. Iljin K, Wolf M, Edgren H, Gupta S, Kilpinen S, Skotheim RI, et al. TMPRSS2 fusions with oncogenic ETS factors in prostate cancer involve unbalanced genomic rearrangements and are associated with HDAC1 and epigenetic reprogramming. Cancer Res. 2006;66(21):10242–6 (Epub 2006/11/03. eng).

    Article  CAS  PubMed  Google Scholar 

  24. Mukhopadhyay S, Ramadass AS, Akoulitchev A, Gordon S. Formation of distinct chromatin conformation signatures epigenetically regulate macrophage activation. Int Immunopharmacol. 2014;18(1):7–11 (PubMed PMID: 24211766).

    Article  CAS  PubMed  Google Scholar 

  25. Crutchley JL, Wang XQ, Ferraiuolo MA, Dostie J. Chromatin conformation signatures: ideal human disease biomarkers? Biomark Med. 2010;4(4):611–29.

    Article  CAS  PubMed  Google Scholar 

  26. Pomerantz MM, Ahmadiyeh N, Jia L, Herman P, Verzi MP, Doddapaneni H, et al. The 8q24 cancer risk variant rs6983267 shows long-range interaction with MYC in colorectal cancer. Nat Genet. 2009;41(8):882–4 (Epub 2009/06/30. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xu Z, Wei G, Chepelev I, Zhao K, Felsenfeld G. Mapping of INS promoter interactions reveals its role in long-range regulation of SYT8 transcription. Nat Struct Mol Biol. 2011;18(3):372–8 (Epub 2011/02/22. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Yao X, Zha W, Song W, He H, Huang M, Jazrawi E, et al. Coordinated regulation of IL-4 and IL-13 expression in human T cells: 3C analysis for DNA looping. Biochem Biophys Res Commun. 2012;417(3):996–1001 (Epub 2012/01/10. eng).

    Article  CAS  PubMed  Google Scholar 

  29. Cortese R, Kwan A, Lalonde E, Bryzgunova O, Bondar A, Wu Y, et al. Epigenetic markers of prostate cancer in plasma circulating DNA. Hum Mol Genet. 2012;21(16):3619–31.

    Article  CAS  PubMed  Google Scholar 

  30. Gao P, Xia JH, Sipeky C, Dong XM, Zhang Q, Yang Y, et al. Biology and clinical implications of the 19q13 aggressive prostate cancer susceptibility locus. Cell. 2018;174(3):576–8918.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Dekker J, Rippe K, Dekker M, Kleckner N. Capturing chromosome conformation. Science. 2002;295(5558):1306–11.

    Article  CAS  PubMed  Google Scholar 

  32. Babu D, Fullwood MJ. 3D genome organization in health and disease: emerging opportunities in cancer translational medicine. Nucleus. 2015;6(5):382–93 (Epub 2015/11/11. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Bastonini E, Jeznach M, Field M, Juszczyk K, Corfield E, Dezfouli M, et al. Chromatin barcodes as biomarkers for melanoma. Pigment Cell Melanoma Res. 2014;27(5):788–800.

    Article  CAS  PubMed  Google Scholar 

  34. Jakub JW, Grotz TE, Jordan P, Hunter E, Pittelkow M, Ramadass A, et al. A pilot study of chromosomal aberrations and epigenetic changes in peripheral blood samples to identify patients with melanoma. Melanoma Res. 2015;25(5):406–11.

    Article  PubMed  Google Scholar 

  35. Carini C, Hunter E, Ramadass AS, Green J, Akoulitchev A, McInnes IB, et al. Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis. J Transl Med. 2018;16(1):18 (Epub 2018/01/31. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Salter M, Corfield E, Ramadass A, Grand F, Green J, Westra J, et al. Initial identification of a blood-based chromosome conformation signature for aiding in the diagnosis of amyotrophic lateral sclerosis. EBioMedicine. 2018;33:169–84.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Carini C, Hunter E, Scottish Early Rheumatoid Arthritis Inception cohort I, Ramadass AS, Green J, Akoulitchev A, et al. Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis. J Translat Med. 2018;16(1):18.

  38. Schroder FH, Hugosson J, Roobol MJ, Tammela TL, Ciatto S, Nelen V, et al. Screening and prostate-cancer mortality in a randomized European study. N Engl J Med. 2009;360(13):1320–8.

    Article  PubMed  Google Scholar 

  39. Hugosson J, Roobol MJ, Mansson M, Tammela TLJ, Zappa M, Nelen V, et al. A 16-yr follow-up of the european randomized study of screening for prostate cancer. Eur Urol. 2019;76(1):43–51.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Andriole GL, Crawford ED, Grubb RL 3rd, Buys SS, Chia D, Church TR, et al. Mortality results from a randomized prostate-cancer screening trial. N Engl J Med. 2009;360(13):1310–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Sun YW, Chen KM, Imamura Kawasawa Y, Salzberg AC, Cooper TK, Caruso C, et al. Hypomethylated Fgf3 is a potential biomarker for early detection of oral cancer in mice treated with the tobacco carcinogen dibenzo[def, p]chrysene. PLoS ONE. 2017;12(10):0186873.

    Google Scholar 

  42. Rodgers JJ, McClure R, Epis MR, Cohen RJ, Leedman PJ, Harvey JM, et al. ETS1 induces transforming growth factor beta signaling and promotes epithelial-to-mesenchymal transition in prostate cancer cells. J Cell Biochem. 2019;120(1):848–60 (Epub 2018/08/31. eng).

    Article  CAS  PubMed  Google Scholar 

  43. Odqvist L, Sanchez-Beato M, Montes-Moreno S, Martin-Sanchez E, Pajares R, Sanchez-Verde L, et al. NIK controls classical and alternative NF-kappaB activation and is necessary for the survival of human T-cell lymphoma cells. Clin Cancer Res. 2013;19(9):2319–30 (Epub 2013/03/29. eng).

    Article  CAS  PubMed  Google Scholar 

  44. Jung JU, Ravi S, Lee DW, McFadden K, Kamradt ML, Toussaint LG, et al. NIK/MAP3K14 regulates mitochondrial dynamics and trafficking to promote cell invasion. Curr Biol. 2016;26(24):3288–302 (Epub 2016/11/28. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Tomlins SA, Mehra R, Rhodes DR, Cao X, Wang L, Dhanasekaran SM, et al. Integrative molecular concept modeling of prostate cancer progression. Nat Gene. 2007;39(1):41–51 (Epub 2006/12/19. eng).

    Article  CAS  Google Scholar 

  46. Tiwari M, Sharma LK, Vanegas D, Callaway DA, Bai Y, Lechleiter JD, et al. A nonapoptotic role for CASP2/caspase 2: modulation of autophagy. Autophagy. 2014;10(6):1054–70 (Epub 2014/06/01. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Lopez-Garcia C, Sansregret L, Domingo E, McGranahan N, Hobor S, Birkbak NJ, et al. BCL9L dysfunction impairs caspase-2 expression permitting aneuploidy tolerance in colorectal cancer. Cancer Cell. 2017;31(1):79–93 (Epub 2017/01/11. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Puccini J, Dorstyn L, Kumar S. Caspase-2 as a tumour suppressor. Cell Death Different. 2013;20(9):1133–9 (Epub 2013/07/03. eng).

    Article  CAS  Google Scholar 

  49. Adamo P, Ladomery MR. The oncogene ERG: a key factor in prostate cancer. Oncogene. 2016;35(4):403–14.

    Article  CAS  PubMed  Google Scholar 

  50. Yuen H-F, Chan Y-P, Cheung W-L, Wong Y-C, Wang X, Chan K-W. The prognostic significance of BMP-6 signaling in prostate cancer. Modern Pathol. 2008;21(12):1436–43.

    Article  CAS  Google Scholar 

  51. Genitsch V, Zlobec I, Thalmann GN, Fleischmann A. MUC1 is upregulated in advanced prostate cancer and is an independent prognostic factor. Prostate Cancer Prostatic Dis. 2016;19(3):242–7 (Epub 2016/05/12. eng).

    Article  CAS  PubMed  Google Scholar 

  52. Saraon P, Trudel D, Kron K, Dmitromanolakis A, Trachtenberg J, Bapat B, et al. Evaluation and prognostic significance of ACAT1 as a marker of prostate cancer progression. Prostate. 2014;74(4):372–80 (Epub 2013/12/07. eng).

    Article  CAS  PubMed  Google Scholar 

  53. Sun J, Hsu FC, Turner AR, Zheng SL, Chang BL, Liu W, et al. Meta-analysis of association of rare mutations and common sequence variants in the MSR1 gene and prostate cancer risk. Prostate. 2006;66(7):728–37 (Epub 2006/01/21. eng).

    Article  CAS  PubMed  Google Scholar 

  54. Steinmann S, Scheibe K, Erlenbach-Wuensch K, Neufert C, Schneider-Stock R. Death-associated protein kinase: a molecule with functional antagonistic duality and a potential role in inflammatory bowel disease. Int J Oncol. 2015;47(1):5–15 (Epub 2015/05/13. eng).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Neubauer E, Latif M, Krause J, Heumann A, Armbrust M, Luehr C, et al. Up regulation of the steroid hormone synthesis regulator HSD3B2 is linked to early PSA recurrence in prostate cancer. Exp Mol Pathol. 2018;105(1):50–6.

    Article  CAS  PubMed  Google Scholar 

  56. Jennbacken K, Vallbo C, Wang W, Damber JE. Expression of vascular endothelial growth factor C (VEGF-C) and VEGF receptor-3 in human prostate cancer is associated with regional lymph node metastasis. Prostate. 2005;65(2):110–6 (Epub 2005/05/10. eng).

    Article  CAS  PubMed  Google Scholar 

  57. Cecconi F, Gruss P. Apaf1 in developmental apoptosis and cancer: how many ways to die. Cell Mol Life Sci. 2001;58(11):1688–97 (Epub 2001/11/15. eng).

    Article  CAS  PubMed  Google Scholar 

  58. Beckedorff FC, Amaral MS, Deocesano-Pereira C, Verjovski-Almeida S. Long non-coding RNAs and their implications in cancer epigenetics. Biosci Rep. 2013;33(4):00061 (Epub 2013/07/24. eng).

    Article  CAS  Google Scholar 

  59. Lai F, Orom UA, Cesaroni M, Beringer M, Taatjes DJ, Blobel GA, et al. Activating RNAs associate with Mediator to enhance chromatin architecture and transcription. Nature. 2013;494(7438):497–501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Kosaka N, Ochiya T. Unraveling the Mystery of Cancer by Secretory microRNA: Horizontal microRNA Transfer between Living Cells. Front Genet. 2011;2:97 (Epub 2012/02/04. eng).

    PubMed  Google Scholar 

  61. Shen SY, Singhania R, Fehringer G, Chakravarthy A, Roehrl MHA, Chadwick D, et al. Sensitive tumour detection and classification using plasma cell-free DNA methylomes. Nature. 2018;563(7732):579–83.

    Article  CAS  PubMed  Google Scholar 

  62. Nuzzo PV, Berchuck JE, Korthauer K, Spisak S, Nassar AH, Abou Alaiwi S, et al. Detection of renal cell carcinoma using plasma and urine cell-free DNA methylomes. Nat Med. 2020;26(7):1041–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Feinberg AP. The key role of epigenetics in human disease prevention and mitigation. N Engl J Med. 2018;378(14):1323–34.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank members of Oxford BioDynamics; Benjamin Foulkes, Chloe Bird, Diana Jaramillo Mahecha, Emily Corfield, Warren Elvidge, Ryan Powell for the laboratory work.

Funding

This work was funded by Oxford BioDynamics.

Author information

Authors and Affiliations

Authors

Contributions

EH, AR, AA, and DP conceived the study. EH, MS, AR, JG, and AA planned, performed and reviewed experiments and analysed the data. HA, RM, MW, and DP enrolled patients and provided samples and clinical data. BMS, WW provided critical assessment of the manuscript. All authors participated in writing and editing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dmitri Pchejetski.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the University of East Anglia Human Tissue committee and all participants gave informed consent before enrolment.

Consent for publication

Not applicable.

Competing interests

Ewan Hunter, Matthew Salter, Aroul Ramadass, Willem Westra, Jayne Green and Alexandre Akoulitchev are employees of Oxford BioDynamics. A. Akoulitchev and A. Ramadass are company directors. Oxford BioDynamics holds patents on the EpiSwitch™ technology. The remaining authors have no conflict of 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: Table S1.

Clinical characteristics of the prostate cancer patients participated in the study. Table S2. Prostate cancer risk group categories. Table S3. List of 425 prostate cancer-related genomic loci tested in the initial array. Table S4. Comparison of pathology and EpiSwitch™ results. Table S5. Markers for prognostic array stratifications category 1 vs 3 and category 2 vs 3. Top 181 markers produced from the prognostic array. Table S6. Markers for high-risk category 3 vs low-risk category 1 and for high-risk category 3 vs intermediate-risk category 2. Table S7. Comparison of pathology and EpiSwitch™ results for category 3 vs 1 classifier. Table S8. Comparison of pathology and EpiSwitch™ results for category 3 vs 2 classifier.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Alshaker, H., Mills, R., Hunter, E. et al. Chromatin conformation changes in peripheral blood can detect prostate cancer and stratify disease risk groups. J Transl Med 19, 46 (2021). https://doi.org/10.1186/s12967-021-02710-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-021-02710-y

Keywords