Skip to main content
Download PDF

Survival and disease burden analyses of occupational pneumoconiosis during 1958–2021 in Huangshi city, China: a retrospective cohort study

BMC Public Health volume 24, Article number: 1437 (2024) Cite this article

Abstract

Background

Pneumoconiosis, a chronic disease stemming from prolonged inhalation of dust particles, stands as a significant global burden of occupational diseases. This study aims to investigate the survival outcomes of pneumoconiosis patients in Huangshi city, China, while also evaluating the disease burden on afflicted patients.

Methods

Data for this study were sourced from the Huangshi Center for Disease Control and Prevention. Survival analyses of pneumoconiosis patients were conducted employing life tables and the Kaplan-Meier method. The Cox proportional hazards models were deployed to identify factors influencing pneumoconiosis patients’ survival duration. Competing risks models were employed to confirm the validity of the model outcomes. Additionally, in the disease burden assessment, disability-adjusted life years (DALYs) were computed for various demographic groups and time frames.

Results

A total of 5,641 pneumoconiosis cases, diagnosed in Huangshi City, Hubei Province between 1958 and 2021, were incorporated into the cohort analysis. The probability of mortality and the risk ratio increased with advancing age. Notably, the median survival time of stage III pneumoconiosis patients was significantly shorter compared with those in stages I and II. The Cox proportional hazards model and competing risks analyses underscored several significant factors influencing survival time, including dust exposure duration (HR = 1.197, 95% CI: 1.104–1.298), age at first diagnosis (HR = 3.149, 95% CI: 2.961–3.349), presence of silicosis (HR = 1.378, 95% CI: 1.254–1.515), and stage II-III pneumoconiosis (HR = 1.456, 95% CI: 1.148–1.848). Cumulatively, DALYs amounted to 7,974.35 person-years, with an average of 1.41 person-years. The period between 2000 and 2019 witnessed the highest disease burden.

Conclusion

Our findings highlight the urgent need for improved prevention, earlier detection, and more effective management strategies for the occupational pneumoconiosis population. This study not only underscores the persistent issue of pneumoconiosis in industrial environments but also serves as a crucial call to action for policymakers and healthcare providers.

Peer Review reports

Introduction

Pneumoconiosis is a chronic occupational disease caused by long-term inhalation of dust particles, and it is a major public health problem globally [1]. In developing countries where coal is the primary source of energy, millions of workers are exposed to coal dust during their occupational activities. Coal worker pneumoconiosis (CWP) constitutes a significant portion of occupational diseases [2]. As of 2021, there have been reported 1.025 million cases of occupational diseases nationwide, with 915,000 cases attributed to pneumoconiosis, and nearly half of them being CWP. The number of new CWP cases remains alarmingly high each year [3]. Common clinical symptoms of pneumoconiosis include coughing, difficulty breathing, phlegm production and chest pain. Currently, there is no effective treatment for the lung damage caused by pneumoconiosis, and treatment mainly involves symptomatic and rehabilitative care [1]. As a lifelong occupational disease, pneumoconiosis leads to reduced patient lifespan and workability, causing significant health and economic losses [4].

Over the recent decades, China has emerged as one of the nations most impacted by industrial dust exposure on a global scale. Amidst its swift economic expansion and surging demand for mineral resources, an array of smaller mines and smelting operations across the country continues to rely on archaic production methods, suffer from a scarcity of proper protective gear, and grapple with subpar health oversight [5]. Consequently, China stands as one of the most heavily afflicted countries by occupational dust hazards internationally, evidenced by its substantial number of pneumoconiosis sufferers and the considerable portion of the populace subjected to dust exposure [6, 7]. According to national prevalence assessments, there were approximately 525,000 cases of pneumoconiosis in China from 1949 to 1996, increasing to 600,000 cases by 2000 [8]. According to data released by the National Health Commission, in 2019 there were 15,898 new cases of pneumoconiosis in China, which represents a significant decrease compared to previous years [9]. Disability adjusted life years (DALYs), as an indicator of disease burden, can be used to compare the health losses of different regions and populations, to provide a vital basis for the rational allocation of health resources and the evaluation of health policies [10]. In 2017, the global total DALYs attributable to occupational pneumoconiosis was 507,425 person-years, of which China accounted for 48.80% (247,619 person-years) [11]. In a study conducted in Jiangsu Province, the cumulative DALYs caused by pneumoconiosis was 154,500.83 person-years in 1956–2021 [12]. According to Cynthia Lu et al., the estimated prevalence of coal workers’ pneumoconiosis in the United States at the beginning of the 21st century exceeded that of the 1990s [13]. It can be seen that pneumoconiosis has brought significant harm both domestically and internationally that cannot be ignored.

In an effort to prevent pneumoconiosis, the International Labor Organization (ILO) and the World Health Organization (WHO) launched a global program to eliminate silicosis by 2030 [14]. In 2016, China released the National Occupational Disease Control Plan (2016–2020), aiming to promote public health and further protect laborers’ rights. Although the government has made significant efforts, pneumoconiosis remains one of the major social burdens.

Huangshi city, situated in Hubei Province of China, boasts abundant mineral resources, particularly in the fields of metallurgy, coal, and cement materials. The surge in mineral-related industries resulted in a significant rise in pneumoconiosis cases in Huangshi city, making it the primary area for such cases within Hubei Province [15]. Therefore, the prevention and treatment of pneumoconiosis in this area are of paramount importance. We analyzed data on patients diagnosed with occupational pneumoconiosis from 1958 to 2019 in Huangshi city. The aim of this study was to investigate the epidemiologic characteristics and survival status of patients with pneumoconiosis in Huangshi city, as well as to identify the risk factors affecting patient survival. The study further explored the disease burden of pneumoconiosis and its temporal trends, providing an important basis for designing comprehensive prevention and control measures for occupational pneumoconiosis.

Methods

Study setting and population

Huangshi city is composed of six administrative districts. The patient data incorporated in this research were sourced from a range of medical units across these districts. Specifically, the data pertaining to pneumoconiosis patients utilized in this study were systematically gathered from a variety of mining and industrial operations within Huangshi city, encompassing coal mines, metal mining sites, cement production facilities, among others, as detailed in Figure S1.

Inclusion criteria: (1) The patients must have been diagnosed by a qualified institution in Huangshi city, with the diagnosis confirmed by at least three radiologists in accordance with the “Diagnostic Criteria for Pneumoconiosis”; (2) Follow-up duration of at least 2 years, or until death, or up to December 31, 2021, with initial exclusions based on these criteria; (3) Comprehensive records, including data on age at birth and history of environmental exposure.

Cause of death encompassed fatalities resulting from pneumoconiosis or its associated complications, including pulmonary heart disease, tuberculosis, and others. Censored cases were delineated as individuals whose demise stemmed from alternate causes, such as accidents, liver failure, renal and autoimmune diseases, among others, or those who remained alive at the conclusion of the observation period.

Diagnostic criteria of pneumoconiosis

The Occupational Disease Control Section of the Centers for Disease Control and Prevention (CDC) conducts annual follow-up on cases of pneumoconiosis for up-to-date data. Pneumoconiosis diagnosis adheres to diagnostic criteria agreed upon by no less than three radiologists. The diagnostic criteria for pneumoconiosis spanning 1956–2021 are outlined below: “Medical Preventive Measures for Silica Dust Workers”, “GB5906-1986 Dust Lung X-ray Diagnosis Standard and Treatment Principles”, “GB5906-1997 Dust Lung X-ray Diagnosis”, “GBZ 70-2002 Pneumoconiosis Diagnostic Standard 13”, “GBZ 70-2009 Diagnostic Standard 14 and Pneumoconiosis Diagnostic Standard 12 in 2015”. All procedures in this study were approved by the Ethics Committee of Huangshi Disease Prevention and Control Center.

The manifestations of pneumoconiosis on X-ray chest radiographs can be categorized into three stages. Stage I pneumoconiosis is characterized by small opacities with an overall profusion of category 1, distributed over at least 2 lung zones. Stage II pneumoconiosis is indicated by small opacities with an overall profusion of category 2, extending beyond 4 lung zones; or small opacities with an overall profusion of category 3, distributed across 4 lung zones. Stage III pneumoconiosis is defined by one of the following conditions: the presence of large opacities with a long axis of at least 20 mm and a short axis of at least 10 mm; small opacities with an overall profusion of category 3, extending beyond 4 lung zones with coalescence of opacities; or small opacities with an overall profusion of category 3, extending beyond 4 lung zones along with the presence of large opacities [16].

Survival analyses

We expressed the categorical variables as frequencies and percentages (%). The chi-square test was used to compare the differences among categorical variable groups, while Analysis of Variance (ANOVA) was employed to assess differences among groups for continuous variables. Survival time was demarcated as the duration from the first diagnosis to either the date of death or the end of the study. Further scrutinizing survival across temporal cohorts within this study population, we employed a life table methodology. This analysis considered various age brackets at the time of initial diagnosis, with an initiation age of 25 and ≥ 70 years, segmented by five-year intervals.

For the detailed survival analysis, we applied the Kaplan-Meier method for the subgroup evaluations and Cox proportional hazards regression models to examine risk factors. The selection of variables for these models was grounded on findings from prior research, with a preliminary univariate Cox proportional hazards model screening where variables showing P < 0.2 were considered for further analysis. These selected variables were subsequently incorporated into a multivariate Cox proportional hazards model, focusing on pneumoconiosis patient survival time as the dependent variable. This was done with the aim of pinpointing significant risk factors influencing patient survival.

Model discrimination was assessed by utilizing receiver operating characteristic (ROC) curves, substantiated by evaluating the area under the curve (AUC), to reflect the distinct capabilities of the model based on subject work characteristics. Concurrently, calibration of the nomograms for the study cohort underwent graphical assessment through the deployment of calibration plots. The Brier score, serving as a precise metric for calibration curve performance, indicates higher model accuracy with decreasing values. A Brier score ranging from 0 to 0.1 denotes an excellent model, whereas scores from 0.1 to 0.25 are generally recognized as good. Comprehensive analyses were executed using R software version 4.3.1 (R Foundation for Statistical Computing, Vienna, Austria). The P < 0.05 was considered significant.

Estimating DALYs

The DALY metric quantitatively measures disease burden by aggregating years of life lost (YLL) due to premature mortality and years lived with disability or disease (YLD), such that DALY = YLL + YLD. According to the method provided by the World Health Organization (WHO) [17], we chose not to discount YLLs, YLDs, or DALYs for time:

$${\rm{YL}}{{\rm{L}}_{{\rm{(c,s,a,t)}}}}{\rm{ = }}{{\rm{N}}_{{\rm{(c,s,a,t)}}}}{\rm{ \times }}{{\rm{L}}_{{\rm{(s,a)}}}}$$
(1)

where: N(c, s,a, t) is the number of deaths due to the cause c for the given age (a) and sex (s) in year (t). L(s, a) is a standard loss function specifying years of life lost for a death at age (a) for sex (s). The life expectancy at different age groups is based on the abridged life tables released by the Department of Economic and Social Affairs Population Division, with data sources as follows: https://population.un.org/wpp/Download/Standard/Mortality/.

$${\rm{YL}}{{\rm{D}}_{{\rm{(c,s,a,t)}}}}{\rm{ = }}{{\rm{I}}_{{\rm{(c,s,a,t)}}}}{\rm{ \times }}{{\rm{D}}_{{\rm{(c,s,a)}}}}{\rm{ \times }}{{\rm{L}}_{{\rm{(c,s,a,t)}}}}$$
(2)

where: I(c, s,a, t) is the incidence (new cases) for a specific cause (c), sex (s), age group (a), and time period (t). D(c, s, a) is the disability weight (DW) assigned to the condition for a specific cause (c), sex (s), and age group (a). The DW for pneumoconiosis is deemed equivalent to chronic obstructive pulmonary disease (COPD), with values of 0.019 for stage I, 0.225 for stage II, and 0.408 for stage III [11]. L(c, s,a, t) is the average duration of the condition (in years) for a specific cause (c), sex (s), age group (a), and time period (t).

Result

Cohort assembly and profiling of study articipants

We collected data on a cohort of 8,756 pneumoconiosis patients in Huangshi city, Hubei Province, from 1958 to 2021. In constructing the cohort, we first eliminated patients with incomplete data, especially those instances lacking critical birth details and exposure history (N = 2,673). Subsequently, we removed 207 cases that did not meet the minimum follow-up requirement of two years, as well as 235 patients who were excluded due to data inconsistencies, such as contradictory or illogical information. The remaining data of eligible patients were then incorporated into the definitive study cohort (N = 5,641) (Fig. 1).

Fig. 1
figure 1

Flow diagram of the study population

Full size image

Demographic and clinical characteristics of pneumoconiosis patients

This study encompassed a cohort of 5,641 patients diagnosed with pneumoconiosis from Huangshi city, spanning from 1958 to 2021. Of these, a total of 2,388 fatalities were recorded. The period between 2000 and 2019 witnessed the greatest surge in both new cases and deaths, accounting for 3,807 and 2,317 respectively. Within the cohort, the prevalence of stage I pneumoconiosis notably high (92.34%), and the elderly population, aged 65 and above, represented 14.47% of this group. A temporal shift in the cohort’s entry point was observed to influence both the average duration of dust exposure and the age at initial diagnosis. The cohort members who joined between 2000 and 2019 presented the highest figures in these categories, averaging 23.53 years for dust exposure and 57.74 years for age at first diagnosis. (Table S1).

As delineated in Table 1, the patients within this pneumoconiosis study predominantly fell into stage I of the disease. Across the three delineated stages of pneumoconiosis, the frequency of male patients was substantially greater than that of females. When examining the types of pneumoconiosis among patients, coal workers’ pneumoconiosis was the most prevalent, documented in 4,027 cases (71.4%). This was followed by silicosis, recorded in 1,470 cases (26.1%). The age group from 45 to 59 years old represented the majority at first diagnosis, encompassing 51.9% of the patients. Notably, with 15–24 years of experience in dust exposure made up the largest portion with 2,667 or 47.3% of the total number of patients. All variables were statistically significant (P < 0.001) at all three stages except gender.

Table 1 Demographic and occupational characteristics
Full size table

Survival analyses of pneumoconiosis patients

Table 2 presents the probability death and hazard ratios for pneumoconiosis across different age groups at first diagnosis, split into 5-year increments. An upward trend in mortality risk was noted (0.002–0.122), with hazard ratios also rising with age (0.001–0.28). Table S2 breaks down case-fatality rates and median survival by stage (I, II, III) and age subgroup, revealing shorter survival among older patients and median survival times being 18, 14, and 9 years across stages, respectively. Using the Kaplan-Meier method, estimated survival of the 5,641 pneumoconiosis patients in Huangshi city from 1958 to 2021. Survival differences across the three pneumoconiosis types were not significant (P = 0.121). Similarly, there was no notable survival difference when comparing stages I and II-III at first diagnosis (P = 0.283). However, patients with ≥ 25 years of dust exposure had significantly lower survival compared to others (P < 0.05). Furthermore, survival varied significantly with age at diagnosis: those ≤ 44 years had the best prognosis, while survival decreased for those 45–49 years, and further for those ≥ 60 years (P < 0.05) (Fig. 2).

Table 2 Cumulative survival rate and probability of death in 5,641 patients with pneumoconiosis: grouping by age of first diagnosis
Full size table
Fig. 2
figure 2

Survival curves of 5,641 pneumoconiosis patients across various cohorts in Huangshi city, Hubei Province. The cohorts were categorized based on the following criteria: (A) the type of pneumoconiosis; (B) the stage of pneumoconiosis; © the duration of dust exposure; (D) the age at initial diagnosis. CWP, coal workers’ pneumoconiosis

Full size image

Cox proportional hazards analyses of pneumoconiosis patients

The model examined six factors for pneumoconiosis survival: dust exposure duration, age at first diagnosis, stage of disease, type of pneumoconiosis, medical insurance, and industry. Through the univariate cox model results, variables with P < 0.2 were included in the multivariate cox model, which included four variables: dust exposure years, age at first diagnosis, stage of disease, and type of pneumoconiosis. The multivariate cox model identified significant risk factors affecting survival time with statistical significance (P < 0.05): Dust exposure years (HR = 1.197, 95% CI: 1.104–1.298), age at first diagnosis (HR = 3.149, 95% CI: 2.961–3.349), silicosis (HR = 1.378, 95% CI: 1.254–1.515), other types of pneumoconiosis (1.456, 95% CI: 1.148–1.848) (Table 3).

Table 3 Cox regression analysis in 5641 patients with pneumoconiosis
Full size table

Utilizing the multivariate model, nomogram plots were created, and risk factor scores were aggregated to predict 15-years and 30-years survival probability; older age at diagnosis was especially predictive (Fig. 3). ROC and calibration curves for 15-years and 30-years survival showed good model discrimination with AUC (0.851 and 0.831). Calibration additionally showed a close concordance between predicted and observed survival, with Brier scores of 0.096 at 15 years and 0.161 at 30 years, confirming the model’s calibration (Fig. 4).

Fig. 3
figure 3

The nomogram of multivariate Cox proportional hazards model

Full size image
Fig. 4
figure 4

The ROC curve and calibration plot of the nomogram. (A) and (B) indicate the ROC curve and AUC of the nomogram in predicting 15- and 30-years survival. The calibration plot for predicting 15-years (C) and 30-years (D) survival. Actual rate of survival is shown on the y-axis, and the nomogram-predicted probability of survival is shown on the x-axis. ROC, receiver operating characteristic; AUC, area under the ROC curve

Full size image

Competing risks models of pneumoconiosis patients

The study employed the competing risks models to assess varying causes of death and conducted sensitivity analyses to confirm the reliability of the multivariate Cox model findings. The analysis indicated that factors such as age at first diagnosis, dust exposure time, silicosis, and stages II-III consistently emerged as potential risk factors in pneumoconiosis, evidencing the model’s stability (Table S3).

Disease burden of pneumoconiosis

Between 1958 and 2021 in Huangshi city, the total DALYs were 7,974.35 person-years, with YLDs of 3,753.37 person-years and YLLs of 4,220.98 person-years. Table 4 details the disease burden across different pneumoconiosis features. In terms of total DALYs, the following were observed for each pneumoconiosis condition: silicosis accounted for 1,863.00 person-years, CWP for 5,894.22 person-years, cement pneumoconiosis for 192.51 person-years, and other types of pneumoconiosis for 24.62 person-years. CWP represented the highest burden in total DALYs. However, when assessing the average burden per individual case (average DALYs), cement pneumoconiosis was found to be higher than the other types.

Table 4 Disease burden of pneumoconiosis patients with different characteristics
Full size table

Patients first diagnosed at age ≤ 44 years had the highest average DALYs (1.64 person-years), YLDs (0.99 person-years), and YLLs (1.31 person-years) among the age groups. While stage I posted the highest accumulated DALYs at 5,732.19 person-years, the average burden per case for stage III surpassed both stages I and II. Moreover, the elevated average YLLs for stage III, at 6.21 person-years, underscores the grave severity of pneumoconiosis in its advanced stages.

While comparing the disease burden over various epochs, it becomes apparent that the period from 2000 to 2019 witnessed the most substantial impact, characterized by the peak in both overall and average DALYs. This surge is attributed to the significant number of individuals diagnosed within this timeframe. Moreover, the decade spanning 1980 to 1999 saw a notably higher average in YLLs when contrasted with other periods (Table S4).

Discussion

Pneumoconiosis poses a serious threat to the health of workers and places a significant burden on society, affecting millions of workers worldwide who engage in hazardous occupations [18]. It is one of the common occupational diseases in China. The primary findings of the retrospective cohort study on occupational pneumoconiosis conducted in Huangshi city, China from 1958 to 2021 reveal significant insights into the evolution of this occupational disease. The study indicates a distinct trend of decreasing survival rates among patients diagnosed with pneumoconiosis, reflecting prolonged exposure to hazardous materials without adequate protective measures. Additionally, the burden of disease, measured in terms of morbidity and mortality, has shown a substantial increase, highlighting the severe impact on workers’ health and the strain on healthcare systems. The data emphasize the critical need for improved occupational safety standards and health monitoring practices to prevent the onset and progression of pneumoconiosis among workers in industrial settings.

Previous studies indicated that the regional distribution of occupational pneumoconiosis is closely tied to the structure of industry [19]. The three most common types of pneumoconiosis included asbestosis, silicosis and CWP. Specifically, silicosis is more prevalent than asbestosis in developing countries, whereas in developed nations, the situation is reversed [20, 21]. According to China’s national occupational disease report covering the years 2015 to 2016, CWP and silicosis were listed as the top two reported new cases of pneumoconiosis [19]. Huangshi city has abundant coal resources, so most of the pneumoconiosis patients in this study were coal miners, followed by silicosis patients. This is different from the distribution of the main types of pneumoconiosis in Guangdong Province [22], Xinjiang Uygur Autonomous Region [23], and Tianjin [24]. For example, silicosis ranks first in Guangdong Province, followed by welders’ pneumoconiosis and other types of pneumoconiosis. Zhejiang Province mainly has non-ferrous and non-metallic mineral resources, so there are more cases of silicosis. Stage I (92.34%) still accounts for the majority of pneumoconiosis patients, followed by stage II and III, which is consistent with the characteristics of pneumoconiosis reported in Shanghai [5] and Tianjin [24].

Compared with stages I and II, the survival time of patients at stage III is much shorter, which is also consistent with previous studies [25,26,27]. This is because pneumoconiosis is a chronic disease with no effective cure, and its progression is accompanied by other complications and fibrosis of the lungs, leading to damage to lung function and affecting patient survival time [28]. It can be seen that prevention and control of disease progression is of great significance for improving the quality of life and survival time in patients with pneumoconiosis.

This study grouped pneumoconiosis patients with different characteristics and constructed their survival curves. In terms of different types of pneumoconiosis, this study found no significant differences in the survival curves among the three categories, which is consistent with a study conducted in Jiangsu Province [12]. However, we did not observe statistically significant differences in survival among different stages of pneumoconiosis in this study, which is inconsistent with a previous study conducted in Zhejiang Province [18]. This discrepancy in findings may be attributable to the smaller number of stage III patients in our study and the relatively similar survival status between stage II and stage I patients. Both the grouping of dust exposure years and the grouping of age at first diagnosis showed that the survival curves of the high-exposure group were lower than those of the low-exposure group, and the survival curves of the elderly group were lower than those of the young group. These results are consistent with previous studies, suggesting that prolonged exposure to dust and delayed diagnosis may both negatively impact the survival time of pneumoconiosis patients.

The Cox proportional hazards model analysis results of this study show that dust exposure time, first diagnosis age, silicosis, other types of pneumoconiosis, and stage II-III pneumoconiosis are important risk factors that affect patient survival time, which is similar to the results of other studies [29]. Dust exposure can influence patient mortality. This may be due to the cumulative effects of dust exposure from extended periods of working in industries with high dust levels. Relevant studies have also shown that even dust with relatively small damage to the human body can cause significant damage to patients’ health when the cumulative dust exposure is large [30]. The age of first diagnosis also affects pneumoconiosis patients’ survival to a certain extent. Patients diagnosed with pneumoconiosis at a younger age have better immunity, which is more conducive to preventing disease progression [31]. Patients with silicosis had a higher hazard ratio than those with CWP, which may be due to the higher content of free silica, which is more harmful to health. Relevant studies also show that patients with silicosis have a lower survival rate than other types of pneumoconiosis patients [32]. In addition, a diagnosis of stage II-III significantly impacts survival time compared to patients in stage I of the disease. Stage II-III pneumoconiosis represents a more advanced stage of disease progression, characterized by further pathological changes in the lungs as the stage of pneumoconiosis increases. The exacerbation of pulmonary fibrosis and more severe impairment of lung function have a direct impact on patients’ quality of life and overall survival [12]. Therefore, it is imperative to conduct regular health check-ups for the early detection of pneumoconiosis and to facilitate timely interventions. Doing so can prevent disease progression and enhance the treatment outcomes and prognosis for patients.

The cumulative DALYs of Huangshi city’s pneumoconiosis patients included in this study from 1958 to 2021 were 7,974.35 person-years, of which the cumulative YLDs were 3,753.37 person-years, and the cumulative YLLs were 4,220.98 person-years. The main contributor to the burden of pneumoconiosis in Huangshi city is the disability resulting from premature death, which contradicts the findings of previous studies [33]. This discrepancy can be attributed to the different weighting of pneumoconiosis disability in this study, considering various stages, as opposed to previous studies that used a fixed value, resulting in lower YLDs estimates. The mean YLLs in this study were 1.77 person-years, which was lower than the previous mean YLLs of 6.35 person-years in the Chinese pneumoconiosis population from 1990 to 2017. The mean YLDs were 0.67 person-years, lower than the 6.22 person-years in the previous study [34]. The average DALYs were higher in cement pneumoconiosis than in other types of pneumoconiosis, but the mean YLLs were higher in silicosis and coal miners’ pneumoconiosis than in cement pneumoconiosis. This may be due to the fact that silicosis and coal miners’ pneumoconiosis cause more severe damage to lung function, leading to premature death. In contrast, patients with cement pneumoconiosis may have slower disease progression, resulting in a higher YLDs. As the age of diagnosis increases, the average DALYs and average YLDs decrease. Early diagnosis of pneumoconiosis correlates with greater loss of healthy life years, as indicated by higher YLDs at younger ages. Conversely, older patients often exhibit weaker immune responses and overall health, heightening mortality risks. Furthermore, a diagnosis of pneumoconiosis at a later stage was found to be associated with increased average YLLs, which supports the COX proportional hazards model indicating that a more advanced initial stage reduces survival time in pneumoconiosis patients. Late diagnosis leads to higher mean YLLs due to disease progression and greater lung damage, resulting in a decreased life expectancy.

This study found that the highest disease burden, as measured by cumulative DALYs, occurred between 2000 and 2019, likely due to enhanced pneumoconiosis screening efforts, resulting in increased diagnosis of new cases. Advances in medical technology extended patient survival and decreased mortality, explaining the lower average YLLs in the 2000–2019 period compared to 1980–1999.

The study faces several limitations. Its retrospective nature introduces potential bias, and the absence of data on disease stage, patients’ economic status, and treatment costs impedes a full economic burden assessment. Besides, the lengthy time frame of the data and shifting testing criteria over different periods may impact the study’s findings.

Conclusion

The historical retrospective cohort study conducted in Huangshi city, China, covering the years 1958–2021, has provided critical insights into the survival trends and disease burden of occupational pneumoconiosis. The findings reveal a significant impact on the affected populations, underscoring the urgency for enhanced preventive measures, early detection, and effective management strategies. This study not only highlights the persistent challenge of pneumoconiosis in industrial settings but also serves as a crucial call to action for policymakers, healthcare providers, and stakeholders to mitigate occupational health risks and improve the quality of life for individuals suffering from this debilitating condition. Furthermore, it underscores the importance of ongoing research and tailored interventions to address the evolving landscape of occupational health hazards.

Data availability

The datasets used during the current study available from the corresponding author on reasonable request (tonygqin@ntu.edu.cn).

References

  1. Qi XM, Luo Y, Song MY, Liu Y, Shu T, Liu Y, Pang JL, Wang J, Wang C. Pneumoconiosis: current status and future prospects. Chin Med J. 2021;134(8):898–907.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Han L, Han R, Ji X, Wang T, Yang J, Yuan J, Wu Q, Zhu B, Zhang H, Ding B, et al. Prevalence characteristics of coal workers’ pneumoconiosis (CWP) in a state-owned mine in Eastern China. Int J Environ Res Public Health. 2015;12(7):7856–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. National Bureau of Statistics. China statistical yearbook. 2022. [http://www.nhc.gov.cn/mohwsbwstjxxzx/tjtjnj/202305/6ef68aac6bd14c1eb9375e01a0faa1fb.shtml]

  4. Bukovitz B, Meiman J, Anderson H, Brooks EG. Silicosis: diagnosis and medicolegal implications. J Forensic Sci. 2019;64(5):1389–98.

    Article  PubMed  Google Scholar 

  5. Mo J, Wang L, Au W, Su M. Prevalence of coal workers’ pneumoconiosis in China: a systematic analysis of 2001–2011 studies. Int J Hyg Environ Health. 2014;217(1):46–51.

    Article  PubMed  Google Scholar 

  6. Wang XT, Zhang G. [Prevalence of pneumoconiosis in a coal enterprise from 2002 to 2013]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua laodong weisheng zhiyebing zazhi = Chinese. J Industrial Hygiene Occup Dis. 2016;34(9):681–2.

    CAS  Google Scholar 

  7. Zhao JQ, Li JG, Zhao CX. Prevalence of pneumoconiosis among young adults aged 24–44 years in a heavily industrialized Province of China. J Occup Health. 2019;61(1):73–81.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Chair SY, Chan JYW, Law BMH, Waye MMY, Chien WT. Genetic susceptibility in pneumoconiosis in China: a systematic review. Int Arch Occup Environ Health. 2023;96(1):45–56.

    Article  PubMed  Google Scholar 

  9. National Health Commission. Statistical bulletin of China’s health development. 2020 [http://www.nhc.gov.cn/guihuaxxs/s10748/202006/ebfe31f24cc145b198dd730603ec4442.shtml].

  10. Jp NA, Imanaka M, Suganuma N. Japanese workplace health management in pneumoconiosis prevention. J Occup Health. 2017;59(2):91–103.

    Article  PubMed  Google Scholar 

  11. Global regional. National disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: a systematic analysis for the global burden of Disease Study 2017. Lancet (London England). 2018;392(10159):1859–922.

    Article  Google Scholar 

  12. Duan Z, Zhou L, Wang T, Han L, Zhang J. Survival and disease burden analysis of occupational pneumoconiosis from 1956 to 2021 in Jiangsu Province. J Occup Environ Med. 2023;65(5):407–12.

    Article  PubMed  Google Scholar 

  13. Lu C, Dasgupta P, Cameron J, Fritschi L, Baade P. A systematic review and meta-analysis on international studies of prevalence, mortality and survival due to coal mine dust lung disease. PLoS ONE. 2021;16(8):e0255617.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Reilly MJ, Timmer SJ, Rosenman KD. The burden of silicosis in Michigan: 1988–2016. Annals Am Thorac Soc. 2018;15(12):1404–10.

    Article  Google Scholar 

  15. Yao Y, Wei T, Zhang H, Xie Y, Gu P, Yao Y, Xiong X, Peng Z, Zhen Z, Liu S et al. Characteristics of diagnosed and death cases of pneumoconiosis in Hubei Province, China, 1949–2019. Int J Environ Res Public Health. 2022;19(23).

  16. Diagnostic criteria for pneumoconiosis. In. Edited by Health CMo; 2009.

  17. Murray CJ, Ezzati M, Flaxman AD, Lim S, Lozano R, Michaud C, Naghavi M, Salomon JA, Shibuya K, Vos T, et al. GBD 2010: design, definitions, and metrics. Lancet (London England). 2012;380(9859):2063–6.

    Article  PubMed  Google Scholar 

  18. Zou H, Shi Z, Zhang Y, Zhou J, Fang X, Zhang Y, Hu Y, Lou X, Zhou L. Epidemiological characteristics and survival analysis on patients with occupational pneumoconiosis in Zhejiang Province from 1987 to 2019. Front Public Health. 2022;10:1006391.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Zhang M, Wang D, Zheng YD, Du XY, Chen SY. [Analyses on the characteristics and the trends of pneumoconiosis notified between 1997 and 2009, in China]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua laodong weisheng zhiyebing zazhi = Chinese. J Industrial Hygiene Occup Dis. 2013;31(5):321–34.

    Google Scholar 

  20. Cullinan P, Reid P. Pneumoconiosis. Prim Care Respir J. 2013;22(2):249–52.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Cheng BW, Su M. [International incidence trend of coal workers’ pneumoconiosis and silicosis]. Zhonghua Lao Dong Wei Sheng Zhi Ye Bing Za Zhi. 2019;37(1):75–8.

    CAS  PubMed  Google Scholar 

  22. Chen M, Wang H, Zhang J, Yu C, Liu W, Xu Y. Distribution of asbestos enterprises and asbestosis cases - China, 1997–2019. China CDC Wkly. 2020;2(18):305–9.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Mohemaiti P, Patigul N, Reshalaiti A, Honglamijiang AS. Characteristics of 3687 pneumoconiosis cases between 1985 and 2006 in Xinjiang Uyghur’s Autonomous Region, China. Occup Environ Med. 2011;68(11):860.

    Article  PubMed  Google Scholar 

  24. Zeng Q, Deng LH, Li XL, Gao Y, Liu J, Yang XY, Tang HJ, Wang X. [Epidemiological characteristics of new pneumoconiosis in different industries in Tianjin]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua Laodong Weisheng Zhiyebing Zazhi = Chin J Industrial Hygiene Occup Dis. 2021;39(2):118–22.

    CAS  Google Scholar 

  25. Xie LZ, Zhou L, Zhao Y, Liu J, Wang W, Zhang W, Lu MX, Yang YX, Zhou ZW, Zhuang JY, et al. [Follow-up and retrospective investigation of patients with pneumoconiosis in Jiangsu Province, China]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua laodong weisheng zhiyebing zazhi = Chinese. J Industrial Hygiene Occup Dis. 2020;38(4):251–5.

    CAS  Google Scholar 

  26. Cui JW, Hu JA, Zhu L. [Survival analysis of patients with pneumoconiosis from 1963 to 2014 in Yueyang]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua laodong weisheng zhiyebing zazhi = Chinese. J Industrial Hygiene Occup Dis. 2016;34(5):360–2.

    CAS  Google Scholar 

  27. Penteado JO, Peres TG, Ramires PF, de Lima Brum R, da Silva Freitas L, Volcão LM, Dos Santos M. Da Silva Júnior FMR: Trends in pneumoconiosis in Brazil, 1979–2019. Occup Med (Lond). 2022;72(6):386–93.

    Article  PubMed  Google Scholar 

  28. Okumura E, Kawashita I, Ishida T. Computerized classification of pneumoconiosis on digital chest radiography artificial neural network with three stages. J Digit Imaging. 2017;30(4):413–26.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Matrat M, Radoï L, Févotte J, Guida F, Cénée S, Cyr D, Sanchez M, Menvielle G, Schmaus A, Marrer E, et al. Occupational exposure to wood dust and risk of lung cancer: the ICARE study. Occup Environ Med. 2019;76(12):901–7.

    Article  PubMed  Google Scholar 

  30. Mao JJ, Chen HL, Li CH, Lu JW, Gu YY, Feng J, et al. Population impact of fine particulate matter on tuberculosis risk in China: a causal inference. BMC Public Health. 2023;23(1):2285.

  31. Tong Y, Kong YY, Bian H, Zheng JZ, Wu YJ, Zhang Y. [Survival and disease burden trend analysis of occupational pneumoconiosis from 1963 to 2020 in Shizuishan City]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua laodong weisheng zhiyebing zazhi = Chinese. J Industrial Hygiene Occup Dis. 2022;40(5):341–7.

    CAS  Google Scholar 

  32. Han B, Yan B, Zhao N, Zhang J, Sun J, Lei X, Li C, Liu H, Chen J. The influence of the functional capacity on subjective well-being and quality of life of patients with silicosis. Aging Ment Health. 2013;17(6):707–13.

    Article  PubMed  Google Scholar 

  33. Cui FT, Shen FH, Chang CF, Xu J, Tang GY, Jiao GL, Gao W, Xu XH, Ding XP. [Analysis of the burden of coal worker’s pneumoconiosis disease in a mining group]. Zhonghua Lao Dong Wei Sheng Zhi ye bing Za Zhi = Zhonghua Laodong Weisheng Zhiyebing Zazhi = Chin J Industrial Hygiene Occup Dis. 2020;38(4):282–5.

    CAS  Google Scholar 

  34. Global burden of chronic respiratory diseases and risk factors, 1990–2019: an update from the global burden of disease study 2019. EClinicalMedicine. 2023;59:101936.

Download references

Funding

This research was supported by grants from the Ministry of Science and Technology of China to XZ (grant 2022YFC2304901), from the National Natural Science Foundation of China to GQ (grant 81370520), from the Science and Technology Support Program of Jiangsu Province, China to GQ (grant BE2015655), and from the Jiangsu Provincial Research Hospital to GQ (grant YJXYY202204-YSA03).

Author information

Author notes

  1. Hai-Lian Chen, Chun-Hu Li and Pei-Yao Zhai contributed equally to this work.

Authors and Affiliations

  1. Joint Division of Clinical Epidemiology, Affiliated Hospital of Nantong University, School of Public Health of Nantong University, Nantong, Jiangsu, China

    Hai-Lian Chen, Chun-Hu Li & Gang Qin

  2. Huangshi Center for Disease Control and Prevention, Huangshi, Hubei, China

    Hai-Lian Chen, Xue Qiao & Zu-Shu Qian

  3. Department of Infectious Diseases, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, Jiangsu, China

    Pei-Yao Zhai

  4. Department of Epidemiology and Biostatistics, School of Public Health of Nantong University, Nantong, Jiangsu, China

    Xun Zhuang & Yu-Long Lian

  5. National Key Clinical Construction Specialty-Department of Respiratory and Critical Care Medicine, Affiliated Hospital of Nantong University, Nantong, Jiangsu, China

    Jian Feng & Gang Qin

Authors
  1. Hai-Lian Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chun-Hu Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Pei-Yao Zhai
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xun Zhuang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yu-Long Lian
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xue Qiao
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jian Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Zu-Shu Qian
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Gang Qin
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Data curation: Hai-Lian Chen, Chun-Hu Li, Xue Qiao. Formal analysis: Chun-Hu Li, Zu-Shu Qian, Hai-Lian Chen and Xun Zhuang. Funding acquisition: Gang Qin, Xun Zhuang. Methodology: Hai-Lian Chen, Chun-Hu Li, Pei-Yao Zhai. Validation: Hai-Lian Chen, Yu-Long Lian. Visualization: Chun-Hu Li, Pei-Yao Zhai. Writing original draft: Hai-Lian Chen, Chun-Hu Li. Writing review & editing: Gang Qin and Jian Feng.

Corresponding authors

Correspondence to Jian Feng, Zu-Shu Qian or Gang Qin.

Ethics declarations

Ethics approval and consent to participate

The Ethics Committee of the Huangshi Center for Disease Control and Prevention has approved the study (Tdb2130). The Ethics Committee granted a waiver of the informed consent form. The rationale for the exemption is that the study is retrospective and analysis using anonymized data.

Consent for publication

Not applicable.

Competing interests

Gang Qin is an editorial board member of BMC Public Health. The remaining authors have no conflicts of interest to declare.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, HL., Li, CH., Zhai, PY. et al. Survival and disease burden analyses of occupational pneumoconiosis during 1958–2021 in Huangshi city, China: a retrospective cohort study. BMC Public Health 24, 1437 (2024). https://doi.org/10.1186/s12889-024-18847-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-024-18847-6

Keywords

BMC Public Health

ISSN: 1471-2458

Contact us