Is area-level socioeconomic deprivation associated with mortality due to circulatory system diseases in Poland?
BMC Public Health volume 23, Article number: 7 (2023)
Socioeconomic deprivation (SED) is known to influence cardiovascular health. However, studies analyzing the relationship between deprivation and circulatory system diseases (CSD) in Central and Eastern Europe are limited. This study aimed to assess the relationship between SED and mortality due to CSD at a population level in 66 sub-regions of Poland.
The 2010–2014 data regarding mortality and SED components were obtained from the Central Statistical Office. An area-based SED index was calculated based on the higher education rates, employment structure, wages, unemployment, and poverty. The dynamics of changes in mortality due to CSD was expressed by the number of deaths prevented or postponed (DPP) in terciles of the SED index. The associations between the mortality from CSD and SED index were analyzed using multivariate Poisson regression models and generalized estimating equations.
Among men, the percentage of DPP in 2014 was 13.1% for CSD, 23.4% for ischemic heart disease (IHD), and 21.4% for cerebrovascular diseases (CeVD). In the case of women, the proportion of DPP was 12.8, 25.6, and 21.6%, respectively. More deprived sub-regions experienced a greater decrease in CSD-related mortality than less deprived sub-regions. The disparity in mortality reduction between more deprived and less deprived sub-regions was even more pronounced for women. After adjusting for smoking prevalence, average BMI, population density, and changes in mortality over time, it was found that the SED index over the 2010–2014 time period was significantly associated with CSD- and IHD-related mortality for men (respectively 5.3 and 19.5% expected mortality increase per 1-unit increase of SED index), and with IHD- and CeVD-related mortality for women (respectively 30.3 and 23.0% expected mortality increase per 1-unit increase of SED index).
Significant differences in mortality changes due to CSD in Poland could be observed in relation to socioeconomic deprivation, resulting in reduced health inequalities. To reduce CSD mortality, more comprehensive preventive measures, including approaches addressing the socioeconomic factors, mainly poverty, education and employment, are needed, particularly in less urbanized areas.
Circulatory system diseases (CSD), which originate from atherosclerotic disease and include both the ischemic heart disease (IHD) and cerebrovascular diseases (CeVD), have been reported for many years the most common causes of death in the developed countries. In nonfatal cases, long-term treatment is required, which limits an individual’s capacity to work and consequently incurs a high cost in terms of both treatment and loss of professional productivity, particularly in those over 50 years of age [1, 2].
The burden of CSD is unequally distributed across Europe, with higher mortality in the central and eastern regions [3, 4]. This higher burden reported in post-communist countries has been somewhat inherited due to several historical reasons associated with the socioeconomic system of that time. Central and Eastern European countries are characterized by a lower gross domestic product per capita, and thereby, a lower expenditure on health care, and also limited access to health care and prevention services. Furthermore, the health care gap between the Eastern and Western populations may result from, or could be explained by, lifestyle and differences in socioeconomic influences . The difference in mortality due to CSD depends on socioeconomic position (SEP), which is determined mainly by education, type of occupation, and income. This indicates the need to include SEP in the risk assessment of CSD along with the classical factors: smoking, high blood pressure, hypercholesterolemia, overweight and obesity, sedentary lifestyle and diabetes [6, 7]. Low SEP may contribute to stress-related biological risk factors associated with the development of cardiovascular diseases (hypertension, proinflammatory cytokines, diabetes, obesity) .
Furthermore, it may lead to poorer health as a result of harmful health behaviors, such as tobacco smoking, an unhealthy diet, and lower physical activity [9, 10]. In addition to individual factors, CSD-related mortality could be linked with socioeconomic deprivation (SED) at the area level [11,12,13,14,15]. Deprivation is a state of observable and demonstrable disadvantages, which can be of social or material nature, compared to the local community or the wider society to which an individual, family, or group belongs . The area-related deprivation serves as a proxy for the SEP of the people who live in those regions, and characterizes the living environment which can influence an individual’s prospects for education, profession, income, or access to resources, and hence, health outcomes [7, 17]. Furthermore, neighborhood factors could have an impact on health through a lack of investments or inequality in resource redistribution in society. Moreover, chronic exposure to psychosocial stress, caused by poverty areas, is linked to the development of multisystemic adverse biological changes, i.e., allostatic burden, which results in cardiovascular diseases . Studies indicate that factors associated with area SED could significantly affect short-term survival, while individual factors may considerably impact long-term survival .
The majority of information on the relationship between SED and cardiovascular disease-related mortality comes from studies conducted in Western Europe and the United States [13,14,15, 19, 20]. At the same time, much fewer data are available from the Central and Eastern European countries [21, 22]. Recent studies conducted in Poland, Russia, and the Czech Republic have confirmed a strong association between the risk of death from CSD and psycho-socioeconomic factors [23, 24]. However, there is a need for a more extensive assessment of the factors determining health inequalities among people in the Central European countries, including Poland. The political transformation that occurred in 1989 in Poland led to significant changes in educational and occupational structures [25, 26]. Despite the considerable increase in SEP, some groups adapted to the new market economy at a slower rate. The initial post-transition period in the 1990s brought adverse consequences to those entering the labor market. In particular, young people with low levels of education, whose poorer material status contributed to poorer health outcomes, were less resistant to these changes . Furthermore, economic changes have clearly reduced employment in the heavy industry and national farm holdings [27, 28], leading to economic and social inequalities in the area of residence, i.e., impoverishment and social exclusion, and consequently excessive alcohol consumption due to changes in availability and affordability [25, 29]. These educational and economic inequalities that occurred in the early years of the transformation may now be affecting CSD mortality, including premature mortality under the age of 65 . Thus far, analyses have primarily focused on changes in CSD mortality caused by changes in exposure to some of classical risk factors (cholesterol level, blood pressure), lifestyle changes, and health care services improvement [31, 32]. No study has investigated the socioeconomic factors at the population level. Territorial variations in CSD mortality in Poland may likely be explained by SED differences and their changes over time.
This study aimed to assess the relationship between SED and mortality due to CSD at a population level in 66 sub-regions of Poland.
An ecological study was performed using the data on SED, lifestyle, urbanization, and mortality due to CSD collected from 66 administrative sub-regions of Poland. The sub-regions were defined according to the Nomenclature des Unités Territoriales Statistiques (NUTS-3) in 2006 . The mean size of the sub-regions was 4738 km2, with a wide range of variation between eight urban-type sub-regions (from 262 to 517 km2) and 58 remaining sub-regions (from 878 to 12,090 km2). Average number of population per sub-region (in year 2010) was 583,786, and ranged from 279,491 to 1,700,112. Descriptive statistics for the sub-regions were weighted to the population of sub-regions using weights as follows:
population x is the population of sub-region x,
population total is the total population of Poland, and.
n is the number of sub-regions.
The sum of weights for all sub-regions equals the number of sub-regions.
Annual data for mortality, all components of the SED index, and population density in 2010–2014 were obtained from the Central Statistical Office. In some analyses - where explicitly stated - the data were averaged over this period. Otherwise - like in longitudinal statistical models, separate data for each year were used. Information on average population BMI and proportion of smokers come from the panel Social Diagnosis Survey conducted repeatedly since 2000 on the same sample of households . For these variables, data from the 2011 survey was used and treated as static for the whole 5-year period between 2010 and 2014.
Data on deaths during 2010–2014 was used to analyze mortality. According to the International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10), the causes for the deaths were encoded as CSD (codes: I00–I99), including IHD (I20–I25) and CeVD (I60–I69).
The age-standardized number of deaths was used for statistical analyses. A direct standardization method was utilized for the analyses, considering the demographic structure in 5-year age groups (from 15 to 19 years to 80–84 years, and above 84 years old) separately for each of the 66 sub-regions of Poland . European Standard Population – 2013 edition was assumed as the standard population, and standardized mortality coefficients were determined accordingly . The coefficients of standardized mortality were calculated as 95% confidence intervals (CIs).
The dynamics of changes in mortality rates in 2010–2014 were assessed by calculating the deaths prevented or postponed (DPP) index . The index value was expressed in percentage and estimated as the difference between the expected number of deaths in 2014, assuming that the mortality rate in 2014 was the same as in 2010, and the actual number of deaths in 2014 relative to the expected number of deaths in 2014:
Area-level SED index
Five SED indicators were selected from the available variables in the system of routine statistical records from the Central Statistical Office. The selected indicators of SED were as follows: people with higher education, people employed in industry and construction, average monthly salary, unemployment rate, and people on social support due to poverty. These variables were chosen as they showed a strong correlation with all-cause mortality . According to published literature data, the SED variables are potential predictors of the geographical distribution of CSD mortality [39,40,41,42,43,44]. Education not only increases the awareness of beneficial health behaviors through medical and prophylactic care and treatment recommendations, but is also one of the essential determinants of employment [39, 40]. The professional class, characterized by a balance between effort and reward, is associated with health. Similarly, people employed in industry and construction are manual workers and exhibit worse health outcomes . Income determines an individual’s access to better-quality goods, services and housing to stay healthy . Thus, unemployment is linked with disease development due to reduced access to medical care and screening programs resulting from the lack of health insurance . Poverty associated with a lack of social support leads to chronic stress and adverse health behaviors [18, 44]. To calculate the SED index, its component variables were standardized via linear transformation such that the expected value of the variables was equal to 0 and the standard deviation was equal to 1. In addition, the sign was reversed for destimulants (proportion of people with higher education and average monthly salary). Finally, the SED index was assumed to be the arithmetic mean of the transformed components and calculated from the averaged values of component variables separately for the years 2010–2014. Therefore, the proposed SED index might be considered the mean of Z-Scores of its components, taking into account inverted scales for destimulant variables. The calculated index values were used to explain the differences in mortality across sub-regions of Poland. In addition, based on the 5-year averaged values of the SED index, the 66 sub-regions were divided into three groups with equal number of areas and with lower and upper tercile of averaged SED index as cut-off points:
the group with the lowest index values and thus the lowest deprivation level (less deprived),
the middle group,
the group with the highest index values (more deprived).
This classification was used to emphasize the inequalities in mortality for descriptive purposes.
Other explanatory variables
Tobacco consumption and body mass index (BMI) were calculated for the 66 sub-regions based on the results of the Social Diagnosis Survey conducted in 2011. The survey involved 12,386 households selected based on the two-stage stratified sampling method—first at the voivodeship level, followed by sampling based on residential location (large cities, small cities, and villages) . The study sample comprised 36,753 persons aged 16 years and over. The participation rate of the study was 72% (i.e., 26,453 persons were finally examined). To eliminate the socioeconomic differences in urbanization between the sub-regions in the model, population density was included as a standardizing factor [19, 45].
The relationship between mortality and the SED index was determined using the multivariate Poisson regression model for count data. As this model reflected the mortality coefficient per 100,000 residents rather than the number of deaths, the age-standardized number of deaths in individual sub-regions was taken as a dependent variable, with the number of people in sub-regions serving as the offset. The generalized linear models with Poisson distribution as probability distribution and logarithm as a link function were used for calculations.
Due to repeated measurements in the same statistical units (mortality in sub-regions in subsequent years during 2010–2014), generalized estimating equations (GEE) were used to obtain generalized linear models for the correlated data . An exchangeable structure was assumed for the working correlation matrix in GEE models.
The Poisson models were presented with exponentially transformed regression coefficients. Therefore, they can be interpreted as the expected relative change of the dependent variable (number of deaths standardized for the age) calculated for a one-unit increase of the independent variable, i.e., relative risk. Therefore, the reported coefficients can also be interpreted as a relative change in mortality coefficients, and not just death numbers. Furthermore, 95% CIs were calculated for the coefficients of the regression model and their p-values of corresponding Wald’s tests.
Variance Inflation Factors were calculated to check the presence of a potential multicollinearity problem in multivariable models. In no case did the VIF exceed the value of 4 taken as a threshold, and the maximum value observed was 3.84.
The regression models developed for the study included three causes of death (CSD, IHD, CeVD) for men and women. Furthermore, as independent variables, calendar year (for assessing the time effect on the dependent variable) and the aforementioned index - level of socioeconomic deprivation in the sub-regions - were included in the models. In addition, the models included lifestyle factors, such as mean BMI and the average percent of smokers in the sub-region (for the female and male population, respectively) and population density in the sub-region using logistic transformation at base 2.
Simple, univariate relationships between the SED index and mortality have been assessed with Spearman’s nonparametric correlation coefficients.
Statistical calculations were performed in the IBM® SPSS® Statistics for Windows statistical package, version 20.0 (IBM Corporation, Armonk, NY, USA). The level of statistical significance was assumed at α = 0.05.
During 2010–2014, a total of 868,418 deaths due to CSD were recorded in Poland, which included 405,235 deaths among men (46.7%). Of the overall death number, the deaths due to CSD in men accounted for 40.6%, while the death proportion in women reached 51.2%. Among CSD, IHD (32.0%) and CeVD (24.6%) were the leading causes of death.
The descriptive statistics for socioeconomic variables, lifestyle, and urbanization, by three sub-groups of Polish sub-regions (classified based on the SED index values), are presented in Table 1. The mean values of the index variables in sub-groups reflect their role as stimulants or destimulants—in the sub-regions characterized by the highest SED index (a more deprived), high values were attained by stimulants and low values by destimulants. Therefore, this group of sub-regions had the highest percentage of unemployment rate (18.4%), percentage of people who are on social support due to poverty (11.1%) and people employed in industry and construction (31.1.%), as well as the lowest percentage of people with university education (12.9%) and people with average monthly salary (3108 PLN). Opposite characteristics occurred in the sub-regions with the lowest SED index values (less deprived), which had the lowest unemployment rate (8.9%) and the lowest percentage of people on social support due to poverty (5.3%). On the other hand, these sub-regions had the highest percentage of people with a university education (22.5%) and the highest average monthly salary (3916 PLN). The precise distributions of SED components are shown on the maps presented in Supplemental Fig. 1.
A substantial difference in population density was noted between the three groups of sub-regions. In the more deprived sub-regions, it was on average 89 people/km2. In contrast in the less deprived sub-regions, it was over 13 times higher—1125 people/km2, which indicates that low SED index values are mainly linked with highly urbanized areas. The percentage of smoking men was similar in all the three groups of sub-regions, ranging from 32.3 to 35.8%, and in women from 20.9 to 22.4%, whereas a higher proportion of smoking men was found in the tercile representing those more deprived, while for women the highest proportion of smokers was observed in less deprived group. No significant differences were found in nutritional behaviors. The mean BMI of male and female populations was very similar in all three groups of sub-regions (26.4–26.6 kg/m2 for males and 25.0–25.4 kg/m2 for females).
The territorial distribution of sub-regions by SED index is presented in Fig. 1. The less deprived sub-regions with the lowest SED index value primarily had large urban agglomerations.
The highest mortality due to CSD was in the more deprived sub-regions and the middle group, and the lowest in those less deprived. At the same time, it can be observed that the dynamics of change in mortality expressed as a percentage of DPP varied among deprivation-defined tercile groups of sub-regions. The percentage of avoided deaths increased significantly with an increase in the SED index (p < 0.001). In the case of men, the IHD-related mortality was the lowest in the more deprived sub-regions. In addition, an increase in the percentage of DPP was observed in these sub-regions (p < 0.001). The relationship between the deaths due to CeVD in men and the SED index was unambiguous. The highest CeVD-related mortality was observed in the sub-regions characterized by an average level of socioeconomic deprivation. However, a statistically significant linear increase in the percentage of DPP was observed in these sub-regions (p < 0.01) (Table 2).
The relationship between the SED index and CSD-related mortality in women was similar to that in men. The age-standardized CSD mortality was the lowest in the sub-regions with the lowest deprivation, while the highest occurred in the middle group sub-regions. The percentage of prevented deaths increased significantly with the SED index value (p < 0.001). Similarly, mortality due to IHD was the lowest in the more deprived sub-regions and the highest (in 2014) in the less deprived sub-regions. The percentage of DPP increased significantly (p < 0.001) with an increase in the SED index value. Similar to men, mortality due to CeVD in women was the highest in middle sub-regions and the lowest in less deprived sub-regions. No significant differences were observed in the percentage of prevented deaths (p = 0.11). The differences were found at the same level in all the three groups of sub-regions studied (Table 3).
Univariate relationships between the SED index and age-standardized mortality in 66 sub-regions of Poland averaged over the 2010–2014 period are presented in Fig. 2. Statistically significant positive correlations were observed for CSD mortality both in males (r = 0.46) and females (r = 0.38) and also for CeVD in males (r = 0.31), however no significant correlations for IHD were observed.
As univariate analysis neglects other potentially important variables, the relationships between CSD mortality and time, SED, lifestyle, and population density are presented in Table 4. As the model is more detailed and adjusts for additional sources of variability, its results are not the same. The relationship between the SED and increased mortality was found for deaths due to CSD and IHD in men and deaths due to IHD and CeVD in women. For both women and men, a statistically significant decrease was observed in mortality due to CSD (2.9% for men and 2.8% for women annually), IHD (6.6% for men and 7.7% for women annually), and CeVD (6.1% for men and 6.6% for women annually). A higher mean BMI was associated with higher mortality due to CSD among women, but no significant relationship was observed between BMI and mortality due to IHD and CeVD. In the case of men, a relationship was noted between smoking and mortality due to CSD. Low SED index (by and large related to urban agglomerations) was found to be related to lower mortality. In comparison, a high population density was associated with greater IHD mortality in women and men and CeVD mortality only in women.
The study showed a considerable geographic variation in CSD mortality and SED in the 66 sub-regions of Poland. For both men and women, the largest mortality was found in the more deprived and middle group sub-regions, and the lowest in the less deprived sub-regions. However, the more deprived and middle group sub-regions benefited more from decreased CSD mortality in men and women. This indicates a reduction in health inequalities, which was also observed in other European countries [47, 48]. In women, the differences in the number of DPP between the SED groups were similar to that in men, except in the case of CeVD for which a nonsignificant trend was observed. The relationship between SED and CSD mortality was significant in men after adjusting for mean BMI, mean smoking rate, and population density. In women, a significant relationship was observed with mean BMI. However, the relationship between CSD mortality and SED was insignificant after adjusting for covariates, although the β coefficients were similar to that of men.
Interpretation of results
Our study indicates that the variation of CSD mortality in Poland’s sub-regions can be explained at least partially by the differences in SED. The results comply with the studies in which synthetic SED indexes were used to assess the relationship between deprivation and mortality due to CSD, IHD and CeVD [13,14,15, 19,20,21,22, 49]. It has been suggested that deprivation significantly influences the mortality due to CSD in men, but not women [13, 50]. We did not show a significant relationship in women, however the effect of SED is similar in size to that in model for men, which may indicate insufficient statistical power of the study.
While the impact of lifestyle on cardiovascular health is well known, it should be noted that there are gender differences in behavioral risk factors. For example, smoking prevalence traditionally affected more men than women . In our study, CSD mortality was shown to be associated with tobacco smoking only in men. This could reflect the higher smoking rates in the group of more deprived sub-regions compared with the less deprived sub-regions. According to the MORGAM study, in Poland, smoking is less prevalent among women with low and secondary education than women with higher education, particularly in less urbanized areas . A study indicated that in Eastern European countries, economic development and social and cultural processes associated with gender empowerment affect the differences in smoking between educated and uneducated women. It appears that the education related different patterns for smoking in Eastern European countries may interfere with the relationship between SED and mortality due to CSD in women, which may have influenced the lack of significance of the association [52, 53]. In our study, another important behavioral factor was BMI, which was associated with mortality due to CSD only in women. Although the prevalence of BMI varies by gender, obesity is more common in men below 45 years old, but at later ages above 45 years old, it predominates in women . Some studies suggest that the risk of developing heart disease in women is increased by obesity coexisting with some classical CSD factors, i.e. diabetes, hypertension and hypercholesterolemia. Compared to men, these factors predispose women to heart disease to a greater extent .
It is worth noting the results obtained in assessing the association between SED and IHD mortality. A significant relationship was not found in univariate analysis, but the inclusion of additional explanatory variables nevertheless showed the existence of such a relationship. This may be related to the problem of omitted-variable bias: in the analysis ignoring the additional variables, the univariable model necessarily imputed their effect to the SED scale. By adding the omitted variables, the SED scale no longer assumes their partial effect, but reflects its “true” effect, which turns out to be statistically significant. Furthermore, this study confirmed the relationship between SED and mortality from IHD, with increased deprivation more strongly increasing mortality in women than in men (by 30.3% vs. 19.5% per SED scale unit). This significant result in both sexes may indicate similar risk factors for IHD, which may be explained by the fact that deprivation co-occurs with stress-related biological risk factors such as hypertension and diabetes [8, 18]. However, for the higher score in women, significant hormonal changes during menopause may play a role, which may also aggravate preexisting risk factors . As in other studies [57, 58], we found the relationship between SED and mortality due to CeVD only in women. An explanation for this meaningful relationship may be psychosocial factors, which play a more significant role in women; for instance, depression may mediate the relationship between low SED and stroke . However, the unclear result in the male population may be associated with a lower risk of stroke among unskilled manual workers compared with high-grade civil servants and executives .
Our results showed inequalities in mortality across sub-regions ranked by deprivation. The lowest CSD mortality was found in the less deprived sub-regions, which predominantly include large urban agglomerations and small cities. This observation may be associated with the fact that these sub-regions have access to different resources, including educational infrastructure, services, and job opportunities . Furthermore, these sub-regions offer a more favorable environment, such as access to gyms and shops selling healthy foods, as well as health care services, which may contribute to better health outcomes in their residents [45, 60, 61]. The more deprived sub-regions are located in the eastern and northwest parts of Poland and are characterized by a low population density. These are considered less attractive to investors, which affects economic development. This confirms, among others, the fact that the areas of many sub-regions with more deprived sub-regions overlap with those areas in which national farm holdings had been liquidated and more harmful effects of the economic transformation are experienced. These areas lack support programs targeting social groups deprived of earlier forms of employment, and are consequently linked with poor health outcomes [26, 27].
Furthermore, the area encompassing sub-regions with a high deprivation in the western part of the country was associated with a more clear trend toward job migrations to work legally abroad. The migrations may exacerbate the optimal development of these areas due to the outflow of human capital and are also associated with a prevalence of social exclusion and poverty, leading altogether to health hazards [26, 62]. Some analyses indicate that the eastern and northwestern parts of the country have high total mortality due to CSD of the age group 25–64 years . The issue of worse cardiovascular health in areas where national farm holdings were eliminated and areas with intensive labor migration is an interesting topic that needs further research.
For IHD, mortality was the highest in the less deprived sub-region group. In the case of CeVD, the highest mortality was found in middle sub-regions, which may suggest the occurrence of other specific factors that are not included in the analysis . Indeed, for both IHD mortality and CeVD, positive relationship with the SED index has been confirmed using multivariable models that account for additional sources of mortality variation. It should also be emphasized that a rapid decrease in CSD mortality was noted in more deprived sub-regions with high SED index values. In the more deprived and middle group sub-regions, the percentage of DPP was higher compared to the less deprived sub-regions.
While the mortality due to CSD and CeVD was the lowest in the less deprived sub-regions, the reduction of mortality between 2010 and 2014 in these sub-regions was smaller compared to the sub-regions with the high SED index. This could be explained by the fact that in the best-developed—highly urbanized—sub-regions, mortality (due to better availability of cardiologic care services and invasive cardiology procedures) had been reduced in earlier years and the scope for improvement of cardiovascular health was narrower . Simultaneously, modern prevention and therapeutic methods are increasingly becoming accessible in more deprived and middle sub-regions. For example, life-saving invasive cardiology procedures have been available after 2000 not only in academic centers but also in district hospitals, which might have resulted in reduced mortality even in smaller centers with a higher SED index. This is also supported by the results of other analyses , according to which the reduction of mortality due to CSD among people with higher education was particularly pronounced between 1991 and 1993 and 2001–2003, whereas during the period 2001–2003 and 2010–2012, the reduction was considerably lower. Noticeably, this mortality reduction was observed at a lower rate in people with low education compared to those with higher education.
Our results have shown that mortality from CSD, IHD and CeVD has decreased, and similar trends prevail in most European Union (EU) countries [1, 4, 5, 63]. It is also worth noting the progressive improvement in the health status of the Polish population, as evidenced by the systematically decreasing mortality rate of both younger and older people [63, 66]. At the same time, however, it must be said, against the background of the EU countries in the analyzed period 2010–2014, the mortality rates due to CSD in Poland were 60% higher compared to the average level for the EU, but already comparing to the old EU-15 countries were higher by 90%. For example, the level of mortality due to CSD in Poland is three times higher compared to France and the Netherlands. In contrast, against the background of the EU-13 countries that joined the EU after 2003, the level of mortality due to CSD in Poland is relatively low, i.e. the rates are 15% lower. Among Poland’s neighboring countries, the level of mortality due to CSD in Poland was 50% higher compared to Germany and similar to the level of mortality in the Czech Republic and Slovakia . The existence of significant disparities in the level of mortality due to CSD between Poland and EU countries, especially the old EU-15, may indicate that there are opportunities for further mortality reduction. The stratification using terciles allowed to identify the residence-related deprivation, which may be necessary at intervention strategies and activities addressing improving cardiovascular health. In the policy process aimed to alleviate health inequalities, the activities should focus not only on poverty reduction throughout income redistribution. It is also critical ensuring equality of health opportunity in the entire population, through education, employment, improved working conditions and preventive care . Decreasing health inequalities ought to be a political and social priority, given that they deteriorate economic productivity and the potential for sustainable and inclusive growth.
Strengths and limitations
To our knowledge, this study is the first to use the SED index to assess the relationship with mortality due to CSD at the population level in Poland. The study was performed considering the whole population of the country. In the 66 sub-regions, a large variation was found in CSD mortality and SED. Furthermore, the sub-regions represented all the characteristics that are typical for Poland. The synthetic SED index enabled an approximate estimation of the singular variables (education, structure in employment, salary, unemployment, and poverty). A database concerning sub-regions defined based on NUTS-3 classification, which is used in the EU member states, was used for the first time in this study . A unique strength of this study was the comparison of mortality due to CSD and the time-related changes using the DPP index in three different environments regarding deprivation level, thus, contrary to other studies reporting mortality trends in administrative areas of the country . Our results showed a potent effect of health inequalities. Therefore, they may contribute to limited literature field dealing with the associations between area-related deprivation and mortality from CSD in Central and Eastern European countries [21, 22]. Similar and comparable socioeconomic levels in specific countries of this region may facilitate cross-comparisons and allow consistently conducting research across the populations.
The results of the study should be interpreted in light of certain limitations. The ecological design does not allow addressing the causality of the relationships. As sub-regions were considered as statistical units in this study, instead of individual persons, it was possible to investigate the inequalities between them, while inequalities within them remain unexplored. Epidemiological analyses for geographical areas lead to the best results when statistical units are populations of small size [69, 70]. It allows for better homogeneity and decreases the problem of averaging of different populations within a single geographical area. Such averaging results in attenuation of studied effects in statistical models. In this paper, population sizes in sub-regions were relatively large, so it was expectable that real effects might be significantly stronger and so their detection could be difficult, if not impossible. However, we were able to confirm statistically significant associations of the SED index with mortality from each of the three causes of death analyzed, although not always for both sexes. One explanation for the difficulty in demonstrating these relationships in each case may be the aforementioned high level data aggregation at the sub-regional level and the blurring of the relationships studied. In contrast to mortality and SED, information on covariates such as education, smoking, and BMI was based on one-point observation (Census 2011, Social Diagnosis Survey 2011), which was the only available for the studied sub-regions within the study period. It is unlikely that these characteristics changed much within the observation time. The relationship between SED and mortality due to CSD could be confounded by the sub-region differences in the exposure to the other uncontrolled factors; therefore, residual confounding should be considered. For example, we did not utilize stress, an important variable, because of the limited accessibility of such data. However, stress is regarded as a significant mediator linking the associations between deprivation and mortality due to cardiovascular diseases . Noticeably, some components of the SED (e.g. percent of people on social support due to poverty) may have, at least partly, reflected chronic stress . Thus, the role of stress should be investigated in future research. Another limitation of the study is the quality of the data on deaths associated with the differences in the reliability of death-cause coding by physicians who filled out death certificates. Our results showed that the highest mortality due to IHD was observed in the less deprived (most urbanized) sub-regions. This may be partially explained by the discrepancies in death certification encoding described earlier in Poland . Such an issue was even reported in countries with highly advanced health information systems . However, territorial differences in death-cause coding could rather contribute to the greater impact of random variability on our results, and it is less likely that the occurrence of a systematic error would explain the observed relationships. The latter is supported by the overall consistency of relationships found for IHD, CeVD, and CSD.
Significant differences in mortality changes due to CSD in Poland could be observed in relation to socioeconomic deprivation, resulting in reduced health inequality. To reduce CSD mortality, more comprehensive preventive measures, including approaches addressing the socioeconomic factors, mainly poverty, education and employment, are needed, particularly in less urbanized areas. Further, well-designed studies should consider the contribution of psychosocial factors and their effect on cardiovascular health.
Availability of data and materials
The datasets analyzed during the current studies are publicly available in Central Statistical Office, Poland (https://stat.gov.pl).
body mass index
circulatory system diseases
deaths prevented or postponed
generalized estimating eqs.
International Statistical Classification of Diseases and Related Health Problems
ischemic heart disease
Nomenclature des Unités Territoriales Statistiques
Gańczak M, Miazgowski T, Kożybska M, Kotwas A, Korzeń M, Rudnicki B, et al. Changes in disease burden in Poland between 1990–2017 in comparison with other central European countries: a systematic analysis for the global burden of disease study 2017. PLoS One. 2020;15(3):e0226766. https://doi.org/10.1371/journal.pone.0226766.
Genowska A, Fryc J, Pinkas J, Jamiołkowski J, Szafraniec K, Szpak A, et al. Social costs of loss in productivity-related absenteeism in Poland. Int J Occup Med Environ Health. 2017;30(6):917–32 https://doi.org/10.13075/ijomeh.1896.01123.
Naghavi M, Abajobir AA, Abbafati C, Abbas KM, Abd-Allah F, Abera SF, et al. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980-2016: a systematic analysis for the global burden of disease study 2016. Lancet. 2017;390(10100):1151–210. https://doi.org/10.1016/S0140-6736(17)32152-9.
Nichols M, Townsend N, Scarborough P, Rayner M. Trends in age-specific coronary heart disease mortality in the European Union over three decades: 1980–2009. Eur Heart J. 2013;34(39):3017–27. https://doi.org/10.1093/eurheartj/eht159.
Hartley A, Marshall DC, Salciccioli JD, Sikkel MB, Maruthappu M, Shalhoub J. Trends in mortality from ischemic heart disease and cerebrovascular disease in Europe: 1980 to 2009. Circulation. 2016;133(20):1916–26 https://doi.org.10.1161/CIRCULATIONAHA.115.018931.
Kaplan GA, Keil JE. Socioeconomic factors and cardiovascular disease: a review of the literature. Circulation. 1993;88(Pt 1):1973–98. https://doi.org/10.1161/01.cir.88.4.1973.
Galobardes B, Shaw M, Lawlor DA, Lynch JW, Davey SG. Indicators of socioeconomic position (part 2). J Epidemiol Community Health. 2006;60(2):95–101. https://doi.org/10.1136/jech.2004.028092.
Steptoe A, Hiltl TJ, Dowd JB, Hamer M. Socioeconomic status and central adiposity as determinants of stress-related biological responses relevant to cardiovascular disease risk. Brain Behav Immun. 2019;77:16–24. https://doi.org/10.1016/j.bbi.2018.11.019.
Zujko ME, Waśkiewicz A, Drygas W, Cicha-Mikołajczyk A, Zujko K, Szcześniewska D, et al. Dietary habits and dietary antioxidant intake are related to socioeconomic status in polish adults: a nationwide study. Nutrients. 2020;12(2):pii: E518. https://doi.org/10.3390/nu12020518.
Algren MH, Ekholm O, Nielsen L, Ersbøll AK, Bak CK, Andersen PT. Associations between perceived stress, socioeconomic status, and health-risk behaviour in deprived neighbourhoods in Denmark: a cross-sectional study. BMC Public Health. 2018;18(1):250. https://doi.org/10.1186/s12889-018-5170-x.
Diez Roux AV, Merkin SS, Arnett D, Chambless L, Massing M, Nieto FJ, et al. Neighborhood of residence and incidence of coronary heart disease. N Engl J Med. 2001;345(2):99–106. https://doi.org/10.1056/NEJM200107123450205.
Coady SA, Johnson NJ, Hakes JK, Sorlie PD. Individual education, area income, and mortality and recurrence of myocardial infarction in a Medicare cohort: the National Longitudinal Mortality Study. BMC Public Health. 2014;14:705. https://doi.org/10.1186/1471-2458-14-705.
Machón M, Aldasoro E, Martínez-Camblor P, Calvo M, Basterretxea M, Audicana C, et al. Socioeconomic differences in incidence and relative survival after a first acute myocardial infarction in the Basque Country. Spain Gac Sanit. 2012;26(1):16–23. https://doi.org/10.1016/j.gaceta.2011.06.013.
Ramsay SE, Morris RW, Whincup PH, Subramanian SV, Papacosta AO, Lennon LT, et al. The influence of neighbourhood-level socioeconomic deprivation on cardiovascular disease mortality in older age: longitudinal multilevel analyses from a cohort of older British men. J Epidemiol Community Health. 2015;69(12):1224–31. https://doi.org/10.1136/jech-2015-205542.
Moissl AP, Delgado GE, Krämer BK, März W, Kleber ME, Grammer TB. Area-based socioeconomic status and mortality: the Ludwigshafen risk and cardiovascular health study. Clin Res Cardiol. 2020;109(1):103–14. https://doi.org/10.1007/s00392-019-01494-y.
Townsend P. Deprivation. J. Soc Policy. 1987;16:125–46. https://doi.org/10.1017/S0047279400020341.
Pampalon R, Raymond G. A deprivation index for health and welfare planning in Quebec. Chronic Dis Canada. 2000;21(3):104–13.
Robinette JW, Charles ST, Almeida DM, Gruenewald TL. Neighborhood features and physiological risk: an examination of allostatic load. Health Place. 2016;41:110–8. https://doi.org/10.1016/j.healthplace.2016.08.003.
Rey G, Jougla E, Fouillet A, Hémon D. Ecological association between a deprivation index and mortality in France over the period 1997–2001: variations with spatial scale, degree of urbanicity, age, gender and cause of death. BMC Public Health 2009 ;9:33. https://doi.org/10.1186/1471-2458-9-33..
Singh GK, Siahpush M, Azuine RE, Williams SD. Increasing area deprivation and socioeconomic inequalities in heart disease, stroke, and cardiovascular disease mortality among working age populations, United States, 1969-2011. Int J MCH AIDS. 2015;3(2):119–33.
Juhász A, Nagy C, Páldy A, Beale L. Development of a deprivation index and its relation to premature mortality due to diseases of the circulatory system in Hungary, 1998–2004. Soc Sci Med. 2010;70(9):1342–9. https://doi.org/10.1016/j.socscimed.2010.01.024.
Slachtová H, Tomásková H, Splíchalová A, Polaufová P, Fejtková P. Czech socio-economic deprivation index and its correlation with mortality data. Int J Public Health. 2009;54(4):267–73. https://doi.org/10.1007/s00038-008-7092-3.
Kozela M, Pająk A, Micek A, Besala A, Kubinova R, Malyutina S, et al. Impact of perceived control on all-cause and cardiovascular disease mortality in three urban populations of central and Eastern Europe: the HAPIEE study. J Epidemiol Community Health. 2017;71(8):771–8. https://doi.org/10.1136/jech-2017-208992.
Tillmann T, Pikhart H, Peasey A, Kubinova R, Pajak A, Tamosiunas A, et al. Psychosocial and socioeconomic determinants of cardiovascular mortality in Eastern Europe: a multicentre prospective cohort study. PLoS Med 2017;14(12):e1002459. https://doi.org/https://doi.org/10.1371/journal.pmed.1002459.
Zatoński WA. Epidemiological analysis of health situation development in Europe and its causes until 1990. Ann Agric Environ Med. 2011;18(2):194–202.
Beblo M, Golinowska S, Laurer C, Pietka K, Sowa A. Poverty dynamics in Poland. Selected quantitative analyses. Warsaw: Center for Social and Economic Research; 2002.
Walsh D, Taulbut M, Hanlon P. The aftershock of deindustrialization - trends in mortality in Scotland and other parts of post-industrial Europe. Eur J Pub Health. 2010;20(1):58–64. https://doi.org/10.1093/eurpub/ckp063.
Halamska M. Transformacja wsi: 1989-2009. Zmienny rytm modernizacji. Studia Regionalne i Lokalne. 2011;44(2):5–25.
Wojtyniak B, Moskalewicz J, Stokwiszewski J, Rabczenko D. Gender-specific mortality associated with alcohol consumption in Poland in transition. Addiction. 2005;100(12):1779–89. https://doi.org/10.1111/j.1360-0443.2005.01247.x.
Leinsalu M, Stirbu I, Vågerö D, Kalediene R, Kovacs K, Wojtyniak B, et al. Educational inequalities in mortality in four eastern European countries: divergence in trends during the post-communist transition from 1990 to 2000. Int J Epidemiol. 2009;38(2):512–25. https://doi.org/10.1093/ije/dyn248.
Zatonski WA, McMichael AJ, Powles JW. Ecological study of reasons for sharp decline in mortality from ischaemic heart disease in Poland since 1991. BMJ. 1998;316(7137):1047–51. https://doi.org/10.1136/bmj.316.7137.1047.
Bandosz P, O'Flaherty M, Drygas W, Rutkowski M, Koziarek J, Wyrzykowski B, et al. Decline in mortality from coronary heart disease in Poland after socioeconomic transformation: modelling study. BMJ. 2012;344:d8136. https://doi.org/10.1136/bmj.d8136.
European Commission. Regions in the European Union. Nomenclature of territorial units for statistics. NUTS 2006/EU-27. Luxembourg: 2007. https://ec.europa.eu/eurostat/web/products-manuals-and-guidelines/-/ks-ra-07-020. Accessed 15 Nov 2020.
The Council for Social Monitoring. Social Diagnosis Report. Warsaw; 2011. http://www.diagnoza.com/index-en.html. Accessed 15 Nov 2020.
Morris JA, Gardner MJ. Calculating confidence intervals for relative risks (odds ratios) and standardized ratios and rates. BMJ. 1988;296(6632):1313–6. https://doi.org/10.1136/bmj.296.6632.1313.
Eurostat. Revision of the European Standard Population. Report of Eurostat’s task force. Luxembourg; 2013. https://ec.europa.eu/eurostat/web/products-manuals-and-guidelines/-/ks-ra-13-028. Accessed 15 Mar 2022.
Hotchkiss JW, Davies CA, Dundas R, Hawkins N, Jhund PS, Scholes S, et al. Explaining trends in Scottish coronary heart disease mortality between 2000 and 2010 using IMPACTSEC model: retrospective analysis using routine data. BMJ. 2014;348:g1088. https://doi.org/10.1136/bmj.g1088.
Polak M, Genowska A, Szafraniec K, Fryc J, Jamiołkowski J, Pająk A. Area-based socio-economic inequalities in mortality from lung cancer and respiratory diseases. Int J Environ Res Public Health. 2019;16(10):pii: E1791. https://doi.org/10.3390/ijerph16101791.
Veronesi G, Ferrario MM, Kuulasmaa K, Bobak M, Chambless LE, Salomaa V, et al. Educational class inequalities in the incidence of coronary heart disease in Europe. Heart. 2016;102(12):958–65. https://doi.org/10.1136/heartjnl-2015-308909.
Ricceri F, Sacerdote C, Giraudo MT, Fasanelli F, Lenzo G, Galli M, et al. The association between educational level and cardiovascular and cerebrovascular diseases within the EPICOR study: new evidence for an old inequality problem. PLoS One. 2016;11(10):e0164130. https://doi.org/10.1371/journal.pone.0164130.
Bertuccio P, Alicandro G, Sebastiani G, Zengarini N, Costa G, La Vecchia C, et al. Mortality by occupation-based social class in Italy from 2012 to 2014. Int J Public Health. 2018;63(7):865–74. https://doi.org/10.1007/s00038-018-1149-8.
Kawachi I, Adler NE, Dow WH. Money, schooling, and health: mechanisms and causal evidence. Ann N Y Acad Sci. 2010;1186:56–68. https://doi.org/10.1111/j.1749-6632.2009.05340.x.
Meneton P, Kesse-Guyot E, Méjean C, Fezeu L, Galan P, Hercberg S, et al. Unemployment is associated with high cardiovascular event rate and increased all-cause mortality in middle-aged socially privileged individuals. Int Arch Occup Environ Health. 2015;88(6):707–16. https://doi.org/10.1007/s00420-014-0997-7.
Al-Turk B, Harris C, Nelson G, Smotherman C, Palacio C, House J. Poverty, a risk factor overlooked: a cross sectional cohort study comparing poverty rate and cardiovascular disease outcomes in the state of Florida. J Investig Med. 2018;66(3):693–5. https://doi.org/10.1136/jim-2017-000621.
Borrell C, Pons-Vigués M, Morrison J, Díez E. Factors and processes influencing health inequalities in urban areas. J Epidemiol Community Health. 2013;67(5):389–91. https://doi.org/10.1136/jech-2012-202014.
Zeger S, Liang K. Longitudinal data analysis for discrete and continuous outcomes. Biometrics. 1986;42(1):121–30.
Di Girolamo C, Nusselder WJ, Bopp M, Brønnum-Hansen H, Costa G, Kovács K, et al. Progress in reducing inequalities in cardiovascular disease mortality in Europe. Heart. 2020;106(1):40–9 https://doi.org.10.1136/heartjnl-2019-315129.
Leyland AH, Dundas R. Declining cardiovascular mortality masks unpalatable inequalities. Heart. 2020;106(1):6–7 https://doi.org.10.1136/heartjnl-2019-315708.
Grau AJ, Dienlin S, Bartig D, Maier W, Buggle F, Becher H. Regional deprivation, stroke incidence, and stroke care—an analysis of billing and quality assurance data from the German state of Rhineland-Palatinate. Dtsch Arztebl Int. 2021;118(23):397–402. https://doi.org/10.3238/arztebl.m2021.0149.
Kondo N, Saito M, Hikichi H, Aida J, Ojima T, Kondo K, et al. Relative deprivation in income and mortality by leading causes among older Japanese men and women: AGES cohort study. J Epidemiol Community Health. 2015;69(7):680–5. https://doi.org/10.1136/jech-2014-205103.
Peters SA, Van der Schouw YT, Woodward M, Huxley RR. Sex differences in smoking-related risk of vascular disease and all-cause mortality. Curr Cardiovasc Risk Rep. 2015;7:473–9. https://doi.org/10.1007/s12170-013-0339-4.
Ferrario MM, Veronesi G, Chambless LE, Tunstall-Pedoe H, Kuulasmaa K, Salomaa V, et al. The contribution of educational class in improving accuracy of cardiovascular risk prediction across European regions: the MORGAM project cohort component. Heart. 2014;100(15):1179–87. https://doi.org/10.1136/heartjnl-2013-304664.
Schaap MM, Kunst AE, Leinsalu M, Regidor E, Espelt A, Ekholm O, et al. Female ever-smoking, education, emancipation and economic development in 19 European countries. Soc Sci Med. 2009;68(7):1271–8. https://doi.org/10.1016/j.socscimed.2009.01.007.
Kenchaiah S, Gaziano JM, Vasan RS. Impact of obesity on the risk of heart failure and survival after the onset of heart failure. Med Clin North Am. 2004;88(5):1273–94. https://doi.org/10.1016/j.mcna.2004.04.011.
Broni EK, Ndumele CE, Echouffo-Tcheugui JB, Kalyani RR, Bennett WL, Michos ED. The diabetes-cardiovascular connection in women: understanding the known risks, outcomes, and implications for care. Curr Diab Rep. 2022;22(1):11–25. https://doi.org/10.1007/s11892-021-01444-x.
Gheisari F, Emami M, Raeisi Shahraki H, Samipour S, Nematollahi P. The role of gender in the importance of risk factors for coronary artery disease. Cardiol Res Pract. 2020;2020:6527820. https://doi.org/10.1155/2020/6527820.
Yang L, Korhonen K, Moustgaard H, Silventoinen K, Martikainen P. Pre-existing depression predicts survival in cardiovascular disease and cancer. J Epidemiol Community Health. 2018;72(7):617–22 https://doi.10.1136/jech-2017-210206.
Novak M, Torén K, Lappas G, Kok WG, Jern C, Wilhelmsen L, et al. Occupational status and incidences of ischemic and hemorrhagic stroke in Swedish men: a population-based 35-year prospective follow-up study. Eur J Epidemiol. 2013;28(8):697–704. https://doi.org/10.1007/s10654-013-9833-8.
Pujades-Rodriguez M, Timmis A, Stogiannis D, Rapsomaniki E, Denaxas S, Shah A, et al. Socioeconomic deprivation and the incidence of 12 cardiovascular diseases in 1.9 million women and men: implications for risk prediction and prevention. PLoS One. 2014;9(8):e104671. https://doi.org/10.1371/journal.pone.0104671.
Koopman C, Vaartjes I, van Dis I, Verschuren WM, Engelfriet P, Heintjes EM, et al. Explaining the decline in coronary heart disease mortality in the Netherlands between 1997 and 2007. PLoS One. 2016;11(12):e0166139. https://doi.org/10.1371/journal.pone.0166139.
Lagström H, Halonen JI, Kawachi I, Stenholm S, Pentti J, Suominen S, et al. Neighborhood socioeconomic status and adherence to dietary recommendations among Finnish adults: a retrospective follow-up study. Health Place. 2019;55:43–50. https://doi.org/10.1016/j.healthplace.2018.10.007.
Rethy LB, McCabe ME, Kershaw KN, Ahmad FS, Lagu T, Pool LR, et al. Neighborhood poverty and incident heart failure: an analysis of electronic health records from 2005 to 2018. J Gen Intern Med. 2021;36(12):3719–27 https://doi.org.10.1007/s11606-021-06785-7.
Wojtyniak B, Goryński P. red. Warsaw: Health situation of the Polish population and its determinants. NIPH-NRI; 2020. [In Polish]
National Health Fund. Ischemic stroke of the brain. Warsaw: NHF; 2019. [In Polish]
Strzelecki Z, Szymborski J. Zachorowalność i umieralność na choroby układu krążenia a sytuacja demograficzna Polski. Warszawa: Rządowa Rada Ludnościowa; 2015.
Genowska A, Zalewska M, Jamiołkowski J, Stepaniak U, Szpak A, Maciorkowska E, et al. Inequalities in mortality of infants under one year of age according to foetal causes and maternal age in rural and urban areas in Poland, 2004-2013. Ann Agric Environ Med. 2016;23(2):285–91. https://doi.org/10.5604/12321966.1203892.
Eurostat: Causes of death statistics. https://ec.europa.eu/eurostat/statistics-explained/index.php?title Accessed 15 Sep 2022.
Marmot M. The health gap: the challenge of an unequal world: the argument. Int J Epidemiol. 2017;46(4):1312–8 https://doi.org.10.1093/ije/dyx163.
Allik M, Leyland A, Travassos Ichihara MY, Dundas R. Creating small-area deprivation indices: a guide for stages and options. J Epidemiol Community Health. 2020;74(1):20–5 https://doi.org.10.1136/jech-2019-213255.
Krieger N, Chen JT, Waterman PD, Soobader MJ, Subramanian SV, Carson R. Geocoding and monitoring of US socioeconomic inequalities in mortality and cancer incidence: does the choice of area-based measure and geographic level matter?: the public health disparities geocoding project. Am J Epidemiol. 2002;156(5):471–82 https://doi.org.10.1093/aje/kwf068.
Tawakol A, Osborne MT, Wang Y, et al. Stress-associated neurobiological pathway linking socioeconomic disparities to cardiovascular disease. J Am Coll Cardiol. 2019;73(25):3243–55 https://doi.org.10.1016/j.jacc.2019.04.042.
Wojtyniak B, Jankowski K, Zdrojewski T, Opolski G. Regional differences in determining cardiovascular diseases as the cause of death in Poland: time for change. Kardiol Pol. 2012;70(7):695–701.
Mikkelsen L, Iburg KM, Adair T, Fürst T, Hegnauer M, von der Lippe E, et al. Assessing the quality of cause of death data in six high-income countries: Australia, Canada, Denmark, Germany, Japan and Switzerland. Int J Public Health. 2020;65(1):17–28. https://doi.org/10.1007/s00038-019-01325-x.
The study was financed by the Polish National Science Centre (grant number 2016/23/D/HS4/01887). Author who received the funding: AG. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The data used in this study are aggregated secondary data with no personal identifying information that can be linked to study participants. Therefore no ethics approval was required for this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jamiołkowski, J., Genowska, A. & Pająk, A. Is area-level socioeconomic deprivation associated with mortality due to circulatory system diseases in Poland?. BMC Public Health 23, 7 (2023). https://doi.org/10.1186/s12889-022-14914-y