Global Burden of Sand Flea Disease Estimated at 575,000 Disability Years

A comprehensive prevalence-based epidemiological modelling investigation has delivered the first quantitative global burden estimate for tungiasis, calculating that this neglected tropical ectoparasitic skin infection causes up to 379,645 base disability-adjusted life-years annually across endemic countries, rising to 575,298 disability-adjusted life-years when comorbid conditions are incorporated. Published in BMC Global and Public Health, the study focuses strictly on years lived with disability across endemic regions in Sub-Saharan Africa and Latin America, establishing that the disease exerts an annual loss of 24.2 base disability-adjusted life-years per 100,000 population, or 36.4 per 100,000 when comorbidities are accounted for. Researchers led by Alistair Antonopoulos evaluated transmission dynamics and health decrements across 52 simulated weekly time points, establishing that the health loss generated by sand flea infestation in endemic zones is broadly comparable in absolute magnitude to major neglected conditions such as soil-transmitted helminth infections and echinococcosis, while exceeding the estimated global burdens attributable to toxocariasis and foodborne fascioliasis.

Tungiasis is a zoonotic parasitic dermatosis triggered when adult female sand fleas belonging to the genus Tunga, primarily Tunga penetrans and Tunga trimamillata, burrow into the epidermis of human and animal hosts. Infestation occurs predominantly on the feet when bare skin contacts infested soil or domestic flooring. Once embedded, the female flea maintains continuous contact with the external environment through a microscopic opening measuring between 250 and 500 micrometers on its terminal abdominal segment, enabling respiration, excretion, copulation, and environmental egg shedding over its lifespan. As internal eggs mature within the flea abdomen, the parasite undergoes hypertrophy, expanding up to 2000-fold in volume within approximately one week. This rapid enlargement generates severe mechanical compression against adjacent cutaneous tissues, provoking intense local inflammation, violent itching, and debilitating pain that disrupts sleep patterns and daily locomotion. The inflammatory process is further exacerbated by endosymbiotic Wolbachia bacteria that release lipopolysaccharides into host tissue upon parasite death, often producing secondary ulcerations, suppuration, lymphangitis, gangrene, and severe permanent tissue disfigurement.

The health burden documented by the research team is concentrated overwhelmingly among socioeconomically disadvantaged communities lacking adequate housing, concrete flooring, and municipal sanitation infrastructure. Epidemiological niche models indicate that approximately 1.1 billion individuals reside within zones exhibiting climatic and ecological conditions conducive to transmission across Sub-Saharan Africa and Latin America. In Sub-Saharan Africa alone, an estimated 668 million individuals live in transmission-prone regions, while approximately 450 million reside in at-risk environments across Latin America. Within these broad geographical envelopes, human exposure is sharply magnified in rural agricultural villages and overcrowded urban informal settlements. Reported community infection rates vary dramatically based on study design and local environmental factors; while Kenya recorded a national population prevalence of 1.35 percent, pooled meta-analyses of focal studies across targeted vulnerable African communities have reported average prevalence rates of 33.4 percent, with school-aged children demonstrating pooled prevalence rates reaching 37.9 percent and certain hyperendemic communities in Uganda exhibiting prevalence rates as high as 84.6 percent.

To convert heterogeneous epidemiological data into standardized disease burden estimates, the authors formulated a specialized prevalence-based disability-adjusted life-year framework focused exclusively on non-fatal morbidity. Mortality was intentionally excluded from the mathematical framework because fatal secondary complications such as severe tetanus, sepsis, and viral infections cannot be reliably parameterized or directly attributed from existing surveillance records. The analytical duration parameter was fixed at one week, with individual health decrements tracked over a 52-week annual simulation cycle to capture dynamic temporal fluctuations in symptom intensity as parasites mature, get manually extracted, and reinfest hosts. Weekly disease burden calculations multiplied the calculated prevalent cases by assigned disability weights and time duration, utilizing the standard Global Burden of Disease multiplicative rule to combine primary physical disfigurement with secondary functional comorbidities. The underlying stochastic simulation completed 1000 Monte Carlo iterations to establish robust mean values and 95 percent uncertainty intervals for national, regional, and global disease burden outputs.

Generating dependable background prevalence ranges required the implementation of a Bayesian hierarchical Beta-Binomial framework capable of handling extensive data gaps and site-level overdispersion. Implemented in Stan through the brms package in R using the No-U-Turn Sampler, the hierarchical structure evaluated observed case numbers as Binomial draws governed by a precision parameter to capture excess variance between localized primary studies. The model simultaneously structured two analytical levels: a regional macro-tier encompassing West and Central Africa, East and Southern Africa, South America, and Central America with the Caribbean, alongside a site-specific tier accounting for variations in study design and community setting. To correct for the inherent bias of focal surveys conducted in heavily affected communities, the authors derived a setting correction factor using national survey data from Kenya, translating localized prevalence measurements into conservative national-level equivalents across data-sparse jurisdictions through four independent sampling chains of 4000 iterations that confirmed statistical convergence.

Because the Global Burden of Disease study lacks formal, disease-specific disability weights for sand flea infection, the research team established proxy health state categories matching validated clinical manifestations from related neglected conditions. Mild infection, defined by fewer than five embedded parasites without functional impairment, was assigned the disability weight for mild onchocerciasis disfigurement at 0.011. Moderate to severe uncomplicated infection, defined by 6 to 30 embedded fleas impairing daily function, received the scabies disfigurement weight of 0.027. Severe infection, characterized by more than 30 embedded parasites causing marked agony and functional restriction, received the severe onchocerciasis skin disease weight of 0.188. To capture the full spectrum of systemic morbidity, the researchers added specialized comorbidity weights to severe cases, incorporating post-dengue chronic fatigue syndrome at 0.219 for persistent insomnia and exhaustion, Guinea worm disease impairment at 0.079 for compromised lower-limb mobility, and mild cognitive impairment at 0.031 for children aged 5 to 15 experiencing educational disruption and concentration deficits.

The demographic structure of tungiasis transmission was modelled by partitioning populations into four standardized age categories based on empirical household surveys: children aged 0 to 4 representing 26 percent of total cases, children aged 5 to 14 accounting for 34 percent, adults aged 15 to 59 comprising 17 percent, and older individuals aged 60 and above representing 23 percent. High-risk age categories were assigned probabilistic ranges for developing moderate or severe disease, utilizing uniform distributions between 0.33 and 0.55 for individuals under 15, and between 0.36 and 0.60 for adults aged 60 and older, while working-age adults received a fixed moderate-to-severe proportion of 0.30. Within the pool of individuals experiencing moderate-to-severe disease, the authors designated 20 percent as severe cases harboring more than 30 embedded fleas, aligning with observed field survey frequencies of 8 to 15.1 percent in hyperendemic Brazilian and Kenyan communities. Furthermore, seasonal variations were incorporated via a uniform sigmoidal multiplier reflecting peak dry-season transmission at month three and reduced wet-season rates at month nine.

The resulting geographical distribution reveals an intense spatial concentration of tungiasis burden across Sub-Saharan Africa, which accounts for nine of the ten most heavily burdened nations globally. Nigeria experienced the highest absolute national toll, with an estimated 64,879 base disability-adjusted life-years, rising to 98,348 when comorbidities were incorporated into the analysis. Ethiopia followed with 41,006 base and 62,116 comorbidity-adjusted life-years, while the Democratic Republic of the Congo registered 37,505 base and 56,835 combined life-years. Tanzania and Uganda sustained 20,565 and 14,762 base disability years respectively, reaching 31,166 and 22,367 years under full comorbidity models. Brazil represented the sole non-African nation appearing within the top ten absolute burden rankings, generating 14,171 base disability years and 21,477 years with comorbidities. Across Central America and the Caribbean, absolute burden was led by Guatemala with 4877 base and 7391 combined life-years, followed by Haiti with 3289 base and 4985 combined life-years.

When national figures were adjusted for total underlying population size, the rankings shifted substantially, highlighting profound per capita suffering across rural nations. Burkina Faso exhibited the highest population-adjusted rate globally at 37.9 base and 57.5 comorbidity-adjusted disability years per 100,000 individuals, followed closely by Chad at 37.8 base and 57.3 combined years. The Democratic Republic of the Congo recorded 35.5 base and 53.7 combined years, South Sudan recorded 34.3 base and 52.0 combined years, and the Central African Republic registered 34.2 base and 51.8 combined years per 100,000 residents. In East and Southern Africa, high per capita tolls were observed in Ethiopia at 31.9 base and 48.3 combined years, Burundi at 31.5 base and 47.7 combined years, and Malawi, Tanzania, and Rwanda each recording approximately 30.9 base and 46.8 combined years per 100,000 residents. Conversely, populous nations like Brazil dropped sharply in per capita rankings to 6.7 base and 10.2 combined years per 100,000 people due to broader urbanization and higher access to sanitation.

A global probabilistic sensitivity analysis utilizing standardized regression coefficients across 1000 simultaneous parameter draws revealed that model variance is overwhelmingly driven by the structural classification of disease severity. The moderate-to-severe case split accounted for the largest individual share of variance with a standardized regression coefficient of 0.88, followed by regional prevalence estimates with a coefficient of 0.24. All remaining model parameters, including individual disability weight bounds, age distribution structures, insomnia durations, and seasonal amplitudes, each contributed less than 15 percent of output variance. The overall linear regression model demonstrated exceptional explanatory power with a coefficient of determination of 0.976. Additionally, an isolated evaluation of seventeen countries situated within transmission niches but excluded due to negligible reported indigenous cases—including Argentina, Chile, Mexico, Uruguay, and North African territories—indicated that their inclusion would have added 32,103 base and 48,644 combined disability-adjusted life-years, representing an 8.5 percent expansion of the global total.

Comparative analyses contextualize the health toll of tungiasis against better-funded neglected tropical diseases. The total comorbidity-adjusted burden of 575,298 disability-adjusted life-years approaches that of Ascaris lumbricoides at 754,000 years, although remaining below the broader soil-transmitted helminth complex at 1.38 million years and schistosomiasis at 1.64 to 2.19 million years. However, when evaluated on a per capita basis across endemic populations, the rate of 36.4 disability-adjusted life-years per 100,000 people substantially exceeds recent global estimates for soil-transmitted helminths at 18.84 per 100,000 and schistosomiasis at 21.895 per 100,000. This stark per capita intensity reflects the acute physical trauma of ectoparasitic ulceration, intense nocturnal pruritus, and loss of physical mobility, contrasting with the slower, chronic systemic presentations characteristic of intestinal helminthiases. The authors noted potential diagnostic overlap in East Africa with podoconiosis, where between 10 and 20 percent of patients in certain surveyed communities exhibit dual affliction, requiring future field differentiation.

The investigators emphasized that their model constitutes an essential preliminary quantification intended to guide future surveillance rather than a definitive epidemiological census. Key operational limitations include the necessary reliance on a single nationwide survey anchor from Kenya, the absence of pathogen transmission models integrating environmental and domestic animal reservoirs, and the application of standardized seasonal multipliers that do not reflect localized multi-modal rainfall regimes. To refine future burden assessments, the authors recommend conducting rigorous national prevalence surveys across five sentinel endemic nations spanning diverse transmission contexts: Ethiopia, Nigeria, Brazil, Burundi, and Madagascar. Furthermore, the findings advocate formal expert elicitation within the Global Burden of Disease architecture to establish bespoke disability weights, alongside coordinated One Health intervention strategies linking the World Health Organization, World Organisation for Animal Health, and United Nations Environment Programme to deliver low-cost, definitive dermatological treatments to underserved populations.

Subject of Research: Global burden of disease estimation for human tungiasis using disability-adjusted life-years modelling

Article Title: Estimation of the global burden of tungiasis in endemic countries using a disability adjusted life-years modelling approach

Article References: Antonopoulos, A., Elson, L., Krücken, J., Thielecke, M., Mutebi, F., Yapi, R. B., Fillinger, U., Feldmeier, H., Nyangacha, R., Kuispond, N.-R. S., Laken, E. K., Carlos, R. S. A., Charlier, J., & von Samson-Himmelstjerna, G. (2026). Estimation of the global burden of tungiasis in endemic countries using a disability adjusted life-years modelling approach. BMC Global and Public Health, 4(1), Article 85. https://doi.org/10.1186/s44263-026-00318-2

Image Credits: AI Generated

DOI: 10.1186/s44263-026-00318-2

Keywords: tungiasis, neglected tropical diseases, disability adjusted life years, public health, epidemiology, sub saharan africa, tunga penetrans, global burden of disease, Estimation, global, burden, endemic