Potential functions and mechanisms of lysine crotonylation modification (Kcr) in tumorigenesis and lymphatic metastasis of papillary thyroid cancer (PTC)

Objectives: To examine the putative functions and mechanisms of lysine crotonylation (Kcr) during the development and progression of papillary thyroid cancer (PTC). Methods: Samples of thyroid cancer tissues were collected and subjected to liquid chromatography-tandem mass spectrometry. Crotonylated differentially expressed proteins (DEPs) and differentially expressed Kcr sites (DEKSs) were analyzed by Motif, dynamic expression model analysis (Mfuzz), subcellular localization, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation, Go Ontology (GO) annotation, and protein-protein interaction analysis (PPI). Validation was performed by immunohistochemistry (IHC). Results: A total of 262 crotonylated DEPs and 702 DEKSs were quantitated. First, for the tumor/normal comparison, a dynamic expression model analysis (Mfuzz) of the DEKSs revealed that clusters 1, 3, and 4 increased with the progression of thyroid cancer; however, cluster 6 showed a dramatic increase during the transition from N0-tumor to N1-tumor. Furthermore, based on GO annotation, KEGG, and PPI, the crotonylated DEPs were primarily enriched in the PI3K-Akt signaling pathway, Cell cycle, and Hippo signaling pathway. Of note, crosstalk between the proteome and Kcr proteome suggested a differential changing trend, which was enriched in Thyroid hormone synthesis, Pyruvate metabolism, TCA cycle, Cell cycle, and Apoptosis pathways. Similarly, for the LNM comparison group, the DEKSs and related DEPs were primarily enriched in Hydrogen peroxide catabolic process and Tight junction pathway. Finally, according to The Cancer Genome Atlas Program (TCGA) database, the differential expression of Kcr DEPs were associated with the prognosis of thyroid cancer, indicating the prognostic significance of these proteins. Moreover, based on the clinical validation of 47 additional samples, Kcr was highly expressed in thyroid tumor tissues compared with normal tissue (t = 9.792, P < 0.001). In addition, a positive correlation was observed between Kcr and N-cadherin (r = 0.5710, P = 0.0015). Moreover, N-cadherin expression was higher in the relatively high Kcr expression group (χ2 = 18.966, P < 0.001). Conclusions: Higher Kcr expression was correlated with thyroid tumorigenesis and lymphatic metastasis, which may regulate thyroid cancer progression by Pyruvate metabolism, TCA cycle, Cell cycle, and other pathways.