Prevalence burden of zIDPs
Globally, there were an estimated 151 376 744 schistosomiasis cases (95% UI 109 062 891 to 198 666 395), 633 404 CE cases (95% UI 517 477 to 782 468), 4 357 190 cysticercosis cases (95% UI 3 150 315 to 5 716 430), and 44 466 329 (95% UI 40 017 218 to 50 034 921) foodborne trematodiases cases, respectively. The age-standardized prevalence rate was 1 914.3 (95% UI 1378.9 to 2510.8) per 100 000 population for schistosomiasis, 7.7 (95% UI 6.3 to 9.5) for CE, 51.3 (95% UI 37.2 to 67.3) for cysticercosis, and 526.7 (95% UI 473.7 to 593.3) for foodborne trematodiases. From 1990 to 2019, the age-standardized prevalence rates of schistosomiasis, cysticercosis, and foodborne trematodiases showed significant global declines, averaging -22.6% (95% CI -0.33 to 0.13), -25.1% (95% CI -0.29 to 0.20), and -50.1% (95% CI -0.53 to 0.47) per year, respectively. Conversely, the prevalence of cystic echinococcosis remained stable over the same period, with a global average change of 4.4% (95% CI -0.03 to 0.10) per year. In addition, the number of cases and the age-standardized prevalence rate of schistosomiasis were markedly higher than those of CE, cysticercosis and foodborne trematodiases (Table 1).
In GBD 2021, the Western Sub-Saharan Africa reported the highest schistosomiasis cases (61 534 853, 95% UI 44 035 020 to 81 338 005), followed by that of Eastern Sub-Saharan Africa (46 949 044, 95% UI 33 903 446 to 61 048 520) and East Asia (11 459 581, 95% UI 8 953 886 to 15 054 894, Fig 1A). Western Sub-Saharan Africa also had the highest age-standardized prevalence rate of schistosomiasis at 13 433.6 per 100 000 population (95% UI 9 651.8 to 17 812.1), followed by Eastern Sub-Saharan Africa (11 617.6, 95% UI 8 558.00 to 15 093.1) and Central Sub-Saharan Africa (7 947.5, 95% UI 5 702.6 to 11 151.9). Furthermore, the age-standardized prevalence rates of schistosomiasis decreased across all GBD regions from 1990 to 2021, except in the Caribbean. The region with the highest average decrease rate was North Africa and the Middle East (-0.83, 95% CI -0.88 to 0.78, Table 1).
For CE, the highest number of cases was observed in South Asia (199 794, 95% UI 161 370 to 248 813), East Asia (114 589, 95% UI 86 308 to 147 199) and North Africa and Middle East (94 754, 95% UI 77 603 to 117 793, Fig 1B). The regions with the highest age-standardized prevalence rates were Southern Latin America (38.2 per 100 000 population, 95% UI 32.0 to 45.2), Central Asia (23.0, 95% UI 18.7 to 28.2) and Eastern Sub-Saharan Africa (17.2, 95% UI 13.8 to 21.7, Table 1). Additionally, six out of 21 GBD regions showed increased trends in age-standardized prevalence rate from 1990 to 2021, with East Asia experiencing the highest average annual increase of 153.6% (95% CI 1.3 to 1.8, Table 1).
In 2021, the regions with the highest number of cysticercosis cases were South AsiaS (755 683, 95% UI 551 993 to 972 079), East Asia (725 883, 95% UI 494 719 to 1 018 980) and Central Latin America (562 969, 95% UI 410 577 to 708 799) (Table 1, Fig 1C). Central Latin America (216.4 per 100 000 population, 95% UI 160.3 to 271.4), Southern Sub-Saharan Africa (208.3, 95% UI 155.3 to 258.4) and Central Sub-Saharan Africa (191.4, 95% UI 144.1 to 235.0) had the highest age-standardized prevalence rate. From 1990 to 2021, the age-standardized prevalence rate of cysticercosis decreased in all GBD regions except High-income North America, with the largest decrease observed in North Africa and Middle East (average annual change -51.9%, 95% CI -0.561 to 0.478, Table 1).
In 2021, the regions with the highest number of foodborne trematodiases cases were East Asia (33 317 223, 95% UI 29 251 039 to 38 353 602), Southeast Asia (4 708 038, 95% UI 4 280 246 to 5 159 478) and High-income Asia Pacific (1 974 178, 95% UI 1 750 147 to 2 202 860, Table 1). The regions with the highest age-standardized prevalence rates were Andean Latin America (2 606.1, 95% UI 2 127.3-3 171), East Asia (1 862.4, 95% UI 1 640.7-2 161.6), and High-income Asia Pacific (821.5, 95% UI 731.3 to 915.2, Table 1). From 1990 to 2021, an increase in the age-standardized prevalence rate of foodborne trematodiases was observed in High-income Asia Pacific (average change 11.0% per year, 95% CI 0·01 to 0·22) and Eastern Europe (0.4%, 95% CI -0.02 to 0.03).
DALYs burden of zIDPs
In 2021, the global number of DALYs (Disability-Adjusted Life Years) due to schistosomiasis was 1 746 333 (95% UI 1 038 122 to 2,984,204), with an age-standardized rate of 21.9 DALYs (95% UI 12.9 to 37.3) per 100 000 population. The age-standardized DALY rate for schistosomiasis significantly decreased from 1990 to 2021, with an average annual change of -0.41 (95% CI -0.52 to -0.32) globally (Table 2). Similarly, in 2021, the number of DALYs due to cystic echinococcosis (CE) was 105 072 (95% UI 78 967 to 133 309), with an age-standardized rate of 1.3 DALYs (95% UI 1.0 to 1.7) per 100 000 population, showing an average annual decrease of -0.72 (95% CI -0.8 to -0.6) from 1990 to 2021. For cysticercosis, the number of DALYs in 2021 was 1 235 939 (95% UI 787 769 to 1 808 363), with an age-standardized rate of 14.6 (95% UI 9.3 to 21.3), and an average annual decline of -0.39 (95% CI -0.47 to -0.29) from 1990 to 2021. Regarding foodborne trematodiases, the number of DALYs in 2021 was 998 029 (95% UI 569 766 to 1 638 112), with an age-standardized rate of 11.8 (95% UI 6.7 to 19.5). The age-standardized DALY rate for foodborne trematodiases also declined, with an average annual decrease of -0.54 (95% CI -0.60 to -0.45) from 1990 to 2021 (Table 2). Overall, the DALY burden (both the number of DALYs and the age-standardized DALY rate) of schistosomiasis was higher than that of cystic echinococcosis, cysticercosis, and foodborne trematodiases.
In 2021, regions of sub-Saharan Africa (central, eastern, and western) had the highest number of DALYs and age-standardized DALY rates (per 100,000 population) for schistosomiasis (Fig. 2A). All schistosomiasis-endemic GBD regions exhibited a decline in the age-standardized DALY rate, with the largest decreases occurring in North Africa and the Middle East (average annual change -0.82, 95% CI -0.85 to -0.78).
For CE, the highest number of DALYs was observed in South Asia (34 006, 95% UI 25 804 to 43 190), Eastern Sub-Saharan Africa (17 337, 95% UI 12 042 to 23 167), and North Africa and the Middle East (17 319, 95% UI 13 165 to 21 974, Fig 2B). The highest age-standardized DALY rates were found in Eastern Sub-Saharan Africa (4.2, 95% UI 3.0 to 5.4), North Africa and the Middle East (3.0, 95% UI 2.3 to 3.7), and Southern Latin America (3.0, 95% UI 2.0 to 4.2). All GBD regions showed obvious decreases in the age-standardized DALY rate, with the highest average annual change in Western Sub-Saharan Africa (-0.81, 95% CI -0.89 to -0.69, Table 2).
In 2021, the highest number of DALYs due to cysticercosis occurred in East Asia (176 228, 95% UI 101 465 to 284 693), Central Latin America (157 382, 95% UI 98 523 to 228 880), and Eastern Sub-Saharan Africa (94 863, 95% UI 64 603 to 129 428, Fig 2C). The highest age-standardized DALY rates were observed in Central Sub-Saharan Africa (72.0, 95% UI 46.4 to 103.7), Southern Sub-Saharan Africa (66.8, 95% UI 42.7 to 94.3), and Central Latin America (60.3, 95% UI 37.9 to 87.3, Table 2). Most regions (except High-income North America) showed decreases in the age-standardized DALY rate for cysticercosis, with the largest decreases in North Africa and the Middle East (-0.65, 95% CI -0.74 to -0.52), East Asia (-0.54, 95% CI -0.67 to -0.40), and South Asia (-0.55, 95% CI -0.63 to -0.43, Table 2).
For foodborne trematodiases, East Asia had the highest number of DALYs (768 297, 95% UI 383 883 to 1 367 826), significantly higher than other regions. Andean Latin America had the highest age-standardized DALY rate (56.1, 95% UI 19.0 to 106.1), followed by East Asia (42.6, 95% UI 21.3 to 77.0, Table 2). The increase in the age-standardized DALY rate was only found in High-income Asia Pacific (average change 0.170 per year, 95% CI 0.042 to 0.319) and Eastern Europe (0.005, 95% CI -0.040 to 0.048), while other regions showed a decreasing trend (Table 2).
The burden of zIDPs by age
For schistosomiasis, the 10-29 age groups had the highest number of cases. The age-standardized prevalence rate showed an upward trend until the 15-19 age group, after which it presented a downward trend with increasing age (Fig 3A). Additionally, the number of male schistosomiasis cases exceeded that of females in all age groups. The maximum number of DALYs for schistosomiasis was found in the 15-34 age groups. The age-standardized DALY rate exhibited an upward trend until the 20-24 age group, and then a downward trend with increasing age (Fig. 3E).
For CE, the 30-49 age groups had the highest number of cases, with the peak of the age-standardized prevalence rate also occurring in these age groups. Interestingly, the number of male cases was less than that of females across all age groups (Fig 3B). The <5 years age group had the highest number of DALYs (both females and males each had more than 5,000). The age-standardized DALY rate did not change much across age groups but showed a significant decrease from the <5 years age group to the 15-19 age group (Fig 3F).
For cysticercosis, the maximum number of cases were concentrated in the 55-74 age groups, while the age-standardized prevalence rate increased with age. Additionally, male cases of cysticercosis were fewer than female cases across all age groups (Fig 3C). The number of DALYs was concentrated in the 20-84 age groups, with the age-standardized DALY rate showing a continuous upward trend with age. Furthermore, the number of DALYs in females was higher than in males across all age groups (Fig 3G).
For foodborne trematodiases, the 45-59 age groups had the highest number of cases and DALYs. The age-standardized prevalence rate and DALY rate both initially rose (in the age groups <5 to 69) and then declined (in those over 69 years old). Additionally, the number of cases and DALYs in males across almost all age groups (except those over 80) were higher than in females (Fig 3D, H).
The burden of zIDPs by SDI
In the five categories of SDI regions (high, middle-high, middle, low-middle, low), the age-standardized prevalence rate of schistosomiasis showed an upward trend in the high SDI region and a downward trend in the other regions from 1990 to 2021. The age-standardized prevalence rate of foodborne trematodiases exhibited a downward trend in the middle SDI region, but rebounded in 2015 in the high and middle-high SDI regions. The age-standardized prevalence rates of CE and cysticercosis remained stable from 1990 to 2021 (Fig 4A).
The age-standardized DALY rate of CE showed no significant change in the five SDI regions from 1990 to 2021. The age-standardized DALY rate of schistosomiasis showed a downward trend in the low, low-middle, and middle SDI regions, and remained stable in the high-middle and high SDI regions. The age-standardized DALY rate of cysticercosis showed a slight downward trend in the low-middle, middle, and high-middle regions. The age-standardized DALY rate of foodborne trematodiases showed a significant downward trend only in the middle SDI region (Fig 4B).
The relation between burden of zIDPs and SDI
The age-standardized prevalence rates of the zIDPs exhibited three distinct patterns with increasing SDI at the national level. Schistosomiasis showed a linear decline in prevalence. The prevalence of cysticercosis was divided into two zones: increasing (<0.6) and decreasing (>0.6), with an SDI threshold of 0.6. The prevalence of CE and foodborne trematodiases showed no significant changes. The highest age-standardized prevalence rates were observed in Mauritius for schistosomiasis (30.8 per 100 000 population, 95% UI: 21.6 to 41.2), in the Republic of Moldova for CE (107.4, 95% UI: 86.2 to 131.6), in Mexico for cysticercosis (281.4, 95% UI: 206.3 to 352.0), and in the Lao People's Democratic Republic for foodborne trematodiases (4 028.8, 95% UI: 4 585.0 to 494.3, Fig 5A, C, E, G).
At the regional level, Western, Central, and Southern Sub-Saharan Africa had higher-than-expected age-standardized prevalence rates of schistosomiasis despite an overall decline. The age-standardized prevalence rates of CE showed a curve that first decreases, then rises, and then decreases again, with higher-than-expected rates in Eastern Sub-Saharan Africa, Central Asia, and Southern Latin America. For cysticercosis, the prevalence decreased with increasing SDI, and Central and Southern Sub-Saharan Africa, Andean and Central Latin America, and High-income North America had higher-than-expected rates. The prevalence rates of foodborne trematodiases showed an initial rise and then a decrease at an SDI of 0.5, with higher-than-expected rates in Latin America, East Asia, Southeast Asia, Eastern Europe, and High-income Asia Pacific (Fig 5B, D, G, H).
At the national level and regional level, the age-standardized DALY of the zIDPs decreased with increases in SDI except foodborne trematodiases (Fig 6). The African countries have higher age-standardized DALYs rates compared to the countries in other regions, and the Central African Republic had the highest rate for schistosomiasis (401.6 per 100 000 population, 95% UI: 289.8 to 540.6), South Sudan for CE (12.3, 95% UI: 8.4 to 17.1), and Zimbabwe for cysticercosis (98.7, 95% UI: 60.2 to 146.5). The Lao People's Democratic Republic had the highest age-standardized DALYs rate of foodborne trematodiases (Fig 6 A, C, E, G).
Sub-Saharan Africa, North Africa and Middle East and Tropical Latin America had higher-than-expected age-standardized DALY rates of scistosomiasis based on the SDI between 1990 and 2021. South Asia, Central Asia, North Africa and Middle East, Southern Latin America and Eastern Europe had higher-than-expected age-standardized DALY rates of CE. For cysticercosis, Central and Southern sub-Saharan Africa, Andean and Central Latin America and High-income North America had higher-than-expected age-standardized DALYs rates. For foodborne trematodiases, Latin America, East Asia, Southeast Asia, Eastern Europe, High−income Asia Pacific had higher-than-expected age-standardized DALY rates (Fig 6B, D, G, H).