Burden of silicosis based on the Global Burden of Disease Study 2021: trend analysis of incidence, mortality, and disability-adjusted life years, and projections for the next 30 years
Original Article

Burden of silicosis based on the Global Burden of Disease Study 2021: trend analysis of incidence, mortality, and disability-adjusted life years, and projections for the next 30 years

Xinxin Zhang1# ORCID logo, Luna Zhao1#, Maolang He1#, Xin Huang2, Dong Liu1

1Department of Pulmonary and Critical Care Medicine, The First Affiliated Hospital of Shihezi University, Shihezi, China; 2Shihezi University School of Medicine, Shihezi, China

Contributions: (I) Conception and design: X Zhang, L Zhao, M He; (II) Administrative support: D Liu; (III) Provision of study materials or patients: X Zhang, L Zhao, M He, X Huang; (IV) Collection and assembly of data: X Zhang, L Zhao, M He; (V) Data analysis and interpretation: X Zhang, L Zhao, M He; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Dong Liu, MS. Department of Pulmonary and Critical Care Medicine, The First Affiliated Hospital of Shihezi University, No. 107 North Second Road Hongshan Street, Shihezi 832000, China. Email: 2322800100@qq.com.

Background: Silicosis, an occupational disease caused by chronic silica exposure, has a high global burden and limited treatment options. This study analyzed the epidemiological trends and future projections of silicosis based on data from the Global Burden of Disease (GBD), aiming to provide data support for public health interventions.

Methods: We extracted and analyzed the data on the incidence, mortality, and disability-adjusted life years (DALYs) of silicosis, as well as the age-standardized rate (ASR) of silicosis from the GBD Study 2021. Using these data, we describe the trends in five dimensions: global, regional, national, age, and sex. We used Joinpoint regression software (V.5.2.0) to calculate the average annual percent changes (AAPCs) in the ASRs from 1990 to 2021. Silicosis trends from 2022 to 2050 were predicted using Bayesian Age-Period-Cohort (BAPC) and Autoregressive Integrated Moving Average (ARIMA) models.

Results: From 1990 to 2021, the global incidence of silicosis, number of deaths, and DALYs showed an upward trend. However, the corresponding ASRs all showed decreasing trends, with AAPC values of −1.1% (−1.1% to −1.0%), −2.5% (−2.7% to −2.3%), and −2.5% (−2.7% to −2.3%), respectively. However, the burden of silicosis varied significantly across countries and regions, with China, South Africa, and Chile having a silicosis age-standardized incidence rate (ASIR) and age-standardized DALYs rate (ASDR) well above the global average. In addition, the ASIR and ASDR of silicosis were generally higher among men. There were also differences between regions at different socioeconomic levels, with the DALYs burden of silicosis being lowest among males and females in low social demographic index (SDI) regions and highest among males in medium and high SDI regions. The BAPC model projected a gradual decrease in the silicosis burden from 2022 to 2050.

Conclusions: Although the disease burden of silicosis showed a decreasing global trend from 1990 to 2021, it is still a global public health concern. Effective preventive and curative measures should be taken to address the challenges posed by silicosis and to protect the lives and health of workers.

Keywords: Silicosis; Global Burden of Disease (GBD); age-standardized rate (ASR); Joinpoint regression; prediction


Submitted Aug 19, 2024. Accepted for publication Dec 20, 2024. Published online Feb 27, 2025.

doi: 10.21037/jtd-24-1341


Highlight box

Key findings

• This study observed a decreasing global burden of silicosis. However, the trends differed across countries and regions. In addition, the silicosis age-standardized incidence rate (ASIR) and age-standardized disability-adjusted life years (DALYs) rate (ASDR) were higher in males.

What is known and what is new?

• Previous studies have reported various health metrics of silicosis in different countries and regions worldwide.

• This cross-sectional study provides the most up-to-date estimates of the incidence, number of deaths, and DALYs resulting from silicosis and their temporal trends based on data from 204 countries and regions, with significant differences by sex, region, country, age, and sociodemographic index. In addition, the global ASIR, age-standardized mortality rate, and ASDR of silicosis were projected for the next 30 years.

What is the implication, and what should change now?

• Silicosis continues to pose a public health issue, and it is essential to enhance health interventions in certain regions to reduce the damage inflicted by this condition.


Introduction

Silicosis is a common occupational disease. It is a pneumoconiosis interstitial lung disease that is progressive and irreversible, and it is mainly caused by long-term exposure to high concentrations of free silica (SiO2) (1). Depending on the duration of exposure to SiO2, silicosis can be categorized as acute, accelerated, or chronic, with the chronic type being more common and the accelerated type being relatively rare. Patients commonly experience symptoms such as chest discomfort, difficulty breathing, mild fever, and nocturnal sweating. As the disease progresses, silicosis can lead to a variety of complications, with respiratory failure being the leading cause of death (2,3). A study by Faubry et al. found that other lung diseases occurring after SiO2 exposure, such as lung cancer, tuberculosis, emphysema, and interstitial pulmonary fibrosis, increase the morbidity and mortality of silicosis (4). Despite the high morbidity and mortality of silicosis, the available treatment options are still very limited (5).

Globally, more than 2,000 people are diagnosed with silicosis each year, and there are more than 12,900 deaths annually, while the disability-adjusted life years (DALYs) for silicosis reached 650,000 in 2019 (6). Therefore, silicosis remains a global public health problem. It is estimated that around 230 million people are exposed to SiO2 annually, of which 40.5 million work in the mining industry (7). A recent study showed that around 2 million people in the United States and Europe are exposed to SiO2 each (8). Workers in South Africa are exposed to SiO2 at concentrations ranging from 9% to 39% (9), and a study in southern Spain showed that from 2009 to 2018, a total of 106 people were diagnosed with silicosis (10). Furthermore, according to the Global Burden of Disease (GBD) Study 2019, approximately 12,900 deaths worldwide were attributed to silicosis, and the burden of silicosis was more pronounced in middle and middle–high regions (11).

There is still a lack of up-to-date literature on the global dynamics of silicosis analysis and surveillance. Therefore, we analyzed the overall trends in silicosis between 1990 and 2021 using data from the GBD database to provide more up-to-date information on its incidence and mortality. Moreover, the latest data on the DALYs of silicosis were predicted, and its trends for the next 30 years were projected and expressed as the average annual percent change (AAPC). We aimed to assess the global status of silicosis and its regional distribution characteristics, and to provide a scientific basis and rational recommendations for the development of effective prevention and control measures. We present this article in accordance with the STROBE reporting checklist (available at https://jtd.amegroups.com/article/view/10.21037/jtd-24-1341/rc).


Methods

Data source

Data from the GBD Study 2021 were analyzed secondarily. The study used existing data to assess the health costs of 371 diseases and injuries in 204 countries or regions by age and sex (12). The overall burden on health and socioeconomic factors that results in poor health, disability, and early mortality is known as the disease burden, and studying the consequences of the disease burden promotes the development of public health policy and the rational allocation of medical resources. The Bayesian meta-regression model (DisMod-MR 2.1), as the core estimation tool of the GBD, uses a mixed-effects method to analyze and integrate data across different age groups, time periods, and geographical regions. The model has the ability to correct data bias and fuse data from different sources to produce comprehensive and consistent estimates of the disease burden and disease trends (13). Estimates of silicosis are based primarily on systematic reviews, hospital records, and claims data reports. Demographic methods can then be applied to impute and supplement any missing data (14).

The sociodemographic index (SDI) is a comprehensive scale for assessing regional social and economic development. It considers factors such as the total fertility rate, the average number of years of education for residents aged ≥15 years, and the per capita income distribution. The SDI value ranges from 0 to 1. Based on the SDI value, 204 countries and regions were classified into five different development levels: low, low–middle, middle, middle–high, and high SDI regions (15). Based on data from the Global Health Data Exchange, we collected estimates of the silicosis incidence, mortality, and DALYs, as well as their 95% uncertainty intervals, based on the GBD 2021 data, and we used the International Classification of Diseases (ICD)-11 code CA60.0 to quantify the disease burden of silicosis (16). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). Since this study is based on public data and does not involve any individual information, ethical approval is not required.

Statistical analysis

To minimize the interference caused by factors such as age and sex, and to facilitate meaningful comparisons across populations, we used the age-standardized rates (ASRs) to express the disease incidence, mortality, and DALYs. In this study, the AAPC was used to describe the average temporal trends in the age-standardized incidence rate (ASIR), age-standardized mortality rate (ASMR), and age-standardized DALYs rate (ASDR) from 1990 to 2021. Joinpoint regression software (version 5.2.0) was utilized to compute the AAPC values. Traditional regression models are mainly applicable to the evaluation of trend fitting on a global scale and are not able to show the characteristics of localized changes over the period under study. The Joinpoint regression model is suitable for describing trends in single disease morbidity or mortality rates by dividing the trend line into several segments and describing these segments in a continuous linear fashion, and it does not strictly require a trend to be present in the data series. We utilized the grid search method (GSM) to calculate all possible join points, selecting the one with the smallest mean squared error (MSE) as the optimal inflection point. Building on this, the optimal number of join points was determined using the Monte Carlo permutation test, allowing for a maximum of 5 join points and a minimum of 0 (17). The AAPC and its 95% confidence interval (CI) were used to express these results. In short, if the AAPC along with its 95% CI was >0, this indicated that the ASRs showed an increasing trend. Conversely, if the AAPC was <0, the ASRs were considered to show a decreasing trend. P<0.05 was considered statistically significant.

The Bayesian Age-Period-Cohort (BAPC) model integrates age, period, and cohort effects to analyze and predict disease trends, accounting for demographic changes (18). This model offers a nuanced approach to projecting disease trends by considering the complex interplay of demographic shifts over time (19). The BAPC model handles data uncertainty and complexity better than traditional models (20). A previous study has been conducted to project the incidence and prevalence of silicosis globally, but no studies have been conducted to project mortality and DALYs (21). Therefore, this study used the BAPC model to predict the future trends in silicosis over the next 30 years, aiming to provide guidance for the prevention and treatment of this condition. The BAPC and integrated nested Laplace approximation packages were used for BAPC model forecasting and visualization. To validate the stability of the prediction results, we performed a sensitivity analysis using the autoregressive integrated moving average (ARIMA) model. All statistical analyses were performed using R software (version 4.4.1).


Results

Global and regional burden of silicosis and disease trends

We analyzed the status of silicosis and the trends in the ASRs in 204 different countries and regions globally from 1990 to 2021. Globally, a decrease was observed in the ASIR, ASMR, and ASDR of silicosis (Figure 1A-1C). Similarly, there was an overall decreasing trend in AAPC for ASIR, ASMR, and ASDR (Figure 2A-2C). In 1990, the global ASIR of silicosis was 0.58/105, and it decreased to 0.41/105 in 2021, which was an overall decrease of 28.67%. The AAPC was −1.1% (−1.1% to −1.0%) (Table 1). In 2021, the medium SDI regions had the highest ASIR of 0.58/105, while the low SDI regions had the lowest ASIR of 0.19/105 (Table 1). The regional trend in ASIR showed that East Asia had the highest ASIR from 1990 to 2021 (Figure 3A). In terms of the ASMR, the global ASMR was 0.26/105 in 1990 and 0.12/105 in 2021, a decrease of 53.61%, with an AAPC of −2.5% (−2.7% to −2.3%) (Table 2). Unlike the ASIR, the lowest ASMR in 2021 was in the high SDI regions (0.05/105), while the highest ASMR was in medium–high SDI regions (0.17/105) (Table 2). In terms of the regional dominance of the ASMR from 1990 to 2021, East Asia ranked first (Figure 3B). In 1990, the global ASDR was 6.58/105, and in 2021, the global ASDR was 3.03/105, a decrease of 53.88% from 1990 to 2021, with an AAPC of −2.5% (−2.7% to −2.3%) (Table 3). In terms of the ASDR by SDI, in 2021, medium and high SDI regions had the highest ASDR of 4.19/105, while high SDI regions had the lowest ASDR of 1.16/105 (Table 3). The regional trend in the ASDR from 1990 to 2021 showed that East Asia had the highest ASDR (Figure 3C). Finally, the spatial distribution of the ASRs for silicosis in 2021 showed that the global ASIR was 0.41/105, ASMR was 0.12/105, and ASDR was 3.03/105. The ASIR of China was 1.11/105, ASMR was 0.31/105, and ASDR was 8.35/105, much higher than the global values (Figure 4A-4C).

Figure 1 Changes in the ASIR (A), ASMR (B), and ASDR (C) of silicosis among different SDI regions. ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years rate; SDI, sociodemographic index.
Figure 2 Trend of ASRs by Joinpoint Regression (5 joinpoints), 1990–2021. (A) ASIR. (B) ASMR. (C) ASDR. The predictive mean is shown as solid line and the vertical dashed line indicates where prediction started. ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years rate.

Table 1

The number of incidence and the ASIR due to silicosis in 1990 and 2021, and changes from 1990 to 2021

Item 1990 2021 1990–2021
Number ASIR, per 100,000 Number ASIR, per 100,000 Percentage changes, % AAPCs, %
Global 23,625.82
(20,019.68–27,406.80)
0.58
(0.49–0.67)
35,482.90
(30,527.07–40,382.75)
0.41
(0.36–0.47)
−28.67
(−31.59 to −25.51)
−1.1*
(−1.1 to −1.0)
Sex
   Male 21,454.11
(18,203.61–24,932.16)
1.16
(0.98–1.34)
31,676.19
(27,323.11–36,033.31)
0.79
(0.69–0.90)
−31.69
(−34.81 to −28.54)
−1.2*
(−1.3 to −1.2)
   Female 2,171.71
(1,729.65–2,683.69)
0.10
(0.08–0.12)
3,806.71
(3,107.51–4,595.61)
0.08
(0.07–0.10)
−12.73
(−18.01 to −7.19)
−0.4*
(−0.5 to −0.4)
SDI
   Low 4,349.49
(3,767.46–4,939.51)
0.21
(0.18–0.25)
4,153.57
(3,689.31–4,682.23)
0.19
(0.16–0.21)
−12.12
(−16.50 to −7.48)
−0.4*
(−0.4 to −0.4)
   Low–middle 503.53
(417.32–587.69)
0.27
(0.23–0.31)
1,042.07
(894.91–1,198.47)
0.21
(0.18–0.24)
−20.36
(−23.85 to −16.63)
−0.7*
(−0.8 to −0.7)
   Middle 9,145.47
(7,617.68–10,851.26)
0.83
(0.69–0.98)
15,711.06
(13,343.88–18,006.32)
0.58
(0.50–0.66)
−29.79
(−33.49 to −25.72)
−1.1*
(−1.2 to −1.1)
   High–middle 7,853.07
(6,635.70–9,148.38)
0.78
(0.66–0.91)
11,346.82
(9,740.09–12,960.67)
0.58
(0.50–0.67)
−25.43
(−29.00 to −21.58)
−0.9*
(−1.0 to −0.9)
   High 1,757.88
(1,479.50–2,048.98)
0.40
(0.35–0.45)
3,218.64
(2,769.45–3,702.39)
0.21
(0.19–0.24)
−46.93
(−49.30 to −44.39)
−2.0*
(−2.1 to −2.0)

*, P<0.05. The data in parentheses are 95% uncertainty intervals. ASIR, age-standardized incidence rate; AAPC, average annual percent change; SDI, sociodemographic index.

Figure 3 Age-standardized rates of silicosis globally and in 21 regions categorized by SDI from 1990 to 2021. (A) ASIR. (B) ASMR. (C) ASDR. For each region, the dots from left to right depict the estimated values for each year from 1990 to 2021. The blue line shows the expected incidence, mortality, and DALYs rate based on the SDI alone. ASR, age-standardized rate; ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years rate; SDI, sociodemographic index.

Table 2

The number of mortality and the ASMR due to silicosis in 1990 and 2021, and changes from 1990 to 2021

Item 1990 2021 1990–2021
Number ASMR per 100,000 Number ASMR per 100,000 Percentage changes, % AAPCs, %
Global 9,973.21
(8,310.09–11,730.91)
0.26
(0.22–0.31)
10,223.28
(8,196.29–12,177.60)
0.12
(0.10–0.14)
−53.61
(−63.40 to −41.49)
−2.5*
(−2.7 to −2.3)
Sex
   Male 9,555.08
(7,965.81–11,249.90)
0.58
(0.49–0.68)
9,562.19
(7,525.32–11,561.93)
0.26
(0.20–0.31)
−55.74
(−65.29 to −44.34)
−2.7*
(−2.8 to −2.5)
   Female 418.12
(243.84–655.81)
0.02
(0.01–0.03)
661.09
(359.25–1,032.18)
0.01
(0.01–0.02)
−29.99
(−51.16−5.92)
−1.1*
(−1.3 to −0.9)
SDI
   Low 354.69
(155.84–562.98)
0.18
(0.08–0.29)
592.84
(244.43–945.28)
0.14
(0.06–0.22)
−25.13
(−43.25 to −2.56)
−0.9*
(−1.0 to −0.8)
   Low-middle 545.19
(302.53–888.66)
0.10
(0.05–0.17)
1,023.71
(603.13–1,435.17)
0.08
(0.05–0.11)
−22.27
(−39.44−6.37)
−0.8*
(−1.0 to −0.7)
   Middle 3,160.49
(2,384.66–4,051.17)
0.31
(0.23–0.40)
4,060.48
(3,074.60–5,195.66)
0.16
(0.12–0.20)
−48.83
(−65.24 to −23.87)
−2.2*
(−2.5 to −2.0)
   High-middle 3,551.58
(3,022.91–4,138.04)
0.38
(0.32–0.44)
3,270.27
(2,575.78–4,013.34)
0.17
(0.13–0.21)
−55.09
(−66.37 to −41.60)
−2.6*
(−2.9 to −2.4)
   High 2,353.89
(2,194.90–2,519.66)
0.21
(0.19–0.22)
1,273.50
(1,104.86–1,425.45)
0.05
(0.05–0.06)
−73.80
(−76.70 to −70.45)
−4.2*
(−4.4 to −4.0)

*, P<0.05. The data in parentheses are 95% uncertainty intervals. ASMR, age-standardized mortality rate; AAPC, average annual percent change; SDI, sociodemographic index.

Table 3

The number of DALYs and the ASDR due to silicosis in 1990 and 2021, and changes from 1990 to 2021

Item 1990 2021 1990–2021
Number ASDR, per 100,000 Number ASDR, per 100,000 Percentage changes, % AAPCs, %
Global 269,198.23
(221,980.78–321,864.02)
6.58
(5.45–7.83)
261,261.69
(214,143.27–311,316.62)
3.03
(2.50–3.61)
−53.88
(−63.81 to −41.41)
−2.5*
(−2.7 to −2.3)
Sex
   Male 257,364.95
(213,398.86–308,451.17)
13.75
(11.51–16.34)
243,482.11
(196,132.64–292,163.76)
6.07
(4.94–7.26)
−55.83
(−65.46 to −43.68)
−2.7*
(−2.9 to −2.5)
   Female 11,833.28
(7,552.04–17,301.33)
0.54
(0.35–0.79)
17,779.58
(11,019.65–25,805.03)
0.39
(0.24–0.57)
−27.47
(−46.66−4.25)
−1.0*
(−1.1 to −0.9)
SDI
   Low 9,138.33
(4,169.61–14,457.95)
4.06
(1.83–6.40)
15,155.16
(6,666.68–23,716.53)
2.98
(1.28–4.71)
−26.60
(−44.61 to −5.01)
−1.0*
(−1.1 to −0.9)
   Low-middle 15,189.70
(9,076.66–23,626.12)
2.41
(1.41–3.81)
26,572.77
(17,038.91–35,993.91)
1.82
(1.14–2.49)
−24.56
(−39.48 to −1.75)
−0.9*
(−1.0 to −0.7)
   Middle 101,900.60
(78,559.80–130,356.66)
8.81
(6.80–11.19)
114,545.12
(91,337.78–144,214.91)
4.19
(3.36–5.28)
−52.45
(−66.58 to −32.15)
−2.4*
(−2.6 to −2.2)
   High-middle 91,108.16
(77,484.42–107,832.50)
9.07
(7.74–10.68)
80,800.40
(65,283.09–100,215.10)
4.19
(3.39–5.20)
−53.76
(−65.06 to −39.34)
−2.5*
(−2.7 to −2.3)
   High 51,680.00
(48,357.52–55,548.59)
4.63
(4.33–4.99)
24,130.01
(21,341.72–27,378.55)
1.16
(1.03–1.32)
−74.87
(−77.61 to −71.42)
−4.3*
(−4.6 to −4.0)

*, P<0.05. The data in parentheses are 95% uncertainty intervals. DALYs, disability adjusted life years; ASDR, age-standardized disability-adjusted life years rate; AAPC, average annual percent change; SDI, sociodemographic index.

Figure 4 Age-standardized rates for silicosis in 204 countries and territories in 2021. A, ASIR. B, ASMR. C, ASDR. Abbreviations: ASIR: age-standardized incidence rate; ASMR: age-standardized mortality rate; ASDR: age-standardized disability-adjusted life years rate.

Country-specific burden of silicosis in 2021

Globally speaking, the ASRs linked to silicosis in 2021 varied greatly among nations. In Kuwait, the ASIR was less than 0.004/105, while in China, the ASIR was 1.11/105. Of all countries, 38 had an ASIR below 0.05/105, while seven countries (Kiribati, Paraguay, Portugal, Monaco, Chile, Democratic People’s Republic of Korea, and China) had an ASIR above 0.50/105 (Figure 4A). Antigua and Barbuda had the lowest ASMRs. Sao Tome and Principe had the highest ASMRs (0.53/105), followed by Chile, Mali, Eswatini, China, Lesotho, and Guinea-Bissau, all with ASMRs exceeding 0.30/105 (Figure 4B). As shown by the ASDR, Kuwait had the lowest rate of silicosis, while Sao Tome and Principe had the highest, followed by Eswatini, Mali, and Lesotho (Figure 4C).

Burden of silicosis by age and sex

In all age groups, the incidence of silicosis and the ASIRs were greater in males than in females. The number of incident cases in both sexes increased with age until it began to decline in the 55–59-year age group, while it began to increase again in the 60–64-year age group, gradually declining again in the 65–69-year age group. The 70–74-year age group had the greatest number of incident cases in both sexes. The ASIR in males increased with age until it began to decline in the 90–94-year age group. In females, the ASIR increased with age (Figure 5A). In 2021, the highest ASIR among females was in the middle SDI region, while the highest ASIR among males was in the middle–high SDI region. The lowest ASIR among both males and females was in the low SDI region (Figure 1A).

Figure 5 Age-specific numbers and rates of silicosis by sex in 2021. (A) Crude incidence rate. (B) Crude mortality rate. (C) Crude DALYs rate. DALYs, disability adjusted life years; UI, uncertainty interval.

In all age groups, males exhibited a greater death count and a greater number of DALYs, along with a higher ASMR and ASDR, than females. The number of deaths and DALYs among males and females increased with age. The number of deaths among males did not begin to decline until the age of 75–79 years, and the number of DALYs did not begin to decline until the age of 55–59 years, after which it began to rise again at the age of 60–64 years and gradually declined in the 70–74-year age group. The number of female deaths did not begin to decline until the age of 80–84 years, and the number of DALYs did not begin to decline until the age of 70–74 years. The ASMR and ASDR increased with age in males. The ASMR began to decline after the age of 90–94 years, and the ASDR began to decline after the age of 85–89 years. The ASMR and ASDR in females increased with age (Figure 5B,5C). In 2021, females had the highest ASMR and ASDR in low SDI regions, while males had the highest ASMR and ASDR in middle–high SDI regions. Conversely, the lowest ASMR and ASDR were observed in both sexes in high SDI regions (Figure 1B,1C).

Predictions for silicosis over the next 30 years

The BAPC model was used to observe and predict the trends in the ASIR, ASMR, and ASDR of silicosis from 1990 to 2050. Figure 6 shows that the burden of silicosis is projected to gradually decrease from 2021 to 2050. In 2021, the global ASIR, ASMR, and ASDR for silicosis were 0.44/105, 0.17/105, and 4.26/105, respectively. By 2050, the global ASIR, ASMR, and ASDR of silicosis are projected to decrease to 0.44/105, 0.07/105, and 1.95/105, respectively. Also, we used the ARIMA model to make the predictions. Figure 7 shows that by 2050, the global ASIR, ASMR and ASDR for silicosis are projected to decrease to 0.37/105, 0.00/105 and 1.02/105, respectively.

Figure 6 Trends in observed and predicted ASRs of silicosis from 1990 to 2050 by Bayesian Age-Period-Cohort model. (A) ASIR. (B) ASMR. (C) ASDR. ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disabilityadjusted life years rate.
Figure 7 Trends in observed and predicted ASRs of silicosis from 1990 to 2050 by Autoregressive Integrated Moving Average model. (A) ASIR. (B) ASMR. (C) ASDR. ASIR, age-standardized incidence rate; ASMR, age-standardized mortality rate; ASDR, age-standardized disability-adjusted life years rate.

Discussion

Silicosis is an important component of the global pneumoconiosis burden. This study reported the temporal trends in the ASIR, ASMR, and ASDR for silicosis from 1990 to 2021 across various regions, countries, age groups, and sexes. Since 1990, there has been an upward trend in the worldwide incidence of silicosis and in the number of deaths from silicosis, whereas the total number of DALYs has shown a downward trend. The global ASIR, ASMR, and ASDR of silicosis demonstrated downward trends. Across all age groups, males had a higher ASIR, ASMR, and ASDR than females. These results provide reliable and comprehensive estimates that may contribute to reducing the global burden of this disease.

The incidence of silicosis, along with silicosis-related deaths and DALYs, as well as the ASIR, ASMR, and ASDR, were greater in males than in females. This may be because industries such as road construction, demolition work, and concrete manufacturing pose a higher risk of silicosis development, and most people working in these industries are male (10). The number of deaths in males slightly increased, while the number of deaths in females significantly increased in 2021 compared to in 1990. This may be explained as follows. First, in some traditional industries, such as gold mining, males and females have the same risk of silicosis exposure (22). Second, a variety of new industries, such as jeans sandblasting, artificial stone production, and glass manufacturing, pose an equal risk of exposure to silicosis to both males and females (23-25). Workers in these industries, regardless of their sex, need to be adequately educated and informed about the potential risks of silicosis. In addition, medical personnel must be sensitized to the occupational history of silicosis in both sexes, which is essential for the early prevention and treatment of this condition. An inadequate understanding of the occupational history of silicosis in females leads to difficulties in early prevention, thus increasing the number of deaths (26). These factors may explain the rise in the number of female deaths and remind us to pay more attention to females in the fight against silicosis. The highest incidence, number of deaths, and DALYs among males and females were in the 65–69-, 75–79-, and 70–74-year age groups, respectively. Chronic silicosis has a latency period of up to several decades, which causes many patients to begin experiencing symptoms only after retirement or as they enter old age (27). Some studies have shown that silicosis mainly affects older workers, and that the 65–79-year age group has the highest risk of developing silicosis (28,29), which is consistent with our findings. This is also a warning that we need to strengthen the health management of silicosis in older workers and conduct regular health checks to detect potential health problems at an early stage.

Detailed studies of silicosis, a traditional occupational disease, were published more than a century ago, and effective occupational health interventions became available (30). The present study found that in 1990, the global incidence of silicosis was 24,000, and the number of deaths was 0.9 million, whereas by 2021, the incidence and number of deaths increased to 35,000 and 10,000, respectively. However, the ASIR, ASMR, and ASDR showed an overall decreasing trend, indicating that positive and effective measures have been taken globally to prevent and control silicosis. The global burden of silicosis is mainly borne by medium SDI and medium–high SDI regions. These regions have relatively well-developed industrial systems and have workplaces with crystalline silica occupational hazards, such as the mining and stone processing industries (31,32). Our study also found that the ASIR for silicosis is lowest in low SDI regions. The main reasons for this are closely related to the level of economic development, such as the lack of diagnostic techniques for silicosis, and an inadequate silicosis registration system (33,34).

Silicosis represents a type of pneumoconiosis that results from prolonged exposure to elevated levels of SiO2 (35). SiO2 is present in many industries, such as the mining, metallurgy, and porcelain industries (36). Recently, several studies have reported other workplaces that pose a risk of SiO2 exposure, such as artificial stone production for indoor building materials, which contain high concentrations of SiO2 ranging from 85% to 93% (37,38). The results of one study showed that young workers using artificial stone became an important emerging population for silicosis (39). In addition, exposure to working environments where manufactured stone is processed is the main reason for the recent increase in the number of silicosis cases, especially in Spain, Australia, and Israel (8). Cohort studies from Australia have found that artificial stone has become a new cause of silicosis, with a prevalence rate of 28.2%, which has led to artificial stone being banned by the Australian government (40,41). Despite the long history of silicosis prevention and treatment, attention needs to be paid to the emergence of new exposure conditions, such as artificial stone production, denim sandblasting, and jewelry manufacturing (10,42). A retrospective cohort study found that the level of silicosis-associated lung function decline in workers using artificial stone was more than five-times that of workers using non-artificial stone (43). The above studies suggest that although the pathogenic form of silicosis is constantly changing, the main cause of silicosis has never changed; in essence, the cause is always SiO2.

Currently, there is a great global effort to prevent and control silicosis, but it is still a public health problem and a priority target requiring prevention and control in developing countries. A report from India noted a high local prevalence rate, mainly identified through regular medical check-ups for workers and the development of a national health program to prevent and control silicosis (44). Meanwhile, educational and training interventions have been used in some countries and regions to improve dust control and the use of respiratory protective equipment (3). Socioeconomic development also affects the occurrence of silicosis, with a lack of control of SiO2 in some areas and workplaces not equipped with appropriate protective equipment (45). A study from Zambia found that mining workers had basic knowledge of SiO2 protection, but their employers did not provide appropriate protective equipment (46). The emergence of new exposure scenarios, such as denim garment sandblasting and jewelry production, means that some workers may be unaware of the hazards, resulting in prolonged exposure (47). Meanwhile, the present study on the spatial distribution characteristics of the ASRs of silicosis in 204 countries and regions around the world found that developing countries, such as China, Chile, Paraguay, and South Africa, had higher ASIRs, ASMRs, and ASDRs than the respective rates at the global level. The study showed that the proportion of cases of silicosis diagnosed in a mobile mine in Zimbabwe was 21% (48). Furthermore, the onset of silicosis may contribute to the emergence of additional respiratory conditions, including tuberculosis and lung cancer (49,50). A systematic evaluation and meta-analysis by Ehrlich et al. 2021 showed that silicosis significantly increased the risk of developing tuberculosis, and that this increase in risk was already evident at milder levels of silicosis (51). Another study showed that silicosis leads to the development of lung cancer and exacerbates the formation of pulmonary fibrosis (52). Although progress has been made in the global fight against silicosis, it remains a major public health problem in developing countries, and there is a need to reduce its incidence through better education, training, regulation, and provision of appropriate protective equipment. A better understanding of its association with other respiratory diseases is also needed.

Limitations

This study has several key limitations that should be considered. First, the data on silicosis were mainly obtained from the GBD 2021 database, which is a global database that collects data from many different sources. This means that there may be inconsistencies among the data, as they may be affected by the collection and reporting methods used in different regions, which could lead to biases in the study results. Second, because some symptoms of silicosis may be similar to the symptoms of other pulmonary diseases, this could have led to the true number of cases being underestimated. Finally, the data on silicosis relied on calculations from surveillance data, which may be subject to lag. Although the GBD provides valuable global health data, we must be cautious when interpreting the results of studies on silicosis and consider the potential impact that these limitations may have on data interpretation and policymaking.


Conclusions

Silicosis continues to be a significant public health issue, and it presents a potentially serious health burden on a global scale. Although the ASRs exhibited a downward trend globally from 1990 to 2021, the number of incident cases, mortality, and DALYs increased significantly in females. Therefore, effective preventive and curative measures should be implemented to address the challenges posed by silicosis, increase the attention paid to silicosis in females, reduce its global burden, and protect the lives and health of workers.


Acknowledgments

The authors thank the GBD Project collaborators for sharing their data publicly, allowing relevant analysis and interpretation.


Footnote

Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://jtd.amegroups.com/article/view/10.21037/jtd-24-1341/rc

Peer Review File: Available at https://jtd.amegroups.com/article/view/10.21037/jtd-24-1341/prf

Funding: The research reported in this publication was funded by Corps Guiding Science and Technology Program Projects (grant number 2023ZD019).

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jtd.amegroups.com/article/view/10.21037/jtd-24-1341/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). As this study is based on public data and does not involve any individual information, ethical approval was not required.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Barnes H, Goh NSL, Leong TL, et al. Silica-associated lung disease: An old-world exposure in modern industries. Respirology 2019;24:1165-75. [Crossref] [PubMed]
  2. Cottin V, Spagnolo P, Bonniaud P, et al. Mortality and Respiratory-Related Hospitalizations in Idiopathic Pulmonary Fibrosis Not Treated With Antifibrotics. Front Med (Lausanne) 2021;8:802989. [Crossref] [PubMed]
  3. Anlimah F, Gopaldasani V, MacPhail C, et al. A systematic review of the effectiveness of dust control measures adopted to reduce workplace exposure. Environ Sci Pollut Res Int 2023;30:54407-28. [Crossref] [PubMed]
  4. Faubry C, Grassion L, Raherison-Semjen C. Silicosis and hemoptysis. Rev Mal Respir 2021;38:300-3. [Crossref] [PubMed]
  5. Li T, Yang X, Xu H, et al. Early Identification, Accurate Diagnosis, and Treatment of Silicosis. Can Respir J 2022;2022:3769134. [Crossref] [PubMed]
  6. Rupani MP. A mixed-methods study on impact of silicosis on tuberculosis treatment outcomes and need for TB-silicosis collaborative activities in India. Sci Rep 2023;13:2785. [Crossref] [PubMed]
  7. Moyo D, Ncube R, Kavenga F, et al. The Triple Burden of Tuberculosis, Human Immunodeficiency Virus and Silicosis among Artisanal and Small-Scale Miners in Zimbabwe. Int J Environ Res Public Health 2022;19:13822. [Crossref] [PubMed]
  8. Raanan R, Zack O, Ruben M, et al. Occupational Silica Exposure and Dose-Response for Related Disorders-Silicosis, Pulmonary TB, AIDs and Renal Diseases: Results of a 15-Year Israeli Surveillance. Int J Environ Res Public Health 2022;19:15010. [Crossref] [PubMed]
  9. Brouwer DH, Rees D. Can the South African Milestones for Reducing Exposure to Respirable Crystalline Silica and Silicosis be Achieved and Reliably Monitored? Front Public Health 2020;8:107. [Crossref] [PubMed]
  10. León-Jiménez A, Hidalgo-Molina A, Conde-Sánchez MÁ, et al. Artificial Stone Silicosis: Rapid Progression Following Exposure Cessation. Chest 2020;158:1060-8. [Crossref] [PubMed]
  11. Chen S, Liu M, Xie F. Global and national burden and trends of mortality and disability-adjusted life years for silicosis, from 1990 to 2019: results from the Global Burden of Disease study 2019. BMC Pulm Med 2022;22:240. [Crossref] [PubMed]
  12. Murray CJL. GBD 2021 Collaborators. Findings from the Global Burden of Disease Study 2021. Lancet 2024;403:2259-62. [Crossref] [PubMed]
  13. Global, regional, and national incidence and mortality burden of non-COVID-19 lower respiratory infections and aetiologies, 1990-2021: a systematic analysis from the Global Burden of Disease Study 2021. Lancet Infect Dis 2024;24:974-1002. [Crossref] [PubMed]
  14. Yang M, Wang D, Gan S, et al. Increasing incidence of asbestosis worldwide, 1990-2017: results from the Global Burden of Disease study 2017. Thorax 2020;75:798-800. [Crossref] [PubMed]
  15. Global burden and strength of evidence for 88 risk factors in 204 countries and 811 subnational locations, 1990-2021: a systematic analysis for the Global Burden of Disease Study 2021. Lancet 2024;403:2162-203. [Crossref] [PubMed]
  16. The Lancet. ICD-11. Lancet 2019;393:2275. [Crossref] [PubMed]
  17. Bai Z, Han J, An J, et al. The global, regional, and national patterns of change in the burden of congenital birth defects, 1990-2021: an analysis of the global burden of disease study 2021 and forecast to 2040. EClinicalMedicine 2024;77:102873. [Crossref] [PubMed]
  18. Riebler A, Held L. Projecting the future burden of cancer: Bayesian age-period-cohort analysis with integrated nested Laplace approximations. Biom J 2017;59:531-49. [Crossref] [PubMed]
  19. Dong B, Zhao Y, Wang J, et al. Epidemiological analysis of chronic kidney disease from 1990 to 2019 and predictions to 2030 by Bayesian age-period-cohort analysis. Ren Fail 2024;46:2403645. [Crossref] [PubMed]
  20. Kiyoshige E, Ogata S, O'Flaherty M, et al. Projections of future coronary heart disease and stroke mortality in Japan until 2040: a Bayesian age-period-cohort analysis. Lancet Reg Health West Pac 2023;31:100637. [Crossref] [PubMed]
  21. Liu X, Jiang Q, Wu P, et al. Global incidence, prevalence and disease burden of silicosis: 30 years' overview and forecasted trends. BMC Public Health 2023;23:1366. [Crossref] [PubMed]
  22. Ndlovu N, Richards G, Vorajee N, et al. Silicosis and pulmonary tuberculosis in deceased female South African miners. Occup Med (Lond) 2019;69:272-8. [Crossref] [PubMed]
  23. Jones RM. Artificial Stone Silicosis: Need for Improved Controls. Ann Work Expo Health 2022;66:137-8. [Crossref] [PubMed]
  24. Henneberger PK, Weissman DN. Old adversaries in new places. Occup Environ Med 2018;75:1-2. [Crossref] [PubMed]
  25. Fireman EM, Fireman Klein E. Association between silicosis and autoimmune disease. Curr Opin Allergy Clin Immunol 2024;24:45-50. [Crossref] [PubMed]
  26. Kerget B, Araz O, Yilmazel Ucar E, et al. Female workers' silicosis diagnosis delayed due to gender bias. Occup Med (Lond) 2019;69:219-22. [Crossref] [PubMed]
  27. Barber CM, Fishwick D, Carder M, et al. Epidemiology of silicosis: reports from the SWORD scheme in the UK from 1996 to 2017. Occup Environ Med 2019;76:17-21. [Crossref] [PubMed]
  28. Mazurek JM, Wood JM, Schleiff PL, et al. Surveillance for Silicosis Deaths Among Persons Aged 15-44 Years - United States, 1999-2015. MMWR Morb Mortal Wkly Rep 2017;66:747-52. [Crossref] [PubMed]
  29. Yi X, He Y, Zhang Y, et al. Current status, trends, and predictions in the burden of silicosis in 204 countries and territories from 1990 to 2019. Front Public Health 2023;11:1216924. [Crossref] [PubMed]
  30. Ramkissoon C, Gaskin S, Song Y, et al. From Engineered Stone Slab to Silicosis: A Synthesis of Exposure Science and Medical Evidence. Int J Environ Res Public Health 2024;21:683. [Crossref] [PubMed]
  31. Knight D, Ehrlich R, Cois A, et al. Predictors of silicosis and variation in prevalence across mines among employed gold miners in South Africa. BMC Public Health 2020;20:829. [Crossref] [PubMed]
  32. Apte SH, Tan ME, Lutzky VP, et al. Alveolar crystal burden in stone workers with artificial stone silicosis. Respirology 2022;27:437-46. [Crossref] [PubMed]
  33. Mohammad S. Rehabilitation of Silicosis Victims of District Karauli, Rajasthan, India. Indian J Community Med 2019;44:347-51. [Crossref] [PubMed]
  34. Akgün M, Ergan B. Silicosis in Turkey: Is it an Endless Nightmare or is There Still Hope? Turk Thorac J 2018;19:89-93. [Crossref] [PubMed]
  35. Li R, Kang H, Chen S. From Basic Research to Clinical Practice: Considerations for Treatment Drugs for Silicosis. Int J Mol Sci 2023;24:8333. [Crossref] [PubMed]
  36. Jumat MI, Hayati F, Syed Abdul Rahim SS, et al. Occupational lung disease: A narrative review of lung conditions from the workplace. Ann Med Surg (Lond) 2021;64:102245. [Crossref] [PubMed]
  37. Mandler WK, Qi C, Qian Y. Hazardous dusts from the fabrication of countertop: a review. Arch Environ Occup Health 2023;78:118-26. [Crossref] [PubMed]
  38. Hoy RF, Jeebhay MF, Cavalin C, et al. Current global perspectives on silicosis-Convergence of old and newly emergent hazards. Respirology 2022;27:387-98. [Crossref] [PubMed]
  39. Hua JT, Zell-Baran L, Go LHT, et al. Demographic, exposure and clinical characteristics in a multinational registry of engineered stone workers with silicosis. Occup Environ Med 2022;79:586-93. [Crossref] [PubMed]
  40. Hoy RF, Dimitriadis C, Abramson M, et al. Prevalence and risk factors for silicosis among a large cohort of stone benchtop industry workers. Occup Environ Med 2023;80:439-46. [Crossref] [PubMed]
  41. Kirby T. Australia bans engineered stone to prevent silicosis. Lancet Respir Med 2024;12:e18. [Crossref] [PubMed]
  42. León-Jiménez A, Mánuel JM, García-Rojo M, et al. Compositional and structural analysis of engineered stones and inorganic particles in silicotic nodules of exposed workers. Part Fibre Toxicol 2021;18:41. [Crossref] [PubMed]
  43. Quan H, Wu W, Yang G, et al. Risk Factors of Silicosis Progression: A Retrospective Cohort Study in China. Front Med (Lausanne) 2022;9:832052. [Crossref] [PubMed]
  44. Rupani MP. Challenges and opportunities for silicosis prevention and control: need for a national health program on silicosis in India. J Occup Med Toxicol 2023;18:11. [Crossref] [PubMed]
  45. Handra CM, Gurzu IL, Chirila M, et al. Silicosis: New Challenges from an Old Inflammatory and Fibrotic Disease. Front Biosci (Landmark Ed) 2023;28:96. [Crossref] [PubMed]
  46. Sifanu M, Taylor TK, Kalebaila KK, et al. Knowledge, Attitude, Behavior Practices and Compliance of Workers Exposed to Respirable Dust in a Zambian Copper Mine. Int J Environ Res Public Health 2023;20:6785. [Crossref] [PubMed]
  47. Murgia N, Akgun M, Blanc PD, et al. Issue 3-The occupational burden of respiratory diseases, an update. Pulmonology 2024;S2531-0437(24)00045-X.
  48. Moyo D, Kavenga F, Moyo F, et al. Health Screening Strategies for Artisanal and Small-Scale Miners for Tuberculosis, Human Immunodeficiency Virus and Silicosis: A Case of the USAID-Supported Kunda Nqob'iTB Project in Zimbabwe. Int J Environ Res Public Health 2024;21:70. [Crossref] [PubMed]
  49. Nishida C, Yatera K. The Impact of Ambient Environmental and Occupational Pollution on Respiratory Diseases. Int J Environ Res Public Health 2022;19:2788. [Crossref] [PubMed]
  50. Jamshidi P, Danaei B, Arbabi M, et al. Silicosis and tuberculosis: A systematic review and meta-analysis. Pulmonology 2025;31:2416791. [Crossref] [PubMed]
  51. Ehrlich R, Akugizibwe P, Siegfried N, et al. The association between silica exposure, silicosis and tuberculosis: a systematic review and meta-analysis. BMC Public Health 2021;21:953. [Crossref] [PubMed]
  52. Krabbe J, Steffens KM, Drießen S, et al. Lung cancer risk and occupational pulmonary fibrosis: systematic review and meta-analysis. Eur Respir Rev 2024;33:230224. [Crossref] [PubMed]
Cite this article as: Zhang X, Zhao L, He M, Huang X, Liu D. Burden of silicosis based on the Global Burden of Disease Study 2021: trend analysis of incidence, mortality, and disability-adjusted life years, and projections for the next 30 years. J Thorac Dis 2025;17(2):872-886. doi: 10.21037/jtd-24-1341

Download Citation