Abstract
Objectives
Evidence from previous documents points to the involvement of the keratin 15 (KRT15) gene in the modulation of signaling networks governing cancer cell death, survival, proliferation, migration, invasion, and metastasis. Nonetheless, its relevance in pan-cancer studies and the precise molecular mechanisms involved in lung cancer remain poorly understood. To comprehensively investigate the clinical relevance of the KRT15 gene in human lung cancer and a diverse array of cancers.
Methods
A comprehensive investigation of the clinical relevance of the KRT15 gene in a diverse array of human tumors was conducted. In light of this, the study also examined the possible link between the KRT15 gene and tumor immunogenic features. Based on the outcomes of pan-cancer analysis, we selected lung adenocarcinoma (LUAD) as the specific tumor type for an in-depth investigation into KRT15-induced signaling pathways and intercellular communications contributing to tumor progression.
Results
According to our study, KRT15 may hold significance as a newly identified biomarker, potentially contributing to both prognostic evaluation and immunotherapeutic targeting across various cancer types. Significantly, KRT15 was hypothesized to function as a guiding marker gene, holding potential for clinical prognostication and personalized tumor-specific therapies in LUAD. Suppression of KRT15 significantly impaired lung cancer cells’ growth, migration, invasion, and survival.
Conclusions
KRT15 possesses the potential to be an innovative biomarker and therapeutic target, playing a role in predicting the prognosis and treatment response of LUAD patients.
Introduction
As one of the most widespread malignancies globally, lung cancer has become the foremost cause of both fatalities and illness [1]. The early signs of lung cancer are relatively harmless, contributing to the challenge of detecting the disease at its initial stages [2], 3]. Advanced lung cancer is characterized by a high metastatic rate and an overall poor prognosis [4], [5], [6]. Lung cancer continues to have a high overall 5-year mortality rate, estimated at 78 % [7], 8]. Approximately 15 % of lung cancer cases are classified as small cell lung carcinoma (SCLC), with the remaining 85 % falling under non-small cell lung carcinoma (NSCLC) [9]. Within the category of NSCLC, there exist three subtypes: large cell carcinoma, lung squamous cell carcinoma (LUSC), and adenocarcinoma (LUAD) [10], [11], [12]. LUAD accounts for approximately 40 % of all lung cancers and is characterized as a malignant epithelial tumor with adenoid differentiation and/or mucous secretion [13]. Despite improvements in diagnosis and therapy in recent years, the overall survival (OS) outlook for individuals with LUAD remains suboptimal [13]. Currently, the complete pathogenesis of LUAD remains unclear. Consequently, the limited understanding of the fundamental processes involved in LUAD hinders the improvement of treatment outcomes. Thus, the paramount objective is to make significant progress in understanding tumor initiation and identifying novel biological markers to optimize prognosis.
Keratin is one of a family of structural fibrous proteins also known as scleroproteins. Ranging between 40 and 76 kDa, keratins (KRTs) are classified as structural proteins [2], 14]. As the primary building blocks of intermediate filaments (IFs) in the intracytoplasmic cytoskeleton of epithelial and endothelial cells, they play a major role in stress protection, cell integrity and structure, and protein targeting [15]. By anchoring into electron-dense desmosomal plaques and forming interconnections with other IFs, KRTs play a pivotal role in providing tensile strength to cellular components [15]. Due to their atypical expression patterns observed in different cancer types, KRTs serve as valuable biomarkers for discriminating between diagnoses and evaluating metastatic status [14]. In light of recent investigations, it is now evident that KRTs expressed in cancer cells possess a dual function; apart from their role as epithelial marker proteins, they act as mediators, engaging in interactions with diverse proteins to modulate signaling networks governing cell death [16], survival [9], proliferation [17], migration [18], invasion [18], and metastasis [19]. Chromosomally located at 17q21.2, keratin 15 (KRT15), commonly referred to as cytokeratin 15, is classified as one of the members in the KRT gene family [20]. KRT15 has been previously documented as a marker specific to skin stem cells within the hair follicle bulge [21], and it is known to have a vital function in regulating epidermal homeostasis [22]. A close link between KRT15 and tumorigenesis has also been demonstrated in recent research [23]. The upregulation of KRT15 was detected in squamous-cell carcinoma samples [24] and has been associated with an unfavorable prognosis in colorectal cancer [25]. Additionally, increased KRT15 expression is noted across various cancers including breast cancer [26], urothelial cell carcinomas [27], ameloblastoma [28], and hepatocellular carcinoma [29], while decreased expression is observed in oral squamous neoplasms [30], prostate tumors [31], and gastric cancer [32]. Further studies have linked KRT15 to breast cancer progression and its potential role as an independent prognostic factor [26], 33]. Additionally, differential expression in gastric cancer, associated with prognosis, was reported by Zhang et al. [32]. In addition, there is evidence to suggest that stem cells with KRT15-positive surface markers could become a source of specific cancer types under certain conditions [34]. Despite extensive research, the precise involvements of KRT15 in LUAD have remained largely unknown.
For the purpose of this study, we acquired RNA-sequencing (RNA-seq) data for LUAD and pertinent clinical patient data from TCGA and GEO databases to explore the potential prognostic value of KRT15. Utilizing these data, we conducted differential gene expression and functional enrichment analyses to pinpoint the genes and functional pathways exhibiting significant associations with KRT15 expression. Subsequently, a comprehensive analysis of tumor immune infiltration was conducted to investigate the potential correlation between KRT15 expression and the presence of immune cells. Furthermore, we conducted an assessment of the connection of KRT15 to clinicopathological and demographic parameters, followed by the analyses of survival and clinical subgroup prognosis, to ascertain the prognostic relevance of KRT15 in LUAD. Lastly, to strengthen the results, biological verification was carried out in vitro utilizing NCI-H1573 and PC-14 lung cancer cell lines. The findings from integrative bioinformatics and in vitro verifications suggest that KRT15 possesses the potential to be an innovative biomarker and therapeutic target, playing a role in predicting the prognosis and treatment response of LUAD patients.
Materials and methods
Analysis of KRT15 in pan-cancer
We conducted a comparison of KRT15 expression differences between the normal and tumor tissue samples, relying on data from The Cancer Genome Atlas (TCGA) and the Cancer Cell Line Encyclopedia (CCLE) databases, as well as the normal tissues sourced from the Genotype-Tissue Expression (GTEx) database. For the analysis of these data, we employed an online software HOME for Researchers (https://www.home-for-researchers.com/static/index.html#/). Additionally, to further clarify the expression and the distribution of KRT15 in lung cancer tissues at the translational level, the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/) was also applied.
Evaluation of the survival outcomes
For the evaluation of the OS and progression-free survival (PFS), the clinical data extracted from the TCGA database was utilized. In addition, to display the connections of KRT15 expression to the prognoses of patients with different kinds of tumors, the forest plots were employed. Lastly, we carried out the Kaplan-Meier (KM) curves to calculate and analyze the survival outcomes of the patients with different tumors.
Assessment of immune infiltration, immune checkpoint, and the status of microsatellite instability (MSI) and tumor mutation burden (TMB)
To further evaluate the immune infiltration, immune checkpoint, and the status of MSI and TMB, the TCGA database was utilized to obtain the clinical data and relevant RNA-sequencing expression profiles of KRT15. Then, the relationships of KRT15 expression with the infiltration of immune cells in different cancer patients were investigated. Next, the connections of KRT15 expression to the genes associated with immune checkpoints, such as SIGLEC15, IDO1, CD274, HAVCR2, PDCD1, CTLA4, LAG3, and PDCD1LG2, were examined. Additionally, based on the previous publications [35], 36], the TMB and MSI were subsequently analyzed. All of the analysis was carried out using the online software HOME for Researchers.
The analysis of KRT15 expression in lung cancer using single-cell RNA sequencing (scRNA-seq) data
For the analysis of the scRNA-seq to explore the expression and co-expression of KRT15 with other LUAD-associated genes, the Cancer SEA (http://biocc.hrbmu.edu.cn/CancerSEA/) online tool was applied.
Cells and cell culture
Human NCI-H1573 and PC-14 lung cancer cells were provided by the China Center for Type Culture Collection (CCTCC). Fetal bovine serum (FBS, 10 %)- and penicillin-streptomycin (1 %)-contained RPMI 1640 medium was utilized to culture the cells in a humid chamber containing 5 % CO2 at 37 °C.
Measurement of cell growth using CCK-8 assay
Transfected NCI-H1573 and PC-14 cells were seeded in a 96-well cell culture plate at a density of 3000 cells per well, with 0.1 mL of medium in each well, and incubated at 37 °C. At 2-h intervals following the incubation with a 10 % CCK-8 solution, OD values at 450 nm were recorded every 24 h for a total of 6 days.
Silencing of KRT15 expression using siRNA transfection
Genechem (Shanghai, China) provided the siRNA-KRT15 (siKRT15) and siRNA-control (siNC) plasmids. The Silencer siRNA construction kit (Qiagen, The Netherlands) was utilized to form siRNA transfection complexes, which were then incubated with the cells along with RNAiFect reagents (Qiagen, The Netherlands). The cells were subjected to a 48-h transfection period, after which they were used in the experiments.
Detection of target genes using RT-qPCR
Extraction of total RNA was carried out using the RNeasy Mini Kit (QIAGEN). The template for reverse transcription (RT) was 1 μg of total RNA, and M-MLV Reverse Transcriptase (Promega) was used for the process. Using the LightCycler 480 machine (Roche) and LightCycler 480 SYBR Green I Master Mix (Roche), quantitative real-time PCR (qPCR) reactions were conducted in triplicate. Normalization of cDNA expression was achieved by comparing it to ACTB, and each reaction was replicated with at least three independent biological samples.
Western blot measurement
Transfected cells were used to extract total cell protein using a cell protein extraction kit. The extracted proteins were separated using a 10 % SDS-PAGE gel and subsequently transferred to a PVDF membrane (Millipore, USA). Then, sealing the membrane with 5 % skim milk took place at room temperature for 1 h. Overnight incubation with the primary antibody was followed by treatment with secondary antibodies and development.
Evaluation of cell migration and invasion using transwell assay
To assess migration and invasion, transwell experiments were conducted using a transwell chamber with 8 μm pore size (Corning), with options for either uncoated or pre-coated conditions using Matrigel Matrix (BD Bioscience, USA). Following a 24-h transfection of NCI-H1573 and PC-14 cells, trypsin digestion was carried out to harvest the cells by centrifugation. The cells were then resuspended in FBS-free medium at a final density of 5 × 104 cells/mL. Afterward, the upper chamber received 250 µL of cell suspension, and the lower chamber was supplemented with 700 µL of medium containing 10 % FBS. After the incubation period of 36 h, the cells in the lower chamber were fixed with 4 % formaldehyde for 20 min and stained with 0.5 % crystal violet (Beyotime) for 15 min. Several randomly designated areas were observed under a light microscope (Olympus) to count the number of cells.
Statistical analysis
Mean ± standard deviation (S.D.) was calculated and presented by employing GraphPad Prism (GraphPad Prism Software 9.3.1, San Diego, CA) and SPSS (IBM SPSS version 23.0.2) software to process the data. Two independent samples were presented as mean ± SD and assessed by Student’s t-test, one-way ANOVA or two-way ANOVA with Tukey post hoc test for multiple comparisons. Additionally, a comparison between LUAD tissue and precancerous tissue was performed using paired t-tests. One-way ANOVA was applied for the comparison of data among multiple groups. Significance was set at p<0.05, p<0.01, p<0.001, and p<0.0001, and each experiment was repeated at least three times.
Results
Evaluation of the expression and prognosis of KRT15 in pan-cancer
Initially, to comprehensively study the expression pattern of KRT15 in various cancer types, we resorted to the TIMER2.0 software for data analysis, utilizing the TCGA database. Analysis of the data demonstrated a significant upregulation of KRT15 expression in UCEC, THCA, STAD, PCPG, PAAD, LUSC, LUAD, KIRP, KIRC, ESCA, COAD, CHOL, and CESC (Figure 1A and B). In contrast, BRCA, HNSC, KICH, and PRAD exhibited a down-regulated level of KRT15 expression (Figure 1A and B). Next, a comprehensive examination was undertaken to explore the plausible relationship between KRT15 expression and the OS rate across all 31 distinct cancer categories available in the TCGA repository. Based on the forest plots, a noteworthy correlation was observed between elevated KRT15 expression and unfavorable overall survival (OS) in BRCA, KIRC, PAAD, LUAD, and UCEC (Figure 1C). Furthermore, it was observed that heightened KRT15 expression correlated significantly with shorter PFS in CESC, KIRC, LGG, LUAD, LUSC, and PAAD (Figure 1D). Taken together, the findings presented above provide collective evidence of elevated KRT15 expression in numerous malignancies, indicating its plausible oncogenic function in these specific cancer contexts.

Evaluation of the expression and prognosis of KRT15 in pan-cancer (A, B) representative plots depicting the expression of KRT15 in multiple cancer types from TCGA utilizing TIMER2.0 (A) and GEPIA (B). *, **, and *** represent the p-values less than 0.05, 0.01, and 0.001 respectively. The samples collected from healthy individuals were used as the normal controls. (C, D) representative forest plots representing the connections of KRT15 expression to OS (C) and PFS (D) in multiple cancer types from TCGA.
Evaluation of the connection between KRT15 expression and immune infiltration, immune checkpoint, TMB, and MSI status in pan-cancer
Initially, an investigation was conducted to explore the plausible association between KRT15 expression and the extent of immune infiltration across various types of cancers. In order to examine the connection between KRT15 expression and immune-infiltrating cells, the researcher employed the TIMER approach. A noteworthy association between KRT15 and infiltrating CD8+ T cells, CD4+ T cells, neutrophils, myeloid DCs, macrophages, and B cells was observed across multiple cancer types, including BLCA, COAD, ESCA, LUSC, MESO, READ, PRAD, STAD, TGCT, and THYM, as depicted by the data (Figure 2A). Following that, an investigation was carried out to examine the association between KRT15 expression and the expression levels of key immune checkpoints (TIGIT, SIGLEC15, PDCD1LG2, PDCD1, LAG3, HAVCR2, CTLA-4, and CD274) in pan-cancer. Based on the findings, a substantial association was evident between KRT15 expression level and one or more of these markers, except for UVM, READ, LAML, GBM, and BLCA (Figure 2B). Furthermore, an investigation was conducted to assess the correlation between KRT15 expression level and TMB. In PAAD, STAD, and READ, we identified a positive correlation between KRT15 expression level and TMB, whereas in PRAD, COAD, ESCA, and BRCA, a negative correlation was observed (Figure 2C). Lastly, a comprehensive analysis of the relationship between KRT15 expression level and MSI status was conducted, revealing a compelling positive link between KRT15 expression level and MSI status in TGCT, UCEC, LUSC, STAD, and SARC (Figure 2D).

Evaluation of the connection between KRT15 expression and immune infiltration, immune checkpoint, TMB, and MSI status in pan-cancer (A, B) representative heat maps depicting the relationships of KRT15 expression with the infiltrated immune cells (A) and key molecules of immune checkpoints (B) in multiple cancer types. *, **, and *** represent the p-values less than 0.05, 0.01, and 0.001 respectively. (C, D) Representative forest plots showing the associations of KRT15 expression with the TMB (C) and MSI status (D) in multiple cancer types.
Evaluation of the connection between KRT15 expression to clinicopathological characters in LUAD
Next, a thorough examination was conducted to explore the connections of KRT15 expression to clinicopathological characters. A considerable upregulation of KRT15 expression levels was detected in LUAD samples in contrast to normal samples (Figure 3A). Additionally, in comparison to the normal samples, the expression of KRT15 was significantly upregulated in both male and female LUAD patients (Figure 3B), the LUAD patients with different races (Figure 3C), and the patients have different smoking histories (Figure 3D). However, no significant differences in KRT15 expression were observed between the patients with different genders (Figure 3B), races (Figure 3C), and smoking histories (Figure 3D). Furthermore, our findings also observed a substantial link of the KRT15 expression to the tumor stages (Figure 3E), tumor size (Figure 3F), node metastasis (Figure 3G), and distant metastasis (Figure 3H). A high level of KRT15 expression in LUAD patients represented a large tumor size (Figure 3F), high N metastasis (Figure 3G), and high distant metastasis (Figure 3H).

Evaluation of the connection between KRT15 expression to clinicopathological characters in LUAD (A) demonstration of KRT15 expression levels is provided for LUAD tissues and matched normal tissues by a violin plot. (B, C) Representative violin plots showing the relationship of KRT15 expression in LUAD patients with different genders (B) or races (C). The samples collected from healthy individuals (tumor free) were used as the normal controls. (D) The connections of KRT15 expression to LUAD patients’ smoking history were shown. The LUAD patients without smoking were used as the normal controls. (E) Representative violin plot showing the expression KRT15 in LUAD patients with different stages. (F–H) Representative violin plots representing the correlations of KRT15 expression with the tumor size (F), node metastasis (G), and distant metastasis (H). **** represents the p-value less than 0.0001. The samples collected from healthy individuals (tumor free) were used as the normal controls.
Relationship of KRT15 expression level with prognostic in LUAD
The prognostic significance of KRT15 in patients with different cancers was explored first through an analysis of data extracted from the TCGA database. By adopting the median value of KRT15 as the dividing threshold, the cancer patients were allocated into high- and low-expression categories. The results indicated a significant relationship between raised KRT15 expression levels and unfavorable OS rates in patients with KIRC (Figure 4A), PAAD (Figure 4B), SKCM (Figure 4C), and THYM (Figure 4D), along with unfavorable PFS rates in patients with KIRC (Figure 4E) and PAAD (Figure 4F). Additionally, we also corroborated that the KRT15 expression exhibited an inverse correlation with both OS (Figure 4G–I) and PFS (Figure 4H) in lung cancer based on data from the GEO database, such as GSE30219 (Figure 4G and H) and GSE31210 (Figure 4I).

Relationship of KRT15 expression level with prognostic in LUAD (A–D) representative KM plots showing the associations of KRT15 expression with the OS of patients with KIRC (A), PAAD (B), SKCM (C), and THYM (D). (E, F) Representative KM plots depicting the associations of KRT15 expression with the PFS of patients with KIRC (E) and PAAD (F). (G, H) The associations of KRT15 expression with the OS (G) and PFS (H) of lung cancer patients from GSE30219 database. (I) The associations of KRT15 expression with the OS (I) of lung cancer patients from GSE31210 database.
Single-cell analysis of KRT15 expression profile and the relationship with functional state in lung cancer
In light of the complex composition of tumor cells, single-cell transcriptome sequencing emerges as a vital tool for the comprehensive evaluation of diverse cancer cells, stromal cells, endothelial cells, and immune cells. Utilizing the TISCH website (http://tisch.comp-genomics.org/home/), the distribution of KRT15 expression in tumor cells was evaluated. Analysis of the GEO databases EMTAB6149 (Figure 5A–C), GSE117570 (Figure 5E–G), and GSE148071 (Figure 5I–K) demonstrated the predominant distribution of KRT15 expression in malignant cells, as opposed to its expression in immune cells, stromal cells, and other cell types. Meanwhile, obvious co-distribution of KRT15 expression with EGFR, a key molecular during the initiation, development, and progression of lung cancer, was also observed in EMTAB6149 (Figure 5C and D), GSE117570 (Figure 5G and H), and GSE148071 (Figure 5K and L) databases. Next, by accessing the CANCER SEA website (http://biocc.hrbmu.edu.cn/CancerSEA/), we validated KRT15 expression in single cells and investigated its relationship with the functional state of tumors. By evaluating data from GSE69405 (Figure 5M) and E-MTAB-6653 (Figure 5N-P), we observed a pronounced association between KRT15 expression levels and vital cellular processes, including the cell cycle, DNA damage, DNA repair, invasion, and DNA damage response.

Analysis of KRT15 expression with single-cell sequencing data and its association with the functional state in LUAD (A) representative single cell plot showing the different cell types of lung cancer samples from GEO database EMTAB6149. (B, C) Representative violin plot (B) and single cell plot (C) depicting the distribution of KRT15 in different cell types of lung cancer samples from GEO database EMTAB6149. (D) The distribution of EGFR in the different cell types of lung cancer sample from GEO database EMTAB6149. (E) Representative single cell plot showing the different cell types of lung cancer samples from GEO database GSE117570. (F, G) Representative violin plot (F) and single cell plot (G) depicting the distribution of KRT15 in different cell types of lung cancer samples from GEO database GSE117570. (H) The distribution of EGFR in the different cell types of lung cancer sample from GEO database GSE117570. (I) Representative single cell plot showing the different cell types of lung cancer samples from GEO database GSE148071. (J, K) Representative violin plot (J) and single cell plot (K) depicting the distribution of KRT15 in different cell types of lung cancer samples from GEO database GSE148071. (L) The distribution of EGFR in the different cell types of lung cancer sample from GEO database GSE148071. (M−P) The relationships of KRT15 expression with indicated functional status in the lung cancer patients from GSE69405 (M) and E-MTAB-6653 (N–P) databases.
Evaluation of KRT15 expression in lung cancer tissue and cells
Utilizing the CCLE database, we detected a significant rise in KRT15 mRNA expression levels in LUAD cell lines (Figure 6A). Then, our analysis of the HPA database revealed elevated KRT15 protein expression levels in LUAD samples, with predominant distribution observed in the cytoplasmic and membranous regions (Figure 6B). Next, through the RT-qPCR assay, we ascertained elevated KRT15 mRNA expression levels in LUAD as compared to normal lung samples, supporting our findings (Figure 6C and D). Additionally, we also confirmed that LUAD samples exhibited dramatically elevated KRT15 expression at the translational level in comparison to the samples from health individuals (Figure 6E and F).

Evaluation of KRT15 expression in lung cancer tissue and cells (A) representative forest plot showing the transcription of KRT15 in different LUAD cell lines obtained from CCLE analysis. (B) The protein of KRT15 at the translation level was analyzed by IHC method using The human protein atlas database. (C, D) Representative bar graph (C) and estimation plot (D) depicting the transcription of KRT15 between tumor and non-cancer tissues from LUAD patients. The adjacent noncancerous tissues were used as the normal controls. (E, F) Representative blot (E) and bar graph (F) representing the translation of KRT15 between tumor and non-cancer tissues from LUAD patients. The adjacent noncancerous tissues were used as the normal controls. ** and **** represent p-values less than 0.01 or 0.0001 respectively.
Silencing of KRT15 suppresses the growth, metastasis, and survival of lung cancer cells
The results obtained above indicate the potential involvement of KRT15 in cell proliferation, metastasis, and cell survival. To further elucidate the functional significance of KRT15 in lung cancer cells, a series of biological function assays were carried out. As shown in Figure 7A, a targeted siRNA directed against the KRT15 sequence led to a significant decrease in KRT15 expression in both NCI-H1573 and PC-14 cells (Figure 7A). The CCK-8 assay demonstrated significant suppression of the growth in NCI-H1573 (Figure 7B) and PC-14 (Figure 7C) cells upon down-modulation of KRT15 (Figure 7B and C). Additionally, inhibition of KRT15 resulted in a significant decrease in migrative (Figure 7D and E) and invasive (Figure 7F and G) capacities of NCI-H1573 (Figure 7D–F) and PC-14 (Figure 7E–G) cells, as evidenced by the transwell assay (Figure 7D–G). Finally, the apoptotic status of NCI-H1573 and PC-14 cells after KRT15 silencing was monitored and the results demonstrated that, in comparison to the control NCI-H1573 and PC-14 cells (siNC), those KRT15-deficient NCI-H1573 (Figure 7H) and PC-14 (Figure 7I) cells exhibited dramatically increased cell apoptosis, suggesting the promotion of KRT15 in the survival of lung cancer cells.

Silencing of KRT15 suppresses the growth, metastasis, and survival of lung cancer cells (A) representative bar graph showing the transcription of KRT15 in human NCI-H1573 and PC-14 lung cancer cells after indicated treatments. (B, C) The growth of human NCI-H1573 (B) and PC-14 (C) lung cancer cells with or without KRT15 deficiency was monitored by CCK-8 method. (D–G) Representative images (up panels) and bar graphs (low panels) representing the migration (D, E) and invasion (F, G) of human NCI-H1573 (D, F) and PC-14 (E, G) lung cancer cells with or without KRT15 silencing. (H, I) Representative histograms (left panel) and bar graphs (right panel) depicting the apoptotic rates of human NCI-H1573 (H) and PC-14 (I) lung cancer cells with or without KRT15 silencing. *, **, and *** represent the p-values less than 0.05, 0.01, and 0.001 respectively. The cells transfected with the siRNA-control (siNC) plasmids were used as the controls.
Discussion
The convergence of rapid medical technological advancements and the tireless dedication of oncologists in research have led to a more scientific, universally applicable, and personalized approach in the diagnosis, monitoring, and treatment of oncological conditions. As a type I keratin, KRT15 exhibits significant expression in the basal layer and stratified epidermis of human skin [37]. Dysregulated expression of KRT15 is observed in individuals diagnosed with a range of cancer types, such as breast invasive carcinoma [20], colorectal cancer [25], renal cell carcinoma [23], and endometrial cancer [38]. According to a previous investigation, KRT15 exhibits upregulation specifically in a subgroup of invasive ureteric and urinary bladder cancers [27]. Additionally, the presence of lymphovascular invasion and higher International Federation of Gynecology and Obstetrics (FIGO) stage in endometrial cancer patients is significantly associated with the KRT15 protein, while KRT15 mRNA exhibits correlation solely with unfavorable FIGO stage [27]. Furthermore, in patients afflicted with esophageal carcinoma, KRT15 is linked to the advanced T stage and TNM stage, along with the occurrence of lymph node metastasis [39]. However, no research has been conducted to explore the clinical role and regulatory mechanism of KRT15 in lung cancer patients.
This pioneering research presents the first comprehensive exploration of KRT15 in pan-cancer, offering a novel and in-depth perspective. Among the most prevalent malignancies worldwide, NSCLC is prominent, with LUAD contributing to 40 % of all NSCLC cases. Based on our analysis, a significant relationship between KRT15 expression levels and the OS and PFS of LUAD patients was observed. This prompted us to undertake a meticulous investigation to further elucidate the intricate biological function of KRT15 in the context of LUAD. Our study aimed to unravel the crucial involvement of the KRT15 gene in the advancement of various malignancies. In order to fulfill our aim, we conducted an exhaustive search of data encompassing a broad range of cancer types and included a substantial number of participants. In the preliminary stage, we observed the disrupted expression of KRT15 in a range of malignancies. Through the implementation of the Cox regression model, we discerned unique associations between KRT15 expression and varied survival markers, suggesting its potential as a promising prognostic component in specific cancers.
Immunotherapy is currently a highly favored and extensively discussed topic in the field of oncology. The presence of tumor-infiltrating lymphocytes, the assessment of TMB, and the evaluation of MSI status are pivotal factors influencing the response to immunotherapy and dictating the clinical outcomes in cancer [40]. According to our investigation, KRT15 consistently elicited activation in all six diverse immune cell types (neutrophils, CD8+ T cells, macrophages, dendritic cells, CD4+ T cells, and B cells) across the majority of cancer types. Across approximately 20 types of cancers, a significant association was observed between the key immune checkpoints (CD274, CTLA-4, HAVCR2, LAG3, PDCD1, PDCD1LG2, SIGLEC15, and TIGIT) and KRT15 expression, excluding UVM, READ, LAML, GBM, and BLCA. The identification of particular gene mutations may enable the anticipation of patients’ prognoses and treatment outcomes. Increased somatic TMB and MSI have demonstrated a positive correlation with enhanced efficacy of immunotherapy and favorable OS across various cancer histology. Mutations in mismatch repair genes and the compromised function of these genes were identified as the underlying factors contributing to MSI. Herein, we identified a positive correlation between KRT15 expression level and TMB in PAAD, STAD, and READ, whereas a negative correlation in PRAD, COAD, ESCA, and BRCA. These comprehensive findings suggest the potential role of KRT15 in promoting tumor growth by upregulating MSI and TMB via the modulation of genes associated with mismatch repair. Tumors harboring high TMB and MSI-H status are characterized by prominent immune cell infiltration, making them more amenable to potential immune checkpoint inhibitor responses [41], [42], [43]. In view of this understanding, these high KRT15-expressed cancer types could offer potential advantages for immune checkpoint inhibitor treatments. The use of single-cell transcriptome sequencing and KEGG pathway analysis highlighted a noteworthy connection between KRT15 expression and cell proliferation in NSCLC. Functional assays further affirmed the inhibitory impact of KRT15-deficiency on lung cancer cell growth and its influence on cell growth.
According to a previous study, in colorectal cancer cells, KRT15 exerts regulatory control over the β-catenin/matrix metalloproteinase-7 (MMP-7) pathway, thereby facilitating the migration and invasion [44]. Additionally, a connection of KRT15 expression with cell apoptosis was also observed [38]. In line with these publications, our results also demonstrated that KRT15 expression is closely connected to metastasis and survival, silencing of KRT15 could dramatically suppress the migration, invasion, as well as survival of lung cancer cells.
Additionally, in our study, we employed three major datasets, The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and the Human Protein Atlas (HPA), to comprehensively investigate the clinical relevance of the KRT15 gene across a diverse array of cancer, with a focused analysis on LUAD. Each dataset was chosen for its unique strengths and contributions to cancer research. TCGA provided a robust foundation with its comprehensive genomic and transcriptomic profiles, enabling us to perform a detailed pan-cancer analysis and identify KRT15’s roles across various cancer types. GEO supplemented our findings by offering additional genomic datasets, which allowed us to validate the functional implications of KRT15 in cancer-specific signaling pathways and tumor behaviors. The HPA was instrumental in corroborating protein expression profiles, linking genetic alterations observed in TCGA and GEO with protein-level evidence, thereby strengthening our hypothesis about KRT15’s involvement in tumor progression.
To ensure the consistency and compatibility of our results across these datasets, we applied rigorous data preprocessing and normalization techniques. This approach mitigated technical variances and aligned the data for comparative analysis. Notably, the consistency of our findings across these distinct datasets was reinforced through cross-validation techniques, where key results related to KRT15’s role as a biomarker and potential therapeutic target were robustly supported. These efforts ensured that our study’s conclusions were based on reliable and reproducible evidence, thus enhancing the overall credibility of our findings regarding the significance of KRT15 in cancer prognosis and therapy.
In this study, while we demonstrated the potential role of KRT15 as a biomarker and therapeutic target across a diverse spectrum of cancers, we acknowledge several limitations. The heterogeneity of cancer types included in our pan-cancer analysis, although valuable for broad biomarker discovery, introduces variability that may affect the generalizability of our findings to individual cancer types. Each cancer type possesses unique genetic, epigenetic, and microenvironmental characteristics that could influence the expression and functional impacts of KRT15. Consequently, the translational relevance of our results may be constrained by this diversity. Furthermore, our reliance on in vitro models for validating the functions of KRT15, while essential for initial mechanistic explorations, may not fully capture the complex interactions and behaviors present within the in vivo tumor milieu. Such models lack the comprehensive representation of immune system interactions and the tumor microenvironment, which are crucial for assessing the therapeutic potential of targets like KRT15. Future investigations employing in vivo models and patient-derived samples will be critical to overcome these limitations and enhance the clinical applicability of our findings.
In conclusion, our study is the inaugural endeavor to undertake a thorough analysis of KRT15 in pan-cancer. Our research unveiled an elevation in KRT15 expression across various malignancies, demonstrating an adverse correlation with OS. Our validation experiment demonstrated a considerable increase in KRT15 expression in LUAD samples when compared to the corresponding non-tumor samples. The implications of this finding highlight KRT15 as a potential biological marker for LUAD. Furthermore, its involvement in LUAD advancement is closely linked to its role in promoting cell growth, metastasis, and survival of LUAD cells.
Funding source: Yichang Natural Science Foundation
-
Research ethics: Not applicable.
-
Informed consent: Not applicable.
-
Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.
-
Use of Large Language Models, AI and Machine Learning Tools: None declared.
-
Conflict of interests: Authors state no conflict of interest.
-
Research funding: Mechanism and clinical significance of lipid PDIM regulation of Glactin-3 (Gal-3) - mediated inflammatory response and granulomatous formation in the cell wall of Mycobacterium tuberculosis. Jin Zhu. Yichang Natural Science Foundation A24-2-063.
-
Data availability: Not applicable.
References
1. Thandra, KC, Barsouk, A, Saginala, K, Aluru, JS, Barsouk, A. Epidemiology of lung cancer. Contemp Oncol 2021;25:45–52. https://doi.org/10.5114/wo.2021.103829.Search in Google Scholar PubMed PubMed Central
2. Bragulla, HH, Homberger, DG. Structure and functions of keratin proteins in simple, stratified, keratinized and cornified epithelia. J Anat 2009;214:516–59. https://doi.org/10.1111/j.1469-7580.2009.01066.x.Search in Google Scholar PubMed PubMed Central
3. Midthun, DE. Early detection of lung cancer. F1000Res 2016;5:739. https://doi.org/10.12688/f1000research.7313.1.Search in Google Scholar PubMed PubMed Central
4. Nooreldeen, R, Bach, H. Current and future development in lung cancer diagnosis. Int J Mol Sci 2021;22:8661. https://doi.org/10.3390/ijms22168661.Search in Google Scholar PubMed PubMed Central
5. Rudin, CM, Brambilla, E, Faivre-Finn, C, Sage, J. Small-cell lung cancer. Nat Rev Dis Prim 2021;7:3. https://doi.org/10.1038/s41572-020-00235-0.Search in Google Scholar PubMed PubMed Central
6. Schabath, MB, Cote, ML. Cancer progress and priorities: lung cancer. Cancer Epidemiol Biomarkers Prev 2019;28:1563–79. https://doi.org/10.1158/1055-9965.epi-19-0221.Search in Google Scholar PubMed PubMed Central
7. Miller, KD, Nogueira, L, Mariotto, AB, Rowland, JH, Yabroff, KR, Alfano, CM, et al.. Cancer treatment and survivorship statistics, 2019. CA Cancer J Clin 2019;69:363–85. https://doi.org/10.3322/caac.21565.Search in Google Scholar PubMed
8. Siegel, RL, Miller, KD, Fuchs, HE, Jemal, A. Cancer statistics, 2021. CA Cancer J Clin 2021;71:7–33. https://doi.org/10.3322/caac.21654.Search in Google Scholar PubMed
9. Li, Y, Zhang, M, Ying, W. Keratin 1 plays significant roles in maintaining the survival and oxidative stress state of B16-F10 melanoma cell lines. Int J Physiol Pathophysiol Pharmacol 2021;13:12–6.Search in Google Scholar
10. Chen, JW, Dhahbi, J. Lung adenocarcinoma and lung squamous cell carcinoma cancer classification, biomarker identification, and gene expression analysis using overlapping feature selection methods. Sci Rep 2021;11:13323. https://doi.org/10.1038/s41598-021-92725-8.Search in Google Scholar PubMed PubMed Central
11. Nicholson, AG, Tsao, MS, Beasley, MB, Borczuk, AC, Brambilla, E, Cooper, WA, et al.. The 2021 WHO classification of lung tumors: impact of advances since 2015. J Thorac Oncol 2022;17:362–87. https://doi.org/10.1016/j.jtho.2021.11.003.Search in Google Scholar PubMed
12. Wang, W, Liu, H, Li, G. What’s the difference between lung adenocarcinoma and lung squamous cell carcinoma? Evidence from a retrospective analysis in a cohort of Chinese patients. Front Endocrinol 2022;13:947443. https://doi.org/10.3389/fendo.2022.947443.Search in Google Scholar PubMed PubMed Central
13. Myers, DJ, Wallen, JM. Lung adenocarcinoma. Treasure Island (FL): StatPearls Publishing; 2023.Search in Google Scholar
14. Takan, I, Karakülah, G, Louka, A, Pavlopoulou, A. “In the light of evolution:” keratins as exceptional tumor biomarkers. PeerJ 2023;11:e15099. https://doi.org/10.7717/peerj.15099.Search in Google Scholar PubMed PubMed Central
15. Moll, R, Divo, M, Langbein, L. The human keratins: biology and pathology. Histochem Cell Biol 2008;129:705–33. https://doi.org/10.1007/s00418-008-0435-6.Search in Google Scholar PubMed PubMed Central
16. Ku, NO, Strnad, P, Bantel, H, Omary, MB. Keratins: biomarkers and modulators of apoptotic and necrotic cell death in the liver. Hepatology 2016;64:966–76. https://doi.org/10.1002/hep.28493.Search in Google Scholar PubMed PubMed Central
17. Dmello, C, Srivastava, SS, Tiwari, R, Chaudhari, PR, Sawant, S, Vaidya, MM. Multifaceted role of keratins in epithelial cell differentiation and transformation. J Biosci 2019;44:33. https://doi.org/10.1007/s12038-019-9864-8.Search in Google Scholar
18. Yoon, S, Leube, RE. Keratin intermediate filaments: intermediaries of epithelial cell migration. Essays Biochem 2019;63:521–33. https://doi.org/10.1042/ebc20190017.Search in Google Scholar
19. Elazezy, M, Schwentesius, S, Stegat, L, Wikman, H, Werner, S, Mansour, WY, et al.. Emerging insights into keratin 16 expression during metastatic progression of breast cancer. Cancers 2021;13:3869. https://doi.org/10.3390/cancers13153869.Search in Google Scholar PubMed PubMed Central
20. Zhong, P, Shu, R, Wu, H, Liu, Z, Shen, X, Hu, Y. Low KRT15 expression is associated with poor prognosis in patients with breast invasive carcinoma. Exp Ther Med 2021;21:305. https://doi.org/10.3892/etm.2021.9736.Search in Google Scholar PubMed PubMed Central
21. Liu, Y, Lyle, S, Yang, Z, Cotsarelis, G. Keratin 15 promoter targets putative epithelial stem cells in the hair follicle bulge. J Invest Dermatol 2003;121:963–8. https://doi.org/10.1046/j.1523-1747.2003.12600.x.Search in Google Scholar PubMed
22. Nagosa, S, Leesch, F, Putin, D, Bhattacharya, S, Altshuler, A, Serror, L, et al.. microRNA-184 induces a commitment switch to epidermal differentiation. Stem Cell Rep 2017;9:1991–2004. https://doi.org/10.1016/j.stemcr.2017.10.030.Search in Google Scholar PubMed PubMed Central
23. Zhang, W, Chen, P, Li, Z, Zhang, R, Zhang, J. Clinical implication of keratin-15 quantification for renal cell carcinoma management: its dysregulation and association with clinicopathologic characteristics and prognostication. Tohoku J Exp Med 2023;260:99–107. https://doi.org/10.1620/tjem.2023.j017.Search in Google Scholar
24. Sanchez-Palencia, A, Gomez-Morales, M, Gomez-Capilla, JA, Pedraza, V, Boyero, L, Rosell, R, et al.. Gene expression profiling reveals novel biomarkers in nonsmall cell lung cancer. Int J Cancer 2011;129:355–64. https://doi.org/10.1002/ijc.25704.Search in Google Scholar PubMed
25. Rao, X, Wang, J, Song, HM, Deng, B, Li, JG. KRT15 overexpression predicts poor prognosis in colorectal cancer. Neoplasma 2020;67:410–4. https://doi.org/10.4149/neo_2019_190531n475.Search in Google Scholar PubMed
26. Folgueira, MA, Brentani, H, Katayama, ML, Patrão, DF, Carraro, DM, Mourão, NM, et al.. Gene expression profiling of clinical stages II and III breast cancer. Braz J Med Biol Res 2006;39:1101–13. https://doi.org/10.1590/s0100-879x2006000800013.Search in Google Scholar PubMed
27. Tai, G, Ranjzad, P, Marriage, F, Rehman, S, Denley, H, Dixon, J, et al.. Cytokeratin 15 marks basal epithelia in developing ureters and is upregulated in a subset of urothelial cell carcinomas. PLoS One 2013;8:e81167. https://doi.org/10.1371/journal.pone.0081167.Search in Google Scholar PubMed PubMed Central
28. Pal, SK, Sakamoto, K, Aragaki, T, Akashi, T, Yamaguchi, A. The expression profiles of acidic epithelial keratins in ameloblastoma. Oral Surg Oral Med Oral Pathol Oral Radiol 2013;115:523–31. https://doi.org/10.1016/j.oooo.2013.01.017.Search in Google Scholar PubMed
29. Zekri, AN, El-Sisi, ER, Abdallah, ZF, Ismail, A, Barakat, A. Gene expression profiling of circulating CD133+ cells of hepatocellular carcinoma patients associated with HCV infection. J Egypt Natl Canc Inst 2017;29:19–24. https://doi.org/10.1016/j.jnci.2016.12.002.Search in Google Scholar PubMed
30. Khanom, R, Sakamoto, K, Pal, SK, Shimada, Y, Morita, K, Omura, K, et al.. Expression of basal cell keratin 15 and keratin 19 in oral squamous neoplasms represents diverse pathophysiologies. Histol Histopathol 2012;27:949–59. https://doi.org/10.14670/HH-27.949.Search in Google Scholar PubMed
31. Shan, M, Xia, Q, Yan, D, Zhu, Y, Zhang, X, Zhang, G, et al.. Molecular analyses of prostate tumors for diagnosis of malignancy on fine-needle aspiration biopsies. Oncotarget 2017;8:104761–71. https://doi.org/10.18632/oncotarget.22289.Search in Google Scholar PubMed PubMed Central
32. Zhang, C, Liang, Y, Ma, MH, Wu, KZ, Dai, DQ, et al.. KRT15, INHBA, MATN3, and AGT are aberrantly methylated and differentially expressed in gastric cancer and associated with prognosis. Pathol Res Pract 2019;215:893–9. https://doi.org/10.1016/j.prp.2019.01.034.Search in Google Scholar PubMed
33. Chong, LY, Cheok, PY, Tan, WJ, Thike, AA, Allen, G, Ang, MK, et al.. Keratin 15, transcobalamin I and homeobox gene Hox-B13 expression in breast phyllodes tumors: novel markers in biological classification. Breast Cancer Res Treat 2012;132:143–51. https://doi.org/10.1007/s10549-011-1555-6.Search in Google Scholar PubMed
34. Giroux, V, Lento, AA, Islam, M, Pitarresi, JR, Kharbanda, A, Hamilton, KE, et al.. Long-lived keratin15+ esophageal progenitor cells contribute to homeostasis and regeneration. J Clin Invest 2017;127:2378–91. https://doi.org/10.1172/jci88941.Search in Google Scholar
35. Bedognetti, D, Ceccarelli, M, Galluzzi, L, Lu, R, Palucka, K, Samayoa, J, et al.. Toward a comprehensive view of cancer immune responsiveness: a synopsis from the SITC workshop. J Immunother Cancer 2019;7:131. https://doi.org/10.1186/s40425-019-0602-4.Search in Google Scholar PubMed PubMed Central
36. Bonneville, R, Krook, MA, Kautto, EA, Miya, J, Wing, MR, Chen, HZ, et al.. Landscape of microsatellite instability across 39 cancer types. JCO Precis Oncol 2017;1:1. https://doi.org/10.1200/po.17.00073.Search in Google Scholar
37. Waseem, A, Dogan, B, Tidman, N, Alam, Y, Purkis, P, Jackson, S, et al.. Keratin 15 expression in stratified epithelia: downregulation in activated keratinocytes. J Invest Dermatol 1999;112:362–9. https://doi.org/10.1046/j.1523-1747.1999.00535.x.Search in Google Scholar PubMed
38. Yang, H, Li, A, Li, A, Zhao, F, Zhang, T. Upregulated keratin 15 links to the occurrence of lymphovascular invasion, stromal cervical invasion as well as unfavorable survival profile in endometrial cancer patients. Med (Baltim) 2022;101:e29686. https://doi.org/10.1097/md.0000000000029686.Search in Google Scholar PubMed PubMed Central
39. Lin, JB, Feng, Z, Qiu, ML, Luo, RG, Li, X, Liu, B. KRT 15 as a prognostic biomarker is highly expressed in esophageal carcinoma. Future Oncol 2020;16:1903–9. https://doi.org/10.2217/fon-2019-0603.Search in Google Scholar PubMed
40. Keshinro, A, Vanderbilt, C, Kim, JK, Firat, C, Chen, CT, Yaeger, R, et al.. Tumor-infiltrating lymphocytes, tumor mutational burden, and genetic alterations in microsatellite unstable, microsatellite stable, or mutant POLE/POLD1 colon cancer. JCO Precis Oncol 2021;5:817. https://doi.org/10.1200/PO.20.00456.Search in Google Scholar PubMed PubMed Central
41. Lin, A, Xu, W, Luo, P, Zhang, J. Mutations status of chemokine signaling pathway predict prognosis of immune checkpoint inhibitors in colon adenocarcinoma. Front Pharmacol 2021;12:721181. https://doi.org/10.3389/fphar.2021.721181.Search in Google Scholar PubMed PubMed Central
42. Picard, E, Verschoor, CP, Ma, GW, Pawelec, G. Relationships between immune landscapes, genetic subtypes and responses to immunotherapy in colorectal cancer. Front Immunol 2020;11:369. https://doi.org/10.3389/fimmu.2020.00369.Search in Google Scholar PubMed PubMed Central
43. Schrock, AB, Ouyang, C, Sandhu, J, Sokol, E, Jin, D, Ross, JS, et al.. Tumor mutational burden is predictive of response to immune checkpoint inhibitors in MSI-high metastatic colorectal cancer. Ann Oncol 2019;30:1096–103. https://doi.org/10.1093/annonc/mdz134.Search in Google Scholar PubMed
44. Chen, W, Miao, C. KRT15 promotes colorectal cancer cell migration and invasion through β-catenin/MMP-7 signaling pathway. Med Oncol 2022;39:68. https://doi.org/10.1007/s12032-021-01619-2.Search in Google Scholar PubMed
© 2025 the author(s), published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.
Articles in the same Issue
- Frontmatter
- Review
- Interaction of fetuin-A with obesity related insulin resistance and diabetes mellitus
- Research Articles
- Preanalytical errors in pediatric blood sampling: a systematic review of common challenges and risks
- The promotive role of reticulocalbin 3 (RCN3) in the pathogenesis of keloid via TGFβ1/Smad2/Smad7 signaling pathway in vitro
- Investigation into drug resistance to cisplatin in cancer stem cell-enriched population in non-small cell lung cancer
- The cytotoxic and antiproliferative effect of Polygala saponin XLIV on the human colorectal carcinoma cell line
- Development of HEK293T cell reference materials for β-thalassemia genetic testing using prime editing
- Investigation of propofol, fentanyl, and midazolam-related toxicity and the protective effect of midazolam on THLE-2 cell lines
- Leucine-rich α-2-glycoprotein 1 can be a novel angiogenic mediator in autosomal dominant polycystic kidney disease
- Asiaticoside reverses the inhibition effect of miR-184 on proliferation, migration and AKT phosphorylation of HTR-8/Svneo cells
- Perioperative D-dimer levels and head and neck cancer surgery: a prospective observational study
- Comprehensive analysis of KRT15 in pan-cancer and verification in lung adenocarcinoma
- Age and sex-dependent effects of nitrosative stress on the osmotic fragility of human red blood cells
- Evaluation of oxidative stress parameters in older patients with urinary incontinence
- Determination of reference change values for thyroid-related biomarkers: TSH, fT3, fT4, Tg, Anti-Tg, and Anti-TPO
- An evaluation of serum boron level in pregnancies with severe pre-eclampsia
- Involvement of SIRT3/FOXO1 and TLR4/NF-κB/NLRP3 pathways in protective effects of Δ(9)-tetrahydrocannabinol on diabetic cardiomyopathy in rats
- Evaluating the potential therapeutic effect of Rosa damascena Mill. essential oil on acetic acid-induced ulcerative colitis in rats
- Assessment of blood collection competency in phlebotomists: a survey study
Articles in the same Issue
- Frontmatter
- Review
- Interaction of fetuin-A with obesity related insulin resistance and diabetes mellitus
- Research Articles
- Preanalytical errors in pediatric blood sampling: a systematic review of common challenges and risks
- The promotive role of reticulocalbin 3 (RCN3) in the pathogenesis of keloid via TGFβ1/Smad2/Smad7 signaling pathway in vitro
- Investigation into drug resistance to cisplatin in cancer stem cell-enriched population in non-small cell lung cancer
- The cytotoxic and antiproliferative effect of Polygala saponin XLIV on the human colorectal carcinoma cell line
- Development of HEK293T cell reference materials for β-thalassemia genetic testing using prime editing
- Investigation of propofol, fentanyl, and midazolam-related toxicity and the protective effect of midazolam on THLE-2 cell lines
- Leucine-rich α-2-glycoprotein 1 can be a novel angiogenic mediator in autosomal dominant polycystic kidney disease
- Asiaticoside reverses the inhibition effect of miR-184 on proliferation, migration and AKT phosphorylation of HTR-8/Svneo cells
- Perioperative D-dimer levels and head and neck cancer surgery: a prospective observational study
- Comprehensive analysis of KRT15 in pan-cancer and verification in lung adenocarcinoma
- Age and sex-dependent effects of nitrosative stress on the osmotic fragility of human red blood cells
- Evaluation of oxidative stress parameters in older patients with urinary incontinence
- Determination of reference change values for thyroid-related biomarkers: TSH, fT3, fT4, Tg, Anti-Tg, and Anti-TPO
- An evaluation of serum boron level in pregnancies with severe pre-eclampsia
- Involvement of SIRT3/FOXO1 and TLR4/NF-κB/NLRP3 pathways in protective effects of Δ(9)-tetrahydrocannabinol on diabetic cardiomyopathy in rats
- Evaluating the potential therapeutic effect of Rosa damascena Mill. essential oil on acetic acid-induced ulcerative colitis in rats
- Assessment of blood collection competency in phlebotomists: a survey study