CDCA3 Predicts Poor Prognosis and Affects CD8+ T Cell Infiltration in Renal Cell Carcinoma

Background Cell division cycle associated 3 (CDCA3) mediates the ubiquitination WEE1 kinase at G2/M phase. However, its contribution to cancer immunity remains uncertain. Methods We first evaluated the effect of CDCA3 on the prognosis of patients with renal cell carcinoma (RCC). The results of bioinformatics analysis were verified by the tissue microarray, immunofluorescence (IF) staining, CCK-8 assay, colony formation, cell cycle, and Western blot. Results Bioinformatics analysis predicated CDCA3 was an independent predictor of poor prognosis in RCC and was associated with poor TNM stage and grade. CDCA3 was related to the infiltration of CD8+ T cells and Tregs. Tissue microarray demonstrated that CDCA3 was strongly associated with poor prognosis and positively relevant to CD8+ T infiltration. In vitro experiments showed that exgenomic interference of CDCA3 could attenuate cellular proliferation, arrest cell cycle, and blockade accumulation of CDK4, Bub3, and Cdc20 in mitosis process. Conclusion CDCA3 presents as a good biomarker candidate to predict the prognosis of RCC patients and potentiates the immune tumor microenvironment (TME) of RCC.


Introduction
Renal cell carcinoma (RCC) is a malignancy from the kidney epithelium and the mobility has steadily increased globally in recent years [1]. The first-line antiangiogenic therapies such as tyrosine kinase inhibitors (TKI) have presented the certain effect for RCC patients, however, the response is discontinued in short time for the majorities [2]. Immune checkpoint inhibitors (ICIs) usher a new time of cancer therapeutic strategies via sparking anticancer immunity [3]. CD8 + T cells serve as an essential effector and partially relevant to the effect of ICI [4]. Traditionally, RCC is considered as an immunogenic cancer, and immunotherapy has shown a certain effect of RCC [5,6]. In clinical practices, we observe the effect of ICIs is diversified, however, scholars fail to find a good candidate to predicate the response and adverse effects (AEs) of ICI in RCC treatment. The biomarkers will also help identify subgroups that respond to immunotherapy and avoid severe AEs.
Cell division malfunctions trigger tumor development and antitumor immune response [7]. CDCA3 has been shown to be a poor prognostic factor for renal papillary cell carcinoma, nonsmall cell lung cancer, etc. [8][9][10]. Scholars reveal that CDCA3 was upregulated in RCC and promote tumor progression and sunitinib resistance [11] via activating the NF-κB/cyclin D1 signaling axis [12]. There are data indicating that CDCA3 can serve as an important biomarker to evaluate the therapeutic sensitivity of TKI and therefore it would be appropriate to underline this aspect also in light of the possible associations of immunological therapies and TKI in various types of malignant tumors [13]. Moreover it should be very interesting to test the role of CDCA3 as a predictive biomarker of toxicity related to a prolonged use of these novel agents in combination therapy of RCC [14]. However, the immune impact of CDCA3 has also not been well reported.
In this paper, we try to evaluate the predicable performance of CDCA3 in RCC and figure out the attribution of     3 Journal of Oncology diagram and calibration curve to better interpret the prognostic significance of CDCA3. Moreover, fold change = 2 was used to compare the differences of gene expression among different expression groups of CDCA3, and a heat map of differentially expressed genes was drawn to show the expression trend in different groups. Finally, considering that CDCA3 can be used as an oncogene to affect the progression of tumor, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis on the upregulated genes of CDCA3 in different expression groups to identify CDCA3 functional pathway localization in tumors.

Correlation between Tumor
Immune Cell Infiltration and CDCA3 Gene Expression. Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBER-SORT) algorithm was used to estimate the infiltration proportion of 22 kinds of immune cells in normal kidney and RCC samples to describe the profile of immune cell infiltration in RCC. The abundance of immune cells infiltration and the expression of 8 important immune checkpoints (CD274, CTLA4, HAVCR2, LAG3, PDCD1, PDCD1LG2, TIGIT, and SIGLEC15) among different CDCA3 expression groups were compared. Finally, we also analyzed the correlation between CDCA3 expression with tumor mutation burden (TMB) and MSI. Tumor Immune Single-cell Hub (TISCH, http://tisch .comp-genomics.org/) is a scRNA-seq database focusing on TME. We obtained the relationship between CDCA3 and RCC TME at single-cell level in TISCH.

Cell Culture and Transfection of Lentivirus.
Caki-1 and 786-O were purchased from the Type Culture Collection (Chinese Academy of Sciences, Shanghai, China). Cells were cultured in RPMI-1640 medium (HyClone, USA) with 10% fetal bovine serum (Gibco, Grand Island, NY, USA). The culture was maintained in a humidified incubator with 37°C, 5% CO 2 . CDCA3 knockdown lentivirus was designed by Obio Technology Corp (Shanghai, China). Then, Caki-1 and 786-O were transfected with the lentivirus, according to the manufacturer's instructions. Two days later, puromycin was added for screening. Knockdown efficiencies of CDCA8 were assessed by Western blot.

Cell Proliferation Assay.
The cells were seeded in 96-well plates (1,000 cells/well) and cultured for 1, 2, and 3 days. After adding 10 μl CCK-8 (Dojindo, Japan) to each well and incubating at 37°C for 2 h, the absorbance at 450 nm was measured by the Rayto-6000 system (Rayto, China). 2.8. Tissue Microarray. The RCC tissue microarray was purchased from Outdo (Shanghai, China) and contains 150 RCC tissues and 30 paired paracancer tissues along with their survival, clinical information, etc. Samples were collected from the National Human Genetic Resources Sharing Service platform (2005DKA21300). All points on the chip were detected by Immunofluorescence (IF). The expression of CDCA3, CD8, CD4, CD68, FOXP3, and PD-1 was detected by intensity and positive number of IF. We divided 150 RCC patients into two groups based on the optimal CDCA3 cut-off value and plotted survival curves to identify their prognostic significance. Finally, we analyzed the correlation between CDCA3 and CD8, CD4, FOXP3, CD68, and PD-1.  5 Journal of Oncology 2.9. Immunofluorescence Staining. Tissue microarray were deparaffinized by graded alcohol and then washed three times with phosphate-buffered saline (PBS), permeabilized with 0.4% Triton X-100 for 30 min, and blocked with goat serum working liquid (Wuhan Boster Biological Technology, Wuhan, China) for 2 hours after antigen retrieval. The sections were then incubated overnight with mixed primary antibodies at 4°C, washed in PBS to remove unbound primary antibodies, and incubated with secondary antibodies in the dark at room temperature (RT) for 1 hour. The sections were counterstained with 4 ′ , 6 diamidino-2phenylindole (Sigma-Aldrich) for 5 minutes and washed with PBS. The primary antibodies included CDCA3 (Proteintech, 15594-1-AP). The fluorophore-conjugated secondary antibodies used were goat anti-rabbit Alexa Fluor 488 (1: 500; Abbkine, Wuhan, China) and goat anti-mouse Alexa Fluor 549 (1: 500; Abbkine, Wuhan, China). Images were captured by confocal laser scanning microscopy (Nikon A1 + R, Japan). The fluorescence intensity was analyzed by using the ImageJ software.
2.10. Statistical Analysis. In this study, R (version 4.0.2), GraphPad Prism 8, and SPSS 20.0 software were used to analyze the data. Survival, survminer, timeROC, rms, Limma, ggplot2, pheatmap, and ClusterProfiler R package were used in this study. The significance of differences between groups was assessed by the student T test. Chisquare test was used for categorical variables, and Wilcoxon test was used for continuous data. Survival differences were calculated using Kaplan-Meier and logarithmic rank tests.

CDCA3 Is Related to Immune Infiltration.
Further, we conducted tissue microarray to try to prove the above results. Figure 3(a) shows that we performed IF staining in RCC tissue microarray. There was a significant positive correlation between CDCA3 and CD8 ( Figure 3(b)). However, our study did not observe the correlation between CDCA3 and CD4, FOXP3, CD68, and PD-1 (Figures 3(c)-3(f)). In conclusion, CDCA3 was closely related to tumor immune cells infiltration and antitumor immunity. And CDCA3 may be important for RCC risk stratification and immunotherapy guidance.

Identification of Molecular
Mechanism of CDCA3. The distribution of different genes among patients with different CDCA3 expression groups was shown in the volcano map (Figure 4(a)). The heat map showed the expression trend     Journal of Oncology of 50 upregulated genes and 50 downregulated genes with the greatest difference (Figure 4(b)). KEGG enrichment analysis showed that the related pathways were mainly concentrated in p53 signal pathway, TGF-β signal pathway, NF-κB signal pathway, and JAK-STAT signal pathway. (Figure 4(c)). GO enrichment analysis showed that its biological function was mainly enriched in spindle organization, regulation of sister chromatid segregation, and nuclear division (Figure 4(d)). These results suggest that CDCA3 mainly affects cell cycle in RCC and may regulate antitumor immune response through NF-κB axis and other important immune-related pathways.

CDCA3 Knockdown Attenuated RCC Cell Proliferation and Arrested Cell Cycle.
To further understand the effect of CDCA3 on the biological behavior of RCC, we constructed CDCA3-knockdown cell lines for functional experiments. Lentiviruses carrying CDCA3 shRNA were used to obtain CDCA3-knockdown Caki-1 and 786-O. The Western blot results showed that the expression of CDCA3 was significantly decreased in the RNAi group, indicating that the CDCA3-knockdown cell lines were successfully constructed ( Figure 5(a)). Meanwhile, the expression of CDK4, BUB3, and Cdc20 was decreased ( Figure 5(a)), which indicating cell cycle arrest. The CCK-8 assay showed that CDCA3 knockdown remarkably attenuated the cell proliferation ( Figures 5(b) and 5(c)). The ability of colony formation was notably impaired after knockdown of CDCA3 gene ( Figure 5(d)). The flow cytometric indicated that CDCA3 knockdown cause G1, S, and G2/M phase arrest ( Figure 5(e)). In general, CDCA3 expression affects cell cycle operation and cell proliferation.

Discussion
Our study suggested that CDCA3 can independently predict prognosis and affect tumor progression in RCC. CDCA3 may also be involved in the regulation of immune-related pathways, and stimulated the infiltration of immune cells, such as CD8 + T cells and Tregs. Importantly, we verified our results in vitro.
Scholars revealed that CDCA3 influence many tumor progression and treatment through a variety of pathways and is associated with poorer prognosis [8,15]. Our results also showed consistency. Patients with high CDCA3 expression had significantly worse survival and clinical stage, which was confirmed by our results and public databases. One important reason is that dysregulation of cell cycle is the basis of abnormal proliferation of tumor cells. We also confirmed that downregulation of CDCA3 blocked the G2/ M phase of cells and reduced cell proliferation ability. This directly proved that the functional localization of CDCA3 was a key regulatory protein in the cell cycle, and its abnormal expression can affect tumor progression and prognosis.
Infiltrating immune cells directly affect the occurrence, development, and treatment of tumors. It has been reported that CDCA3 is closely related to immune infiltration in hepatocellular carcinoma [16]. Our results showed that CDCA3 affected tumor infiltration of various immune cells, including CD8 + T cells in endogenous and exogenous data. Previous studies have shown that CD8 + T cells can recognize tumor-specific antigens and played a role in tumor control [17]. The high density of tumor infiltrating CD8 + T cells has been proved to be associated with a good prognosis of most cancers [18], but the infiltration of CD8 + T cells in   RCC was associated with a poor prognosis [19], this is consistent with our survival outcomes. Since immune cells hardly express CDCA3, antitumor therapy targeting CDCA3 may not cause damage to immune cells, which is a potential treatment. Our study was firstly proved that CDCA3 may be involved in the regulation of immune cell infiltration and tumorigenesis in RCC. But more importantly, the specific pathway through which CDCA3 affects immune infiltration needs further study.
As we know, immune checkpoint is a key molecule in tumor immune escape pathway. There were a lot of evidences showed that immune checkpoints were related to the benefit degree of ICIs treatment, which can be used as biomarkers for ICIs treatment [20][21][22]. Our results showed that patients with high expression of CDCA3 also expressed high levels of CTLA4 and PD-1. This initially showed that there was a close relationship between CDCA3 and immune checkpoints and further suggested that CDCA3 may participate in the immune pathway of RCC by regulating immune regulatory factors, which may be a potential target for immunotherapy. Moreover, findings suggested that TMB may predict clinical response to ICIs [23]. The neoantigen produced by TMB may be an important reason for stimulating antitumor response. In our study, we found that there was a positive correlation between CDCA3 and TMB, also suggesting that CDCA3 may predict the benefits of immunotherapy.
In summary, reactive TME is the key to immunotherapy, and CDCA3 helps to evaluate this phenomenon. Furthermore, enrichment analysis was performed to evaluate the actual molecular mechanism of CDCA3 in RCC. CDCA3 has been suggested to influence the NF-κB pathway to mediate tumor progression [12]. Our results supported this point. CDCA3 is also involved in P53 and TGF pathways. NF-κB is involved in the regulation of inflammation and innate immunity in tumor development. P53 also plays an important role in immune system. P53 mutation in cancer triggers B cell antibody response and CD8 + killing T cell response [24]. TGF-β can inhibit the proliferation, activation, and effector function of T cells. In addition, TGF-β further enhances immunosuppression in TME by promoting Tregs differentiation and destroying T cell immunity [25]. These evidences suggest that CDCA3 has a reasonable influence on TME, but more specific studies are needed to uncover the regulatory mechanisms. TME is recognized as a complex dynamic ecosystem, which is composed of malignant tumor cells, various infiltration immune cells, fibroblasts, and a variety of cytokines. In this ecosystem, immune response plays an important role in tumorigenesis and development. RCC has always been regarded as an immunogenic malignant tumor [26][27][28], and it is usually insensitive to chemotherapy and radiotherapy. Immunotherapy is regarded as another therapeutic target in addition to chemotherapy and radiotherapy [29]. Clinicians are focusing on immunotherapy to create a new era of RCC treatment, trying to break through the traditional barrier [30]. The first thing to use immunotherapy is to evaluate the immune status, which is the premise of personalized treatment. Therefore, find a biomarker that can better indi-cate the immune status and curative effect of patients, which provides an important reference for immunotherapy of RCC. When we focus on CDCA3, the problem seems to become transparent. CDCA3 has the potential to evaluate prognosis and TME and helps to hierarchically label patients at high risk. Then apply medical intervention in advance, select appropriate treatment strategies, and improve the prognosis. However, our study has its limitations. First, the specific mechanism of CDCA3 on CD8 + T cells and its influence on immunotherapy of renal cell carcinoma need to be further explored; second, we have proved that CDCA3 can block the cell cycle, but there is no further study on the biological mechanism.

Conclusion
CDCA3 can be used as an oncogene to affect the prognosis of RCC patients. Downregulation of CDCA3 causes cell stagnation in G2/M phase, promotes cell apoptosis, and reduces proliferation ability. More importantly, the immunological implications of CDCA3 have also been preliminarily evaluated. CDCA3 may participate in the regulation of immune infiltration in tumor microenvironment by affecting the expression of many immune regulatory factors and TMB, which is expected to provide valuable reference for clinical ICIs treatment. Overall, CDCA3 can be used as a biomarker to evaluate prognosis and CD8 + T cell infiltration in RCC. Targeted therapy against CDCA3 is a promising new therapeutic modality, and focusing on it may help to improve the management of therapeutic resistance in the combination of ICI and TKI, but this needs further research to confirm.

Data Availability
All public data access addresses are visible in the manuscript. Data archiving will be made available on reasonable request and all of the authors are responsible to the data.