Establishment of an optimized CTC detection model consisting of EpCAM,MUC1 and WT1 in epithelial ovarian cancer and its correlation with clinical characteristics

Tongxia Wang ,Yan Gao ,Xi Wang ,Junrui Tian ,Yuan Li ,Bo Yu ,Cuiyu Huang,Hui Li,Huamao Liang,David M.Irwin,Huanran Tan,Hongyan Guo

1Department of Obstetrics and Gynecology,Peking University Third Hospital,Beijing 100191,China;2Department of Pharmacology,School of Basic Medical Sciences,Peking University Health Science Center,Beijing 100191,China;3Department of Laboratory Medicine and Pathobiology,University of Toronto,Toronto M5S1A8,Canada

Abstract Objective:Emerging studies have demonstrated the promising clinical value of circulating tumor cells (CTCs)for diagnosis,disease assessment,treatment monitoring and prognosis in epithelial ovarian cancer.However,the clinical application of CTC remains restricted due to diverse detection techniques with variable sensitivity and specificity and a lack of common standards.Methods:We enrolled 160 patients with epithelial ovarian cancer as the experimental group,and 90 patients including 50 patients with benign ovarian tumor and 40 healthy females as the control group.We enriched CTCs with immunomagnetic beads targeting two epithelial cell surface antigens (EpCAM and MUC1),and used multiple reverse transcription-polymerase chain reaction (RT-PCR) detecting three markers (EpCAM,MUC1 and WT1)for quantification.And then we used a binary logistic regression analysis and focused on EpCAM,MUC1 and WT1 to establish an optimized CTC detection model.Results:The sensitivity and specificity of the optimized model is 79.4% and 92.2%,respectively.The specificity of the CTC detection model is significantly higher than CA125 (92.2% vs.82.2%,P=0.044),and the detection rate of CTCs was higher than the positive rate of CA125 (74.5% vs. 58.2%,P=0.069) in early-stage patients (stage I and II).The detection rate of CTCs was significantly higher in patients with ascitic volume ≥500 mL,suboptimal cytoreductive surgery and elevated serum CA125 level after 2 courses of chemotherapy (P＜0.05).The detection rate of CTCEpCAM+ and CTCMUC1+ was significantly higher in chemo-resistant patients (26.3% vs.11.9%;26.4% vs.13.4%,P＜0.05).The median progression-free survival time for CTCMUC1+ patients trended to be longer than CTCMUC1? patients,and overall survival was shorter in CTCMUC1+ patients (P=0.043).Conclusions:Our study presents an optimized detection model for CTCs,which consists of the expression levels of three markers (EpCAM,MUC1 and WT1).In comparison with CA125,our model has high specificity and demonstrates better diagnostic values,especially for early-stage ovarian cancer.Detection of CTCEpCAM+ and CTCMUC1+ had predictive value for chemotherapy resistance,and the detection of CTCMUC1+ suggested poor prognosis.

Keywords:Circulating tumor cells;epithelial ovarian cancer;optimized detection model;diagnosis and prognosis

Introduction

Epithelial ovarian cancer (EOC) is the most lethal gynecologic malignancy,with a 5-year survival rate of less than 35% (1,2).Carbohydrate antigen 125 (CA125) is the most commonly used serological marker for the clinical diagnosis of EOC,but its diagnosis rate for early-stage disease is only 50% (3).In addition,its specificity is also quite poor,as elevated serum CA125 levels can also be seen in inflammatory diseases and endometriosis (4).Therefore,there is a need to explore additional detection methods that might improve early diagnosis of ovarian cancer.

Detection of circulating tumor cells (CTCs),as a form of liquid biopsy,is a non-invasive tumor detection method that has attracted extensive attention in recent years (5) due to its advantages such as easy access,minimal invasiveness and good patient compliance.CTCs are cells that are shed into the bloodstream from a primary tumor,metastases,or recurrent tumors,and therefore possess tumor-specific characteristics such as large size,low density,less deformation,and a high nucleoplasmic ratio (N/C) (6),thus,making them good targets for liquid biopsy.The clinical value of CTCs in ovarian,breast,prostate and colorectal cancers has been confirmed in several studies and meta-analyses (5,7-13),with nearly thirty articles concerning the clinical application of CTCs to ovarian cancer being published in the last decade.Due to the variety of detection techniques used,the reported detection rates in these publications are variable,ranging from 12%to 100% (14-22).Currently,a common standard for CTC detection has not been established.

Here,through a review of relevant research databases and publications,we identified three markers (EpCAM,MUC1 and WT1) that are significantly highly expressed in ovarian cancer cells and then used these markers in our CTC detection model.In our previous study (12),we explored the diagnostic value of these markers in CTCs of ovarian cancer,however we only enrolled healthy females as controls,which lead to a sample limitation,as the actual clinical need is to differentiate ovarian cancer from benign ovarian tumors.Therefore,in this study we enrolled EOC patients as the experimental group,and benign ovarian tumor patients as well as healthy females,as the control group.We aim to analyze the diagnostic value of the above-mentioned markers in the detection of CTCs in EOC and the correlation of levels of these markers with clinical features,and to establish an optimized CTC detection model with high sensitivity and specificity to assist in clinical applications.

Materials and methods

Study population

This study included patients with EOC,patients with benign ovarian tumors,and healthy females who were admitted to the Department of Gynecology,Peking University Third Hospital between September 2013 and December 2020.Peripheral venous blood (5 mL) was drawn from all patients before initial treatment,with serum CA125 levels being measured at this time.All EOC patients were pathologically verified,with tumor stage determined according to the International Federation of Gynecology and Obstetrics (FIGO) staging,and patients overlapped with other malignant tumors will be excluded.Initial treatment included advanced chemotherapy or full staging of ovarian cancer surgery and tumor cytoreductive surgery,and postoperative adjuvant chemotherapy mainly based on platinum.Patients with benign ovarian tumors refer to patients with benign tumors that have been pathologically proven to be of ovarian origin,including ovarian chocolate cysts,teratomas,and serous cystadenomas,etc.The healthy female population refers to adult female patients who have no obvious abnormalities in the regular physical examination every year and whose CA125 is in the normal range.Clinical data including medical history,physical examination,laboratory examination and imaging evaluation were obtained from all enrolled patients.This study was approved by the Peking University Biomedical Ethics Committee (No.IRB00001052-16022).All patients signed the informed consent form before participation.

Study methods

We enriched CTCs with immunomagnetic beads targeting to two epithelial cell surface antigens (EpCAM and MUC1),which are not expressed on the surface of blood cells,and then used multiple reverse transcriptionpolymerase chain reaction (RT-PCR) to detect three markers (EpCAM,MUC1 and WT1) to quantify the enriched CTCs.The binding of two antibody-labeled immunomagnetic beads (EpCAM and MUC1) to the EOC cell line SKOV3 is shown inSupplementary Figure S1.The primer sequences of three genes are shown inSupplementary Table S1.A DNA 1000 LabChips (Agilent Technologies,Santa Clara,USA) was used to detect the amplified sequences using an Agilent 2100 Bioanalyzer(Agilent Technologies,Santa Clara,USA) according to the manufacturer’s instructions.The Agilent 2100 bioanalyzer has a list of Ladder with known DNA concentration,and polymerase chain reaction (PCR) amplification product concentration of each marker gene is obtained by comparing with Ladder.

Then we focused on EpCAM,MUC1 and WT1 to establish an optimized detection model.We used a binary logistic regression analysis,where we setYas a binary dependent variable,EpCAM,MUC1 and WT1 are independent variables.The logistic regression model equation is linearly expressed as:

PCEpCAM:PCR concentration of EpCAM;PCMUC1:PCR concentration of MUC1;PCWT1:PCR concentration of WT1.

Criteria for determining treatment efficacy and chemosensitivity are according to the 2021 National Comprehensive Cancer Network (NCCN) clinical practice guidelines (23).Chemo-resistance was defined as recurrence within less than 6 months when achieving complete recovery for an initial treatment,or disease progression during chemotherapy.Chemosensitive patients were defined as recurrence after more than 6 months when achieving complete recovery for an initial treatment (23).

Prognostic criteria were assessed as progression-free survival (PFS) and overall survival (OS) of patients with EOC.During follow-up,the evidence of recurrence or progression was based on the imaging or pathology confirmation on biopsy or based on serum CA125 levels according to Response Evaluation Criteriain Solid Tumors(RECIST) Version 1.1 and adhering to the Gynecological Cancer Intergroup (GCIG) definition as greater than,or equal to,two times the upper limit of the reference range on two occasions at least 1 week apart (24).OS was calculated from date of treatment (surgery or chemotherapy) to date of death or date of last follow-up (end of follow-up).PFS was defined as the time from treatment(surgery or chemotherapy) to physical,biological or radiological evidence of disease progression,or death as a result from any cause (23).

Statistical analysis

Data were statistically analyzed using IBM SPSS statistical software (Version 23.0;IBM Corp.,New York,USA).Binary logistic regression was used to build an optimized CTC detection model.PCR concentration was compared using a Mann-Whitney U test. Receiver operating characteristic (ROC) curve analysis was used to evaluate the ability of different PCR marker product concentrations to distinguish EOC from non-cancerous control groups.These results were expressed as area under the ROC curve(AUC),with 95% CI and P values.Sensitivity,specificity and Youden index of the diagnostic results were calculated.The threshold corresponding to the maximum value of Youden index was taken as the best threshold for clinical diagnosis.Count data were expressed as n (%) andχ2test was used.Multivariate Cox regressions were used on OS analysis,and P＜0.05 was considered statistically significant.The survival curve was drawn by Graphpad Prism 8.4.2(GraphPad Software,Inc.,San Diego,CA,USA) for the prognostic correlation analysis.

Results

Clinical characteristics of enrolled patients

A total of 160 patients with EOC were enrolled in the experimental group,whereas 90 patients were enrolled in the non-cancerous control group,including 50 patients with benign ovarian tumor (i.e.,teratoma,ovarian chocolate cysts and serous cystadenoma) and 40 healthy females.The general clinical characteristics of the patients are presented inTable 1.

Table 1 Patients’ clinical information

Expression level differences of three marker (EpCAM,MUC1 and WT1) in EOC patients vs.non-cancerous controls

mRNA abundance of EpCAM,MUC1 and WT1 was quantified using multiplex RT-PCR.Comparison of the mRNA abundance levels of EpCAM,MUC1 and WT1 in EOC patients and non-cancerous controls showed that all of the three markers were significantly more abundant inEOC (P＜0.001,Figure 1).

We plotted ROC curves (Supplementary Figure S2) to find the sensitivity and specificity of single marker(EpCAM,MUC1 and WT1) in distinguishing EOC from non-cancerous controls.The thresholds,sensitivity and specificity of each single marker are shown inTable 2.The results of this analysis showed that each marker gene had a high specificity (92.1%?95.5%),but relatively low sensitivity (56.1%?62.6%) (Table 2).

Establishment of an optimized detection model for CTCs that distinguishes EOC from non-cancerous controls

To assess the contributions of three marker genes more comprehensively,we used a binary logistic regression analysis.The regression coefficients of EpCAM,MUC1 and WT1 are shown inSupplementary Table S2.The regression equation of the optimized CTC detection model is:

PCEpCAM:PCR concentration of EpCAM;PCMUC1:PCR concentration of MUC1;PCWT1:PCR concentration of WT1.

As shown inFigure 2,the optimal threshold forYis 0.617,the corresponding sensitivity is 79.4%,and the specificity is 92.2%.So we defineY≥0.617 as positive for CTC detection.And we also performed ten-fold cross validation on the dataset (Supplementary Table S3).

Figure 1 Expression of three markers in EOC patients vs. non-cancerous controls.mRNA abundance levels are quantified by PCR.Expression of (A) EpCAM (P＜0.001);(B) MUC1 (P＜0.001) and (C) WT1 (P＜0.001) in EOC patients vs. non-cancerous controls,respectively.EOC,epithelial ovarian cancer;PCR,polymerase chain reaction;EG,experimental group;CG,control group.

Correlation analysis between CTC detection model and clinical features in EOC patients

In comparison with CA125,the specificity of the optimized model was significantly higher than that of CA125 (92.2%vs.82.2%,P=0.044),while the sensitivity was also slightly higher than that of CA125 (79.4%vs.75.6%,P＞0.05).Especially in early-stage EOC (stage I and II),the detection rate of CTCs was significantly higher than the positive rate for CA125 (74.5%vs.58.2%,P=0.069).But in advancedstage EOC (stage III and IV),the detection rate of CTCs was close to CA125 (81.9%vs.84.8%,P=0.579).

Univariate analysis showed that the detection rate of CTCs was significantly higher in patients with ascitic volume ≥500 mL,suboptimal cytoreductive surgery and elevated serum CA125 levels after 2 courses of chemotherapy (P＜0.05,Table 3).

Application value of CTC optimized detection model and CTCs characterized by single marker genes in predicting chemo-resistance and prognosis in EOC patients.

Among the 160 EOC patients,115 were found to be chemo-sensitive and 28 were chemo-resistant.The results showed that the detection of CTCs through the optimized detection model was not correlated with chemotherapy resistance (P＞0.05,Supplementary Table S4).But the detection rate of CTCEpCAM+(defined as EpCAM PCR concentration ≥0.195 ng/μL) and CTCMUC1+(defined as MUC1 PCR concentration ≥ 0.355 ng/μL) were significantly higher in chemo-resistant EOC patients(P＜0.05,Supplementary Table S4).

Table 2 Detection thresholds,sensitivity and specificity of single marker to distinguish EOC patients from non-cancerous control group

The time distribution of PFS was 4?84 months,with a median PFS of 34 months.The time distribution of OS was 5?84 months,with a median OS of 42 months.While no statistical difference was found,PFS of CTCMUC1+patients trended to be shorter than that of CTCMUC1?patients (P=0.357,Figure 3A).But OS of CTCMUC1+patients was significantly shorter (P=0.043,Figure 3B).In addition,PFS and OS of CTC+(detected by the optimized model),CTCEpCAM+and CTCWT1+patients did not show any significant difference with their corresponding negative patients (P＞0.05,Supplementary Figure S3?5).

Figure S1 Confocal laser imaging of immunomagnetic beads combined with SKOV3 cell line.Immunomagnetic beads were incubated with two monoclonal antibodies (EpCAM and MUC1) without staining;ovarian cancer cell line was SKOV3 cell line,and nucleus was stained with DAPI.

Figure S2 ROC curves for three markers to distinguish EOC patients from the non-cancerous control group.AUC of EpCAM,MUC1 and WT1 is 0.794,0.721 and 0.791,respectively.ROC,receiver operating characteristic;EOC,epithelial ovarian cancer;AUC,area under the curve.

Figure 2 ROC curve for three combined markers to distinguish EOC patients from the non-cancerous control group(AUC=0.893;Y=0.617;specificity=0.922;sensitivity=0.794).ROC,receiver operating characteristic;EOC,epithelial ovarian cancer;AUC,area under the curve.

Figure 3 Comparison of PFS and OS curves in CTCMUC1+ and CTCMUC1? patients.(A) Median PFS of CTCMUC1+ and CTCMUC1?patients was 35.0 and 45.0 months,respectively,but there was no statistical difference (HR=1.242;95% CI:0.787?1.959;P=0.357);(B) OS of CTCMUC1+ patients tended to be significantly shorter than CTCMUC1? patients (HR=1.900;95% CI:1.020?3.540;P=0.043).PFS,progression-free survival;OS,overall survival;CTC,circulating tumor cell;HR,hazard ratio;95% CI,95% confidence interval.

Figure S3 Comparison of PFS and OS between CTC+and CTC? EOC patients.(A) Median PFS of CTC+and CTC? patients was 28.5 and 48.0 months,respectively,but there was no statistical difference (HR=1.137;95% CI:0.638?2.025;P=0.386);(B) OS of CTC+and CTC? patients was not statistically different (HR=0.838;95% CI:0.388?1.808;P=0.652).PFS,progression-free survival;OS,overall survival;CTC,circulating tumor cell;HR,hazard ratio;95% CI,95% confidence interval.

Figure S4 Comparison of PFS and OS curves between CTCEpCAM+ and CTCEpCAM? patients.(A) PFS of CTCEpCAM+ and CTCEpCAM?patients was not statistically different (HR=1.049;95% CI:0.659?1.670;P=0.841);(B) OS of CTCEpCAM+ and CTCEpCAM? patients was not statistically different (HR=0.960;95% CI:0.505?1.823;P=0.899).PFS,progression-free survival;OS,overall survival;CTC,circulating tumor cell;HR,hazard ratio;95% CI,95% confidence interval.

Figure S5 Comparison of PFS and OS curves between CTCWT1+ and CTCWT1? patients.(A) Median PFS of CTCWT1+ and CTCWT1?patients was 42.0 and 35.0 months,respectively,but there was no statistical difference (HR=0.869;95% CI:0.546?1.384;P=0.308);(B) OS of CTCWT1+ and CTCWT1? patients was not statistically different (HR=0.735;95% CI:0.391?1.381;P=0.201).PFS,progression-free survival;OS,overall survival;CTC,circulating tumor cell;HR,hazard ratio;95% CI,95% confidence interval.

Multivariate Cox regression analysis was performed on the effects of CTCMUC1+,staging,the level of CA125,ascites volume,cytoreductive surgery,lymph node metastasis,and chemotherapy sensitivity on the prognosis of patients.The results showed that advanced staging(P=0.006) and chemotherapy resistance (P＜0.001) were independent risk factors affecting PFS of patients(Supplementary Table S5).The volume of ascites ≥500 mL(P=0.006),suboptimal cytoreductive surgery (P=0.046) and chemotherapy resistance (P＜0.001) were independent risk factors affecting OS of patients (Supplementary Table S6).

Discussion

CTCs have become popular as it is expected that they will become screening markers for early diagnosis,treatment monitoring,chemotherapy reaction and prognostic evaluation (12,25).A large amount of research concerning ovarian cancer CTCs has been reported in the recent 5 years (Table 4).Though promising evidence for their clinical application in assessing patient prognosis,real-time monitoring of treatment effects and disease recurrence has been reported (12,25),many researchers still focus on continuous optimization of CTC detection technologies.

Classically,the Cell Search? system is a positive enrichment method using an EpCAM targeted immunomagnetic bead to separate CTCs from peripheral blood samples (43),and is the only FDA-approved CTC detection technology that can be used for metastatic breast,colorectal and prostate cancers (11,40).However,the biggest limitation of this detection method is that EpCAM is a marker of epithelial cells and not a tumor-specific marker (44).When tumor cells undergo epithelialmesenchymal metastasis (EMT),CTCs expressing mesenchymal-derived markers become undetectable,resulting in the loss of CTC detection (45,46).Besides,we reviewed the researches of CTC detection methods for EOC over the past 5 years,a series of new detection methods have emerged,such as streptavidin-biotin cascade amplification effect coupled with folate receptor binding method (25),invasiveness CTC subgroup labeling method(18,29) and a high-throughput microfluidic technology system designed by organically combining physical and biological characteristics (37,39).Also,based on the antigen-antibody binding enrichment methods similar to CellSearch system,studies reported the use of immunomagnetic beads coupled with two or more antibodies(including epithelial markers:EpCAM,MUC-1,CK,HER2 and ESA;interstitial markers:Vimentin and Ncadherin).The above mentioned detection methods significantly increased the detection rates by 58%?95%(Table 4).

Table 3 Univariate analysis of factors affecting CTC detection in EOC

Table 4 Overview of studies on detection methods for CTCs in ovarian cancer.

Table 4(continued)

Our study also used an antigen-antibody binding method to screen and separate CTCs with immunomagnetic beads,which was then combined with multiple RT-PCR to detect CTCs in the peripheral blood of EOC patients.To compensate for the decrease in detection rate caused by weak expression of a single target antigen,three molecules(EpCAM,MUC1 and WT1),which are significantly more abundant in ovarian cancer cells,based on a literature review,were selected as gene markers to establish an optimized CTC detection model.Early in 2018,our laboratory successfully used this method to detect CTCs in the peripheral blood of EOC patients (12).If PCR concentration of any of the six marker (EpCAM,HER2,MUC1,WT1,P16 and PAX8) was greater than or equal to 0.3 ng/μL,it was defined as positive for CTCs.Based on this detection method,we detected CTCs in 90% (98/109)of the newly diagnosed EOC patients,demonstrating the feasibility of this detection method.However,our previous study only enrolled healthy females in the control group which inevitably becomes a limitation,as benign ovarian tumor is the first need to be differentiated from EOC.Our pre-experiment showed that,when compared with benign ovarian tumor patients,only EpCAM,MUC1 and WT1 transcripts were significantly more abundant in EOC patients (P＜0.05).This result suggests that the other threemarker (HER2,P16 and PAX8) fail to differentiate benign ovarian tumor from EOC.Therefore,to establish optimized detection model,we expanded the sample sizes and focused on EpCAM,MUC1 and WT1 to achieve optimal differentiation between the 160 EOC and 90 noncancerous patients.

Similar to our results,relevant studies demonstrated that EpCAM,MUC1 and WT1 are significantly more abundant in ovarian cancer (32,47-50).Among them,EpCAM is a type I transmembrane glycoprotein located on epithelial cells (51) and is a tumor-associated antigen that is highly expressed in ovarian cancer and is also the most widely used marker for ovarian cancer CTC detection technology (47,48).Zhaoet al.(52) used the Gene Expression Omnibus (GEO) database to analyze differentially expressed genes in ovarian cancer and found that the abundance of EpCAM transcripts in ovarian cancer differs significantly from normal tissues (P=0.0048),and especially is high in ovarian cancers with recurrence,metastasis,invasiveness and chemotherapy resistance (53).This supports EpCAM role in the occurrence and metastasis of ovarian cancer.Mucin 1 (MUC1) gene encodes a membrane-bound O-glycosylated protein that plays an essential role in forming protective mucous barriers on epithelial surfaces (54),which is overexpressed in 90%?100% of serous carcinomas (55-57).The role ofMUC1 has been verified in mouse models where MUC1 expression is necessary for the occurrence of ovarian cancer(56,57).The study by Houet al.(49) confirmed that the expression level of MUC1 in malignant EOC tissues was significantly higher than that in borderline,benign and normal ovarian tissues (P＜0.05).WT1(Wilms tumor gene)is located in the nuclei and encodes a tumor-associated antigen (58).In epithelial ovarian tumors,WT1expression has been detected,and studies showed that ovarian cancers with high expression of WT1 are more aggressive (50) and suggestive of poor prognosis (59,60).

Based on single indicators for CTC detection,our results showed that the specificity of each gene marker was high(92.1%?95.5%),but their sensitivity was relatively low(56.1%?62.6%),which suggests that we need to combine multiple markers to construct an optimized model to increase the sensitivity of detection while ensuring a high specificity.However,the combined application of these three markers for the detection of CTCs in ovarian cancer has been rarely reported.We used a binary logistic regression analysis to set up the regression equation of the optimized CTC detection model.When the threshold was set as 0.617 (Y≥0.617 as positive for CTCs),the sensitivity was 79.4%,and the specificity was 92.2%.Compared with the application of a single index,the sensitivity of the optimized detection model was improved without any evident decrease in specificity.Similar to our results,Abreuet al.(42) published results from their collection of peripheral blood CTCs from 38 patients with advanced ovarian high-grade serous cancer and also established an optimized gene panel for CTC RT-qPCR detection.Their model was composed of GAPDH,TIMP1,CK19 and MUC1 (AUC:0.938,P＜0.01;95% CI:0.868?1),and its detection sensitivity was 63%,which is lower than our optimized model,but had a better specificity of 100%.

Then we compared the diagnostic value of our CTC detection model with serum CA125 level.CA125 is widely used as a clinical serological marker,which is commonly found in the serum of EOC patients (61),so it can assist the diagnosis of EOC and can be used to monitor the progress of the disease (62).However,the sensitivity of CA125 in early-staged EOC is quite low,resulting in an early diagnosis rate of only 50% (3).Also,the specificity of CA125 is rather poor,as inflammation,endometriosis and other benign diseases can also present with elevated serum CA125 levels (4).Therefore,relying solely on serum CA125 levels to assist the diagnosis of ovarian cancer,monitor treatment effects and recurrence has limitations.The detection of CTCs can to some extents make up for the limitations of CA125.Studies (12,19,35) have shown that,compared with CA125,CTCs have a higher detection rate for early diagnosis.Our study also confirms this point of view.In early-stage EOC (stage I and II),the detection rate of our CTCs detection model was significantly higher than the positive rate of CA125 (74.5%vs.58.2%,P=0.069),and also showed a higher specificity than that of CA125 (92.2%vs.82.2%,P＜0.05).This suggests the advantage of our CTC detection model in the early diagnosis of EOC and the feasibility to combine CTCs and CA125 to improve the diagnostic rate of EOC.Moreover,univariate analysis showed that the detection rate of CTCs was significantly higher in patients with ascitic volume ≥500 mL and in patients with suboptimal cytoreductive surgery.The large amount of ascites and the incomplete cytoreductive surgery all indicate the presence of a large tumor burden and advanced stage,which suggests that the detection of CTCs can reflect the severity of the disease.

A large number of studies have confirmed that CTCs can indirectly reflect the tumor proliferation and chemosensitivity through dynamically monitoring the changes in CTCs counts during chemotherapy,indicating that CTCs can be used as a predictive marker to assess the effect of chemotherapy (63-65).However,results from relevant studies regarding this question have been controversial.Some studies concluded that there was no statistically significant association between CTC detection and chemoresistance (26,66).Although the results of our study showed no statistical difference between CTCs detection and chemo-resistance.In clinical practice,we found that patients with elevated serum CA125 levels after more than 2 courses of postoperative chemotherapy often had poor outcomes and usually ended up with a high incidence of chemoresistance or uncontrolled tumor progression.Our study also found that the detection rate of CTCs in this group of patients was significantly higher,suggesting a correlation between CTCs and chemo-resistance.Besides,an in-depth analysis of the abundance of the different marker genes suggested that CTCs characterized by positive expression of EpCAM or MUC1 have a predictive value for chemoresistance in EOC patients.

Recent studies have been devoted to exploring the predictive value of CTCs as tumor markers for prognosis of ovarian cancer,but the result is yet undetermined(12,18,21,34,36,38,39).We reviewed 20 articles published in the past 10 years,and more than half of the studies concluded that the detection of CTCs in peripheral blood is closely related to the shortening of PFS and OS (Table 4).Three published meta-analysis articles summarize relevant studies that CTCs could predict the prognosis of ovarian cancer patients (67-69),but there are also a few studies that concluded that CTCs are not significantly related to PFS and OS (17,26,70).Unfortunately,our study showed that there was no significant correlation between the CTC detection model and the prognosis.However,we found that patients with CTCs characterized by positive expression of MUC1 (CTCMUC1+) had a significantly shorter OS,and a trend for shorter PFS than patients without,which suggests that the detection of CTCMUC1+may indicate poor prognosis.This finding indicates a correlation between different CTCs’ molecular phenotypes and clinical features.The rapid development of single-cell sequencing has paved the way for the molecular characterization of CTCs (71).Our following studies are needed to resolve the roles of the molecular phenotypic characteristics of CTCs,and explore the mechanisms for ovarian cancer metastasis,recurrence and chemotherapy resistance,and potentially guide individualized treatment of ovarian cancer patients.

Table S1 Primer sequences of three genes

Table S2 Regression coefficients

Table S3 Ten-fold cross validation results

Table S4 Detection of CTCs through optimized detection model and CTCs characterized by single marker gene (EpCAM,MUC1 and WT1) in EOC patients with different levels of chemo-sensitivity

Table S5 Multivariate Cox regression analysis on PFS in patients with EOC

Table S6 Multivariate Cox regression analysis on OS in patients with EOC

There are some limitations of our study that should be improved in subsequent studies.The detection method of our study was not to capture intact CTCs from the blood,but instead to obtain mRNA after cell lysis from enriched cells and to evaluate expression of marker genes by RTPCR.In addition,since the follow-up period of most patients was short,future studies should extend the followup period and expand sample size to validate the above findings.

With the continuous optimization of CTCs detection technology,with increased sensitivity and specificity,it may become possible to isolate viable and intact CTCs.Development of CTCsin vitrocell culture technology,the construction of biomimetic models,and the combined application of genetic assessment and bioinformatics analysis should help reveal the molecular biological characteristics of primary ovarian cancer,metastases and recurrences.This should lead to improve understanding of the mechanisms of ovarian cancer metastasis and chemotherapy resistance,identify new therapeutic targets,and guide the individualized treatment of ovarian cancer patients.

Conclusions

This study used immunomagnetic beads targeting EpCAM and MUC1 for CTC separation,then quantified the expression of EpCAM,MUC1 and WT1 in peripheral blood for CTC detection,and further used these three markers to establish an optimized detection model for CTCs.Compared with a single marker gene,the optimized detection model had higher sensitivity and specificity,and had better diagnostic value for early-staged ovarian cancer than CA125.However,for predicting drug resistance and prognosis,positive CTCs characterized by a single marker appeared more advantageous than the optimized detection model.Detection of CTCEpCAM+and CTCMUC1+has predictive value for chemotherapy resistance and the detection of CTCMUC1+may indicate poor prognosis.

Acknowledgements

None.

Footnote

Conflicts of Interest:The authors have no conflicts of interest to declare.