1. Introduction
Detection of potentially small populations of non-haematopoetic circulating tumour cells (CTCs) within whole blood in cancer patients represents a significant technical challenge [
1,
2]. Isolating CTCs is generally attempted by using physical features such as size and enrichment techniques to confirm the identity of epithelial rather than haematopoietic cells. This has relied on surface proteins such as the Epithelial Cell Adhesion Molecule (EpCAM) to confirm/refute their origin as cancer cells.
The most widely used method for CTC analysis is the CELLSEARCH
® CTC Test, FDA approved in breast, colorectal and prostate cancer, where identification of CTCs broadly rely on EpCAM, despite evidence of considerable heterogeneity of EpCAM expression during the cell cycle, in different cell lines and a lack of EpCAM in cells undergoing epithelial-mesenchymal transition (EMT) [
3,
4,
5]. The CellSearch test has been explored in patients with ovarian cancer in multiple publications with very low yields (e.g., single CTCs in 7.5 mL of whole blood, which may contain 75 million leucocytes and 50 billion erythrocytes) and is thus of limited utility [
6] even in patients with recently diagnosed advanced epithelial ovarian cancer (aEOC) prior to any treatment. Most groups have used EpCAM as the common key target/antigen for the identification and confirmation of epithelial origin, with or without numerous other epithelial markers such as MUC1, MUC16, CK8, CK18, CK19, CA125, DPP4, ERCC1, HER2, MOC-31, PPIC. In clinical practice, immunohistochemistry assays on tissue likely to be of ovarian/fallopian tube/primary peritoneal origin utilise different antigens such as CK7, CK20 and WT1 to confirm the diagnosis of high-grade serous cancer prior to definitive treatment. We sought to explore the use of these antigens to identify CCs in aEOC patients.
The primary aim of this study was to examine whole blood of patients with aEOC using ImageStream™ (Amnis), a high-definition circulating cell (HD-CC) imaging flow cytometry system, without using surface protein-based enrichment techniques. This technology allows multiple antibody applications and presents CCs in sufficient definition to satisfy diagnostic pathology image quality requirements [
7,
8]. Clinical immunohistochemistry cell surface and nuclear immunohistochemical antibody markers, initially pan-cytokeratin (PCK; AE1/AE3), CD45, and DRAQ5™ were used to identify and characterise non-haematopoetic and haematopoetic cells. Subsequently WT1 antibody was used to specifically identify cells of high grade serous (ovarian, fallopian tube/primary peritoneal) origin. Correlations were made with serum CA125 levels during treatment and CC levels in 4 patients. We have expanded on our observations, to include further interrogation of some liquid biopsies by means of gene expression from whole blood.
2. Materials and Methods
2.1. Cell Culture
SKOV3 and MDAH-2774 ovarian cancer cells (ATCC) were grown in DMEM (Dulbecco’s Modified Eagle’s Medium, Gibco, Paisley, UK) supplemented with 10% FBS (fetal bovine serum, Gibco), 1% penicillin/streptomycin (Gibco) and 1% L-glutamine (Gibco); at 37 °C/5% CO2. PEO1 cells were grown in RPMI-1640, 10% FBS, 1% penicillin/streptomycin at 37 °C/5% CO2.
Healthy volunteer (control) whole blood (1mL) was spiked individually with 2 cell lines, SKOV3 and MDAH-2774, to mimic liquid biopsies that were processed as before and analysed using the enumeration features of the ImageStream™ based on size criteria and staining for both EpCAM and pan-cytokeratin (AE1/AE3) antibodies. This procedure was repeated with the introduction of different quantities of SKOV3 and MDAH-2774 cells, to assess the efficiency of cell retrieval and loss when using EpCAM and AE1/AE3 markers. The cell suspension was resuspended in different volumes to give a variety of cell concentrations between 20–200,000 cells/mL and spiked in 1mL blood from a donor.
2.2. Patients
Patients and healthy volunteers included in this report were enrolled in the prospective CICATRIx clinical study where blood samples are collected from healthy volunteers and patients with advanced cancer attending Mount Vernon Cancer Centre (East and North Hertfordshire NHS Trust) for exploratory biomarkers. All volunteers and patients provided written informed consent for participation in the study and for use of their donated tissue (patients only) and blood specimens.
The CICATRIx study: Sample collection study to explore circulating tumour cells, cell free DNA and leucocytes with ImageStream analysis in patients with various cancers. Protocol number RD2016-08 was approved by the West Midlands–South Birmingham Ethics Committee (reference 16/WM/0196).
TRANS-METROBIBF: Sample collection study to investigate the predictive role of Circulating Tumour Cells (CTCs) and their gene expression in relation to outcome in multiply-relapsed ovarian cancer patients. Protocol number: RD2013-01 was approved by the South East Coast–Surrey Ethics Committee (reference 14/LO/0792).
The studies were performed in accordance with the Declaration of Helsinki.
5–10 mL of whole blood was taken from patients at regular intervals with a tissue diagnosis of high-grade serous ovarian cancer (HGSOC), in one of the three groups:
(1) Relapse—patients whose aEOC had relapsed following prior remission, and who required further chemotherapy treatment.
(2) NACT—patients whose primary treatment at diagnosis was Neo-Adjuvant Chemotherapy.
(3) post-PDS—patients who had had Primary Debulking Surgery as their first treatment for aEOC.
39 healthy female volunteers donated 5–10 mL whole blood, which was also screened for cells that would stain with pan-cytokeratin and WT1 antibodies. Whole blood samples were assayed within 6 days of venesection.
2.3. Clinical Cohorts
Three cohorts of patients with aEOC were included in this exploration of CCs (
Figure 1):
(1) In the Relapse cohort, patient samples were taken prior to the first infusion of the course of chemotherapy and at the time of subsequent infusions. This group included 5 patients who participated in the METRO-BIBF clinical trial, a phase II, randomised, placebo-controlled, multicentre, feasibility study of low dose (metronomic) cyclophosphamide with and without nintedanib (BIBF 1120)/placebo (all patients received oral cyclophosphamide). Bloods tests were performed every 6 weeks at clinic review in the METRO-BIBF patients and their samples assayed for pan-cytokeratin, DRAQ5™ and CD45, but not WT1 expression. All other relapsed patients had WT1 expression examined.
(2) In the NACT cohort, patient samples were taken following confirmed histological diagnosis and prior to the start of their chemotherapy and at each subsequent treatment. Many of these patients underwent interval surgery following 3–4 chemotherapy cycles.
(3) In the post-PDS cohort, patient samples were taken after surgery and prior to any adjuvant chemotherapy and at each subsequent chemotherapy cycle.
Initially, a comparative assessment was made of a single sample, 30 mL from each of four patients, decanted into 6–7 individual tubes of different types: EDTA, Cell-Free DNA Blood Collection Tube (Streck), PAXgene Blood DNA Tube (Qiagen), and Cell-Free DNA Collection Tube (Roche) to establish the best vehicle for cell preservation prior to assay. Samples were analysed within 4 h of venesection for Day 1 and then daily for a further 6 days to assess the preservation ability of each type of collection tube. Subsequent patient and “control” samples from volunteers were collected in Roche tubes and assayed within 6 days of venesection.
2.4. CA125 Measurements
Serum CA125 assays were done for all aEOC patients, as part of standard practice, i.e., approximately each time they had chemotherapy treatment.
2.5. Preparing Blood Samples for Imagestream™
One ml of whole blood from each patient was decanted from Roche tubes into a 15 mL Falcon tube and mixed with 9 mL of red blood cell (RBC) lysis buffer (G Biosciences), inverted 10 times and incubated for 10 min with gentle agitation. The solution was then centrifuged at 2500 RPM for 10 min, the supernatant removed and a further 3 mL of RBC lysis buffer added to resuspend the pellet, before being incubated for 10 min at room temperature with gentle agitation. The solution was spun for a further 10 min at 2500 RPM, and the supernatant was aspirated. The pellet was then washed in 1.5 mL of PBS, and the resulting suspension was moved to a 1.5 mL microcentrifuge tube and spun at 3600 RPM for 3 min. Samples were then fixed immediately.
2.6. Fixing Cells
Cell pellets were transferred to a 1.5 mL microcentrifuge tube and resuspended in 1mL of ice cold 4% paraformaldehyde (PFA, Sigma, Gillingham, UK) for 7 min. The cell suspension was centrifuged for 2 min at 3600 RPM and the PFA was removed. Cells were washed in prewarmed PBS twice, centrifuging for 2 min at 3600 RPM between each wash. After the cells were fixed and washed with PBS, the pellet was suspended in 0.5% Triton X in PBS (v/v) and incubated on ice for 10 min, followed by centrifugation at 3600 RPM for 3 min to enable membrane permeabilisation for cytoplasmic/nuclear staining.
2.7. Staining Cells
Cells were incubated in blocking buffer (10% fetal bovine serum, Gibco, in PBS) for 1 h with gentle agitation, centrifuged for 3 min at 3600 RPM and the blocking buffer was removed. The cells were then incubated either in conjugated antibodies or in appropriate primary antibodies (i.e., WT1, AE1/AE3, CD45, CD34) with blocking buffer overnight at 4 °C with gentle agitation. Following primary antibody incubation, cells were centrifuged for 3 min at 3600 RPM and antibody aspirated. The cells were washed in PBS (Gibco) tween (0.2%) to remove any remaining antibody and centrifuged again for 3 min at 3600 RPM. In case of EpCAM, PBS was removed, and the cells were incubated in secondary antibody (Alexa Fluor 488 anti-mouse) diluted in blocking buffer (1:1000 v/v) for 1 h with gentle agitation. From this step onwards, the cells were protected from light as the fluorophore conjugated to the secondary antibody is light sensitive. Further centrifugation for 3 min at 3600 RPM was undertaken and the secondary antibody was removed. Following a wash in PBS tween and centrifugation, the PBS was removed, and the cells were resuspended in 99 μL Accumax (Innovative Cell Technologies, San Diego, CA, USA) to dissociate any cellular aggregates. 1 μL of DRAQ5™ (Biostatus Ltd. Loughborough, UK) nuclear stain was added before visualisation on the ImageStream. All the data files were then analysed on the IDEAS software [
9].
For the differential expression of AE1/AE3 and WT1 in SKOV3 and PEO1 cells, they were visualised under a Leica DM4000 microscope, x60.
2.8. Imagestream™
To positively identify ovarian cancer CCs, pan-cytokeratin (AE1/AE3), CD45 (to exclude lymphocytes), WT1 (a nuclear/cytoplasmic stain specific to high grade serous ovarian/fallopian tube/primary peritoneal cancer) and DRAQ5™ a general nuclear stain was all applied in different channels. EpCAM was applied to a subset of patients’ CCs to assess the value of EpCAM-related enumeration. CCs were identified and quantified using the IDEAs software.
2.9. RNA Extraction/cDNA Synthesis
If RNA extraction was not performed immediately, 0.5 mL of whole blood was added to 1.5 mL of RNAlater® (Life Technologies, Paisley, UK) in an Eppendorf tube, left at room temperature for up to 72 h and stored at −80 °C until further use. Samples stored in RNALater were first defrosted from −20 °C at room temperature then centrifuged at 13,000 RPM for 1 min, followed by aspirating supernatant. RNA from blood samples were extracted using the Ribopure RNA extraction kit (Ambion, Fisher, UK) according to the manufacturer’s protocol. cDNA concentration was normalised using RNA concentrations determined using the NanoDrop (Thermo Scientific, Waltham, MA, USA).
2.10. qPCR TATAA
The CTC GrandPerformance panel (TATAA Biocenter, Göteborg, Sweden) was used for the quantification of ovarian cancer-associated gene expression using reverse transcription quantitative real-time PCR (qPCR). The isolated mRNA was reverse-transcribed into cDNA using GrandScript (TATAA Biocenter), pre-amplified in a multiplex-PCR using GrandMaster PreAmp (TATAA Biocenter) and quantified by single-plex qPCR using GrandMaster mix (TATAA Biocenter). Data analysis and classification of samples was done with the GenEx software (MultiD).
2.11. Statistical Analysis
All statistical tests were performed using GraphPad Prism® Software (GraphPad Software). Statistical analyses were performed using one-way ANOVA followed by Tukey’s Multiple Comparison Test or student’s t-test with significance determined at the level of p < 0.05. P values are indicated in graphs as follows; * p = 0.01–0.05, ** p = 0.001–0.009, and *** p < 0.0009.
4. Discussion
In this study, rather than attempting to ‘enrich’ our samples for circulating tumour cells (CTCs), we chose to examine the whole liquid biopsy following red blood cell lysis. By co-staining with CK and WT1 markers, and CD45 to exclude WBCs, we show that there are abundant non-haematopoietic CCs (CK
+, WT1
+ and CD45
−) in patients with aEOC compared to healthy female controls. Detection of these CCs is reproducible, is reliably undertaken and can be informative. Previous studies may have underestimated the importance of CCs in patients with stage III/IV aEOC, which seem to represent a window in the microenvironment in which ovarian CTCs flourish. We show here that the size of CCs in ovarian cancer is not significantly different from the size of the rarer non-haematological cells found in healthy controls, suggesting that efforts to identify ovarian CTCs based on size alone can be challenging and further studies needed to assess if treatments can impact on cell size. We have also shown, via a 79 gene panel, that 4 genes were found upregulated mirroring changes at tissue level, whereas EpCAM does not appear up-regulated in relapsed aEOC patients. Hence the important EMT derived CTC-population is also likely to be missed by currently used techniques involving EpCAM reliant enrichment [
11].
A strength of flow cytometry is the ability to detect heterogeneity within a patient’s “CC-rich” sample, including sub-populations such as tumour cells and those of haematopoetic origin. However, flow cytometric assays are at risk of over-estimating rare events through non-specific antibody binding. Another false positive problem relates to antibody binding to soluble antigen forms. To mitigate these, we used DRAQ5™ staining for nuclear material to exclude apoptotic cells, cell debris, and platelet aggregates. Repeated gating and fluorescence compensation strategies may also augment inter-laboratory variability. Ideally >1 cell/µL is required to fulfil flow cytometric criteria with a coefficient variation of 10% [
12].
A further limiting factor includes inter and intra-patient variability, described as differences in patient traits, characteristics, genetics and/or pharmacokinetics (including in some instances exposure to environmental influences) [
13]. As highlighted earlier, the exact mechanism of CC shedding in malignancy is not well understood [
5] and composition of CCs may differ from patient to patient, possibly in response to different biological and physiological states. Despite these potential limitations, we have shown that there are significantly lower numbers of CCs in patients with relapsed aEOC compared to newly presenting and immediately post-operative patients. In general, relapsed patients are not as systemically unwell as those at presentation, probably reflecting more timely interventions when we are already aware of the likelihood of relapse. Exploring the composition of the CCs in larger relapse patient cohorts is required to determine the relevance of this finding.
We chose to explore the presence/absence of the Wilm’s Tumour protein, WT1 in CCs as nuclear/cytoplasmic expression of WT1 has 96% specificity and 100% sensitivity as a marker in serous tubal, ovarian, peritoneal cancer [
14,
15]. To the best of our knowledge, no other reports detail WT1 expression on CCs in ovarian cancer and its potential to be a biomarker candidate for the identification of ovarian cancer CTCs. The median number of WT1
+ CCs in our patient cohorts was 2914 per mL whole blood, which is significantly more than the median of 667 per mL found in the healthy volunteers. Most patients in this study had confirmed WT1 positive tissue diagnoses of serous ovarian, tubal, or primary peritoneal cancer. Surprisingly, in two young female volunteers, we detected ≥4000 WT1
+ CCs/mL which, however, is consistent with the presence of circulating endothelial cells found in females during the menstrual phase [
16]. Indeed, in our study, these two controls with high detectable levels WT1 were of menstrual age (31 years). The average age of all healthy volunteers was 40 years (range 21–64 years), with 7 of them not having any WT1
+ cells detected. Future studies using a larger age-matched control cohort will provide a better insight into how reliable this may be in premenopausal women in terms of indicating disease burden.
We found a significantly higher number of cytokeratin-positive CCs amongst our patient population compared to healthy controls. Techniques for concomitant membrane and nuclear staining are not sufficiently optimised yet to show how many of these are also WT1+ to confirm that these are tumour cells. We postulate that a proportion of them might be circulating endothelial cells (CECs). The presence of CECs has been recognised as a potential biomarker of vascular damage in many cancers including ovarian (500 CECs/mL) and renal cell (800 CECs/mL) cancers [
17,
18,
19]. Further work is on-going to confirm/refute this utilizing specific antibodies for endothelial cells (e.g., CD106, CD105, CD34, CD146). Preliminary evidence suggests presence of CD34
+ (marker for CECs) cells in the blood of ovarian cancer patients (
Figure S4). Another reason for the high numbers of CCs described in this study might be the absence of EPCAM enrichment, allows inclusion of cancer cells undergoing epithelial–mesenchymal transition (EMT), and not just those with pure epithelial phenotype.
EMT is a fundamental developmental process during which epithelial cells lose their polarity and junctional architecture to become mesenchymal cells that are motile and contributes in tumour metastasis [
20]. High grade serous ovarian cancer is believed to originate in the fimbrial ends of the fallopian tube where cells express a mesenchymal phenotype [
21]. E-cadherin is highly expressed in primary ovarian carcinomas, but this is lost in the advanced stages; a hallmark of EMT. This is suggestive of EMT facilitating the invasive ovarian cancer phenotype [
22]. In a more recent study, the authors argue that EMT is not complete in epithelial ovarian cancer and tumors can retain both phenotypes [
23]. In terms of EMT events recorded in liquid biopsies, breast cancer CTCs express mesenchymal and epithelial markers, but mesenchymal cells are highly enriched in CTCs and are associated with disease progression [
24]. To the best of our knowledge, there is no data relating to an EMT phenotype amongst CTCs of ovarian / fallopian tube origin, however during EMT transition there is a decrease in the expression of cytokeratins in ovarian cancer cells [
25]. Therefore, we hypothesised that the difference in staining found here could be also attributed to ovarian CCs undergoing EMT. In a proof of principle experiment using SKOV3 cells (which exhibit an intermediate mesenchymal phenotype) and PEO1 cells (which exhibit an epithelial phenotype); we show that more PEO1 than SKOV3 cells stained positive with CK.
Although, generally less sensitive than circulating-free DNA (cfDNA) [
26,
27] whole CC analysis can be performed at DNA, RNA (mRNA or micro-RNA) and protein level, providing further insight into the biology and microenvironment of cancer metastases. Contamination of CTCs with other non-haematopoietic CCs, is little different to examining solid tumour tissue where cancer cells are within a microenvironment of normal and reactive stromal tissue. Our data from whole blood provides novel insights into the differential expression of genes from liquid biopsies but also poses certain limitations. For example, WBCs expression patterns vary with both inflammatory and nutritional states which are difficult to predict to allow accurate assessment of changing CC phenotypes. Future studies should concentrate on assessing gene expression of CCs only. However, current methods rely heavily on EPCAM dependent isolation techniques which will miss cells undergoing EMT transition. Future studies could capitalise on data available from databanks, although the current bioinformatic indexes (e.g., Oncomine, The Cancer Genome Atlas; TCGA) rely mostly on microarray/RNAseq data from tissue rather than from liquid biopsies. Exploration of the CC signatures of larger cohorts of aEOC patients, and correlation with paired diagnostic tissue, is planned to identify changes which will hopefully yield alternative therapeutic strategies or targets that can be harnessed for prognostic purposes. We suggest that this simple method of extracting CC-rich information has the potential for exploring part of the circulating microenvironment within the remit of nearside patient testing and rapid responses required for 21st century treatment.