Abstract

Background. Urolithiasis is common worldwide and can predispose to urinary tract infections and renal failure. We aimed to explore the global, regional, and national burden of urolithiasis between 1990 and 2019, stratified by sex, age, and sociodemographic index (SDI). Methods. From 1990 to 2019, data on the number of incident cases of urolithiasis, associated deaths, and disability-adjusted life years (DALYs) were extracted from the 2019 Global Burden of Disease (GBD) study. The trends for the incidence rate, mortality, and DALYs were evaluated using estimated annual percentage changes (EAPCs). Results. The incidence of urolithiasis increased by 48.57%, from 77.78 million incident cases in 1990 to 115.55 million in 2019, while its age-standardized incidence rate (ASIR) decreased. The ASIR increased slightly in the low SDI regions (EAPC = 0.33; 95% confidence interval [CI]: 0.24–0.43), while ASIRs in other SDI regions decreased. The incidence of urolithiasis by age presented a unimodal distribution, with the peak observed in patients aged between 50 years and 70 years. Urolithiasis-related mortality and DALYs also increased over time. Yet, the age-standardized death rate (ASDR) decreased by 2.05% (95% CI, −2.25% to −1.85%) per year, and the annual age-standardized DALY rate decreased by 1.77% (95% CI, −1.92% to −1.63%). The mortality and DALYs increased with age. The incidence, mortality, and DALYs were greater in males than those in females. The burden of urolithiasis showed obvious differences in its regional distribution over the past three decades. Conclusion. From 1990 to 2019, ASIR, ASDR, and age-standardized DALY rate of urolithiasis have decreased. Yet, particularly significant differences exist in the geographic, age, and sex distribution. Thus, medical resources should be rationally allocated and adjusted according to the geographic and demographic distribution of urolithiasis.

1. Introduction

Urolithiasis is highly prevalent in urology and can present as an asymptomatic condition or as a painful, recurrent disease. Affected by a variety of factors, such as the climate and seasonal temperature variation, dietary habits, water quality, direct occupational exposure, inheritance and genetic constitution, latitude, and comorbidities, the incidence and prevalence of urolithiasis are characterized by a distinct geographical variation [13]. The overall prevalence was 7.54% in mainland China between 1990 and 2016 [4], 5% to 10% in Europe in the decade before 2011, [5], 8.8% in North America between 2007 and 2010 [6], and 5.7% in Iran in 2005 [7].

The incidence and prevalence of urolithiasis have been on the rise. Although rarely life-threatening, urolithiasis can often cause intense pain and adversely affect patients’ quality of life. Moreover, because of a large number of new and recurrent cases, the high rate of surgical intervention, and the advent of novel technology, global health care costs related to the management of stones are relatively high. In 2000, the treatment of urolithiasis did cost up to $5.3 billion per year in the United States alone [8]. Thus, further insight into the global urolithiasis burden is essential for the allocation of limited health resources and formulation of rational policies.

The global burden of disease (GBD) study examined the burden of hundreds of diseases and injuries in 195 countries and territories around the world and provided an opportunity to comprehensively assess different aspects of human health, including the distribution and development trends of urolithiasis [9]. For a better understanding of the trends in urolithiasis burden according to geography, age, sex, and social development index (SDI), we used data from the 2019 GBD to describe the global, regional, and national trends in the incident rate, deaths, and DALYs associated with urolithiasis from 1990 to 2019.

2. Methods

2.1. Study Data

Study data on the urolithiasis burden—the annual incidence rate, death, disability-adjusted life years (DALYs), and their age-standardized rates (ASR)—were extracted from the 2019 GBD study using the global health data exchange (GHDx) query tool (https://ghdx.healthdata.org/gbd-results-tool). We also collected data on sex, age, and SDI to investigate their influence on the urolithiasis disease burden. The SDI, which ranges from 0 to 1, is a comprehensive indicator of social and demographic development. According to the order of their SDI value, 192 countries and regions around the world are divided into low SDI, middle (low-middle, middle, high-middle) SDI, and high SDI countries. The epidemiological trends for urolithiasis were observed at the regional, national, and global levels.

2.2. Statistical Analysis

The trends for urolithiasis incidence and mortality rates were evaluated by calculating the annual age-standardized incidence rate (ASIR), the age-standardized death rate (ASDR), age-standardized DALYs rate, and their respective estimated annual percentage changes (EAPCs). Urolithiasis DALYs were computed as the sum of the years lived with disability and the years of life lost [10]. According to the age group construction of the standard population, the ASRs (per 100,000 population) were calculated using the following formula:(ai refers to the incidence of the ith age group. denotes the number of persons (or weight) in the same age subgroup i of the assigned reference standard population) [11].

EAPCs is a generally well-accepted method to describe ASR using a regression model, and it quantitatively calculates the average annual rate of change of ASR for a specified period [11]. The regression line is used to estimate the natural logarithm of the rates, i.e., y = α + βx + ɛ, where y = ln (ASR), and x = the calendar year. The EAPC calculation formula, 100 × (exp(β) − 1), and its 95% confidence intervals (CI) can also be calculated from the linear regression model [11, 12]. All statistical data were analyzed using R version 3.6.3 (The R Foundation for Statistical Computing, Vienna, Austria), and a two-sided was considered statistically significant.

3. Results

3.1. The Change in the Incidence of Urolithiasis

Globally, the incidence of urolithiasis was 1.16 × 108 (95% uncertainty interval (UI): 0.93–1.40 × 108) in 2019 and 0.78 × 108 (UI: 0.62–0.95 × 108) in 1990. From 1990 to 2019, its incidence increased by 48.57% (Table 1, additional file 1: figure S1A). However, the global ASIR showed a decreasing trend, with an average annual decrease of 0.83% from 1696.18/100,000 persons (95% UI, 1358.11–2078.11) in 1990 to 1394.03/100,000 persons (95% UI: 1126.4–1688.16) in 2019 (EAPC = −0.83; 95% CI: −0.92 to −0.74) (Table 1, additional file 1: figure S1B). The annual ratio of elder (aged >60 years) to younger patients remained relatively stable each year (additional file1: figure S2A). Moreover, ASIR in both male and female patients decreased similarly, and ASIR in male patients was more than that in female patients (figures 1(a) and 1(b), additional file 1: figure S1B). Between 1990 and 2019, the incidence of urolithiasis plotted against age showed a unimodal distribution, with the peak value observed in patients aged between 50 years and 70 years (figures 1(a) and 1(b), additional file1: figures S3S5). The male : female ratio tends to increase with age, increasing to approximately 3 in 2019 (Figure 2). However, urolithiasis tended to occur in the younger population (<60 years) (additional file1: figure S2A)

With respect to the SDI level, the ASIR in the low SDI quintile presented a slightly increasing trend, with an EAPC of 0.33 (95% CI: 0.24–0.43). On the contrary, ASIRs in other quintiles decreased, with the decrease in the high-middle SDI quintile being the most significant (Table 1, figures 3(a)–3(c)). The higher the SDI level, the more the elderly patients among all urolithiasis incidence patients (figure 4(a)). There were similar trends of male : female ratio in all SDI regions (Figure 2). Meanwhile, a significant positive correlation was detected between ASIR and SDI (R = 0.46, ) (additional file1: figure S6A). Furthermore, EAPC was negatively correlated with ASIR (R = −0.34, ), implying that urolithiasis increased more slowly in countries with high incidence than in countries with low incidence (Figure 5(a)).

As for the findings by the GBD regions and countries, the ASIR of urolithiasis increased in 143 countries and 11 regions (the top three countries: Jordan, Romania, and Germany), decreased in 32 countries and 9 regions (the bottom three countries: Poland, China, and Indonesia), and remained stable in 17 countries and 1 region (Albania, Costa Rica et al.). The top three countries with high ASIRs were the Russian Federation (4541.88 per 100,000 people), Ukraine (4282.60 per 100,000 people), and Latvia (4156.67 per 100,000 people). The bottom three countries were Burundi (525.01 per 100,000 people), South Sudan (533.43 per 100,000 people), and Madagascar (535.88 per 100,000 people). The top three countries with EAPC were Jordan (2.10), Romania (2.01), and Germany (2.00), while the bottom three countries were Poland (−3.87), China (−2.8), and Indonesia (−2.17). All the above results are shown in figure 6(a) and additional file 1: tables S1S3, figure S7A. In addition, the majority of countries had a unimodal age distribution (additional file 1: table S4).

3.2. The Change in the Mortality of Urolithiasis

Globally, mortality because of urolithiasis increased by 17.12% from 113.38 × 102 (95%UI: 72.78–137.77) in 1990 to 132.27 × 102 (95% UI: 106.16–162.67) in 2019 (Table 2, additional file 1: figure S1C). Over the past three decades, ASDR decreased significantly with an average annual decrease of 2.05% from 0.3/100,000 persons (95% UI, 0.20–0.37) in 1990 to 0.17/100,000 persons (95% UI: 0.14–0.21) in 2019 (EAPC = −2.05; 95% CI: −2.25 – −1.85) (Table 2). The death cases increased rapidly with age, especially for the elderly (figure 1(b); additional file 1: figures S8S10). ASDR in both male and female patients decreased significantly, and the ASDR of male patients was higher than that of female patients between 1990 and 2019 (Table 2, figure 1(b), additional file 1: figure S1C). The mortality ratio of male and female patients fluctuated with age (additional file 1: figure S11).

With respect to the SDI, ASDRs in the different SDI regions decreased (Table 2, figure 3(b)). EAPC was negatively associated with ASDR (R = −0.44, ). However, a correlation between SDI and EAPC or ASDR was not revealed (figure 5(b), additional file 1: figure S6B). Contrary to the change in ASIR, older patients (>60 years) died from urolithiasis (figure 4(b), additional file1: figure S1A). The number of annual young deaths decreased year by year, however, the number of elderly deaths increased gradually, irrespective of SDI regions, especially in patients >95 years of age (additional file 1: figures S1C & S8S10).

At the level of GBD regions and countries, 69 countries and 7 regions had an increased ASDR, 104 countries and 11 regions had a reduced ASDR, and 19 countries and 3 regions had a stable ASDR. The top three ASDRs were those of Armenia, Kazakhstan, and the Philippines, while the bottom three were those of North Macedonia, Montenegro, and Lebanon. Details are described in figure 6(b), additional file 1: tables S1, S2, & S5, figure S7B. The death cases increased rapidly with age irrespective of sex in the majority of countries (additional file 1: table S7).

3.3. The Change in DALYs of Urolithiasis

On a global level, there were 5,167,310 (95% UI, 3,741,330 to 6,357,170) DALYs in 1990 and 6,043,090 (95% UI, 4,773,540 to 7,451,940) DALYs in 2019, an increase of 16.95%. The age-standardized DALYs rate demonstrated a downward trend with an EAPC of −1.77 (95% CI, −1.92 to −2.21), declining from 11.75/100,000 persons (95% UI, 8.57–14.39) in 1990 to 7.35/100,000 persons (95% UI, 5.82–9.04) in 2019 (Table 3). The age-standardized DALYs rate of males was higher than that of female cases between 1990 and 2019 (Table 3). Investigating from the SDI standpoint, all four SDI regions witnessed a drop in the age-standardized DALYs rate, and DALYs were not associated with SDI (Table 3, additional file 1: figure S6C). Besides, there was no significant correlation between EAPC and SDI (R = −0.016, p = 0.85), and there was a negative association between EAPC and the age-standardized DALYs rate (figure 5(c), R = −0.52, ). The DALYs rate rose in both male and female patients with age (figure 1(c)). Compared with 1990, there was no obvious decrease in DALYs in elderly patients in 2019, and similar trends in the age distribution of mortality were seen in DALYs (additional file 1: figures S12S14). The DALY ratio of male and female patients also fluctuated with age (additional file 1: figure S15).

At the level of GBD regions and countries, 73 countries and 6 regions had an increased age-standardized DALYs rate, 104 countries and 11 regions had a decreased age-standardized DALYs rate, and 15 countries and 4 regions had a stable age-standardized DALYs rate. The top three age-standardized DALYs rates were seen in Armenia, Russian Federation, and the Philippines. The bottom three age-standardized DALYs rates were seen in Cabo Verde, Panama, and El Salvador. Details are displayed in figure 6(c), additional file 1: tables S1, S2, & S6, and additional file 1: figure S7C. Similarly, The DALYs increased rapidly with age, irrespective of sex in almost every country (additional file 1: table S8).

4. Discussion

Currently, urolithiasis remains a major global public health problem and warrants our attention. Several previous epidemiological studies on urolithiasis have mainly focused on individual regions or countries [5, 13, 14], however, global epidemiological data on urolithiasis are lacking. Based on the 2019 GBD study, we described the incidence, mortality, and DALY of urolithiasis at the global, regional, and national levels with their corresponding current trends and survival patterns from 1990 to 2019.

Our findings demonstrated that though the global incidence of urolithiasis increased by 48.57% in 2019 compared with its incidence in 1990, the ASIR actually decreased. There are several explanations for this trend. Zhu et al. reported that the increase in incidence might be attributed to population growth and the aging process of the population [15]. Advances in the detection of both symptomatic and asymptomatic stones by modern imaging methods and the widespread use of CT scans also contributed to the increased incidence of kidney stones in the study by Kittanamongkolchai et al. [16].

We found that urolithiasis incidence was associated with age distribution, which has not changed markedly over the last 30 years, and that the 40 to 60 years age group in both adult males and females was more likely to suffer from urolithiasis. This result is consistent with that of other scholars [8, 17]. The shift to middle age-preponderance of urolithiasis can be attributed to multiple factors, such as changing work, diets, and lifestyles [5, 18, 19]. For example, to prevent or treat osteoporosis, women aged between 50 years and 79 years often consume calcium and vitamin D supplements, which can cause hypercalciuria. Consequently, this population may have a high incidence of urolithiasis [20]. Consistent with previous studies [3, 19, 21, 22], male predominance in the incidence of urolithiasis was also revealed in our study. From 1990 to 2019, the male : female ratio of urolithiasis incidence has gradually increased to approximately 3. This gender disparity may be associated with changes in the diet and an increase in metabolic syndromes, such as diabetes and obesity [3].

EAPC was found to be negatively associated with ASIR, implying that urolithiasis increased more slowly in countries with high incidence than in countries with low incidence. The ASIR also differed between SDI quintiles. ASIR was higher in high SDI, middle-high SDI, middle SDI, and low-middle regions than in low SDI regions. Except in low regions, ASIR had a decreasing global pattern in other quintiles. Interestingly, ASIR of the Russian Federation and its neighboring countries, such as Ukraine and Latvia, was substantially higher than that in others. Meanwhile, African countries—Madagascar, South Sudan, and Burundi—had a lower ASIR. This phenomenon is likely to be because of an increase in meat consumption and the prevalence of diabetes and obesity, all of which are risk factors for urolithiasis [23]. Another study pointed out that the diversity and greater numbers of oxalate-degrading bacteria existing in the gastrointestinal tract played a crucial role in preventing urolithiasis in Black South Africans [24].

Although mortality and DALY in both males and females increased in 2019, the ASDR and age-standardized DALY rate decreased year by year. This decreasing pattern in the ASDR and age-standardized DALY rate of urolithiasis over the last three decades is intertwined with surgical innovations and better treatment guidelines [25]. Surgical treatment might not be the first choice, especially for stones in the inferior pole of the kidney, and pharmacological treatment can be beneficial for spontaneous stone passage. Moreover, it is technological advances that have made it possible to opt for less invasive surgical interventions that are safe, effective, and associated with shorter recovery times and lesser discomfort [26].

The proportion of the elderly in both ASDR and age-standardized DALYs rate of urolithiasis increased. It may be because of their poor tolerance to severe infections and susceptibility to complications, including decreased bone density, cardiovascular disease, and chronic kidney disease [27]. Similarly, there was a demonstrable negative relationship among EAPC, the change in ASDR and age-standardized DALYs rate from 1990 to 2019, and their baseline value in 1990. The higher the ASDR and age-standardized DALYs rate in 1990, the more they changed. This finding may be explained by the fact that countries with higher ASDR and age-standardized DALYs rates were more likely to prioritize urolithiasis prevention programs because of public health concerns. However, in our study, the SDI level did not have an impact on ASDR and age-standardized DALYs rate, which may indicate that better healthcare is less important for a decrease in ASDR and age-standardized DALYs rate. Moreover, as seen in a prior study, we also found a significant regional distribution [15]. For instance, ASDR and age-standardized DALYs rates were very high in Eastern Europe, Central Asia, and Southeast Asia, while extremely low in North Africa, the Middle East, Southern Latin America, and Southern Sub-Saharan Africa.

Despite the GBD study providing high-quality estimates of the global urolithiasis burden, several limitations are inevitable. Firstly, limited data can be obtained from lower SDI countries, and the data we extracted is not representative. Therefore, the data only reveals the general situation of a population in specific regions and countries. Secondly, different countries and regions have different diagnoses, screening standards, and monitoring systems for urolithiasis. Thus, there may be differences in the quality of data. Thirdly, because of the lack of data on risk factors, risk factors for urolithiasis were not assessed.

5. Conclusions

Over the last 30 years, ASIR, ASDR, and age-standardized DALYs rates of urolithiasis have decreased. Yet, significant differences in the geographic, age, and sex distributions were observed. Thus, based on the current epidemiological characteristics of urolithiasis, medical resources should be rationally allocated or adjusted.

Abbreviations

ASR:Age-standardized rate
ASIR:Age-standardized incidence rate
ASDR:Age-standardized death rate
EAPC:Estimated annual percentage changes
CI:Confidence interval
DALYs:Disability-adjusted life years
GBD:Global burden of disease study
SDI:Sociodemographic index
UI:Uncertainty interval
YLD:Years lived with disability
YLL:Years of life lost.

Data Availability

The data used in this study is available from the global health data exchange query tool (https://ghdx.healthdata.org/gbd-results-tool).

Ethical Approval

Not applicable.

Not applicable.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors made substantial contributions to the conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; agreed to submit to the current journal; gave final approval of the version to be published; and agreed to be accountable for all aspects of the work. Xiaoyuan Qian and Junlai Wan contributed equally to this work.

Acknowledgments

This study was financially supported by a grant from the National Natural Science Foundation of China (No. 81900647) and Key Research and Development Program of Hubei Science and Technology Department (No. 2020BCB052).

Supplementary Materials

The supplementary material file 1: the supplementary tables and figures are already cited in our manuscript “Epidemiological trends of urolithiasis at the global, regional, and national levels: a population-based study.” The supplementary material file 2: figures are already cited in the manuscript “Epidemiological trends of urolithiasis at the global, regional, and national levels: a population-based study.” (Supplementary Materials)