Expression and gene regulatory network of S100A16 protein in cervical cancer cells based on data mining

S100A16 protein belongs to the S100 family of calcium-binding proteins, which is widely distributed in human tissues and highly conserved. S100 calcium-binding proteins possess broad biological functions, such as cancer cell proliferation, apoptosis, tumor metastasis, and inflammation (Nat Rev Cancer 15:96–109, 2015). The S100A16 protein was initially isolated from a cell line derived from astrocytoma. The S100A16 protein, consisting of 103 amino acids, is a small acidic protein with a molecular weight of 11,801.4 Da and an isoelectric point (pI) of 6.28 (Biochem Biophys Res Commun 313:237–244, 2004). This protein exhibits high conservation among mammals and is widely expressed in various human tissues (Biochem Biophys Res Commun 322:1111–1122, 2004). Like other S100 proteins, S100A16 contains two EF-hand motifs that form a helix-loop-helix structural domain. The N-terminal domain and the C-terminal domain of S100A16 are connected by a "hinge" linker.S100A16 protein exhibits distinct characteristics that distinguish it from other S100 proteins. A notable feature is the presence of a single functional Ca2 + binding site located in the C-terminal EF-hand, consisting of 12 amino acids per protein monomer (J Biol Chem 281:38905–38917, 2006). In contrast, the N-terminal EF-hand of S100A16 comprises 15 amino acids instead of the typical 14, and it lacks the conserved glutamate residue at the final position. This unique attribute may contribute to the impaired Ca2 + binding capability in the N-terminal region (J Biol Chem 281:38905–38917, 2006). Studies have shown an integral role of S100 calcium-binding proteins in the diagnosis, treatment, and prognosis of certain diseases (Cancers 12:2037, 2020). Abnormal expression of S100A16 protein is implicated in the progression of breast and prostate cancer, but an inhibitor of oral cancer and acute lymphoblastic leukemia tumor cell proliferation (BMC Cancer 15:53, 2015; BMC Cancer 15:631, 2015). Tu et al. (Front Cell Dev Biol 9:645641, 2021) indicate that the overexpression of S100A16 mRNA in cervical cancer(CC) such as cervical squamous cell carcinoma and endocervical adenocarcinoma as compared to the control specimens. Tomiyama N. and co-workers (Oncol Lett 15:9929–9933, 2018) (Tomiyama, N) investigated the role of S100A16 in cancer stem cells using Yumoto cells (a CC cell line),The authors found upregulation of S100A16 in Yumoto cells following sphere formation as compared to monolayer culture.Despite a certain degree of understanding, the exact biological function of S100A16 in CC is still unclear. This article explores the role of S100A16 in CC through a bioinformatics analysis. Referencing the mRNA expression and SNP data of cervical cancer available through The Cancer Genome Atlas (TCGA) database, we analyzed S100A16 and its associated regulatory gene expression network in cervical cancer. We further screened genes co-expressed with S100A16 to hypothesize their function and relationship to the S100A16 cervical cancer phenotype. Our results showed that data mining can effectively elucidate the expression and gene regulatory network of S100A16 in cervical cancer, laying the foundation for further investigations into S100A16 cervical tumorigenesis.