- 1Department of Dermatology, Changde Hospital, Xiangya School of Medicine, Central South University (The First People’s Hospital of Changde City), Changde, Hunan, China
- 2Department of Dermatology, The Tenth Affiliated Hospital of Southern Medical University (Dongguan City People's Hospital), Dongguan, China
Objective: To evaluate the risk of local recurrence in patients with multiple actinic keratoses (AKs), to analyze the influence of host-, lesion-, and treatment-related factors, and to develop and validate a Cox regression model for predicting recurrence risk.
Methods: This retrospective study enrolled 148 patients with multiple AKs during January 2019-February 2022. Baseline characteristics, treatment modalities, and follow-up outcomes were collected. The Kaplan-Meier method was used to estimate cumulative recurrence rates. Cox regression analyses were performed to identify independent risk factors. The interaction effect between lesion count and treatment response was further assessed. A multivariate Cox predictive model was constructed, and its performance was validated using the 24-month calibration curve, Harrell’s C-index, and Brier score. A nomogram was developed for individualized risk prediction.
Results: The cumulative recurrence rates at 12, 24, and 36 months were 22.3%, 35.6%, and 44.7%, respectively. Multivariate analysis identified higher lesion count (11-20: HR = 2.39; > 20: HR = 2.96) and incomplete treatment response (HR = 2.43) as independent risk factors. Immunosuppression and regular sunscreen use were not significant. Although visual analysis suggested elevated risk with more lesions and incomplete response, their interaction term was not statistically significant. The model demonstrated moderate discrimination (C-index = 0.630) and good calibration (Brier score = 0.176). The nomogram enabled individualized risk estimation.
Conclusion: Lesion burden and incomplete treatment response significantly predict recurrence in multiple AKs patients. The developed Cox model and nomogram offer a clinically useful tool for identifying high-risk individuals and optimizing management strategies.
1 Introduction
Actinic keratosis (AK) is a chronic cutaneous photodamage condition induced by long-term ultraviolet (UV) radiation exposure, representing an early stage in the development of squamous cell carcinoma (SCC) or a precancerous lesion. Its histological hallmark is atypical hyperplasia of epidermal keratinocytes, clinically presenting as scaly or crusted patches (1). Epidemiological studies indicate that the prevalence of AK increases significantly with age, particularly in individuals with fair skin, chronic sun exposure, and those who are immunosuppressed (2). With an aging population and cumulative photodamage, AK has become one of the most common precancerous skin lesions worldwide.
Clinically, patients with multiple AKs often present with numerous lesions, and the affected skin areas frequently exhibit a state termed “field cancerization”. This concept posits that within sun-exposed skin, even apparently normal surface areas may harbor underlying genetic mutations and cellular damage, predisposing the field to frequent local recurrences and new lesion development (3). Research indicates that multiple AKs not only signify a higher cumulative burden of photodamage but are also closely associated with an increased risk of local recurrence and progression to SCC (4).
Current treatment modalities for AK are diverse, including cryotherapy, photodynamic therapy (PDT), topical 5-fluorouracil (5-FU), imiquimod, and various combination or sequential regimens (5). Although these methods can effectively clear lesions in the short term, recurrence remains a significant challenge in clinical management. Previous studies suggest that recurrence is influenced not only by the treatment modality but also by multiple factors, including host factors [e.g., immunosuppression, skin phototype, history of non-melanoma skin cancer (NMSC)], lesion characteristics (e.g., number, location, thickness/morphology), and sun protection behaviors (6). However, systematic studies assessing the risk of recurrence in patients with multiple AKs remain limited, particularly regarding long-term follow-up data stratified by lesion number and location characteristics. Therefore, this study, based on a retrospective analysis of clinical data from patients with multiple AKs, aims to systematically evaluate the relationship between host-, lesion-, and treatment-related characteristics and the risk of local recurrence.
2 Materials and methods
2.1 Study population
This single-center retrospective cohort study included 148 AK patients who visited the Department of Dermatology at our hospital and received standard treatment between January 2019 and February 2022. All cases were clinically and dermoscopically diagnosed by at least two senior dermatologists, with pathological confirmation when necessary. The study was approved by the Ethics Committee of Changde Hospital, Xiangya School of Medicine, Central South University (The First People’s Hospital of Changde City) and adhered to the principles of the Declaration of Helsinki. Due to the retrospective nature of the study and the absence of privacy risks, the Ethics Committee waived the requirement for informed consent. All data were de-identified and anonymized prior to analysis.
Inclusion criteria (1): Aged ≥ 18 years (2); Clinically confirmed multiple AKs (≥ 5 lesions per patient) (3); Received systematic and standardized initial treatment at our center (including local laser therapy, PDT, oral retinoids, imiquimod, or combination/sequential regimens) (4); Completed at least 12 months of clinical follow-up post-treatment with trackable recurrence status (5); Possessed complete clinical data providing the main study variables.
Exclusion criteria (1): Lacking key exposure or outcome information (2); Diagnosed with in situ or invasive SCC requiring surgical excision at initial presentation (3); Comorbid severe systemic diseases or conditions affecting follow-up (e.g., active malignancy, severe immunodeficiency) (4); Lost to follow-up or incomplete data preventing determination of recurrence time.
The primary outcome was the time to local recurrence, analyzed using the Cox proportional hazards regression model. The required number of events (E) was calculated using the formula: , where α is the two-sided significance level (set at 0.05, Z1−α/2 = 1.96), β is the type II error rate (set at 0.20, power = 80%), HR is the expected hazard ratio, and p is the proportion in the exposed group. Based on previous research on AK recurrence (7), HRs for various factors typically range from 1.8-2.0. Assuming HR = 2.0 and p = 0.5, approximately 65 events were required. Considering the adjustment for multiple confounders in the multivariate Cox model and Peduzzi’s rule (requiring ≥ 10 events per predictor variable), and planning to include 6–7 primary variables, the total number of events should ideally be ≥ 70. Given reported 1–3 year local recurrence rates of 25%-35% for multiple AKs, a sample size (N) of approximately 233 patients was theoretically needed to achieve 80% power. However, the final cohort comprised 148 patients. With an expected recurrence rate of 30%, this yielded approximately 44 events, providing moderate statistical power (≈ 70%) for analyzing primary risk factors, while subgroup and interaction analyses should be considered exploratory.
2.2 Data collection and variable definitions
Clinical data for all enrolled patients were sourced from the hospital’s electronic medical record system and outpatient follow-up records. Data extraction followed a standardized protocol, with dual independent entry and cross-verification to ensure accuracy. Data extraction covered the period from the initial visit to the last follow-up. Variables collected included:
2.2.1 Host factors
Age, sex, skin phototype (Fitzpatrick classification, grouped into light and dark types), immunosuppression status (organ transplant, long-term immunosuppressive medication), and history of NMSC.
2.2.2 Lesion factors
Number of lesions (categorized as 5–10, 11–20, and >20), primary lesion location (head/scalp, forearm/dorsal hand, trunk or lower limbs), and lesion morphology assessed by clinical examination and dermoscopy. Lesions were classified as thin or non-thin according to the Olsen clinical grading system: thin lesions corresponded to Olsen grade I actinic keratoses, characterized by slightly palpable erythematous macules or patches with minimal hyperkeratosis, whereas non-thin lesions included Olsen grade II–III lesions, presenting as moderately to markedly hyperkeratotic, clearly palpable lesions. When grading discrepancies occurred, consensus was reached by two senior dermatologists.
2.2.3 Treatment-related factors
Initial treatment modality (local laser therapy, PDT, oral retinoids, topical imiquimod, or combination/sequential therapy), treatment course completion (yes/no), and treatment response (complete or partial). Treatment response (8) was defined as: Complete clearance: disappearance of all clinically visible AKs in the treated area; Partial clearance: reduction in lesion count by ≥75%.
2.2.4 Protective factor
Regular sunscreen use (yes/no).
2.2.5 Recurrence definition and determination
Local recurrence was defined as the reappearance of a clinically (with dermoscopic correlation) or pathologically confirmed AK lesion within a 2 cm radius of the completely healed treatment area, assessed at least 1 month post-treatment completion. The recurrence time was calculated from the end of initial treatment to the date of first confirmed recurrence. For patients without recurrence, the last follow-up date or study cutoff date was used as the censoring time. Cumulative recurrence rates at 12, 24, and 36 months were recorded for survival analysis.
2.3 Data quality control
All data were independently verified and entered by two researchers, with discrepancies resolved by a senior dermatologist. The extent and pattern of missing data were assessed before analysis. Variables with more than 5% missing data included regular sunscreen use (6.8%) and lesion morphology (5.4%), while all other variables had missingness below 5%.
Missing values were assumed to be missing at random and were handled using multiple imputation by chained equations (MICE). A total of five imputations were generated, incorporating all variables included in the Cox regression models and the outcome indicator. Estimates from the imputed datasets were pooled according to Rubin’s rules for subsequent analyses.
2.4 Statistical analysis
Data were analyzed using R software (version 4.3.0). A two-sided P-value< 0.05 was considered statistically significant. Continuous variables were tested for normality. Normally distributed data were presented as mean ± standard deviation ( ± s) and compared using the t-test, while non-normally distributed data were presented as median (interquartile range) and compared using the Wilcoxon rank-sum test. Categorical variables were presented as frequencies and percentages, and compared using the χ² test or Fisher’s exact test. Local recurrence was the endpoint event. The Kaplan-Meier method was used to plot recurrence-free survival (RFS) curves and calculate cumulative recurrence rates at 12, 24, and 36 months, with inter-group comparisons using the log-rank test. Cox proportional hazards regression models were used to analyze factors influencing local recurrence. Univariate analysis included all clinical variables, calculating HRs and 95% confidence intervals (CIs). Variables with clinical importance or P< 0.10 in univariate analysis were considered for inclusion in the multivariable model, with careful attention to model parsimony. Subgroup and interaction analyses were performed to assess result consistency across populations. Considering the relatively limited number of recurrence events, the number of covariates included in the multivariable Cox model was restricted based on clinical relevance and events-per-variable considerations to minimize the risk of overfitting. Interaction analyses were performed for exploratory purposes only and were not used for model selection. The proportional hazards assumption was verified using Schoenfeld residual tests. All tests were two-sided.
3 Results
3.1 Baseline characteristics
A total of 148 patients with multiple AKs met the inclusion criteria. The cohort comprised 81 males (54.73%) and 67 females (45.27%), with a mean age of 69.43 ± 7.58 years. During follow-up, 54 patients (36.49%) experienced local recurrence. Patients were divided into recurrence (n = 54) and non-recurrence (n = 94) groups. No significant differences were observed between groups regarding age, sex, skin phototype, immunosuppression status, or history of NMSC (all P > 0.05). Lesion number was closely associated with recurrence: the recurrence group had a higher lesion count (15.76 ± 8.57 vs. 13.03 ± 7.81, P = 0.054). Stratified by lesion number category, the recurrence group had significantly higher proportions of patients with > 20 lesions and 11–20 lesions (P = 0.015). Among treatment-related factors, treatment course completion (P = 0.040) and treatment response (P = 0.002) differed significantly between groups. The recurrence group had a higher proportion of incomplete treatment courses (25.93% vs. 12.77%) and a significantly higher rate of partial clearance (57.41% vs. 30.85%) (Table 1).
3.2 Incidence and temporal distribution of local recurrence
Recurrence events were concentrated in the early post-treatment period. The cumulative recurrence rates at 12, 24, and 36 months were 22.3%, 35.6%, and 44.7%, respectively, indicating a gradually increasing risk over time. The median RFS was not reached (as over half of the patients did not experience recurrence during follow-up) (Table 2). The Kaplan-Meier curve showed that most recurrences occurred within the first 12 months post-treatment, after which the risk increase plateaued (Figure 1).
3.3 Analysis of risk factors for local recurrence
Using local recurrence as the dependent variable, Cox proportional hazards regression models were employed. Multivariate analysis included host-, lesion-, and treatment-related factors prespecified based on previous literature (9–11) and clinical relevance, along with variables with P< 0.10 in univariate analysis. Lesion count was entered into the Cox regression model as a categorical variable, with patients with 5–10 lesions serving as the reference group and those with 11–20 and >20 lesions modeled as separate categories.
In univariate Cox regression analysis, age, lesion number category (11–20 vs. 5–10), and treatment response were significantly associated with recurrence risk (all P< 0.05) (Table 3). Each 1-year increase in age was associated with a 3.6% increase in recurrence risk (HR = 1.036, 95% CI 1.001–1.072, P = 0.045). Patients with 11–20 lesions had a significantly higher risk of recurrence compared with those with 5–10 lesions (HR = 1.765, 95% CI 1.032–3.017, P = 0.038). Incomplete treatment response was also associated with an increased risk of recurrence (HR = 2.378, 95% CI 1.383–4.086, P = 0.002). Regular sunscreen use was not statistically significantly associated with recurrence risk in univariate analysis (HR = 0.636, 95% CI 0.350–1.154, P = 0.136).
In the multivariable Cox regression model, lesion burden remained one of the strongest independent predictors of local recurrence (Table 4). Compared with patients with 5–10 lesions, those with 11–20 lesions (HR = 2.390, 95% CI 1.269–4.503, P = 0.007) and those with more than 20 lesions (HR = 2.957, 95% CI 1.367–6.395, P = 0.006) had a significantly increased risk of recurrence. Incomplete treatment response also remained independently associated with recurrence risk after adjustment (HR = 2.431, 95% CI 1.387–4.263, P = 0.002).
Immunosuppression and regular sunscreen use, which were included in the multivariable model as clinically relevant covariates based on prior evidence, were not independently associated with recurrence risk after adjustment (immunosuppression: HR = 1.684, P = 0.192; regular sunscreen use: HR = 0.661, P = 0.180). These findings indicate that, after accounting for lesion burden and treatment response, regular sunscreen use did not demonstrate an independent protective effect against local recurrence in this cohort.
3.4 Subgroup and interaction effect analysis
To further explore differences in recurrence risk across clinical subgroups, stratified and interaction effect analyses were conducted. Stratified Cox regression by treatment response status showed that lesion count was positively associated with recurrence risk in both the complete clearance (HR = 1.037, 95% CI 0.996-1.080) and partial clearance (HR = 1.042, 95% CI 0.994-1.093) groups, suggesting increased risk with higher lesion counts, albeit with similar magnitudes of increase between groups.
Subsequently, an interaction term (lesion count × treatment response) was included in the Cox model. The interaction term was not statistically significant (B = 0.005, HR = 1.005, 95% CI 0.944-1.070, P = 0.874) (Table 5), indicating no significant interaction between lesion number and treatment response regarding their effect on recurrence risk. This suggests that treatment response did not significantly modify the effect of lesion number.
Table 5. Cox regression analysis of lesion count, treatment response, and their interaction on local recurrence risk.
Despite the non-significant interaction term, heat maps and marginal effect plots were generated to visualize trend variations across different ranges (Figures 2A, B). The heat map indicated that predicted recurrence risk gradually increased with lesion number in the partial clearance group, while the increase was relatively milder in the complete clearance group. The marginal effect plot further confirmed that the risk curve for the partial clearance group was consistently higher than that for the complete clearance group, but the slopes were similar, consistent with the non-significant interaction term. These findings collectively indicate that lesion number is a shared risk factor in both groups, but treatment response did not significantly alter the strength of its association with recurrence risk.
Figure 2. Interaction effect analysis between lesion count and treatment response on local recurrence risk. (A) Heat map of the interaction effect; (B) Marginal effect plot of predicted recurrence risk.
3.5 Model validation and predictive analysis
The discriminative ability and calibration of the multivariate Cox model were further validated. The Harrell’s C-index of the model was 0.640, indicating moderate predictive ability for distinguishing between patients with and without local recurrence. The 24-month calibration curve, constructed using bootstrap resampling (B = 200), showed good agreement between predicted and observed recurrence probabilities (Figure 3A), indicating good calibration. Furthermore, the model’s Brier score was 0.176, suggesting low overall prediction error and supporting the model’s stability and reliability. Based on the multivariate Cox model, a nomogram for predicting local recurrence risk was constructed (Figure 3B). The nomogram integrated key predictors such as immunosuppression status, lesion number category, treatment response, and regular sunscreen use. By summing the points assigned to each variable, the total points corresponded to an individualized probability of recurrence.
Figure 3. Predictive performance and risk visualization of the multivariate Cox model. (A) Calibration curve of the multivariate Cox model; (B) Nomogram for predicting local recurrence risk in patients with multiple AKs.
4 Discussion
AK is a clinical manifestation of field cancerization, characterized by multiple visible and subclinical lesions within a photodamaged field, carrying the potential for recurrence and progression to SCC (12–14). Although various treatments like cryotherapy can effectively clear individual lesions, recurrence is common, and risk varies considerably among patients. Therefore, assessing factors associated with recurrence and developing risk prediction models hold significant clinical value for guiding individualized follow-up and precise treatment.
Based on a real-world cohort of patients with multiple AKs, this study delineated the temporal pattern of local recurrence and identified key influencing factors. The cumulative recurrence rates at 12, 24, and 36 months were 22.3%, 35.6%, and 44.7%, respectively, with most events occurring within 12–24 months post-treatment, consistent with previous reports on the timeline of AK recurrence (15). Relevant European multicenter data also indicate that recurrences typically occur within the first two years post-treatment, often earlier and more concentrated in areas of greater photodamage (16, 17). Our observations further support that AK is a chronic photodamage disease based on field cancerization, where recurrence reflects both the reactivation of residual lesions and the continuous development of new lesions within the compromised skin field. This interpretation is consistent with prospective long-term studies that more comprehensively capture field cancerization, repeated treatment cycles, and cumulative keratinocyte carcinoma burden. Such studies indicate that recurrence over extended follow-up often reflects ongoing field failure rather than true relapse of a single treated lesion (18). Accordingly, the recurrence rates observed in our study likely represent a composite of residual lesion reactivation and early emergence of new lesions within a chronically damaged field.
Importantly, the treatment patterns represented in this cohort (2019–2022) should be interpreted in light of evolving real-world management strategies for actinic keratoses. In recent years, short-course, field-directed therapies with high adherence, such as tirbanibulin 1% ointment, have been increasingly adopted in clinical practice. Large real-world evidence has demonstrated favorable treatment completion rates and sustained effectiveness with tirbanibulin, particularly compared with longer topical regimens that are more susceptible to non-completion in routine care (19). Consequently, predictors of treatment completion and response identified in the present study may differ in magnitude when applied to contemporary cohorts treated predominantly with short-duration, high-adherence field therapies.
Mechanistically, UV radiation (particularly UVB) causes direct DNA damage to epidermal keratinocytes, leading to the accumulation of driver mutations in genes like TP53, NOTCH1, CDKN2A, and FGFR3, forming multiple latent mutant clones. These clones do not progress synchronously but expand gradually, driven by the local microenvironment, inflammatory responses, and declining skin repair capacity, giving rise to clinical and subclinical lesions (20). Previous high-throughput sequencing studies have revealed high similarity in mutational profiles and clonal architecture between AK lesions, photodamaged epidermis, and early cutaneous SCC, suggesting that AK is inherently a chronic, recurrent disease with polyclonal subpopulations (21). This implies that clinical recurrence likely reflects the re-expansion of previously residual or subclinical clones rather than entirely new events. The observed recurrence peak within the first two years in our study aligns with this clonal kinetic model.
In this study, the number of lesions was a core predictor of local recurrence, with higher numbers correlating with increased risk. Previous research indicates that lesion count indirectly reflects the extent of compromised skin and the degree of field cancerization, serving as a quantitative measure of photodamage burden and representing the number and competitive pressure of mutant clones (22). Consequently, even if surface lesions are cleared by initial treatment, deeper or marginal mutant clones may continue to expand, increasing recurrence potential. The marginal effect analysis in our study reinforces this point: the risk-amplifying effect of incomplete treatment response became more pronounced with higher lesion numbers, suggesting that patients with high lesion loads require greater attention to the adequacy of lesion clearance post-treatment.
Incomplete treatment response was another independent high-risk factor for recurrence, with a clear pathological basis. Although modalities like PDT, chemical peels, and cryotherapy effectively target abnormal keratinization, their depth of effect is limited. They may fail to eradicate lesions with hyperkeratosis, associated superficial dermal inflammation, or areas with heterogeneous mutations. The presence of numerous subclinical mutant clones in field-cancerized skin, which prior studies suggest have a higher potential for expansion after superficial therapies (23), explains the significantly higher recurrence risk in patients with incomplete response, consistent with our findings. In addition, lesion- and field-directed pre-conditioning represents a practical but unmeasured determinant of treatment effectiveness in real-world AK management. Randomized evidence has shown that keratolytic pretreatment, such as 10% urea applied before methyl aminolevulinate photodynamic therapy (MAL-PDT), significantly enhances protoporphyrin IX (PpIX) uptake and improves clinical outcomes, particularly in hyperkeratotic AKs (24). The lack of detailed information on pre-conditioning strategies in the present cohort may partly contribute to heterogeneity in treatment response and subsequent recurrence risk. Immunosuppression is a potential risk factor for recurrence but did not reach statistical significance in this cohort. Possible explanations include the low proportion of immunosuppressed patients, limiting statistical power. Furthermore, while sun protection mechanistically reduces cumulative photodamage, its subjective nature and the difficulty in reversing long-term photodamage might explain why regular sunscreen use did not show a significant protective effect in this study.
Regarding model performance, the multivariate Cox model demonstrated moderate discriminative ability (Harrell’s C-index = 0.638), consistent with the multifactorial nature of AK. The calibration curve indicated good predictive accuracy at 24 months, closely matching the observed recurrence rates. This suggests that the model reliably represents the risk structure for recurrence during medium-term follow-up. The nomogram further digitizes the model into a visual scoring system, potentially enabling risk stratification based on lesion number, treatment response, immune status, and sunscreen use, and may assist clinical risk assessment.
In terms of clinical application, the model’s value is threefold: First, it aids in formulating individualized follow-up strategies, suggesting shorter intervals for high-risk patients (e.g., those with > 20 lesions or incomplete response). Second, it can guide treatment selection; for instance, patients with extensive field cancerization might benefit from field-directed therapies with broader coverage and greater depth, such as 5-aminolevulinic acid PDT, imiquimod, or sequential/combination regimens. Third, it can be used for patient education and management, as quantified recurrence risk may improve adherence to sun protection and long-term management.
5 Study significance and limitations
This study underscores the central role of lesion burden and treatment response in AK recurrence, providing plausible clinical and biological explanations. It also proposes a predictive model for recurrence risk that may support individualized follow-up and risk assessment in patients with multiple AKs.
Nevertheless, several limitations should be acknowledged. First, this was a single-center study with a relatively limited sample size, which may introduce selection bias and limit generalizability. In addition, the absence of novel biomarkers, such as genetic mutation profiles or advanced skin imaging parameters (e.g., reflectance confocal microscopy), may have constrained the predictive performance of the model. Future studies should focus on external validation and explore the integration of multi-omics data and artificial intelligence–assisted imaging to enhance risk prediction. Second, although multivariable Cox regression was applied, the relatively limited number of recurrence events may increase the risk of model overfitting, particularly in subgroup, interaction, and nomogram analyses. Accordingly, findings from subgroup and interaction analyses should be interpreted as exploratory and hypothesis-generating rather than confirmatory. Finally, several limitations inherent to retrospective real-world cohorts warrant consideration. Local recurrence was operationally defined as the appearance of an AK lesion within a 2 cm radius of the treated area; in the context of field cancerization, this definition may partially capture newly emerging lesions rather than true local relapse, potentially leading to overestimation of recurrence. In addition, treatment-selection bias (confounding by indication) is unavoidable, as treatment choice and completion are influenced by lesion thickness, field severity, anatomical site, seasonality, and patient preference. Although multivariable adjustment was performed, residual confounding cannot be fully excluded.
6 Conclusion
While our Cox-based nomogram provides clinically relevant risk stratification for local recurrence in patients with multiple AKs, its generalizability should be interpreted with caution. External validation in contemporary cohorts incorporating evolving treatment paradigms and real-world confounders is warranted.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by Changde Hospital, Xiangya School of Medicine, Central South University (The First People’s Hospital of Changde City). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because the retrospective nature of the study and the absence of privacy risks.
Author contributions
ZL: Writing – original draft, Methodology, Visualization, Conceptualization. YT: Writing – review & editing, Data curation, Formal Analysis. PH: Writing – original draft, Data curation, Formal Analysis. JL: Visualization, Validation, Writing – review & editing. ZC: Resources, Writing – original draft, Supervision. MW: Writing – review & editing, Project administration, Software. JZ: Conceptualization, Investigation, Writing – review & editing, Supervision.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was supported by the 2025 Changde Municipal Scientific and Technological Innovation Guiding Plan Projects (No. 2025ZD246).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Li Z, Lu F, Zhou F, Song D, Chang L, Liu W, et al. From actinic keratosis to cutaneous squamous cell carcinoma: the key pathogenesis and treatments. Front Immunol. (2025) 16:1518633. doi: 10.3389/fimmu.2025.1518633
2. Worley B, Harikumar V, Reynolds K, Dirr MA, Christensen RE, Anvery N, et al. Treatment of actinic keratosis: a systematic review. Arch Dermatol Res. (2023) 315:1099–108. doi: 10.1007/s00403-022-02490-5
3. Bernal Masferrer L, Gracia Cazaña T, Bernad Alonso I, Álvarez-Salafranca M, Almenara Blasco M, Gallego Rentero M, et al. Topical immunotherapy for actinic keratosis and field cancerization. Cancers (Basel). (2024) 16(6):1113. doi: 10.3390/cancers16061133
4. Zou DD, Sun YZ, Li XJ, Wu WJ, Xu D, He YT, et al. Single-cell sequencing highlights heterogeneity and Malignant progression in actinic keratosis and cutaneous squamous cell carcinoma. Elife. (2023) 12:e85270. doi: 10.7554/eLife.85270
5. Arcuri D, Ramchatesingh B, Lagacé F, Iannattone L, Netchiporouk E, Lefrançois P, et al. Pharmacological agents used in the prevention and treatment of actinic keratosis: A review. Int J Mol Sci. (2023) 24(5):4989. doi: 10.3390/ijms24054989
6. Stockfleth E and Bégeault N and Delarue A. Intensity of local skin reactions during 5-fluorouracil treatment related to the number of actinic keratosis lesions: A post hoc, exploratory analysis. Dermatol Ther (Heidelb). (2022) 12:467–79. doi: 10.1007/s13555-021-00668-9
7. Dirschka T, Pellacani G, Micali G, Malvehy J, Stratigos AJ, Casari A, et al. A proposed scoring system for assessing the severity of actinic keratosis on the head: actinic keratosis area and severity index. J Eur Acad Dermatol Venereol. (2017) 31:1295–302. doi: 10.1111/jdv.14267
8. Kandolf L, Peris K, Malvehy J, Mosterd K, Heppt MV, Fargnoli MC, et al. European consensus-based interdisciplinary guideline for diagnosis, treatment and prevention of actinic keratoses, epithelial UV-induced dysplasia and field cancerization on behalf of European Association of Dermato-Oncology, European Dermatology Forum, European Academy of Dermatology and Venereology and Union of Medical Specialists (Union Européenne des Médecins Spécialistes). J Eur Acad Dermatol Venereol. (2024) 38:1024–47. doi: 10.1111/jdv.19897
9. Steeb T, Wessely A, Petzold A, Brinker TJ, Schmitz L, Schöffski O, et al. Long-term recurrence rates of actinic keratosis: A systematic review and pooled analysis of randomized controlled trials. J Am Acad Dermatol. (2022) 86:1116–9. doi: 10.1016/j.jaad.2021.04.017
10. Stockfleth E. The importance of treating the field in actinic keratosis. J Eur Acad Dermatol Venereol. (2017) 31 Suppl 2:8–11. doi: 10.1111/jdv.14092
11. Welzel J. From actinic keratosis to squamous cell carcinoma – answers to some open questions. Br J Dermatol. (2019) 180:699–700. doi: 10.1111/bjd.17637
12. Kim YS, Jung SH, Park YM, Park GS, Kim HS, Maeng LS, et al. Targeted deep sequencing reveals genomic alterations of actinic keratosis/cutaneous squamous cell carcinoma in situ and cutaneous squamous cell carcinoma. Exp Dermatol. (2023) 32:447–56. doi: 10.1111/exd.14730
13. Sotiriou E, Kiritsi D, CHaitidis N, Arabatzis M, and Lallas A and Vakirlis E. Daylight photodynamic therapy for actinic keratosis and field cancerization: A narrative review. Cancers (Basel). (2025) 17(6):1050. doi: 10.3390/cancers17061050
14. Paradisi A, Bocchino E, Mannino M, Gualdi G, D’Amore A, Traini DO, et al. The state of the art in the treatment of actinic keratosis and field cancerization: A narrative review. J Pers Med. (2025) 15(9):421. doi: 10.3390/jpm15090421
15. Reinhold U, Philipp-Dormston WG, Dirschka T, Ostendorf R, Aschoff R, Berking C, et al. Long-term follow-up of a randomized, double-blind, phase III, multi-centre study to evaluate the safety and efficacy of field-directed photodynamic therapy (PDT) of mild to moderate actinic keratosis using BF-200 ALA versus placebo and the BF-RhodoLED(R) lamp. J Eur Acad Dermatol Venereol. (2025) 39:1449–59. doi: 10.1111/jdv.20452
16. Agnetta V, Williamson S, Bisbee E, Torres A, Hooey L, Motaparthi K, et al. A retrospective review of histopathologic features associated with increased risk of recurrence of non-melanoma skin cancer after mohs micrographic surgery. J Clin Aesthet Dermatol. (2022) 15:27–9.
17. Kumar SV, Kumar V, and Mishra SK and Sati A. Clinicopathological profile, management and clinical outcomes in recurrent cases of ocular surface squamous neoplasia at a tertiary care centre. Int Ophthalmol. (2023) 43:371–80. doi: 10.1007/s10792-022-02432-0
18. Calzavara-Pinton P, Arisi M, Bettolini L, Rovaris S, Maione V, Soglia S, et al. Management of all patients with multiple actinic keratoses requires lifelong follow-up and repeated treatment cycles: A 12-year prospective observational study of a cohort of 81 patients. Dermatology. (2025), 1–9. doi: 10.1159/000548336
19. Carugno A, Nazzaro G, Argnani L, Bellinato F, Zappia E, Masala MV, et al. Tirbanibulin 1% real-world use in 679 actinic keratoses patients during warm season (SUMM-AK study). J Eur Acad Dermatol Venereol. (2025). doi: 10.1111/jdv.70217
20. Korecka K, Polańska A, and Dańczak-Pazdrowska A and Navarrete-Dechent C. Assessing field cancerization and actinic keratosis using ultraviolet-induced fluorescence dermatoscopy after the application of 5-aminolevulinic acid - An observational study. Photodiagnosis Photodyn Ther. (2024) 46:104056. doi: 10.1016/j.pdpdt.2024.104056
21. Thomson J, Bewicke-Copley F, Anene CA, Gulati A, Nagano A, Purdie K, et al. The genomic landscape of actinic keratosis. J Invest Dermatol. (2021) 141:1664–74.e7. doi: 10.1016/j.jid.2020.12.024
22. Malvehy J, Stratigos AJ, Bagot M, Stockfleth E, and Ezzedine K and Delarue A. Actinic keratosis: Current challenges and unanswered questions. J Eur Acad Dermatol Venereol. (2024) 38 Suppl 5:3–11. doi: 10.1111/jdv.19559
23. Sinclair R, Baker C, Spelman L, Supranowicz M, and MacMahon B. A review of actinic keratosis, skin field cancerisation and the efficacy of topical therapies. Australas J Dermatol. (2021) 62:119–23. doi: 10.1111/ajd.13447
24. Artelli GL, Cattaneo I, Rubelli L, Ariasi C, Bighetti S, Tomasi C, et al. Can urea 10% Promote photosensitizer uptake before MAL-PDT for the treatment of facial actinic keratoses? Results of a randomized clinical trial. Photodermatol Photoimmunol Photomed. (2025) 41:e70058. doi: 10.1111/phpp.70058
Keywords: actinic keratosis, Cox regression, nomogram, prediction model, recurrence
Citation: Liu Z, Tang Y, He P, Long J, Chen Z, Wan M and Zhu J (2026) Clinical characteristics and local recurrence risk in patients with multiple actinic keratoses: a retrospective clinical data analysis. Front. Oncol. 16:1770341. doi: 10.3389/fonc.2026.1770341
Received: 18 December 2025; Accepted: 12 January 2026; Revised: 09 January 2026;
Published: 11 February 2026.
Edited by:
Hao Chang, University of Wisconsin-Madison, United StatesReviewed by:
Tianshun Zhang, University of Minnesota Twin Cities, United StatesStefano Bighetti, University of Brescia, Italy
Copyright © 2026 Liu, Tang, He, Long, Chen, Wan and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianjian Zhu, bGlkb2N1dGVAMTYzLmNvbQ==; Miao Wan, OTcwNjI4MjcxQHFxLmNvbQ==
Zhaogui Liu1