ORIGINAL RESEARCH article

Front. Public Health, 04 June 2025

Sec. Infectious Diseases: Epidemiology and Prevention

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1580221

Worldwide trends and future projections of fungal skin disease burden: a comprehensive analysis from the Global Burden of Diseases study 2021

Deng Li&#x;Deng Li1Siqi Fan&#x;Siqi Fan1Haochen ZhaoHaochen Zhao2Jiayi SongJiayi Song1Linfen GuoLinfen Guo1Wei LiWei Li1Xuewen Xu
Xuewen Xu1*Qingfeng Li,
Qingfeng Li1,3*
  • 1Department of Plastic and Burns Surgery, West China Hospital, Sichuan University, Chengdu, China
  • 2Department of Urology, West China Hospital, Sichuan University, Chengdu, China
  • 3Department of Plastic and Reconstructive Surgery, Shanghai Ninth People's Hospital Affiliated to Shanghai Jiaotong University School of Medicine, Shanghai, China

Background: Fungal skin diseases represent pervasive global health concerns, predominantly arising from dermatophytes, yeasts, and molds.

Objective: This study aimed to estimate the disease burden associated with fungal skin diseases in 2021. Additionally, it sought to analyze trends from 1990 to 2021 and forecast future patterns.

Methods: This observational study first utilized data from the Global Burden of Disease (GBD) database covering the years 1990 to 2021. We specifically used data from GBD 2021 to evaluate the global incidence, prevalence, and disability-adjusted life years (DALYs), disaggregated by age, gender, socio-demographic index (SDI), and GBD regions. Linear regression models were then employed to identify temporal trends, estimating the annual percentage change. Cluster analysis examined disparities across 45 GBD regions. To forecast future disease burden, we applied the age-period-cohort model and the autoregressive integrated moving average model.

Conclusion: In 2021, there were approximately 1.73 billion global cases of fungal skin diseases. Males had higher age-standardized rates for incidence, prevalence, and DALYs compared to females. Age-specific analyses showed that although younger groups experienced the highest incidence rates, ASRs increased with age, especially among older populations. Regionally, low and middle SDI areas faced the greatest burden, with Asia having the highest incidence and Oceania the lowest. Projections suggest significant increases in incidence, prevalence, and DALYs, notably in middle- and low-income regions. These results highlight meaningful spatiotemporal disparities in fungal skin diseases and emphasize the need for strategic allocation of resources to mitigate these challenges and reduce the growing burden across various global populations.

1 Introduction

Fungal infections, specifically dermatomycoses and fungal skin diseases, form a distinct category of skin disorders primarily caused by dermatophytes, yeasts, and molds. As of 2017, these fungal skin conditions became the most prevalent skin disorders worldwide, representing 10.09% of all skin diseases (1, 2). In a 2016 global comparison of 328 diseases and injuries, fungal skin diseases were found to have the fourth highest incidence, with a staggering 2.1 billion cases reported. This prevalence underscores the significant socioeconomic impact they impose (2).

Fungal skin diseases are caused by diverse pathogens with distinct ecological niches. Superficial infections such as tinea corporis in tropical South Asia and Sub-Saharan Africa are predominantly driven by dermatophytes (3), which thrive in humid climates. In contrast, deep infections such as chromoblastomycosis in rural Brazil (4), and mycetoma in Sudan (5) are linked to soil-borne molds, reflecting occupational exposure to contaminated environments. High-income regions face unique challenges: Candida auris outbreaks in ICUs and azole-resistant Aspergillus fumigatus in agricultural settings underscore the role of healthcare and environmental selection pressures (6). Beyond geographic location and climate, other factors such as the host’s immune status (7), the pathogenicity of infectious agents (2, 8), and access to medical services also influence the incidence of fungal infections among populations (9).

Recent studies have highlighted the escalating threat of antifungal resistance. For instance, the emergence of Trichophyton indotineae—a terbinafine-resistant dermatophyte strain originating from South Asia—has led to near-epidemic outbreaks globally (10, 11). Concurrently, the global spread of Candida auris and azole-resistant Aspergillus fumigatus underscores the urgency of monitoring resistance patterns (12). Furthermore, host genetic predispositions (e.g., TLR2/4 polymorphisms) and immunological vulnerabilities (e.g., HIV coinfection or diabetes) have been increasingly recognized as critical modulators of disease susceptibility (7, 8). These factors collectively complicate therapeutic management and necessitate updated epidemiological frameworks.

Advances in diagnostic technologies, such as PCR-based assays and MALDI-TOF mass spectrometry, have improved species-level identification and antifungal susceptibility profiling (13, 14). However, these tools remain underutilized in high-burden regions, highlighting a critical implementation gap. Meanwhile, epidemiological studies have identified hotspots of fungal transmission linked to urbanization, overcrowding, and occupational exposures—factors increasingly relevant in rapidly developing economies (1517). These studies highlight the profound implications these diseases have on public health, notably affecting daily life quality and imposing financial strains. The impact is often measured using disability-adjusted life years (DALYs), which combine the potential years lost to early mortality with years lived with disability. Additionally, the socio-demographic index (SDI) provides an encompassing metric of development by accounting for factors such as income, education, and fertility rates across various regions.

Despite progress in researching fungal skin diseases, several challenges remain unresolved. One major issue is the unequal representation in data collection globally, especially in low- and middle-income countries (18), hindering the accurate assessment of the global disease burden.

This study leverages updated GBD 2021 data to address these gaps, providing a comprehensive analysis of global trends from 1990 to 2021 and projecting disease trajectories through 2046. Previous studies have used Joinpoint regression analysis to review temporal trends in these diseases (19). Building on this work, we elucidate the evolving demographic and geographic patterns of fungal skin diseases by integrating age-period-cohort and ARIMA models, emphasizing the synergistic roles of socioeconomic inequity, environmental factors, and healthcare infrastructure. Our findings aim to inform targeted prevention strategies and advocate for equitable resource allocation in an era of escalating antifungal resistance and climatic volatility.

2 Materials and methods

2.1 Data source

This study utilizes data sourced from the GBD 2021 database. The GBD initiative, which commenced in 1991, has provided comprehensive empirical evaluations of global health statuses over the past 30 years (20, 21). Each iteration of the study has broadened in scope and intricacy, incorporating a wider array of causes and risks, expanding geographical coverage, and refining age-specific analyses. For the 2021 iteration, data from 328,938 unique sources were analyzed, resulting in the generation of over 60.7 billion estimates. These estimates addressed 25 distinct age groups, from birth to 95 years and older, across 204 countries and territories (22). The territories were further categorized into 21 geographic regions and stratified into five Socio-Demographic Index (SDI) categories—low, low-middle, middle, high-middle, and high SDI—as well as into seven broader super-regions. Separate estimates were generated for males, females, and overall populations by an international team of collaborators dedicated to rigorous data evaluation and analysis (23).

This section elaborates on the methodologies for data collection and estimation utilized by the GBD 2021. Detailed descriptions of these methodologies, including data search and selection processes, have been thoroughly documented in previous GBD publications. The procedures for data collection and analysis adhered closely to the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) (24). The study employed DALYs to quantify the burden of diseases, calculated as the cumulative years lost to illness, disability, or premature mortality in a specific demographic.

The specific methodologies of GBD 2021, along with comprehensive analyses and model specifications, are detailed in other publications. The standard life expectancy utilized in this study was derived from the lowest age-specific mortality rates recorded in any country.

2.2 Statistical analysis

To analyze temporal trends in fungal skin diseases globally from 1990 to 2021, we estimated the Estimated Annual Percentage Change (EAPC). The detailed methodology for EAPC calculation has been extensively documented in previous studies. In summary, a regression model was employed to capture the yearly percentage variation in the Age-Standardized Rate (ASR) (25), using the natural logarithm of the fitted rate equation: y = α + βx + ɛ, where y denotes ln (ASR) and x represents the calendar year. The EAPC of ASR is computed as 100 × (exp(β) - 1). These analyses were conducted using the R programming language (R Core Team, version 3.5.3, Vienna, Austria), with statistical significance assessed at a two-sided p-value threshold of < 0.05 (26).

We implemented an Age-Period-Cohort (APC) model to forecast the burden of fungal skin diseases (incidence, prevalence, and DALYs) across age groups from 2022 to 2046. This approach quantifies three temporal drivers: (1) age effects (biological susceptibility and lifetime exposure accumulation), (2) period effects (population-level changes like diagnostic guidelines revisions), and (3) cohort effects (generational risks from early-life antifungal stewardship programmes) (27). The log-linear model was specified as: Y = μ + α(age) + β(period) + γ(cohort) + ε, where α, β, γ represent partial regression coefficients for each temporal dimension, μ is the intercept, and ε accommodates overdispersion in Poisson-distributed outcomes (28).

For sensitivity analysis, we contrasted APC projections with seasonal Autoregressive Integrated Moving Average (ARIMA) Model (1, 2) forecasts (29), with autoregressive (lag = 2) and moving average (lag = 1) terms selected via autocorrelation function minimization (30).

3 Results

3.1 Global disease burden attributable to fungal skin diseases in 2021

In 2021, there were approximately 1,729,224,009 cases of fungal skin diseases (95% uncertainty intervals (UI): 1,562,709,709 to 1,894,666,053) (Table 1). Notably, the incidence rates were similar across both genders, as illustrated in Supplementary Figure S1. However, an examination of prevalence and DALYs revealed that males experienced higher values than females in terms of absolute case numbers and ASRs. Specifically, incidence, prevalence, and DALYs were 1.01, 1.07, and 1.09 times greater in males compared to females. Likewise, the ASIR ASPR and ASDR were 1.06, 1.11, and 1.12 times higher in males, respectively.

Table 1
www.frontiersin.org

Table 1. Temporal trends of fungal skin diseases by gender, SDI and GBD regions, 1990 to 2021.

In terms of age distribution, the global ASRs demonstrated a notable increase across age groups ranging from 25 to 95+, peaking in the 95 + age group with an ASR for DALYs of 232.32 per 100,000 (95% uncertainty interval [UI]: 90.88–497.8) in 2021. Conversely, for individuals from birth to 24 years old, the global age-standardized rates (ASRs) initially rose, peaking in the 5–9 age group, and then began to decrease.

However, when considering the number of cases, a different pattern emerged. The younger age group (5–9 years) exhibited the highest incidence (140,377,507 cases, 95% UI: 103,902,413-182,566,892) and prevalence (63,559,849 cases, 95% UI: 46,563,538-87,068,326), with a gradual decrease observed in older age groups. This trend was also reflected in the DALYs, which were significantly higher among the younger demographic (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. The age distribution of ASRs (A) and numbers (B) of fungal skin diseases-related incidence, prevalence and DALYs in 2021.

At the SDI region level, the high SDI regions exhibited the lowest ASRs for incidence and prevalence in 2021, with values of 16,148 (95% uncertainty interval [UI]: 14,687-17,628) and 5,117 (95% UI: 4,726-5,561) per 100,000, respectively. However, the ASRs increased significantly towards the low SDI regions. A similar pattern was observed for DALYs, with the highest burden noted in the low SDI regions (ASR of 71.8, 95% UI: 29.26–149.06) (Table 1).

The middle SDI region reported the highest number of incidence cases at 503,290,634 (95% UI: 454,024,197-554,675,359), prevalence cases at 173,183,772 (95% UI: 155,963,056-192,835,881), and DALYs cases at 963,136 (95% UI: 393,633-1,990,445). Conversely, the high SDI region had the lowest number of cases, followed by the high-middle SDI region. Notably, the low SDI region had fewer incidence cases but more DALYs and prevalence cases compared to the low-middle SDI region.

Among the 45 GBD regions in 2021, Asia ranked first in terms of fungal skin disease-related incidence (883,292,233, 95% uncertainty interval [UI]: 796,351,242-972,024,863), prevalence (296,194,797, 95% UI: 266,859,176-329,067,591), and DALYs (1,647,941, 95% UI: 675,987-3,379,970). Conversely, Oceania ranked last for incidence (2,329,919, 95% UI: 2,081,626-2,608,322), prevalence (739,478, 95% UI: 652,407-832,720), and DALYs (4,155, 95% UI: 1,687-8,565) (Figure 2).

Figure 2
www.frontiersin.org

Figure 2. GBD Region-specific analysis of fungal skin diseases burden in 2021. ASRs (A) and number of cases (B) of incidence, prevalence, and DALYs across different GBD regions.

Regarding the corresponding ASRs, Western Africa emerged as the top GBD region for incidence (40,345, 95% UI: 36,125-44,592), prevalence (18,033, 95% UI: 16,222-20,064), and DALYs (101, 95% UI: 41–209). In contrast, High-income North America ranked lowest for incidence (8,439, 95% UI: 7,866-8,990), prevalence (2,506.94, 95% UI: 2,407-2,606), and DALYs (13.67, 95% UI: 6–28). Further, Western Sub-Saharan Africa and Eastern Sub-Saharan Africa reported elevated incidence and DALYs compared to most other GBD regions. Additionally, Sub-Saharan Africa – WB and Eastern Sub-Saharan Africa exhibited higher prevalence rates than the majority of other GBD regions.

3.2 Global distribution in 2021

Supplementary Table S1 provides the global distribution of fungal skin diseases varied significantly across different regions and health metrics in 2021. Ethiopia reported the highest ASIR of 47,245 (95% uncertainty interval [UI]: 42,242-52,965), followed by Nigeria and Mali. Additionally, Mali exhibited the highest ASRs for prevalence and DALYs, with values of 23,697 (95% UI: 20,973-26,339) and 133 (95% UI: 54–274), respectively.

In 2021, the United States recorded the lowest ASRs for both incidence (8,401; 95%UI: 7,814–8,916) and disability-adjusted life years (DALYs) (14; 95% UI: 6–28), with Canada and Greenland following. For prevalence, Canada had the lowest ASR (2,490; 95% UI: 2,245–2,782), followed by the United States. Contrastingly, India reported the highest numbers in 2019 for incidence (265,162,118; 95% UI: 235,600,132–294,116,264), prevalence (87,230,762; 95% UI: 77,093,668–98,204,047), and DALYs (485,519; 95% UI: 197,480–994,783), with China and Nigeria next in line. Remarkably, the numbers for incidence, prevalence, and DALYs cases were nearly negligible in Tokelau, followed by Niue.

3.3 Global trends from 1990 to 2021

Supplementary Figure S2 presents the results of a cluster analysis based on the EAPC values of ASIR, ASPR and ASDR of fungal skin diseases from 1990 to 2021. The analysis reveals distinct trends among various regions.

Specifically, regions such as Southern Africa, Western Europe, Southern Sub-Saharan Africa, and North America exhibit significant decreases in ASRs, denoted by purple branches. In contrast, several regions, including Eastern Sub-Saharan Africa, Eastern Africa, Central Africa, and Central Sub-Saharan Africa, demonstrate significant increases, indicated by blue branches. Despite these notable variances, regions like the Americas (including both North and South America) show localized stability or minor decreases, as highlighted in green.

The majority of regions worldwide exhibit an increase ranging from 0 to 200% across all three metrics, as predominantly indicated by the light blue shading in all maps (Supplementary Figure S3). Notable exceptions to this general trend include Eastern Europe and parts of Southern Africa. Specifically, Eastern Europe demonstrates substantial decreases in these metrics, while parts of Southern Africa show increases of up to 600–800%.

Figure 3 and Table 1 illustrates the global distribution of estimated annual percentage changes (EAPCs) for ASIR, ASPR, and ASDR of fungal skin diseases from 1990 to 2021.

Figure 3
www.frontiersin.org

Figure 3. The global distribution of EAPCs for age-standardized incidence, prevalence and DALYs rates of fungal skin diseases from 1990 to 2021. (A) ASIR; (B) ASPR; (C) ASDR.

Figure 3A displays significant regional variability, with notable declines observed in parts of North America, Europe, and East Asia. Conversely, several African and South American nations exhibit increases in ASIR. Figure 3B (ASPR) generally shows a more widespread decline across various regions, including North America, parts of Europe, and regions in Asia, with isolated increases noted in certain African and Latin American countries. In contrast, Figure 3C (ASDR) highlights substantial increases in Eastern Europe and parts of Central Asia, which contrast sharply with decreases observed in other global regions.

3.4 Prediction of future trends of fungal skin diseases

According to projections using the APC model for 2022–2046, both males and females are expected to see a rise in the number of incidence, prevalence, and DALYs cases. This increase is accompanied by a corresponding growth in ASRs across all categories. Specifically, the anticipated growth rates for incidence, prevalence, and DALYs cases are 34.3, 32.0, and 30.1%, respectively. There are distinct gender-specific trends observed, with males consistently displaying higher ASRs across all metrics compared to females.

According to Figure 4, the number of prevalence and DALYs cases is projected to continue being higher in males than in females in the future. Specifically, the number of prevalence and DALYs cases is forecasted to increase by 35.0 and 33.0% for females, and by 29.1 and 27.4% for males, respectively. Notably, females are expected to experience a faster growth rate in the number of prevalence and DALYs cases, outpacing that of males. However, the trend in incidence cases will be more intricate. From 2022 to 2033, more incidence cases are forecasted in males than in females, but this trend is expected to reverse from 2034 to 2046. In general, females are expected to see a faster relative increase in incidence cases by 37.0% compared to males (31.7%) during this timeframe.

Figure 4
www.frontiersin.org

Figure 4. Prediction by age-Period-Cohort (APC) model for global trends in ASR and number of cases of fungal skin diseases from 2022 to 2046. (A) ASIR and number of incidence cases; (B) ASR and number of prevalence cases; (C) ASDR and number of DALYs cases.

In line with the APC model findings, the ARIMA model also predicts an upward trend in incidence, prevalence, and DALYs cases for both genders, as well as corresponding ASRs (Figure 5). Globally, incidence and DALYs cases are estimated to rise by 37.2 and 38.6% from 2022 to 2050. Unlike the APC model’s predictions, the ARIMA model forecasts that the number of incidence cases in females will exceed that of males by 2024. Moreover, prevalence cases are anticipated to increase by 40.4%, predominantly impacting males.

Figure 5
www.frontiersin.org

Figure 5. Prediction by ARIMA Model for global trends in ASR and number of cases of fungal skin diseases from 2022 to 2046. (A) ASIR and number of incidence cases; (B) ASR and number of prevalence cases; (C) ASDR and number of DALYs cases.

4 Discussion

This study comprehensively analyzes the global burden of fungal skin disease from 1990 to 2021, highlighting an increasing trend of worldwide fungal skin disease, accounting for approximately 8% of total DALYs owing to skin disease.

The increasing global burden of fungal skin diseases over the past three decades stems from multifaceted environmental, demographic, and healthcare challenges. Climate change has created favorable conditions for fungal proliferation and adaptation, with rising temperatures and humidity enhancing pathogen survival and resistance mechanisms (31). Demographic shifts, particularly population aging in regions like East Asia, have amplified susceptibility due to age-related immune decline and comorbidities (32). Healthcare disparities persist in low-resource settings, where limited diagnostic capacity and fragmented surveillance systems delay detection and exacerbate transmission (33). Rapid urbanization and population growth have intensified overcrowded living conditions, facilitating fungal spread in densely populated areas (34). Concurrently, expanding immunocompromised populations, driven by HIV/AIDS and immunosuppressive therapies, have heightened vulnerability to severe and opportunistic fungal infections (4). Paradoxically, hygiene practices, including prolonged use of protective equipment during the COVID-19 pandemic, disrupted skin integrity, increasing risks of superficial infections (35). Diagnostic and therapeutic gaps remain critical, as inadequate access to mycological testing and emerging antifungal resistance undermine treatment efficacy (33). Finally, ecological disruptions, such as pollution-induced microbiome alterations and zoonotic spillover risks, underscore evolving interactions between pathogens, hosts, and environments (36).

Fungal skin diseases result from a complex interplay of pathogen virulence, host susceptibility, and environmental factors. Host genetic factors play a crucial role in infection outcomes. Our study shows gender-specific disparities in ASIR, potentially linked to biological factors, such as androgen-influenced sebum composition, and behavioral risks (37), like occupational exposure in males (38). Environmental conditions also influence disease patterns—tropical climates support thermophilic species such as Aspergillus fumigatus (39), whereas overcrowding in low-SDI regions promotes the spread of anthropophilic species like Trichophyton tonsurans (40). Importantly, climate change is altering fungal distributions, for instance, the suitable habitats of Colletotrichum spp. expand northward with rising temperatures in southern China, which may indirectly increase the risk of human contact - based skin infections through agricultural activities (31, 41).

In 2021, males exhibited higher ASIR, ASPR and ASDR compared to females, likely driven by occupational exposures, hormonal influences, and skin physiology, emphasizing the need for gender-sensitive prevention strategies. Age-stratified analyses further identify distinct risk profiles: the older adult (≥80 years) experience disproportionately high burden due to immunosenescence and comorbidities, while children under 9 years old face early-life susceptibility linked to hygiene practices and immature immune defenses. These dual demographic peaks underscore the necessity of age-tailored interventions, ranging from microbiome-modulating probiotics in geriatric care to school-based screening programs in pediatric populations.

The global burden of fungal skin diseases demonstrates distinct regional patterns shaped by climatic and socioeconomic drivers. Although tropical climates in South Asia and Sub-Saharan Africa naturally favor fungal growth through sustained humidity, our analysis identifies poverty-related factors as equally critical. Low SDI regions exhibit the highest ASRs, underscoring how poverty amplifies exposure through overcrowded living conditions, inadequate footwear usage, and limited access to clean water and healthcare (33). South Asia and Sub-Saharan Africa’s surveillance gaps likely mask true disease prevalence, as resource-limited systems prioritize acute infections over dermatological monitoring—a cycle that perpetuates underfunding of skin health programs (42). Cultural practices like shared bathing facilities may further interact with environmental and economic factors to sustain transmission (43). Middle-SDI regions bear the greatest absolute caseload, driven by rapid urbanization concentrating populations in informal settlements with poor sanitation, combined with healthcare systems that enable passive case reporting without effective interventions (44).

While Sub-Saharan Africa is emblematic of the interplay between poverty and fungal burden, the continent’s cultural and ecological diversity demands a nuanced subregional perspective. In general, tinea capitis, tinea corporis, tinea pedis, and onychomycosis are the main manifestations, and immunosuppressed patients (e.g., HIV-infected patients) are prone to develop severe infections. In sub-Saharan Africa, where HIV is highly prevalent, the problem of dual infection is a concern (45), in addition, tinea capitis is high in school-age children in this area (46), and the impact of forest degradation and wood fuel use on sanitary conditions exacerbates the risk of infection (47); North Africa has a high incidence but limited diagnostic resources; In places such as Ethiopia in East Africa, the high infection rate is closely related to hot and humid environments, poor sanitation and animal contact (48). In South Africa, tinea capitis is the most common subtype, followed by tinea corporis and onychomycosis (43), Terbinafine-resistant strains of ringworm have been reported in parts of South Africa, leading to an increase in cases of refractory infection (49).

Effective mitigation requires region-specific strategies that bridge biomedical and structural interventions. In low-SDI settings, scaling up point-of-care diagnostics and community-led hygiene initiatives could reduce transmission. Moreover, leveraging digital health tools for improved disease monitoring and resource allocation can help bridge existing healthcare gaps, particularly in high-prevalence, resource-constrained regions. High-SDI regions can prioritize antifungal stewardship and hospital decontamination protocols to curb nosocomial outbreaks. Emerging therapies, such as antifungal peptides and nanoparticle-based drug delivery systems enhancing drug penetration, offer promise but require validation in real-world settings (50, 51).

This study’s limitations, including the aggregation of superficial and deep fungal infections, highlight the need for pathogen-specific subtyping in future research. Molecular tools like whole-genome sequencing could disentangle the epidemiology of dermatophytes versus invasive molds, as demonstrated in recent clade analyses of C. auris (52). Furthermore, the studies available for comparison across various GBD regions might be limited by the geographic scope, leading to certain populations being over- or under-represented relative to their actual population size. Additionally, the ability to report on the worldwide incidence of fungal diseases might be hindered by factors such as the absence of mandatory reporting for these diseases, suboptimal performance of diagnostic tests, and low levels of clinician awareness in particular areas. Future research should aim to bridge this gap, ensuring more comprehensive data collection in regions heavily burdened by fungal skin diseases. Longitudinal studies assessing climate-driven shifts in fungal biogeography and resistance patterns will be critical to informing adaptive public health policies.

5 Conclusion

This study demonstrates an escalating global burden of fungal skin diseases, characterized by significant projections of increased incidence, prevalence, and DALYs until 2046. These findings emphasize the imperative need for intensified global health strategies, particularly in regions with low SDI and tropical areas. Targeted interventions for older adults and very young individuals, as well as strategic resource allocation, are crucial to mitigate the growing impact of fungal skin diseases.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://vizhub.healthdata.org/gbd-results/.

Author contributions

DL: Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing. SF: Methodology, Writing – original draft. HZ: Conceptualization, Writing – original draft. JS: Writing – original draft. LG: Conceptualization, Writing – review & editing. WL: Conceptualization, Methodology, Writing – review & editing. XX: Conceptualization, Writing – review & editing. QL: Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

The authors gratefully acknowledge the Global Burden of Disease Collaborative Network and the Institute for Health Metrics and Evaluation (IHME) for providing open access to the GBD dataset through their online results tool (https://vizhub.healthdata.org/gbd-results/). The use of these data is governed by IHME's terms of service and the Creative Commons Attribution 4.0 International License (CC-BY 4.0). We extend our sincere gratitude to Dr. Binbin Su for his rigorous statistical review, particularly in validating trend analysis methodologies and refining projection models. We also deeply thank Dr. Jingyi Li for her invaluable dermatological expertise. Dr. Li critically reviewed the clinical and mycological interpretations, significantly enhancing the depth and accuracy of the discussion.

Conflict of interest

The authors declare that the research 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 authors declare that no Gen AI was used in the creation of this manuscript.

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.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1580221/full#supplementary-material

References

1. Urban, K, Chu, S, Scheufele, C, Giesey, RL, Mehrmal, S, Uppal, P, et al. The global, regional, and national burden of fungal skin diseases in 195 countries and territories: a cross-sectional analysis from the global burden of disease study 2017. JAAD Int. (2021) 2:22–7. doi: 10.1016/j.jdin.2020.10.003

PubMed Abstract | Crossref Full Text | Google Scholar

2. Liu, X, Zhang, Y, Hong, Y, and Zhang, H. Global burden of fungal skin diseases: an update from the global burden of diseases study 2019. Mycoses. (2024) 67:e13770. doi: 10.1111/myc.13770

PubMed Abstract | Crossref Full Text | Google Scholar

3. Barac, A, Stjepanovic, M, Krajisnik, S, Stevanovic, G, Paglietti, B, and Milosevic, B. Dermatophytes: update on clinical epidemiology and treatment. Mycopathologia. (2024) 189:101. doi: 10.1007/s11046-024-00909-3

PubMed Abstract | Crossref Full Text | Google Scholar

4. Seas, C, and Legua, P. Mycetoma, chromoblastomycosis and other deep fungal infections: diagnostic and treatment approach. Curr Opin Infect Dis. (2022) 35:379–83. doi: 10.1097/qco.0000000000000870

PubMed Abstract | Crossref Full Text | Google Scholar

5. Litvintseva, AP, Bakhiet, S, Gade, L, Wagner, DD, Bagal, UR, Batra, D, et al. Genomics and metagenomics of Madurella mycetomatis, a causative agent of black grain mycetoma in Sudan. PLoS Negl Trop Dis. (2022) 16:e0010787. doi: 10.1371/journal.pntd.0010787

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kriegl, L, Egger, M, Boyer, J, Hoenigl, M, and Krause, R. New treatment options for critically important WHO fungal priority pathogens. Clin Microbiol Infect. (2024) 31:922–30. doi: 10.1016/j.cmi.2024.03.006

PubMed Abstract | Crossref Full Text | Google Scholar

7. Cifaldi, C, Ursu, GM, D'Alba, I, Paccoud, O, Danion, F, Lanternier, F, et al. Main human inborn errors of immunity leading to fungal infections. Clin Microbiol Infect. (2022) 28:1435–40. doi: 10.1016/j.cmi.2022.06.031

PubMed Abstract | Crossref Full Text | Google Scholar

8. Deng, R, Wang, X, and Li, R. Dermatophyte infection: from fungal pathogenicity to host immune responses. Front Immunol. (2023) 14:1285887. doi: 10.3389/fimmu.2023.1285887

PubMed Abstract | Crossref Full Text | Google Scholar

9. Martinez-Rossi, NM, Peres, NTA, Bitencourt, TA, Martins, MP, and Rossi, A. State-of-the-art dermatophyte infections: epidemiology aspects, pathophysiology, and resistance mechanisms. J Fungi. (2021) 7:629. doi: 10.3390/jof7080629

PubMed Abstract | Crossref Full Text | Google Scholar

10. Nenoff, P, Verma, SB, Vasani, R, Burmester, A, Hipler, UC, Wittig, F, et al. The current Indian epidemic of superficial dermatophytosis due to Trichophyton mentagrophytes-a molecular study. Mycoses. (2019) 62:336–56. doi: 10.1111/myc.12878

PubMed Abstract | Crossref Full Text | Google Scholar

11. Ebert, A, Monod, M, Salamin, K, Burmester, A, Uhrlaß, S, Wiegand, C, et al. Alarming India-wide phenomenon of antifungal resistance in dermatophytes: a multicentre study. Mycoses. (2020) 63:717–28. doi: 10.1111/myc.13091

Crossref Full Text | Google Scholar

12. Lockhart, SR, Chowdhary, A, and Gold, JAW. The rapid emergence of antifungal-resistant human-pathogenic fungi. Nat Rev Microbiol. (2023) 21:818–32. doi: 10.1038/s41579-023-00960-9

PubMed Abstract | Crossref Full Text | Google Scholar

13. Begum, J, Mir, NA, Lingaraju, MC, Buyamayum, B, and Dev, K. Recent advances in the diagnosis of dermatophytosis. J Basic Microbiol. (2020) 60:293–303. doi: 10.1002/jobm.201900675

PubMed Abstract | Crossref Full Text | Google Scholar

14. Ghenciu, LA, Faur, AC, Bolintineanu, SL, Salavat, MC, and Maghiari, AL. Recent advances in diagnosis and treatment approaches in fungal keratitis: a narrative review. Microorganisms. (2024) 12:161. doi: 10.3390/microorganisms12010161

PubMed Abstract | Crossref Full Text | Google Scholar

15. Arenas, R, Sánchez-Cardenas, CD, Ramirez-Hobak, L, Ruíz Arriaga, LF, and Vega Memije, ME. Sporotrichosis: from KOH to molecular biology. J Fungi. (2018) 4:62. doi: 10.3390/jof4020062

PubMed Abstract | Crossref Full Text | Google Scholar

16. Casadevall, A, Kontoyiannis, DP, and Robert, V. On the emergence of Candida auris: climate change, azoles, swamps, and birds. MBio. (2019) 10:19. doi: 10.1128/mBio.01397-19

PubMed Abstract | Crossref Full Text | Google Scholar

17. Misseri, G, Ippolito, M, and Cortegiani, A. Global warming "heating up" the ICU through Candida auris infections: the climate changes theory. Crit Care. (2019) 23:416. doi: 10.1186/s13054-019-2702-4

Crossref Full Text | Google Scholar

18. Weber, AM, Gupta, R, Abdalla, S, Cislaghi, B, Meausoone, V, and Darmstadt, GL. Gender-related data missingness, imbalance and bias in global health surveys. BMJ Glob Health. (2021) 6:405. doi: 10.1136/bmjgh-2021-007405

PubMed Abstract | Crossref Full Text | Google Scholar

19. Qin, Q, Su, J, Liu, J, Chen, R, Wei, W, Yuan, Z, et al. Global, regional, and national burden of fungal skin diseases in 204 countries and territories from 1990 to 2021: an analysis of the global burden of disease study 2021. Mycoses. (2024) 67:e13787. doi: 10.1111/myc.13787

PubMed Abstract | Crossref Full Text | Google Scholar

20. Rudd, KE, Johnson, SC, Agesa, KM, Shackelford, KA, Tsoi, D, Kievlan, DR, et al. Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the global burden of disease study. Lancet. (2020) 395:200–11. doi: 10.1016/s0140-6736(19)32989-7

PubMed Abstract | Crossref Full Text | Google Scholar

21. GBD 2019 Diseases and Injuries Collaborators. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the global burden of disease study 2019. Lancet. (2020) 396:1204–22. doi: 10.1016/s0140-6736(20)30925-9

PubMed Abstract | Crossref Full Text | Google Scholar

22. GBD 2021 Diabetes Collaborators. Global, regional, and national burden of diabetes from 1990 to 2021, with projections of prevalence to 2050: a systematic analysis for the global burden of disease study 2021. Lancet. (2023) 402:203–34. doi: 10.1016/s0140-6736(23)01301-6

PubMed Abstract | Crossref Full Text | Google Scholar

23. GBD 2021 Anaemia Collaborators. Prevalence, years lived with disability, and trends in anaemia burden by severity and cause, 1990-2021: findings from the global burden of disease study 2021. Lancet Haematol. (2023) 10:e713–34. doi: 10.1016/s2352-3026(23)00160-6

PubMed Abstract | Crossref Full Text | Google Scholar

24. Stevens, GA, Alkema, L, Black, RE, Boerma, JT, Collins, GS, Ezzati, M, et al. Guidelines for accurate and transparent health estimates reporting: the GATHER statement. Lancet. (2016) 388:e19–23. doi: 10.1016/s0140-6736(16)30388-9

PubMed Abstract | Crossref Full Text | Google Scholar

25. Hu, W, Fang, L, Ni, R, Zhang, H, and Pan, G. Changing trends in the disease burden of non-melanoma skin cancer globally from 1990 to 2019 and its predicted level in 25 years. BMC Cancer. (2022) 22:836. doi: 10.1186/s12885-022-09940-3

PubMed Abstract | Crossref Full Text | Google Scholar

26. Zhang, W, Zeng, W, Jiang, A, He, Z, Shen, X, Dong, X, et al. Global, regional and national incidence, mortality and disability-adjusted life-years of skin cancers and trend analysis from 1990 to 2019: an analysis of the global burden of disease study 2019. Cancer Med. (2021) 10:4905–22. doi: 10.1002/cam4.4046

PubMed Abstract | Crossref Full Text | Google Scholar

27. He, Q, Mok, TN, Sin, TH, Yin, J, Li, S, Yin, Y, et al. Global, regional, and National Prevalence of gout from 1990 to 2019: age-period-cohort analysis with future burden prediction. JMIR Public Health Surveill. (2023) 9:e45943. doi: 10.2196/45943

PubMed Abstract | Crossref Full Text | Google Scholar

28. Huang, L, Zhang, D, and Liu, M. Global trends in refractive disorders from 1990 to 2021: insights from the global burden of disease study and predictive modeling. Front Public Health. (2025) 13:1449607. doi: 10.3389/fpubh.2025.1449607

PubMed Abstract | Crossref Full Text | Google Scholar

29. Liang, W, Hu, A, Hu, P, Zhu, J, and Wang, Y. Estimating the tuberculosis incidence using a SARIMAX-NNARX hybrid model by integrating meteorological factors in Qinghai Province. China Int J Biometeorol. (2023) 67:55–65. doi: 10.1007/s00484-022-02385-0

PubMed Abstract | Crossref Full Text | Google Scholar

30. Arinaminpathy, N, Reed, C, Biggerstaff, M, Nguyen, AT, Athni, TS, Arnold, BF, et al. Estimating community-wide indirect effects of influenza vaccination: triangulation using mathematical models and bias analysis. Am J Epidemiol. (2024) kwae365. doi: 10.1093/aje/kwae365

PubMed Abstract | Crossref Full Text | Google Scholar

31. Sedik, S, Egger, M, and Hoenigl, M. Climate change and medical mycology. Infect Dis Clin N Am. (2025) 39:1–22. doi: 10.1016/j.idc.2024.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

32. Wang, Y, Tang, S, Ma, T, Zhang, L, Xiao, Y, and Chen, S. Burden of fungal skin diseases in five east Asian countries: an analysis of the global burden of disease study 2021. Mycoses. (2025) 68:e70016. doi: 10.1111/myc.70016

PubMed Abstract | Crossref Full Text | Google Scholar

33. Badiane, AS, Ramarozatovo, LS, Doumbo, SN, Dorkenoo, AM, Mandengue, C, Dunaisk, CM, et al. Diagnostic capacity for cutaneous fungal diseases in the African continent. Int J Dermatol. (2023) 62:1131–41. doi: 10.1111/ijd.16751

PubMed Abstract | Crossref Full Text | Google Scholar

34. Tischner, Z, Páldy, A, Kocsubé, S, Kredics, L, Dobolyi, C, Sebők, R, et al. Survival and growth of microscopic fungi derived from tropical regions under future heat waves in the Pannonian biogeographical region. Fungal Biol. (2022) 126:511–20. doi: 10.1016/j.funbio.2022.04.005

PubMed Abstract | Crossref Full Text | Google Scholar

35. Jabrodini, A, Bahrololuom, SM, Amirinia, F, Kharazi, M, Mohamadi, T, Mehrabanpour, R, et al. Prevalence and comparison of superficial-cutaneous fungal infections before and during the COVID-19 pandemic in shiraz, Iran: a retrospective study (2018-2021). Indian J Microbiol. (2024) 64:1570–6. doi: 10.1007/s12088-023-01174-6

PubMed Abstract | Crossref Full Text | Google Scholar

36. Yan, D, Li, M, Si, W, Ni, S, Liu, X, Chang, Y, et al. Haze exposure changes the skin fungal community and promotes the growth of Talaromyces strains. Microbiol Spectr. (2023) 11:e0118822. doi: 10.1128/spectrum.01188-22

PubMed Abstract | Crossref Full Text | Google Scholar

37. Rao, HH, and McClelland, EE. A new overview of sex Bias in fungal infections. J Fungi. (2024) 10:607. doi: 10.3390/jof10090607

PubMed Abstract | Crossref Full Text | Google Scholar

38. Egger, M, Hoenigl, M, Thompson, GR, Carvalho, A, and Jenks, JD. Let's talk about sex characteristics-as a risk factor for invasive fungal diseases. Mycoses. (2022) 65:599–612. doi: 10.1111/myc.13449

PubMed Abstract | Crossref Full Text | Google Scholar

39. Pinzan, CF, Valero, C, de Castro, PA, da Silva, JL, Earle, K, Liu, H, et al. Aspergillus fumigatus conidial surface-associated proteome reveals factors for fungal evasion and host immunity modulation. Nat Microbiol. (2024) 9:2710–26. doi: 10.1038/s41564-024-01782-y

PubMed Abstract | Crossref Full Text | Google Scholar

40. Chen, XQ, and Yu, J. Global demographic characteristics and pathogen Spectrum of tinea Capitis. Mycopathologia. (2023) 188:433–47. doi: 10.1007/s11046-023-00710-8

PubMed Abstract | Crossref Full Text | Google Scholar

41. Fu, C, Peng, Y, Yang, F, He, Z, Ali, H, and Xu, D. Potentially suitable geographical area for Colletotrichum acutatum under current and future climatic scenarios based on optimized MaxEnt model. Front Microbiol. (2024) 15:1463070. doi: 10.3389/fmicb.2024.1463070

PubMed Abstract | Crossref Full Text | Google Scholar

42. Mohammad, KA, Ismail, HM, Shekhany, KAM, Yashooa, RK, Younus, DA, Abdullah, SK, et al. Fungal disease incidence and prevalence in Iraq - preliminary estimates. J Mycol Med. (2024) 34:101516. doi: 10.1016/j.mycmed.2024.101516

PubMed Abstract | Crossref Full Text | Google Scholar

43. Seudjip, NLJ, Libeko, SK, Ntshila, LK, Bunga, PM, Mvumbi, GL, Kabututu, PZ, et al. Retrospective epidemiology of dermatomycosis in Kinshasa, Democratic Republic of Congo, from 2000 to 2023. Mycoses. (2024) 67:e70010. doi: 10.1111/myc.70010

PubMed Abstract | Crossref Full Text | Google Scholar

44. Xihua, Y, Xiaofeng, Q, Lulu, T, Zhiyu, K, Min, D, and Yi, W. Fungal skin disease incidence, prevalence and disability-adjusted life years in four Asian countries (1990-2019). Mycoses. (2025) 68:e70004. doi: 10.1111/myc.70004

PubMed Abstract | Crossref Full Text | Google Scholar

45. Global variations in pubertal growth spurts in adolescents living with perinatal HIV. AIDS. (2023) 37:1603–15. doi: 10.1097/qad.0000000000003602

PubMed Abstract | Crossref Full Text | Google Scholar

46. Fulgence, KK, Marie, KPC, Akoua, VB, Massafoma, KEG, Etienne, AK, Abibatou, KT, et al. Dermatophytosis and the associated risk factors among primary school children in southern and Central Côte d'Ivoire. Mycoses. (2023) 66:869–75. doi: 10.1111/myc.13624

PubMed Abstract | Crossref Full Text | Google Scholar

47. Sulaiman, C, and Abdul-Rahim, AS. Relationship between wood fuel energy consumption and forest degradation at regional and sub-regional levels of sub-Saharan Africa: the role of control of corruption and government effectiveness. Environ Sci Pollut Res Int. (2022) 29:74512–25. doi: 10.1007/s11356-022-21108-x

PubMed Abstract | Crossref Full Text | Google Scholar

48. Haro, M, Alemayehu, T, and Mikiru, A. Dermatophytosis and its risk factors among children visiting dermatology clinic in Hawassa Sidama, Ethiopia. Sci Rep. (2023) 13:8630. doi: 10.1038/s41598-023-35837-7

PubMed Abstract | Crossref Full Text | Google Scholar

49. Madarasingha, NP, Thabrew, H, Uhrlass, S, Eriyagama, S, Reinal, D, Jayasekera, PI, et al. Dermatophytosis caused by Trichophyton indotineae (Trichophyton mentagrophytes ITS genotype VIII) in Sri Lanka. Am J Trop Med Hyg. (2024) 111:575–7. doi: 10.4269/ajtmh.24-0027

PubMed Abstract | Crossref Full Text | Google Scholar

50. Xin, Y, Quan, L, Zhang, H, and Ao, Q. Emerging polymer-based Nanosystem strategies in the delivery of antifungal drugs. Pharmaceutics. (2023) 15:866. doi: 10.3390/pharmaceutics15071866

PubMed Abstract | Crossref Full Text | Google Scholar

51. Wu, Y, Jiang, W, Cong, Z, Chen, K, She, Y, Zhong, C, et al. An effective strategy to develop potent and selective antifungal agents from cell penetrating peptides in tackling drug-resistant invasive fungal infections. J Med Chem. (2022) 65:7296–311. doi: 10.1021/acs.jmedchem.2c00274

PubMed Abstract | Crossref Full Text | Google Scholar

52. Yang, JX, Ma, GN, Li, YT, Shi, YP, and Liang, GW. Resistance and virulence genes characteristic of a South Asia clade (I) Candida auris strain isolated from blood in Beijing. Clinics. (2024) 79:100497. doi: 10.1016/j.clinsp.2024.100497

PubMed Abstract | Crossref Full Text | Google Scholar

Glossary

GBD - Global Burden of Diseases

DALYs - Disability-Adjusted Life Years

ASR - Age-Standardized Rate

ASIR - Age-Standardized Incidence Rate. This is the number of new cases of a disease in a standardized population, adjusted for age differences. It allows for comparisons of how often a disease occurs between populations with different age structures, usually expressed per 100,000 people

ASPR - Age-Standardized Prevalence Rate. This shows the proportion of people affected by a disease at a specific point in time, adjusted for age. It includes both new and existing cases, and is also measured per 100,000 people, which helps compare how widespread a disease is across populations of varying age structures

ASDR - Age-Standardized DALY Rate. This represents the rate of Disability-Adjusted Life Years (DALYs), which combines the years lost due to early death and the years lived with disability from a disease. It is age-adjusted and expressed per 100,000 people, giving an estimate of the overall burden of the disease in different populations

SDI - Socio-Demographic Index

GATHER - the Guidelines for Accurate and Transparent Health Estimates Reporting

APC - Age-Period-Cohort

ARIMA - Autoregressive Integrated Moving Average. An autoregressive integrated moving average model is a form of regression analysis that gauges the strength of one dependent variable relative to other changing variables. The model’s goal is to predict future securities or financial market moves by examining the differences between values in the series instead of through actual values.

EAPC - Estimated Annual Percentage Change

CI - Confidence Interval

UI - Uncertainty Interval

Keywords: fungal skin diseases, disability-adjusted life-years, Global Burden of Disease, incidence, prevalence, trend, age-period-cohort

Citation: Li D, Fan S, Zhao H, Song J, Guo L, Li W, Xu X and Li Q (2025) Worldwide trends and future projections of fungal skin disease burden: a comprehensive analysis from the Global Burden of Diseases study 2021. Front. Public Health. 13:1580221. doi: 10.3389/fpubh.2025.1580221

Received: 20 February 2025; Accepted: 12 May 2025;
Published: 04 June 2025.

Edited by:

Gloria M. González-González, Autonomous University of Nuevo León, Mexico

Reviewed by:

Saumya Panda, Belle Vue Clinic, India
Yi Xue, National Clinical Research Center for Child Health, China
Olusola Ayanlowo, University of Lagos, Nigeria
María Guadalupe Frías De León, Hospital Regional de Alta Especialidad de Ixtapaluca, Mexico

Copyright © 2025 Li, Fan, Zhao, Song, Guo, Li, Xu and Li. 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: Xuewen Xu, eHd4X3djaHN1QDE2My5jb20=; Qingfeng Li, ZHIubGlxaW5nZmVuZ0B5YWhvby5jb20=

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.