Abstract
Background/Aim: Autophagy plays a crucial role in maintaining cellular homeostasis and has been implicated in the pathogenesis of knee osteoarthritis (OA). However, data on radiographic stage-dependent transcriptional variation of autophagy-related genes in patients with knee OA, particularly using peripheral blood samples, remain limited. The aim of this study was to evaluate whether disease severity was associated with stage-dependent changes in the expression of selected autophagy-related genes within a patient cohort.
Patients and Methods: A total of 200 patients diagnosed with knee OA were included in the study. Disease severity was classified according to the Kellgren-Lawrence radiographic grading system. Peripheral blood samples were collected, and the expression levels of selected autophagy-related genes were analyzed using quantitative real-time polymerase chain reaction [autophagy-related 5 (ATG5), ATG7, unc-51-like kinase 1 (ULK1), microtubule-associated protein 1 light chain 3 beta (LC3B), WD repeat domain phosphoinositide-interacting protein 1 (WIPI1), neighbor of BRCA1 gene 1 (NBR1), forkhead box O3 (FOXO3), transcription factor EB (TFEB)]. Relative gene expression was calculated using the ΔCt method, and comparisons were performed across radiographic stages. Associations between gene expression levels and systemic inflammatory markers were also assessed.
Results: Significant stage-dependent differences were observed in the expression of ULK1, TFEB, WIPI1, and NBR1 (p<0.05), with higher ΔCt values (reduced relative expression) in advanced radiographic stages compared with early-stage disease. In contrast, ATG5, ATG7, LC3B, and FOXO3 expression remained stable across radiographic stages. Furthermore, no significant associations were observed between expression of autophagy-related genes and systemic inflammatory status, as assessed by C-reactive protein levels.
Conclusion: In patients with knee OA, regulatory and early autophagy-related genes exhibit radiographic stage-associated transcriptional alterations in peripheral blood, while expression of core autophagy machinery genes remain relatively stable. These findings suggest that disease severity is associated with selective transcriptional changes in autophagy-related pathways within the OA patient population and support further investigation of stage-dependent molecular patterns in knee OA.
Introduction
Knee osteoarthritis (OA) is a long-term and progressive condition of the joint where the cartilage gradually deteriorates. At the same time, the subchondral bone remodels, there is inflammation of the synovial lining, and the patient experiences pain, eventually leading to loss of function and a decrease in their quality of life (1). The number of people suffering from knee OA increases with age and it is one of the main reasons for disability globally (2). Even though knee OA is clinically very problematic, the molecular pathways that cause the disease and its progression are still not fully understood, and the current treatment options mainly focus on relieving the symptoms rather than changing the course of the (3).
Autophagy is a highly conserved intracellular degradation pathway that facilitates the clearance and reuse of damaged organelles and macromolecules via lysosomal processing. The latter is vital for maintaining cellular homeostasis and is especially crucial under stress scenarios such as nutrient scarcity, oxidative stress, and inflammation (4). In the case of joint tissues, autophagy is one of the mechanisms that the body uses to maintain cartilage and cell survival. Dysregulation of autophagy has been associated with the initiation and the worsening of OA (5, 6). Both experimental and clinical research has demonstrated that autophagic activity is likely to be maintained at normal levels or even upregulated in the early stages of OA. On the contrary, a gradual reduction of autophagy in advanced disease has been reported, which might lead to cell dysfunction and tissue (7).
Autophagy is a multistep process that involves different molecular components that regulate initiation, vesicle nucleation, elongation, cargo recognition, and lysosomal fusion (8). The unc-51-like kinase 1 (ULK1) complex serves as a main initiator of autophagy upon cellular energy and nutrient signaling, whereas WD repeat domain phosphoinositide-interacting protein 1 (WIPI1) and autophagy-related proteins take part in autophagosome formation and membrane dynamics (9, 10). Transcriptional regulators such as transcription factor EB (TFEB) coordinate lysosomal biogenesis and expression of autophagy-related genes, thus linking the environmental stress response to the long-term regulation of autophagic capacity (11). Furthermore, selective autophagy receptors, such as neighbor of BRCA DNA repair-associated 1 gene 1 (NBR1), mediate the specific degradation of certain cellular components and thus indicate that different steps of autophagy may be differentially affected during disease (12).
Previous studies have explored how autophagy affects OA by measuring autophagy-related markers in cartilage tissue, synovium, or experimental models (5, 6, 13). These studies have indeed shed light on the role of autophagy in OA pathophysiology; however, the results reported differ based on disease stage, tissue type, and experimental setting (5). Although changes in autophagy-related gene expression have been reported, a comprehensive understanding of stage-dependent transcriptional regulation in human knee OA remains limited, particularly in studies using minimally invasive and clinically applicable sample sources (14).
Many studies have compared autophagy-related markers between patients with OA and healthy individuals, but such case-control studies do not really investigate whether molecular changes are directly related to the severity of the disease shown by an X-ray. Thus, without a healthy control group, the current study aimed to assess the changes in expression of some autophagy-related genes in an OA depending on the stage of the disease according to Kellgren-Lawrence (K-L) grade (15).
Hence, to investigate the stage-dependent expression patterns of certain genes related to autophagy in patients with knee OA, the study utilized peripheral blood samples. The study targeted genes participating in autophagy initiation and regulation, and in selective autophagy, as well as the essential structural components of the autophagy machinery, aiming to detect changes in gene transcription linked to the severity of the disease. A better knowledge of these molecular patterns can help explain the role of autophagy regulation in the progression of knee OA and may also pave the way for the development of molecular indicators that would complement clinical and radiological.
Patients and Methods
Study population. A total of 200 patients diagnosed with knee OA were included in the study. Patients were recruited from the outpatient clinics of Pamukkale University Hospital based on clinical and radiographic diagnostic criteria. Inclusion criteria comprised age ≥40 years and radiographic evidence of knee OA according to the K-L grading system. Exclusion criteria included inflammatory rheumatic diseases, previous knee surgery, malignancy, active infection, chronic autoimmune disorders, and current use of systemic corticosteroids or immunosuppressive therapy. The present study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of Pamukkale University Faculty of Medicine (approval number: E-60116787-020-805062). Written informed consent was obtained from all participants prior to enrollment.
The diagnosis of OA was established based on clinical evaluation and radiographic findings. Disease severity was classified according to the K-L grading system (stages 1-4) using standard anteroposterior knee radiographs and this was used to stratify patients for stage-dependent analyses. Grade 0 represents a normal joint structure; grade 1 indicates doubtful joint space narrowing and possible osteophytic lipping; grade 2 shows definite osteophyte formation with possible joint space narrowing; grade 3 is characterized by moderate multiple osteophytes, definite joint space narrowing, and subchondral sclerosis; grade 4 demonstrates severe joint space narrowing with large osteophytes, marked sclerosis, and bony deformity (Figure 1).
Kellgren-Lawrence radiographic grading of knee osteoarthritis.
Clinical and laboratory assessment. Demographic characteristics (age and sex) and routine laboratory parameters were recorded for all participants. Peripheral blood samples were collected to assess systemic inflammatory status and gene expression profiles. Complete blood counts were measured using an automated hematology analyzer (Sysmex XN-1000; Sysmex Corporation, Kobe, Japan). C-Reactive protein (CRP) levels were determined by an immunoturbidimetric method using an Abbott Architect c16000 analyzer (Abbott Diagnostics, Abbott Park, IL, USA). Erythrocyte sedimentation rate was measured with an automated system (Alifax TEST1; Alifax S.p.A., Padua, Italy). All laboratory analyses were conducted at the Central Laboratory of Pamukkale University Hospital in accordance with the manufacturer’s instructions and international quality control standards. The laboratory data obtained were used in statistical analyses along with the stage of knee OA and gene expression levels related to autophagy. The overall workflow of the study, including patient selection, clinical evaluation, blood sample collection, and gene expression analysis, is summarized in Figure 2.
Overview of patient evaluation, molecular analysis, and key outcomes in patients with knee osteoarthritis (OA). ATG5: Autophagy-related 5; ATG7: autophagy-related 7; FOXO3: Forkhead box O3; LC3B: microtubule-associated protein 1 light chain 3 beta; NBR1: neighbor of BRCA DNA repair-associated 1 gene 1; TFEB: transcription factor EB; ULK1: unc-51-like kinase 1; WIPI1: WD repeat domain phosphoinositide-interacting protein 1.
Blood sample collection and RNA isolation. Peripheral venous blood samples were obtained between 08:00 and 10:00 a.m. under fasting conditions when possible. Samples intended for gene expression analysis were collected into EDTA-containing tubes and transported to the laboratory within 2 h of collection. Total RNA was isolated from peripheral blood using the QIAamp RNA Blood Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. RNA concentration and purity were assessed spectrophotometrically (A260/A280 ratio), and only samples meeting quality criteria were included. RNA samples were stored at −80°C until complementary DNA (cDNA) synthesis. cDNA synthesis was performed using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA) following the manufacturer’s protocol. Equal amounts of total RNA were used for each reverse transcription reaction to ensure consistency across samples.
Quantitative real-time polymerase chain reaction (qPCR). Gene expression analyses were conducted using qPCR with PowerUp™ SYBR™ Green Master Mix (Applied Biosystems) on an Applied Biosystems 7500 Fast Real-Time PCR System. The autophagy-related genes analyzed were autophagy-related 5 (ATG5), ATG7, ULK1, WIPI1, microtubule-associated protein 1 light chain 3 beta (MAP1LC3B ; hereafter referred to as LC3B) NBR1, forkhead box O3 (FOXO3), and TFEB. Primer sequences are provided in Table I. Primer specificity was confirmed by melting-curve analysis. PCR cycling conditions were as follows: initial denaturation at 95 °C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 s and annealing/extension at 60°C for 60 s. Gene expression levels were normalized using glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and β-actin (ACTB) as reference genes. Relative expression was calculated using the ΔCt method, as no healthy control group was available to serve as a calibrator for ΔΔCt analysis.
Primer sequences for qPCR analysis.
Statistical analysis. Statistical analyses were performed using IBM SPSS Statistics (version 26.0; IBM Corp., Armonk, NY, USA). Data distribution was assessed using histograms, Q-Q plots, and the Shapiro-Wilk test. Continuous variables were expressed as the mean±standard deviation, or median (interquartile range), depending on their distribution. Comparisons of gene expression levels across OA stages were performed using one-way analysis of variance or the Kruskal-Wallis test, as appropriate. Post hoc comparisons were conducted with Tukey HSD or Mann-Whitney U-tests, with Holm-Bonferroni correction applied for multiple comparisons.
Correlations between gene expression and laboratory parameters were evaluated using Pearson or Spearman correlation analyses. Receiver operating characteristic (ROC) curve analysis was performed to assess the ability of selected genes to discriminate early-stage (K-L1-2) from advanced-stage (K-L 3-4) OA. A value of p<0.05 was considered statistically significant.
Results
A total of 200 patients diagnosed with knee OA were included in the study. After classifying patients into four groups according to the K-L radiographic grading system, the distribution of patients across disease stages was uneven, with fewer individuals in stage 1. Demographic characteristics and routine laboratory parameters, including inflammatory markers, are summarized in Table II.
Demographic and clinical characteristics of the study population.
Stage-dependent expression of autophagy-related genes. The expression levels of selected autophagy-related genes were analyzed in peripheral blood samples using quantitative real-time PCR, and ΔCt values were compared across K-L stages 1-4. The distribution of ΔCt values for all analyzed autophagy-related genes according to K-L stage is presented in Figure 3. Eight genes were evaluated: ATG5, ATG7, ULK1, LC3B, TFEB, FOXO3, WIPI1, and NBR1. Statistically significant stage-dependent differences were observed for ULK1, TFEB, WIPI1, and NBR1 expression (Kruskal-Wallis test, p<0.05). Post hoc analyses demonstrated that these differences were primarily driven by comparisons between early-stage OA (K-L stages 1-2) and advanced-stage OA (K-L stage 4). Specifically, ULK1 and TFEB ΔCt values were significantly higher in advanced stages compared with early stages, indicating lower relative mRNA expression in patients with more severe disease. Similarly, WIPI1 and NBR1 expression levels exhibited significant stage-dependent alterations, with reduced expression observed predominantly in advanced-stage OA.
Distribution of ΔCt values for the expression of autophagy-related genes according to the severity of knee osteoarthritis based on Kellgren-Lawrence (K-L) grade. Box plots show expression levels of (A) autophagy-related 5 (ATG5), (B) autophagy-related 7 (ATG7), (C) unc-51-like kinase 1 (ULK1), (D) microtubule-associated protein 1 light chain 3 beta (LC3B), (E) transcription factor EB (TFEB), (F) forkhead box O3 (FOXO3), (G) WD repeat domain phosphoinositide-interacting protein 1 (WIPI1), and (H) neighbor of BRCA DNA repair-associated 1 gene 1 (NBR1). The boxes represent the median and interquartile range; whiskers indicate minimum and maximum values; circles denote outliers. Comparisons between K-L stages were performed using the Kruskal-Wallis test, followed by post hoc pairwise comparisons with Holm-Bonferroni correction where appropriate. Statistically significant stage-dependent differences were observed for ULK1, TFEB, WIPI1, and NBR1. Significantly different at *p<0.05, **p<0.01.
In contrast, no statistically significant differences were detected in the expression levels of ATG5, ATG7, LC3B, or FOXO3 across OA stages (p>0.05), suggesting relative stability of these genes throughout disease progression.
Correlation between autophagy-related gene expression and inflammatory markers. The relationship between autophagy-related gene expression and systemic inflammation was assessed using Spearman correlation analysis with CRP levels. No statistically significant correlations were identified between CRP levels and the expression of any of the analyzed autophagy-related genes (all p>0.05). Although weak correlation tendencies were observed for certain gene-CRP pairs, these did not reach statistical significance and were not indicative of a consistent association. Given the lack of significant findings, correlation results are summarized descriptively, and detailed graphical representations were omitted to maintain clarity and focus on the primary stage-dependent analyses.
ROC curve analysis. ROC curve analysis was conducted to assess the ability of selected autophagy-related genes to discriminate early-stage (K-L grades 1-2) from advanced-stage (K-L grades 3-4) knee OA. Individually, ULK1, TFEB, WIPI1, and NBR1 showed modest discriminatory capacity. A multigene model incorporating these four genes demonstrated improved discrimination compared with single-gene analyses. The corresponding ROC curves and areas under the curve are shown in Figure 4.
Receiver operating characteristic (ROC) curves evaluating the ability of autophagy-related genes to discriminate early-stage (Kellgren-Lawrence grades 1-2) from advanced-stage (Kellgren-Lawrence grades 3-4) knee osteoarthritis. ROC curves are shown for autophagy-related 5 (ATG5), ATG7, forkhead box O3 (FOXO3), microtubule-associated protein 1 light chain 3 beta (LC3B), neighbor of BRCA DNA repair-associated 1 gene 1 (NBR1), transcription factor EB (TFEB), unc-51-like kinase 1 (ULK1), WD repeat domain phosphoinositide-interacting protein 1 (WIPI1), along with a combined multigene model. The corresponding areas under the curve (AUC) and 95% confidence intervals are given. The diagonal line represents the line of no discrimination. Singly, ULK1, TFEB, WIPI1, and NBR1 showed modest discriminatory capacity. A model incorporating all four of these genes demonstrated improved discrimination compared with single-gene analyses.
Figure 5 summarizes the integrative interpretation of the results without implying a direct mechanistic pathway.
Integrative interpretation of the results without implying a direct mechanistic pathway.
Discussion
In this work, we carried out a study on transcriptional changes related to autophagy at different stages in patients with knee OA by using blood samples. The main result is that ULK1, TFEB, WIPI1 and NBR1 genes showed significant stage-dependent changes, concomitant with OA progression, while ATG5, ATG7, LC3B, and FOXO3 did not change significantly at the mRNA expression level with OA stage. This finding indicates that the transcriptional changes observed along OA progression might mainly be concerned with regulation and early-stage components of the autophagy pathway rather than affecting all autophagy-related genes evenly.
Autophagy is a complex, tightly controlled cellular pathway with several stages, and loss of control at different points in this pathway can result in various biological consequences (26). ULK1 is the main autophagy-initiating protein that responds to cellular stress signals, whereas TFEB is the predominant transcription factor that regulates the expression of autophagy- and lysosome-related genes (27, 28). WIPI1 functions in the initial steps of autophagosome formation and membrane remodeling, and it is identified as a component of the core early autophagy machinery (29). On the other hand, NBR1 is a selective autophagy adaptor protein that plays a role in cargo recognition and delivery of ubiquitinated proteins to autophagosomes (30). The finding that the expression of these regulatory and early-process genes varies at different stages and that their transcriptional differences are dependent on the stage supports the idea that the molecular changes associated with OA may mainly impact the upstream regulatory control of autophagy (31, 32).
Our results largely align with previous studies investigating the role of autophagy in OA (31, 32). Cartilage-based and experimental investigations have demonstrated that autophagic activity and the expression of regulatory autophagy-related genes, including ULK1-associated signaling pathways, decline with aging and OA progression (16, 31, 32). In contrast, certain core components of the autophagy machinery appear to remain relatively stable at the transcriptional level. Our findings showing no significant stage-dependent changes in ATG5, ATG7, and LC3B mRNA expression are consistent with this observation. Collectively, both previous experimental data and our results suggest that OA progression may preferentially affect upstream regulatory nodes of the autophagy pathway rather than uniformly altering downstream structural components (26, 31, 32).
TFEB, a major transcription factor that integrates cellular stress responses with metabolic and proteostatic adaptations, regulates transcriptions of genes involved in lysosomal biogenesis and autophagy (11). Through the coordinated lysosomal expression and regulation network, TFEB coordinates numerous autophagy and lysosome genes, playing a central role in the long-term maintenance of cellular homeostasis (33). In the context of OA, studies on experimental OA models and human osteoarthritic cartilage tissues showed that the expression of TFEB and/or nuclear translocation was suppressed (11, 28, 31-33). A decrease in TFEB activity has been linked to lysosomal dysfunction, insufficient autophagic flux, and an increased sensitivity of chondrocytes to oxidative stress, which ultimately leads to the degradation of the cartilage matrix and the loss of cells, thereby accelerating joint degeneration (11, 32).
Recent experimental studies have demonstrated that pharmacological or genetic activation of TFEB enhances autophagic flux, improves lysosomal function, and attenuates inflammatory signaling pathways in cellular and animal models (11, 17, 28). In OA models, TFEB activation has been associated with improved chondrocyte survival and reduced cartilage degeneration (31, 32). These findings suggest that TFEB may represent not only a downstream consequence of degenerative processes but also a potential disease-modifying therapeutic target. Within this context, our study detected stage-dependent TFEB transcriptional changes in peripheral blood, which, as OA severity increases, support systemic impairment of autophagy upper regulatory mechanisms (31). This situation is consistent with the idea that OA progression might be associated with a gradual impact on the regulatory nodal points which are responsible for the initiation and maintenance of autophagy (31-33).
However, it is not realistic to consider the transcriptional profiles obtained from peripheral blood to be a direct reflection of the molecular changes specific to joint tissue. Peripheral TFEB expression is less about intra-articular chondrocytes and more about the combined response of circulating immune cells exposed to systemic inflammatory and metabolic signals (34). Thus, while our findings are consistent with tissue-based studies showing a decrease in TFEB activity in late-stage OA (32), these results should not be considered as direct evidence of intra-articular TFEB dysfunction but rather as an indication of systemic regulatory changes associated with disease stage (24). In the future, holistic studies in which peripheral biomarkers are considered together with joint tissue expression, protein level analyses, and functional autophagy assessments will clarify the real role of TFEB in OA.
WIPI1 is a key regulator of the early stages of autophagy and is recruited to phosphatidylinositol 3-phosphate-enriched membranes during phagophore formation, where it contributes to membrane expansion (29, 35). In recent years, gene expression studies have shown that genes involved in early phases of autophagy can exhibit transcriptional changes in degenerative and aging-related diseases whereas the structural markers representing the late stages of autophagosome (e.g., LC3B) may remain more stable (20-22, 26, 31, 36). This suggests that rather than a complete loss of autophagy, regulatory and initiating mechanisms are selectively.
In this context, the stage-dependent transcriptional changes of WIPI1 identified in our study suggest that the genes involved in the early stages of autophagy are differently regulated during OA progression. However, these findings do not provide direct evidence whether autophagy is functionally upregulated or downregulated; instead of demonstrating functional autophagy activation or suppression, these results indicate that disease severity is associated with transcriptional modulation of genes regulating early autophagy.
Hence, our WIPI1 data do not reveal the presence of a selective autophagy type or a specific autophagic activation state, but rather show that early autophagy control points are transcriptionally affected in OA.
On the other hand, NBR1 is a selective autophagy adaptor protein that directly binds ubiquitinated cargo and facilitates its delivery to autophagosomes for degradation (19, 30). Experimental studies have demonstrated that NBR1 functions in cooperation with LC3 to mediate selective autophagic clearance of protein aggregates and other ubiquitinated substrates (19). NBR1 is particularly involved in protein quality control and cellular homeostasis under chronic stress conditions (37). There are only a few studies directly investigating the role of NBR1 in OA. However, the changes in expression of selective autophagy adaptor proteins have been shown to be associated with chronic inflammation, oxidative stress, and degenerative processes in various disease models (37, 38).
In this context, in our work, the changes in NBR1 mRNA levels depending on the disease stage suggest that the molecular pathways associated with selective autophagy could be involved in the adaptation process in the progression of OA (38). Nevertheless, it is not possible to say that these changes indicate a specific selective autophagy subtype (e.g., mitophagy) or a functional autophagy output directly. Instead, our data demonstrate that NBR1 might be regulated as part of the cellular stress response in OA, but the functional significance of this regulation needs to be uncovered by advanced mechanistic and tissue-based.
Remarkably, no substantial correlations were found between the expression of autophagy-related genes and the systemic inflammatory marker CRP. This suggests that the changes in gene transcription identified in this study were not simply the result of systemic inflammation as measured by CRP. On the other hand, since correlation analyses were carried out on all OA stages combined and immune cell subpopulations were not examined, more comprehensive stage-specific and cell-type-specific studies are.
At the mRNA level, the stage-dependent expression patterns of ULK1, TFEB, WIPI1, and NBR1 observed within this patient cohort suggest that regulatory autophagy-related genes may have potential as complementary molecular indicators of disease severity, particularly when interpreted alongside clinical and radiographic assessments. This interpretation is supported by the modest discriminatory performance observed in ROC analysis (Figure 4). Nevertheless, the absence of a healthy control group, the cross-sectional study design, and the lack of protein-level or functional validation preclude classification of these genes as diagnostic biomarkers. Rather, they should be considered exploratory severity-associated molecular markers requiring further validation.
One of the major features of the present work is the within-cohort design. Since there was no healthy control group, the transcriptional differences found should be seen as severity-associated changes in K-L stages rather than dysregulation compared to a normal baseline. Therefore, this method focused on determining if the gradual development of the disease was related to changes in autophagy gene expression in the OA population. By focusing exclusively on intra-cohort comparisons, this approach minimizes inter-group variability and potential confounding effects that may occur when external control groups are included.
Several limitations should be acknowledged. Firstly, the cross-sectional design of this study allows assessment of gene expression at a single time point and does not provide information about longitudinal molecular changes during OA progression. Secondly, only blood mRNA was analyzed without autophagy estimations in tissues or protein level or functional studies. Further work including repeated sampling over time, analysis of joint tissues, and use of functional tests will shed more light on autophagy regulation in knee OA.
In conclusion, this study has revealed that the genes primarily responsible for the initiation, regulation, and formation of the autophagosome early stages are differentially expressed depending on the stage of the disease in knee OA. At the same time, the genes essential in structural autophagy hardly changed at the mRNA level in peripheral blood. Our findings are consistent with previous reports indicating that dysregulation of upstream regulatory components of the autophagy pathway is associated with OA progression and cartilage degeneration (16, 31, 32). Experimental and mechanistic studies have further emphasized that impairment of key regulatory nodes, including ULK1- and TFEB-mediated pathways, may contribute to disrupted autophagic homeostasis in degenerative conditions (11, 28, 33). Together, these data support the concept that disease severity in OA may be linked to alterations at the regulatory level of the autophagy network rather than to uniform disruption of core structural machinery (26, 31, 32). This offers a platform for future integrative and mechanistic.
Acknowledgements
This study was supported by Strategic Research and Innovation Programme for the Development of the Medical University–Plovdiv (SRIPD-MUP)” Contract No BG-RRP-2.004-0007-C01.
Footnotes
Authors’ Contributions
Conceptualization, Z.Ö. and A.K.; methodology, Z.Ö. and A.K.; software, A.K.; validation, A.K., Z.Ö., M.G. and N.Ö.; formal analysis, A.K.; investigation, Z.Ö., M.G., N.Ö., and A.K.; resources, Z.Ö., and A.K.; data curation, A.K.; writing–original draft preparation, A.K.; writing–review and editing, G.L., M.B., and S.P.; visualization, A.K.; supervision, A.K.; project administration, A.K.; funding acquisition, A.K. All Authors read and agreed to the published version of the manuscript.
Conflicts of Interest
The Authors declare that they have no conflicts of interest.
Artificial Intelligence (AI) Disclosure
No artificial intelligence (AI) tools, including large language models or machine learning software, were used in the preparation, analysis, or presentation of this manuscript.
- Received January 16, 2026.
- Revision received February 28, 2026.
- Accepted March 6, 2026.
- Copyright © 2026 The Author(s). Published by the International Institute of Anticancer Research.
This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.











