PCDH8 is a novel prognostic biomarker in thyroid cancer and promotes cell proliferation and viability

Protocadherin 8 (PCDH8), a calcium-dependent transmembrane protein in the protocadherin family, regulates cell adhesion and signal transduction. While some studies have provided indirect evidence that PCDH8 has cancer-promoting properties, this association is controversial. In particular, its involvement in thyroid cancer (THCA) remains unclear. We aimed to elucidate the role of PCDH8 in THCA using bioinformatic analysis. Subsequently, the results were experimentally validated. The analysis conducted using the R programming language and online web tools explored PCDH8 expression levels, prognostic, and clinical implications, and its relationship with the tumor immune microenvironment in THCA. Furthermore, we examined the association between PCDH8 and co-expressed genes, highlighting their involvement in several biological processes relevant to THCA. The potential of PCDH8 as a therapeutic target for this pathology was also explored. Immunohistochemical (IHC) staining was performed on samples from 98 patients with THCA, and experimental validation was carried out. PCDH8 was significantly elevated in cancer tissues and associated with poor prognosis, several clinical factors, and immune cell and checkpoint abundance. Cox regression and survival analyses, together with Receiver Operating Curves (ROC) indicated that PCDH8 was an independent prognostic factor for THCA. Furthermore, PCDH8 impacts cell viability and proliferation, promoting tumorigenesis. Also, it influences tumor cell sensitivity to various drugs. Thus, PCDH8 might be a potential therapeutic target for THCA. IHC, cell culture, MTT, and colony formation experiments further confirmed our findings. This analysis provided insights into the potential carcinogenic role of PCDH8 in THCA, as it impacts cell viability and proliferation. Thus, PCDH8 might play an important role in its prognosis, immune infiltration, and diagnosis. Supplementary Information The online version contains supplementary material available at 10.1007/s10142-024-01312-3.

Thyroid cancer (THCA), the most common endocrine cancer worldwide is typically treated with conventional modalities such as major surgery, thyroid stimulation, hormone suppression, and 131I therapy which have shown satisfactory efficacy (Verburg et al. 2017).Despite a high fiveyear survival rate exceeding 97% and a ten-year rate of 85%, 10 to 28% of patients experience relapse 10 to 27 years after treatment (Grogan et al. 2013;Schneider and Chen 2013;Riesco-Eizaguirre et al. 2016;Luo et al. 2019a, b).However, the incidence and mortality of THCA have gradually increased in the last decade (Sturgeon et al. 2016).Moreover, THCA is a chronic disease, necessitating long-term healthcare and medication (Jiang et al. 2019).However, for aggressive malignant forms, 5-20% of patients who undergo total thyroidectomy experience local recurrence or distant metastasis (Pellegriti et al. 2004; American Thyroid Association Guidelines Taskforce on Thyroid et al. 2009).Novel biomarkers are needed to enhance understanding and prevention of the disease and provide better follow-up diagnosis, leading to a tailored treatment.PCDH8, with its unique characteristics and functions, may be associated with THCA, though the relationship remains unclear.In our study, we analyzed the Cancer Genome Atlas (TCGA) database to investigate the differential expression of PCDH8 between paired paracancerous thyroid tissues and tumor tissues and used human THCA samples to validate our analysis.We employed enrichment analysis techniques to explore potential mechanisms underlying the involvement of PCDH8 in THCA tumorigenesis.Cox regression and survival analyses were performed to assess its diagnostic and prognostic value.
We investigated the role of PCDH8 in immune infiltration and its functions in the tumor microenvironment.

Public cohort
We extracted THCA mRNA expression data (FPKM) and clinical information from the Cancer Genome Atlas (TCGA) database (https:// portal.gdc.cancer.gov/), including 510 tumor and 58 normal samples.The expression matrix in FPKM format and clinical data are available as Supplementary materials (Supplement 1: RNA_seq_FPKM and Supplement 4: Clinical).Throughout the study, we adhered to TCGA's publication guidelines for the appropriate data use and reporting.

Human THCA samples and IHC
From January 2013 to December 2018, 98 human papillary thyroid carcinoma samples and their adjacent non-tumorous specimens were collected at the First Affiliated Hospital of Xi'an Jiaotong University.IHC staining was performed to examine the expression of PCDH8 as previously described (Sun et al. 2020).The anti-PCDH8 antibody was purchased from Santa Cruz Biotechnology (#sc-377348) and was utilized at a working concentration of 1:50.All samples underwent histopathological validation, confirming their diagnostic accuracy.Patients associated with these samples had not undergone any radiotherapy, chemotherapy, or I 131 treatment before surgery.The study adhered to the ethical standards of the Research Ethics Committee of The First Affiliated Hospital of Xi'an Jiaotong University and the 1964 Helsinki Declaration and its later amendments.Written informed consent was obtained from all patients included in this study.Baseline data for the 98 patients are provided as Supplement 5 Baseline of patients.

PCDH8 expression in pan-cancer and THCA
We examined PCDH8 expression in 33 types of human cancer using THCA data from TCGA and the R package.We used gene set cancer analysis (GSCA) for visualization of pan-cancer analysis.PCDH8 expression levels were compared between THCA tissues and adjacent normal tissues using unpaired and paired t-tests.Differentially expressed genes (DEGs) were identified using the limma software package (Ritchie et al. 2015) with cutoffs of log|FC|> 1.0 and p-value < 0.05.Furthermore, we compared PCDH8 expression between normal and thyroid cancer tissues using the Human Protein Atlas database (HPA) (www.prote inatl as.org) (Ponten et al. 2008).The HPA database uses antibody profiling for precise protein localization, offering access to protein expression profiles in 32 human tissues.

Genetic alterations and the mutation landscape of PCDH8 in THCA
We obtained RNA-sequencing profiles, genetic mutation data, and clinical information for THCA from the TCGA dataset.Mutation data were downloaded and visualized using the maftools package in R software.Details are available in the Supplementary information (Supplement 2).

Clinical factor analysis
We compiled the clinical characteristics of patients with thyroid cancerTHCA and analyzed their correlation with PCDH8 expression level.Additionally, we explored the relationship between PCDH8 and certain clinical tumor markers using the ggplot2 package in R.

Survival and prognosis analyses
The area under the curve (AUC) of receiver operating characteristic (ROC) curves was used to analyze the diagnostic value of PCDH8 in THCA.Utilizing Kaplan-Meier plots, hazard ratios (HR), and log-rank p-values, we analyzed the relationship between PCDH8 and survival in THCA.
For evaluating the association between progression-free interval (PFI) and clinical-pathological factors in THCA, we carried out both univariate and multivariate Cox regression analyses.ROC curves for diagnosis and nomogram were drawn using R packages (pROC and survival).All analyses were based on patients with tumors included in the public cohort.PFI is defined as the duration during which a patient survives without experiencing further deterioration of the disease after receiving treatment.The outcome indicator for PFI is based on the occurrence of either deterioration or death (Supplement 4 Clinical).

Protein-protein (PPI) network, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs
Utilizing transcriptome data of 568 patients from the TCGA-THCA database, we categorized them into high and low PCDH8 expression groups based on the median value.Differential analysis, using the limma package, identified 415 genes (Supplement 3 DEGs-logFC2).Subsequently, we carried out a GO (Gene Ontology 2006) analysis to explore specific biological properties, encompassing biological processes, cellular components, and molecular functions.KEGG enrichment analysis (Kanehisa et al. 2017) was carried out to identify pathways associated with tumor occurrence and development.Results from GO and KEGG analyses were processed using enrichplot, cluster-Profiler, and ggplot2 packages.The PPI network of PCDH8 was established from the STRING database and visualized using GeneMANIA.A statistically significant difference and enrichment were defined as p-values < 0.05 and counts over 6, respectively.

Gene set enrichment analysis (GSEA)
Furthermore, we categorized 568 THCA samples into high and low PCDH8 expression groups based on the median value.GSEA was utilized to highlight significant differences between the two expression groups and identify correlations between gene expression and biological pathways.The calculations were repeated 1000 times, considering results with p-values < 0.05 and q-values < 0.05 as statistically significant.

Immune infiltration analysis
We examined the relationship between PCDH8 expression and immune responses in THCA, focusing on the abundance of tumor-infiltrating immune cells and relevant molecules.Furthermore, we evaluated the impact of PCDH8 on immune checkpoints in THCA utilizing the ggplot2 package in R and web tools.Additionally, immune infiltration levels were evaluated using algorithms such as EPIC and MCP-counter, and the tumor mutation burden (TMB) score was calculated.These analyses collectively offer comprehensive insights into the immune landscape of THCA and its potential implications in tumor development and immunotherapy response.

LinkedOmics and drug sensitivity analysis
This study utilized LinkedOmics, a comprehensive platform for multi-omics data evaluation across various tumor types (Vasaikar et al. 2018).The LinkInterpreter module from LinkedOmics was used to investigate the PCDH8 kinase target.GSCALite (http:// bioin fo.life.hust.edu.cn/ web/ GSCAL ite/), an online tool, facilitated gene cluster analysis in tumors (Liu et al. 2018).We analyzed the relationship between PCDH8 and drug sensitivity expressed as the inhibitory concentration 50 (IC50), using 481 drugs obtained from the Therapeutics Response Portal and applying Spearman correlation.

Assessment of PCDH8 knockdown efficiency in THCA cells
After transfecting PCDH8 shRNA into TPC-1 and FTC-133 cells for 48 h or 72 h, we collected mRNA and protein from these cells to validate the knockdown efficiency of PCDH8 in these cells using quantitative reverse transcription polymerase chain reaction (qRT-PCR) and Western blot, respectively.The utilized protocols have been previously described (Sun et al. 2020).The primers for qRT-PCR were: PCDH8 forward: CAG TCC GAT ACA GCA CCT TC; PCDH8 reverse: GCT TGT GTC ACC CGA TAC TT; GAPDH forward: GGT GTG AAC CAT GAG AAG TATGA; GAPDH reverse: GAG TCC TTC CAC GAT ACC AAAG.Anti-PCDH8 and anti-GAPDH antibodies for Western blot were purchased from Santa Cruz Biotechnology (#sc-377348; 1:500) or Abcam (#ab8245; 1:2000) respectively.

MTT and colony formation assay
The PCDH8 shRNA depletion (#sc-76085) was purchased from Santa Cruz Biotechnology, Inc.Following the designated intervention, THCA cells were used for assessing cell viability and proliferation using MTT assay and clone formation assay, respectively.The utilized protocols have been previously published (Sun et al. 2020).

Statistical analysis
Statistical analysis was carried out using R and GraphPad Prism 5.For correlation analyses, we employed Spearman's test.The results were considered significant for p-values < 0.05.

Pan-cancer analysis and differential expression of PCDH8 in THCA
In our study, we initially assessed the expression of PCDH8 across 32 types of human cancer.As illustrated in Fig. 1A, B, PCDH8 expression varied among these cancer types compared to their corresponding normal tissues.Remarkably, higher expression of PCDH8 was observed in several cancers, such as bladder urothelial carcinoma, breast invasive carcinoma, colon adenocarcinoma, lung squamous cell carcinoma, pancreatic adenocarcinoma, pheochromocytoma, paraganglioma, THCA, esophageal carcinoma, thymoma, and uterine carcinosarcoma.However, PCDH8 expression did not change in cervical squamous cell carcinoma, endocervical adenocarcinoma, cholangiocarcinoma, kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, prostate adenocarcinoma, and uterine corpus endometrial carcinoma.
In terms of THCA, both paired and unpaired t-tests indicated that PCDH8 expression was significantly increased  1C, D).Based on the HPA database, their IHC results showed that the protein expression level of PCDH8 was significantly higher in thyroid cancer tissues than in normal tissues (Fig. 1E, F).Furthermore, we observed distinctive patterns and performance of PCDH8 in THCA, particularly concerning chemokines, major histocompatibility complex, and survival (Fig. 1G, H, I).These findings suggest that PCDH8 may have unique roles and specific implications in THCA, extending beyond its general expression in other cancers.

Clinical, pathological, and prognosis analyses
We analyzed the correlation between PCDH8 and several key clinical factors in THCA.These findings are summarized in Table 1.Our results indicate that the expression of PCDH8 is significantly associated with the type of primary neoplasm focus (Table 1, Fig. 2A), the PFI (Fig. 2B), history of thyroid gland disorders (Fig. 2C), methylation levels (Fig. 2D), age (Fig. 2E), race (Fig. 2F), pathologic stage (Fig. 2G), and thyroglobulin levels (Fig. 2H).However, no significant correlation was observed between PCDH8 and overall survival (OS), possibly due to the generally favorable prognosis of thyroid cancer.
Kaplan-Meier survival curve analysis revealed that patients with THCA and upregulated PCDH8 expression have a significantly poorer PFI than those with lower PCDH8 expression levels (p = 0.016, HR = 2.25 [1.16-4.37])(Fig. 3A).Thus, an increased PCDH8 expression is associated with a worse prognosis in patients with THCA.Furthermore, our investigation revealed a strong correlation between PCDH8 overexpression and an unfavorable prognosis across various clinical subgroups (Fig. 3A, Fig. S1).

Diagnostic and prognostic significance of PCDH8
As indicated by the ROC curve, PCDH8 effectively discerned thyroid tumors from healthy tissue (Fig. 3C, AUC = 0.747, 95% CI: 0.714-0.780),indicating its potential for identifying patients with THCA from healthy individuals.Favorable ROC results, with AUC values of approximately 0.7, were observed in specific subgroups (Fig. 3C).A nomogram and Forest diagram integrating clinical factors with PCDH8 expression were created, illustrating survival probabilities at 1, 3, and 5 years (Fig. 3D, E).
Cox regression analysis was conducted to evaluate the independent predictive value of PCDH8.Given the generally favorable prognosis and low mortality rate of THCA,   2).In multivariate Cox regression, M1 staging (HR = 4.247, p = 0.005) and PCDH8 expression (HR = 2.419, p = 0.035) remained significantly correlated with PFI in patients with THCA (Table 2).Therefore, high expression of PCDH8 can independently predict a poorer prognosis of THCA.

Genetic modification of PCDH8 in THCA
We mapped the genomic location of PCDH8 (Fig. 4A) to provide a comprehensive overview.Subsequently, we analyzed the overall gene mutation profile in patients with THCA, identifying nonsense mutations as the most frequent, with single nucleotide polymorphisms (SNPs) being the predominant type.C > T transitions were the most common nucleotide changes observed.On average, patients with THCA had a median of nine mutations.Notably, high mutation rates were observed in genes such as BRAF, NRAS, and TTN (Fig. 4B).PCDH8 mutations were observed in 1.4% of patients with THCA, most of them being nonsense mutations (Fig. 4C).This suggests a potential PCDH8 impact on the development and progression of the disease in a small percentage of THCA cases.Understanding the genetic landscape of THCA, including the presence of PCDH8 mutations, is of critical importance.Further research is needed to elucidate the functional consequences of these mutations and their potential role in THCA pathogenesis.

GO and KEGG analyses
To investigate the biological functions of PCDH8 in THCA, we identified the top 415 associated DEGs based on |log2 Fold change|> 1 and p < 0.05 thresholds.Subsequently, we constructed a PPI network for these DEGs using the STRING database and Cytoscape software (Fig. 6A).The PPI network provided insights into potential interactions and functional relationships among these genes, revealing possible biological roles of PCDH8 in THCA.
We conducted a GO enrichment analysis.In the biological processes, genes were mainly enriched in humoral immune response, complement activation of the classical pathway, defense response to bacterium, leukocyte migration, and humoral immune response mediated by circulating immunoglobulin.In the cellular component, genes were enriched in immunoglobulin complex, circulating, external side of plasma membrane, and anchored component of membrane.The major molecular functions were cytokine activity, interleukin-1 receptor binding, receptor-ligand activity, and endopeptidase inhibitor activity (Fig. 6B, C, D, E).KEGG pathway analysis revealed enrichment in pathways such as cytokine-cytokine receptor interaction and neuroactive ligand-receptor interaction, closely related to tumor progression (Fig. 6C, D, and E).

Identification of PCDH8-related signaling pathways using GSEA
Based on normalized enrichment scores (NES), we identified several pathways significantly associated with PCDH8 expression, such as cell cycle, p53 signaling, chemokine signaling, antigen processing and presentation, and graft versus host disease pathways (Fig. 6H).These findings suggest that PCDH8 may play a crucial role in the development and progression of THCA by participating in these metabolic pathways.

LinkedOmics analysis
Using the LinkedOmics platform, we conducted a comprehensive analysis of PCDH-8-associated genes exhibiting differential expression in THCA.In total, 19,927 genes correlated with PCDH8.The top 50 genes with significant positive and negative correlations with PCDH8 are illustrated in Fig. 6F.Further analysis revealed that PCDH8 and its neighboring genes were primarily enriched in various biological processes, including cell cycle checkpoint, DNA replication, translational initiation, DNA strand elongation, neutrophil-mediated immunity, and chromosome segregation (Fig. 6G).Moreover, KEGG pathway analysis indicated their involvement in pathways such as ribosome, bladder cancer, small cell lung cancer, interleukin (IL-17) signaling pathway, leishmaniasis, and allograft rejection (Fig. 6G).These findings shed light on the specific mechanisms through which PCDH8 may contribute to the development and progression of THCA.Notably, enrichment in pathways associated with cell viability and proliferation was consistently observed across KEGG, GO, and GSEA analyses, strongly supporting the role of PCDH8 in THCA through its involvement in regulating cell viability and proliferation.

PCDH8 kinase targets
The LinkedOmics database was utilized to identify potential kinase targets of PCDH8 in THCA.The identified network included NUAK family SNF1-like kinase 1 (NUAK1), glycogen synthase kinase 3 beta, mitogen-activated protein kinase-activated protein kinase 2, and polo-like kinase 1 (PLK1).GeneMANIA was used to construct a PPI network and determine specific functions and potential regulatory mechanisms of these kinases and PCDH8 (Fig. 7A).The analysis revealed their involvement in the canonical wingless-type MMTV integration site family (Wnt) signaling pathway, mitotic nuclear division, regulation of mRNA catabolic process, regulation of RNA stability, G2/M transition of the mitotic cell cycle, and microtubule cytoskeleton organization involved in mitosis (Fig. 7A).
Additionally, we assessed the survival and clinical relevance of PCDH8-regulated kinases in THCA.Interestingly, all four kinases exhibited a significant association with poor prognosis in THCA (Fig. 7B, D).Specifically, NUAK1 was deeply involved in pathways such as RTK, PI3K-AKT, and RAS-MAPK, known to promote THCA development (Fig. 7C).These findings support the hypothesis that PCDH8 is involved and regulates THCA progression through downstream signaling pathways involving these kinases.

Impact of PCDH8 on drug sensitivity
Using data from the Cancer Drug Sensitivity Genomics database, we investigated the relationship between PCDH8 expression and the IC50 values of various drugs.Our analysis revealed that reduced PCDH8 expression was significantly associated with the effectiveness of six molecules: AT13387, BRD-K33514849, CD-437, tacedinaline, tivantinib, and tozasertib (Fig. 7E).These findings suggest that PCDH8 is negatively regulated by these molecules, highlighting its potential as a promising therapeutic target for THCA.Further research and development of drugs targeting PCDH8 could improve treatment outcomes for patients with THCA.

THCA samples and IHC
We collected 98 human THCA samples and their adjacent non-tumor samples.Employing IHC staining, we demonstrated that PCDH8 expression was substantially higher in THCA samples compared to the adjacent non-tumorous specimens (Fig. 8, Fig.S2).These results are consistent with the other findings of the present analysis.

MTT and colony formation assay
Given the significant impact of PCDH8 on cell viability and proliferation pathways as indicated by KEGG, GO, and GSEA analyses, we performed MTT and colony formation assays to further explore its mechanism in THCA.Notably, the knockdown of PCDH8 in TPC-1 and FTC-133 cells substantially decreased cell viability and proliferation (Fig. S3, Fig. 9A, B).Thus, PCDH8 promotes tumor cell proliferation and viability in THCA, potentially being an important driver of THCA progression.

Discussion
THCA, along with breast, lung, colorectal, and uterine cancers, is a frequently occurring malignant tumor with an increasing annual incidence (Siegel et al. 2020).Although patients with THCA typically have a favorable prognosis, those who experience recurrence or distant metastasis often face more adverse outcomes.Consequently, there is an essential need to identify patients with THCA early, a necessity for both researchers and clinicians.PCDHs play important roles in tumor cell adhesion, proliferation, apoptosis, and migration, potentially serving as prognostic biomarkers and therapeutic targets in THCA (Liu et al. 2019;Luo et al. 2019b;Vega-Benedetti et al. 2019;Fujiki et al. 2023).Exploring the role of PCDHs in THCA could provide valuable insights for effectively identifying and managing THCA.
PCDH8, a member of the PCDH family, encodes an integral membrane protein believed to contribute to cell adhesion.Its oncogenic role has been investigated in various types of cancer.In esophageal squamous cell carcinoma, PCDH8 is regulated by miR-200c and functions as a regulator of metastasis (Yu et al. 2020).In gastric cancer, PCDH8 promotes invasion and metastasis through enhanced interaction with extracellular matrix receptors, possibly through upregulation of laminin subunit gamma-2 (Lin et al. 2018).In glioma, PCDH8 might function as a biomarker for early detection and prognosis, its expression level being directly associated with glioma progression (Zong et al. 2017).However, the involvement of PCDH8 in THCA remains underexplored.Investigating its role in THCA could provide valuable insights into its potential as a biomarker or therapeutic target for this specific cancer type.
In this study, we conducted a comprehensive analysis of the role of PCDH8 in THCA.Pan-carcinoma analysis revealed a distinct difference in PCDH8 expression between tumor and noncancerous samples in THCA.This was validated through unpaired and paired t-tests using public databases.Additionally, staining results from the HPA database supported these findings, indicating that PCDH8 likely plays a crucial role in THCA development.These findings align with those of Lin et al., revealing that PCDH8 overexpression is associated with malignant features and poor prognoses (Lin et al. 2018).Additionally, we explored the association between PCDH8 expression and clinical factors in THCA, revealing significant correlations with age, primary neoplasm focus type (including multifocal and unifocal), PFI, and a history of thyroid gland disorders.Patients with multifocal THCA were more likely to experience recurrence or persistence compared to those with unifocal THCA.This is consistent with the results of Ting Zhang et al. indicating that multifocality is associated with tumor aggressiveness and plays a crucial role in tumorigenesis and cancer progression (Zhang et al. 2021).While the 2015 update of the ATA guidelines did not consider multifocality as a risk factor, these newer studies have shed light on its significance.
Cox regression analysis identified PCDH8 as an independent prognostic factor for predicting PFI in patients with THCA.The ROC curve analysis demonstrated high sensitivity and specificity, indicating the potential of PCDH8 as a diagnostic and prognostic biomarker in THCA.Further research on the underlying mechanisms and functional roles of PCDH8 in THCA could lead to novel therapeutic strategies for this disease.Tumor microenvironment immune cells and molecules play a crucial role in various cancers (Bai et al., 2023).Herein, we observed a negative correlation between PCDH8 expression and immune infiltration.The degree of immune cell infiltration strongly impacts patient prognosis in solid tumors (Ino et al. 2013;Schneider et al. 2018).Our findings revealed significant associations between PCDH8 and various immune cells, molecules, and checkpoints, including Th1 cells, aDC,, T cells, eosinophils, NK cells, Tgd cells, CD56 bright cells, pDC, cytotoxic cells,, cytotoxic T-lymphocyteassociated protein 4, Bcl-2, tumor growth factor beta 2 (TGFβ, PD-L1, CD276, VSIR and some chemokines.Thus, the PCDH8 gene affects immune cells and related molecules in the tumor microenvironment, and its overexpression is correlated with poorer prognosis in THCA.Specifically, the positive correlation of PCDH8 with negative immune checkpoints (CTLA-4, PD-L1, CD276) and negative regulator cytokines such as TGF-β, known to reduce immune infiltration (Krummel and Allison 1995;Keir et al. 2008), suggests a potential role in immune escape.This mechanism contributes to THCA progression.The expression of PD-L1 in thyroid cancer has been extensively studied in prior research (Ulisse et al. 2019).In summary, PCDH8 may be involved in thyroid cancer by inhibiting immune infiltration and promoting immune escape.
Enrichment analysis between high and low PCDH8 expression groups was carried out to uncover its potential mechanisms in THCA.GO analysis revealed significant enrichment in pathways related to humoral immune response, complement activation of classical pathway, defense response to bacteria, leukocyte migration, humoral immune response mediated by circulating immunoglobulin, immunoglobulin complex, cytokine activity, and interleukin-1 receptor binding.These pathways align with the role of PCDH8 in the tumor microenvironment and immune response.
need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 1
Fig. 1 PCDH8 was differentially expressed in THCA and many other cancers.(A, B) The PCDH8 mRNA expression in tumor and normal thyroid tissues.(C, D) The PCDH8 mRNA expression in paired THCA samples and unpaired THCA samples.(E, F) The PCDH8

Fig. 4
Fig. 4 Investigation of the mutation landscape of PCDH8 (A) The location of PCDH8 in the genome.(B, C) The alteration landscape of PCDH8 in THCA

Fig. 5 ◂Fig. 6
Fig. 5 Immune infiltration analysis.(A) The association between 12 immune cell types and PCDH8 mRNA expression based on ssGSEA in THCA (single sample gene set enrichment analysis).(B) The correlation between PCDH8 expression and immune chemokines in THCA.(C) The association of PCDH8 expression with TGFB2, ICAM1 and BCL-2 of THCA patients.(D) The generally correlation between 24 immune cells and PCDH8 based on Estimate algorithm in THCA.(E) Heat map of single gene co-expression between PCDH8 and some vital immune checkpoints.(F) TMB (tumor mutation burden) of PCDH8 in THCA.(G, H) The association analysis of 3 immune cell types with PCDH8 mRNA expression based on EPI and MCP algorithm.MCP, a method based on marker gene sets to quantify tumor immune cells, fibroblasts, and epithelial cells.EPI, an immune algorithm that provides the relative components of cells.*p < 0.05; **p < 0.01; ***p < 0.001 ◂

Fig. 8
Fig. 8 The expression of PCDH8 in THCA samples and matched adjacent non-tumorous specimens.The magnification of the left is × 0.4, and scale bars = 200 μm.The magnification of the right

Fig. 9
Fig. 9 Depletion of PCDH8 obviously mitigated the proliferation of THCA cells.Knockdown of PCDH8 in THCA cells was achieved by transfection of PCDH8 shRNA, then the cell viability (A) and cell proliferation (B) were assessed by MTT and colony formation assay respectively.**p < 0.01, ***p < 0.001