FAM64A is an androgen receptor-regulated feedback tumor promoter in prostate cancer

Endocrine therapy for prostate cancer (PCa) mainly inhibits androgen receptor (AR) signaling, due to increased androgen synthesis and AR changes, PCa evolved into castration-resistant prostate cancer (CRPC). The function of Family With Sequence Similarity 64 Member A (FAM64A) and its association with prostate cancer has not been reported. In our research, we first reported that FAM64A is up-regulated and positively associated with poor prognosis of patients with prostate cancer (PCa) by TCGA database and immunohistochemistry staining. Moreover, knockdown of FAM64A significantly suppressed the proliferation, migration, invasion, and cell cycle of PCa cells in vitro. Mechanistically, FAM64A expression was increased by dihydrotestosterone (DHT) through direct binding of AR to FAM64A promoter, and notably promoted the proliferation, migration, invasion, and cell cycle of androgen-dependent cell line of PCa. In addition, abnormal expression of FAM64A affects the immune and interferon signaling pathway of PCa cells. In conclusion, FAM64A was up-regulated by AR through directly binding to its specific promoter region to promote the development of PCa, and was associated with the immune mechanism and interferon signaling pathway, which provided a better understanding and a new potential for treating PCa.


INTRODUCTION
Prostate cancer (PCa) is one of the most common malignant tumors in the male urinary tract and its incidence ranks fifth among malignant tumors worldwide and second among male malignancies. According to the statistics, the total number of global cases of PCa in 2016 is 1.4 million, and the total number of deaths is 375,304 [1]. Due to PCa patients have no obvious symptoms in the early stage, the majority of them are diagnosed in the late stage, resulting in a low cure rate and higher mortality rate. For patients with advanced PCa, the current standard firstline treatment is androgen deprivation treatment (ADT), which can effectively reduce tumor burden, improve patient quality of life and prolong overall survival time. However, after ADT treatment with a median time of 18-24 months, almost all patients gradually evolved into castration-resistant prostate cancer (CRPC) [2], and 86% of CRPC eventually turned into metastatic castration-resistant prostate cancer (mCRPC) [3]. Interestingly, activated AR molecules can both enhance and inhibit the expression of genes related to PCa progression. This hormone-driven AR signal is essential for the development, differentiation, and normal function of the prostate. Correspondingly, CRPC continues to express AR and AR regulatory genes after androgen ablation therapy [4].
The human FAM64A gene (Family With Sequence Similarity 64 Member A, also known as RCS1, PIMREG, or CATS) is located on chromosome 17p13 and contains eight exons. The mRNA of human FAM64A is about 1517 bp (NM_001195228.1), which encodes a 248 amino acid protein (NP_001182157.1) [5]. FAM64A plays important biological functions in various cells by accelerating the cell cycle. Initial research found that the mutations of the FAM64A gene could regulate the size of Saccharomyces cerevisiae cells [6]. FAM64A gene was highly expressed in cardiomyocytes of mouse embryonic period, and it was proved that FAM64A could participate in the process of the cell cycle as a regulator of mitosis [7]. The expression of the FAM64A gene is abnormally regulated in a variety of tumor tissues. Current studies have found that FAM64A was remarkably highly expressed in tumor tissues and cells of patients with leukemia [7], breast cancer [8], and pancreatic cancer [9], which may be related to cell cycle disorders, but its specific molecular mechanism is still unclear.
Janus kinase signal transducer and activator of transcription (JAK-STAT) signal transduction were found in almost all immunoregulatory processes, including those involved in tumor cell recognition and tumor-driven immune escape. STAT activation of countless cytokines and growth factors induces driving events such as hematopoiesis, immune adaptation, inflammation, tissue repair, adipogenesis, and apoptosis. Most immune responses triggered by cytokines rely on STATs. This pathway was first discovered in the framework of interferon (IFN)-signaling research. The anti-tumor immune response is largely driven by STAT1 and STAT2 induced type I and II IFN and downstream program-enhanced IFN [10].
In the present study, we found that FAM64A was overexpressed in PCa tissues and cells, and its expression was positively correlated with poor prognosis in PCa patients. FAM64A expression was activated by AR through direct binding of AR to FAM64A promoter, and notably promoted the proliferation, migration, invasion, and cell cycle of PCa cells. In addition, abnormal expression of FAM64A affected the immune mechanism and IFN signaling pathway in PCa. The results indicate that the DHT/ AR-FAM64A-IFN axis may serve as a novel promising diagnostic and therapeutic target for PCa.

FAM64A is upregulated in PCa tissues and cell lines
To investigate the role of FAM64A in the progress of PCa, we first analyzed the expression level of FAM64A in 499 cases of PCa tissue and 52 cases of normal prostate tissue from the TCGA database. It was shown that the FAM64A mRNA level was significantly higher in PCa tissue compared with normal prostate tissue (Fig. 1a). Furthermore, we analyzed the correlation of FAM64A expression and clinicopathological characteristics of 499 patients with PCa. The dysregulation of FAM64A in PCa was associated with the clinical tumor stage, lymph node metastasis, and Gleason score, but not the age of patients with PCa ( Fig. 1b-d and Table 2). In addition, the Kaplan-Meier survival analysis and the Kaplan-Meier disease-free survival analysis showed a poor prognosis for PCa patients displaying high expression levels of FAM64A (Fig. 1e, f).
The above results were also identified by the IHC data. There was a significant increase of FAM64A protein in 97 cases of PCa tissues compared to the adjacent normal prostate tissues (Fig. 1i, j, Table 1). In addition, the expression of FAM64A was related to the patient's Lymph node status, Ki-67, TOPO-II, ERG, and the difference was statistically significant (P < 0.05) ( Table 3). The multivariate Cox-regression analysis showed that the patient's T stage was related to survival time, and the risk coefficients were 0.001 (Table 4). In summary, our results indicated that FAM64A expression was up-regulated in PCa tissues and cells, and was associated with PCa progression. For further study, FAM64A expression levels at both transcriptional and translation levels were detected by qPCR and Western blotting, respectively. It was revealed that FAM64A was overexpressed in PCa cells (PC-3, LNCaP, 22Rv1, and DU145) compared to the human prostatic epithelial cells (P69) (Fig. 1g, h). LNCaP and 22Rv1 cell lines were selected for FAM64A knockdown in subsequent experiments.
FAM64A promotes the proliferation, migration/invasion, and cell cycle of PCa cells To explore the biological function of FAM64A in PCa cells, we transfected FAM64A-specific siRNAs or negative control siRNAs in LNCaP and 22Rv1 cells and analyzed the efficacy of siRNAs knocking down FAM64A. The qPCR and Western blotting results demonstrated that FAM64A expression was decreased in both mRNA and protein levels (Fig. 2a). Further, we detected the role of FAM64A on cell proliferation, migration/ invasion, and apoptosis of PCa cells by CCK-8, wound-healing/transwell-invasion assays, and flow cytometric analysis, respectively. The CCK-8 results showed that the knockdown of FAM64A dramatically inhibited the proliferation of LNCaP and 22Rv1 cells (Fig. 2b). Colony formation data revealed that the growth of LNCaP and 22Rv1 cells was significantly inhibited by siR-FAM64A (Fig. 2c). The Transwell and wound-healing analysis indicated that siR-FAM64A reduced the invasion and migration abilities of LNCaP and 22Rv1 cells strikingly (Fig. 2d-f). The silence of FAM64A had no effect on the apoptosis of PCa cells (Fig. 2g). As shown in Fig. 2h and i, siR-FAM64A cells comprised a significantly lower proportion of cells in the G2-M phase compared to the control. In addition, the protein analysis showed that cycle-related protein cyclin B1 and cyclin D were positively correlated with FAM64A protein expression (Fig. 2j). These results suggested that FAM64A plays a role in promoting tumorigenesis of PCa.
DHT activated the AR to regulate the expression of FAM64A by targeting its promoter sequences To determine if FAM64A is an androgen-regulated gene, we preliminarily searched the NURSA database and found that the FAM64A mRNA level was regulated after treatment of androgen in PCa cells (Table 5). Furthermore, Pearson's correlation analysis revealed that the expression of FAM64A was positively correlated with the expression of AR in prostate cancer (TCGA: P < 0.001, R = 0.864, Fig. 3a). To verify the above result, PCa cells cultured with hormone-free CSS medium were treated with different concentrations of DHT. Then the FAM64A or AR expression was examined by qRT-PCR and Western blotting. In the AR-positive human PCa cell lines LNCaP and 22Rv1, FAM64A and AR were elevated in a dosedependent manner by low concentration of androgen at both the mRNA (Fig. 3b) and protein level (Fig. 3c). Knockdown of AR attenuated the effect of DHT on FAM64A mRNA (Fig. 3d, e) and protein expression (Fig. 3f), confirming the requirement of AR for DHT-mediated promotion of FAM64A.
To identify that FAM64A is a direct target gene of AR, we performed the JASPAR database to predict the androgen response sequence in the promoter region of the FAM64A gene (Fig. 3g). There were four potential androgen response sequences in the promoter region of the FAM64A gene (Fig. 3h). Further, the analysis of ChIP-qPCR data revealed a putative AR-binding site within the promoter region of the FAM64A gene following androgen treatment (Fig. 3i, j). To confirm that these binding sites were responsible for AR regulation, we cloned the FAM64A wild or mutant promoter upstream of luciferase and conducted reporter assays (Fig. 3k, l). As expected, the transcriptional activity of this construct was promoted by androgens in LNCaP and 22Rv1 cells. Together, these data indicate the direct transcriptional control of AR on FAM64A gene.

AR-mediated FAM64A expression regulates cell proliferation, migration, invasion, and cell cycle of PCa cells
To assess the functional significance of AR-mediated FAM64A expression in PCa cells, PCa cell lines LNCaP and 22Rv1 were treated with DHT and/or transfected by FAM64A-specific siRNA.
The qPCR results showed that the expression of FAM64A decreased by siR-FAM64A in PCa cells was reversed after the DHT treatment ( Fig. 4a, b). Further, the knockdown of FAM64A inhibits cell proliferation of PCa cells. Interestingly, exposure to DHT enhanced FAM64A-mediated PCa cell proliferation ( Fig. 4c). At the same time, clonogenic analysis confirmed that DHT promoted the colony formation capability of FAM64A-mediated PCa cells (Fig. 4d). In addition, wound healing and transwell invasion analysis revealed that DHT could reverse the inhibition of knocking down FAM64A on the migration and invasion of PCa cells (Fig. 4e, f). As a chromosomal segregating gene, FAM64A was previously confirmed to promote the cell cycle from G2 to M phase. Here exposure to DHT could accelerate the process of FAM64A-mediated cell cycle G2 to M (Fig. 4g) and the expression of cyclin B1 and cyclin D (Fig. 4h). These results suggest that DHT   activated AR to affect PCa cell proliferation, migration, invasion, and cell cycle by upregulating FAM64A expression.

FAM64A inhibits immunity and IFN-signaling pathway in PCa
To investigate the molecular mechanisms of FAM64A in process of PCa, we introduced transcriptome sequencing. We applied the GEO2R to identify differential expression analysis sequence (DEGs) regulated by FAM64A in PCa cells. 129 DEGs (74 upregulated genes and 55 downregulated genes) were identified in the siR-FAM64A cell samples and siR-NC samples. We set the cut-off criterion as P < 0.05 and [logFC] ≥ 2 (Fig. 5a). The Venn software was used to screen 50 overlapping up-regulated genes and 37 overlapping down-regulated genes (Fig. 5b, Table 6). With the aid of the RRA software package, the 87 DEGs in the integrated microarray analysis were located on the heat map (Fig. 5c). Then, we mapped the differentially expressed genes related to FAM64A to the protein interaction network through the STRING database. Key genes analysis was performed according to the topological properties of the genes in the network (Fig. 5d). The qPCR data verified the top 18 up-regulated genes and top 11 downregulated genes in transcriptome sequencing of siR-NC and siR-FAM64A group PCa cells ( Fig. 5e and Table 6). Further, the Database for Annotation, Visualization, and Integrated Discovery (DAVID) was introduced for GO and KEGGenrichment analysis to study the biological function of DEGs. DEGs are mainly rich in the multi-organism process, immune system process, response to stimulus, negative regulation of the biological process, signaling etc. (Fig. 5f and Table S1). In cellular component terms, the DEGs were mainly involved in the organelles, extracellular region, extracellular region part, membrane-enclosed lumen, and cell part. For molecular functions, the DEGs were mainly enriched in binding receptor regulator activity, molecular function regulator, protein tag, and enzyme regulator activity. Through KEGG pathway analysis of the DEGs, we could find that DEGs is mainly enriched in cell growth and death, cell motility cellular community, transport, and catabolism in cellular processes, and is mainly enriched in environment information processing, signal transduction, signaling molecules, and interaction in membrane transport ( Fig. 5g and Table S2). To explore the functional enrichment of the identified DEGs, we submitted the upregulated genes and downregulated genes to the DAVID for functional pathway analysis including GO ( Fig. 5h and Table S3). The most significantly altered pathway was a response to the virus, a response to type I IFN, a negative regulation of viral life cycle, an innate immune response, etc.
Finally, we downloaded these differential genes through the TCGA database and analyzed the linear correlation with the seq-HTSeq-FPKM-UQ value of FAM64A in PCa patients, and found that the expression levels of SNHG18, IFIT3, OAS1, PMSB9, OAS2, UBE2L6, IGFBP3, ZNF48, and the expression level of FAM64A are linearly correlated and consistent with the sequencing results (Fig. 5i). The dimerization of the IFN receptor initiates the phosphorylation of JAK1, which phosphorylates and activates STAT1 and STAT2 proteins. These proteins form a well-characterized complex IFN-stimulated gene factor 3 (ISGF3) with interferon regulatory factor 9 (IRF9). ISGF3 translocates to the nucleus and binds to IFN-stimulated response elements (ISREs) in the gene promoter, resulting in IFN-stimulated gene transcription [11]. This is a classic pathway related to type I IFN immunity signaling pathway. FAM64A expression was strongly associated with the expression genes of immunity and type I interferon by JAK-STAT-signaling pathway. Compared with the control group, the phosphorylation levels of STAT1 and STAT2 in LNCaP and 22Rv1 cells silenced by FAM64A were reduced, whereas the overall increase in the protein levels of IRF9 (Fig. 5j). Taken together, these data confirm the concept that FAM64A inhibition modulates the expression of IFN and immunologically relevant transcripts by the JAK-STAT-signaling pathway and highlights the complexity of this regulation.

DISCUSSION
The present study is the first to identify FAM64A as an oncogenic factor in process of PCa. In this study, we combine PCa cells, human PCa samples, and bioinformatics analysis to explore the expression, regulation, function, and mechanism of FAM64A in PCa through the DHT/AR axis. FAM64A modulates the expression of IFN and immunologically relevant transcripts by the JAK-STATsignaling pathway.
We found that the expression of FAM64A in PCa cell and tissue samples was higher than that in normal ones, and was related to the patient's T stage, Gleason score, lymph node metastasis, number of positive lymph nodes, overall survival, and disease-free survival by analysis of the TCGA database. In addition, we found that FAM64A promoted not only the proliferation and migration/ invasion of PCa cells but also the cell cycle, which suggests that FAM64A may be closely related to the occurrence and progression of PCa. Several studies have shown that overexpression of FAM64A accelerated the proliferation, migration, and invasion of esophageal cancer cells, breast cancer, pancreatic cancer, lymphoma, and other tumor tissues than normal tissues [7][8][9]12]. Arcangelo et al. reported that FAM64A plays a role in controlling cell proliferation and can be regarded as a sign of proliferation [13]. In glioblastoma cells, when the cells enter the S phase, FAM64A would be progressively up-regulated as the cell cycle progresses, and the FAM64A biological protein level was the highest in the G2 phase [13]. This finding reinforces the idea that FAM64A promotes PCa cell proliferation by promoting the G2 to M phase transition.
As a nuclear transcription factor and a steroid hormone receptor, AR plays an important role in the growth of PCa [14][15][16]. In canonical AR signaling, binding of androgens to AR causes the dimerization and activation of AR. AR dimer translocating into the nucleus binds to the promoters/enhancers of target genes and recruits RNA polymerase-II to initiate gene transcription [17]. In this study, it is showed that AR expression in PCa tissues is positively correlated with FAM64A. And the expression of AR and FAM64A are increased gradually following the increasing concentration of DHT. In addition, Chip-qPCR and luciferase reporter analysis verified that AR can bind to FAM64A promoter and promote its transcription. These results indicate that FAM64A is a new target gene of AR in PCa. In line with previous reports, we found that DHT treatment significantly accelerated cell proliferation, cell migration/invasion, and cell cycle progression. Advanced PCa relies on the internal and external mechanisms of a variety of cancer cells to promote tumorigenesis, evade immune surveillance and restrain IFN signaling, and actively block immune responses [18,19]. Whole-genome sequencing of androgen-binding sites in PCa cell lines also Y. Zhou et al. demonstrated that GO functional enrichment analysis includes death, immune system process, multi-biological process, and virus reproduction [20,21]. Using transcriptome sequencing, IGFBP3, CCL5, TGFb, and CD24 are related to immune surveillance, SNHG18, IFIT3, OAS1, OAS2, and UBE2L6 are related to IFN signaling. SNHG18 inhibits the nuclear and cytoplasmic transport of ENO1 to promote the occurrence and metastasis of glioma [22]. IFIT3, a member of IFN-induced protein with tetratricopeptide repeats gene family mediating a variety of protein-protein and protein-RNA interactions, is involved in various cellular regulations [23]. IFITs are reported as prognostic markers to determine the clinical outcome of many cancers, such as glioblastoma, hepatocellular carcinoma, breast cancer, and pancreatic cancer [24][25][26][27]. OAS1 and OAS2 are composed of IFN-induced antiviral enzymes, which lead to the instability of virus-derived dsRNA with RNase L function, and the enzyme function has been well characterized [28]. Ubiquitin-conjugating enzyme E2 L6 (UBE2L6), an IFNinducible ubiquitin E2, modulates ISGylation in leukemic cells to strongly interfere with the neutrophil differentiation of ATRA-treated APL cells [29]. Studies have shown that CD4+ T cells may activate STAT3 through CCL5 to affect the chemosensitivity of PCa, and the expression of MHC I increases the immunogenicity of PCa tumor cells [30]. Nuclear IGFBP3 is bound to chromatin and regulates the expression of chemokines and cytokines (including CCL5 and TGFb). These factors contribute to the establishment of an immunosuppressive tumor environment by driving the increased number of regulatory T cells in tumors [31]. CD24 was specifically overexpressed in the tumor compartment, orchestrates a novel innate immune checkpoint through interaction with the inhibitory receptor sialic acid-binding Ig-like lectin 10 (Siglec-10) on tumor-associated macrophages in primary ovarian and breast cancer samples [32]. The JAK-STAT pathway has been identified as the key node of many pathways, and mediates a large number of processes required for the occurrence and development of tumor cells. It determines the fate of T lymphocytes, and regulates their recruitment, survival, status, and eventual death [33,34]. IL17 can promote the invasion of PCa which is secreted by T cells by increasing several EMT transcription factors and MMP7. However, the collaboration between type I and type II IFN supports innate and adaptive anti-tumor immune responses, STAT1 protein plays an important role in the interconnection of these pathways [35]. Under the stimulation of type I IFN (such as IFN α and IFN β), transcription factors such as signal converters and transcription activators (such as STAT1 and STAT2) are phosphorylated, and together with IRF9, the ISGF3 trimeric complex is formed. The complex localizes to the nucleus, and its response elements (ISREs) bind to DNA and drive the expression of IFN-stimulating genes (ISGs) [11]. In melanoma and lung cancer models, impaired phosphorylation of STAT1 in response to IFN-γ results in low MHC induction rate and reduced IFN-γ sensitivity due to defective JAK signal transduction [36]. Activation of STAT has been implicated in treatment naïve and advanced PCa, the levels of activated, nuclear STAT proteins were significantly elevated in metastatic CRPCa compared with benign prostatic hyperplasia (BPH) [37]. At the same time, JAK/STAT and AR pathways have synergistic functions in tumor cells, which may be related to the progression of PC [38,39]. We have also proved that the expression of FAM64A can affect the immunity of PCa cell lines and the expression of IFN-relevant transcripts through the JAK-STAT pathway. Therefore, the activation of AR by DHT leads to the high expression of FAM64A which promotes the proliferation and invasion of PCa cell lines, FAM64A could impact the expression of immunity and IFN relevant transcripts by the JAK-STAT pathway in PCa.

CONCLUSION
In summary, FAM64A was up-regulated in both PCa tissues and cell lines, and the overexpression of FAM64A in PCa cells could promote cell proliferation and invasion, and induce cell cycle progression. Moreover, AR activated by DHT bind the promoter sequence of the FAM64A gene to transcriptionally activate FAM64A and further accelerated the function of FAM64A in PCa. In addition, the lever expression of FAM64A might participate in the process of PCa by regulating the expression of immune and IFN genes by the JAK-STAT pathway in PCa. Understanding the mechanism of FAM64A might provide us a new therapeutic target for PCa treatments.
The TCGA data was downloaded to investigate the expression of gene mRNA levels in PCa and normal prostate tissues, and the correlation of FAM64A expression and clinicopathological characteristics of patients with PCa.

Transfection and RNA silencing
PCa cells were transfected with siRNAs using Lipofectamine 3000 (Invitrogen Inc., CA, USA). FAM64A-specific siRNAs (siR-FAM64A), AR siRNAs (siR-AR), or the corresponding negative control siRNAs (siR-NC) were synthesized by GenePharma Co., Ltd (Suzhou, China). The siRNAs' sequences were as followed (actually, we have two pairs of siR-AR and siR-FAM64A, and chose one pair of them, Materials S1, Fig. S1): siR-AR forward: 5′-GAAAAUGAUUGCACUAUUGATT-3′, siR-AR reverse: 5′-UCAAUAGUGCAAUCAUUUCTT-3′; siR-FAM64A forward: 5′-GGCUCAUGCCCACCCAUTT-3′, siR-FAM64A reverse: 5′-AUGGGUGGGCAUGUGAGCCTT-3′; siR-NC forward: 5′-UUCUCCGAACGUGUCACGUTT-3′, siR-NC reverse: 5′-ACGUGACACGUUCGGAGAATT-3′. Fig. 3 DHT activated the AR to regulate the expression of FAM64A by targeting its promoter sequences. a The linear correlation analysis was a positive correlation between AR expression and FAM64A expression in prostate cancer patients by TCGA database. b, c The expression of FAM64A was accompanied by increased DHT (10 nM) induced AR at both transcription and translation levels. d-f The mRNA expression of AR and FAM64A in cells transfected with siRNA-AR and DHT (10 nM) at both transcription and translation levels. g Canonical full-ARE sequences. h The AR-binding sites in FAM64A promoter from ChIP-seq results and the prediction by website JASPAR. i, j DNA immunoprecipitated with AR antibody and IgG (as a nonspecific control) were used as templates for PCR using primers that amplified the ARbinding sites within the FAM64A promoter. k The designed mutants of key residues in the ARE4 T-test. l LNCaP and 22Rv1 cells were cotransfected with reporter plasmids (containing the wide type and mutants of ARE), Renilla luciferase reporter gene, and AR expression plasmid and then harvested to measure luciferase activity treat the data as above. *P < 0.05, **P < 0.01, ***P < 0.001.

Immunohistochemistry staining for FAM64A
Immunohistochemistry staining was performed on Single spot tissue microarrays (TMA) slides (HProA150CS01, Shanghai Outdo Biotech Company) using the FAM64A rabbit antibody (1:20, Thermo Fisher Scientific, Inc. USA) and the rabbit streptavidin-biotin detection system (Beijing Zhongshan Golden Bridge Biotechnology Co., Ltd., Peking, China). The primary PCa specimens, which were excised from patients and confirmed including different Gleason scores from 6 to 10, were selected as PCa group; and the normal tissues next to the tumor were selected as an adjacent normal group. Of note, the result was assessed according to the intensity of staining (0, 1+, 2+, and 3+) and the percentage of positive cells (0 (0-10%), 1 (10-25%), 2 (26-50%), 3 (51-75%), and 4 (76-100%)) by two experienced pathologists. Finally, the product of the fraction of staining intensity and positive rate was used as the final score. Those less than or equal to the median value were considered as the low-expression group, while those higher than the median value were considered as the high-expression group.

Wound healing assay
PCa cells were suspended and seeded into a six-well plate. After the cells attached to the bottom and formed a monolayer, linear wounds were created by scratching the center of the monolayer with a 10 μl pipette tip. The detached cells were carefully washed off with PBS thrice. Then, 2 mL of culture medium without FBS was added to each well, and the wounds were photographed at 0, 48, and 96 h under a microscope (Olympus, Japan). The width of each scratch wound was recorded and analyzed.

Transwell invasion assay
Following transfection, PCa cells were seeded onto the basement membrane matrix (EC matrix, Chemicon, Temecula, CA) present in the insert of a 24-well culture plate. Complete culture medium was added to the lower chamber as a chemoattractant. After a further 48 h, the noninvading cells and EC matrix were gently removed with a cotton swab. Invasive cells located on the lower side of the chamber were fixed with 4% paraformaldehyde and stained with crystal violet, air-dried and photographed with a microscope (Nikon Elipse Ci, Japan).

Cell cycle assay
PCa cells with siR-NC and siR-FAM64A cultured for 48 h were fixed overnight in 75% alcohol at 4°C. The fixed cells were resuspended in PBS containing PI/RNase staining buffer (BD Biosciences, San Diego, CA), and incubated for 15 min at room temperature. Cells were subjected to flow cytometric analysis of DNA content using a flow cytometer (BD Accuri TM C6 Plus, BD Bioscience, CA). The percentages of cell cycle distribution were calculated by Flowjo VX software.

Transcriptome sequencing and data processing
Total RNA was extracted using the Total RNA Extractor (Trizol) kit (B511311, Sangon, China) according to the manufacturer's protocol, and treated with RNase-free DNase I to remove genomic DNA contamination. RNA integrity was evaluated with a 1.0% agarose gel. Thereafter, the quality and quantity of RNA were assessed using a NanoPhotometer ® spectrophotometer (IMPLEN, CA, USA) and an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA). The high-quality RNA samples were subsequently submitted to the Sangon Biotech (Shanghai) Co., Ltd. for library preparation and sequencing.

Statistical analysis
All error bars in graphical data represent mean ± SD. Pearson correlation coefficients were used to determine the correlation between the expressions of two genes. Correlation analysis between FAM64A Fig. 5 The expression level of FAM64A could regulate the relevant genes of immunity and interferon by the JAK-STAT-signaling pathway in PCa. a Red dots represent genes with fold change ≥ 2 and P < 0.05, green dots represent genes with fold change ≤ −2 and P < 0.05, and the other dots represent the rest of genes with no statistically significant change in NC and siRNA-FAM64A cells by Volcano map. FC fold change. b Venn diagram of up-regulated genes and down-regulated genes. c The red color represents up-regulated genes and the green color represents down-regulated genes in Heat maps. d Construction of protein-protein interaction (PPI) network of differential expressed genes (DEGs), Red circles are upregulated genes in the PPI network, while green ones are downregulated. e The q-PCR experiment verifies that the highest log2FoldChange value of 18 up-regulated genes and 13 down-regulated genes which had been found by transcriptome sequencing to siRNA-NC and siRNA-FAM64A PCa cell lines. f, g The biological functions of the integrated DEGs were explored by GO and KEGG enrichment analysis. h The upregulated genes and downregulated genes were subjected to the DAVID for functional pathway analysis by GO enrichment. i By the TCGA database, the linear relationship of the seq-HTSeq-FPKM-UQ value of FAM64A in prostate cancer patients was analyzed. It was found that the expression level of SNHG18, IFIT3, OAS1, PMSB9, OAS2, UBE2L6, IGFBP3, ZNF48 was linearly related to the expression level of FAM64A, consistent with sequencing results. j Western blot results of key proteins expression for IFN and JAK-STATsignaling pathway upon the knockdown of FAM64A.
expression and various clinical-pathological variables was performed by χ 2 test and Fisher's exact test. All statistical analysis was performed using SPASS 19.0 and GraphPad Prism 6 software and P < 0.05 was considered as a statistically significant difference between groups.

DATA AVAILABILITY
The dataset used during this study is available from TCGA and NURSA database.