Discussion
We investigated the association between a municipal-level ADI and premature CVD mortality. More deprived regions are associated with an elevated mortality rate from premature CVD. When stratifying by urbanicity, we found that area deprivation is more associated with the outcomes in metropolitan regions than in other provinces. These associations persisted even after adjusting for the number of hospitals and area size.
Our findings are consistent with those of the existing literature, highlighting the association between regional socioeconomic conditions and CVD mortality despite variations in geographical units, deprivation indices and research designs across different studies. A cross-sectional study in the USA showed a significant association between ADI and cardiovascular premature mortality at the county level, while explaining 44% of county variability in cardiovascular mortality.31 Similarly, a UK cohort study followed up 3924 men aged 60–79 years for 12 years and conducted a multilevel analysis.32 After adjusting for individual socioeconomic status, smoking and drinking status, physical activity and body mass index, significantly higher odds of CVD mortality were observed in the more deprived areas compared with less deprived areas. In France, the ratio of SMR in the most deprived area over that in the least deprived area was 1.27 for all-cause deaths and 1.31 for CVD deaths.33 However, a Japanese study showed a weak negative correlation between municipal socioeconomic conditions and ischaemic heart disease mortality, possibly resulting from regional differences in urbanicity.34
In this context, we stratified the regions into metropolitan cities and other provinces to consider the distinct differences in area characteristics. Municipalities in metropolitan areas are smaller in size, are more densely populated and have easier access to public resources, including healthcare services, than those in other provinces. After stratification, the association between ADI and premature mortality was more substantial in metropolitan areas than in other provinces. Furthermore, the pronounced gradient in metropolitan cities remained evident after adjusting for other geographical factors influencing healthcare service accessibility, which are also key determinants of CVD. The significant association of area deprivation in urban settings can be attributed, in part, to increased exposure to adverse environments. Factors such as limited space for physical activity, increased accessibility to fast but low-quality food and elevated pollution levels may contribute to this phenomenon. For instance, one study found that regions with high unemployment rates were more often exposed to road traffic and that high traffic exposure was significantly associated with a higher risk of coronary calcification.35 These findings imply a greater risk for urban residents in deprived areas. Another study in Korea pointed out that the previous urban advantage in CVD mortality has diminished; instead, the gap between and within cities has worsened.36
Moreover, adjusting for the number of healthcare facilities and the area size weakened the association between ADI and premature CVD mortality in other provinces but slightly increased those in metropolitan areas. This suggests that tailored regional strategies should be implemented to decrease the CVD burden. While socioeconomic gradient is significantly associated with CVD mortality in metropolitan cities, the distribution and accessibility of medical resources can be critical in non-metropolitan areas where healthcare facilities are likely to be more sparsely distributed.37
The association between area deprivation and cardiovascular mortality is multifactorial and cannot be attributed to a single specific mechanism. For instance, the collective resource model elucidated that direct material deprivation at the regional level results in insufficient collective resources such as public services, recreational facilities, social support and healthcare services in deprived areas, thereby affecting residents’ health.38 Wilkinson and Pickett’s psychosocial model, on the other hand, emphasises that the perception of one’s relative social status, whether accurate or not, can lead to increased exposure to unfavourable environments and physiological manifestations of negative emotions, such as chronic stress and a sense of helplessness.39–41 Also, previous studies reported that the association between an unhealthy lifestyle and CVD mortality becomes more prominent with increasing levels of deprivation.42 43
This study had some limitations. First, temporal interpretations should be made with caution because of the cross-sectional design of this study. Although the outcome in this study was CVD death, and reverse causality was less likely, the cumulative effect of area deprivation could not be measured. Living in a deprived area for a certain period may have increased an individual’s risk of CVD mortality. However, an individual facing imminent death may have had substantial medical expenses and moved to a more deprived area. This may also have been reflected in the higher number of deaths in less affluent areas. Second, we combined several national datasets collected in 2020, including the demographic structure based on midyear population registration, deprivation index from the 2020 census data and death information from the 2020 national death statistics. However, each dataset may have been collected during different periods in 2020, possibly leading to the omission of information on population immigration and emigration. Given the global impact of COVID-19 on public health in 2020, it may also influence our results, despite the high participation rate in the Census survey and the relative robustness of death statistics. We performed the sensitivity analyses for all available years and found consistent results. Finally, there may have been residual confounding due to unknown confounders, although we attempted to account for other regional characteristics through stratification and adjustment.
Despite these limitations, our study has implications for further studies and prevention strategies by examining the association between area deprivation and premature CVD deaths. The ecological approach used in this study allows for the consideration of regional variations and psychosocial effects that cannot be measured at the individual income level. As Wilkinson noted, socioeconomic inequalities at the regional level are affected by the difference in the degree of residential segregation between the rich and poor and their perception of deprivation in the broader societal context.40 Furthermore, using the municipality as the unit of analysis is appropriate; the association was blurred when analysing at the provincial level in additional analyses. Aggregating the data into 17 provinces may have led to the omission of the social gradients in CVD mortality that were present at the municipal level. Another approach, a multilevel analysis confirmed a higher CVD mortality in deprived areas, as shown in the primary analysis.
Our findings indicate a significant association between area deprivation and premature CVD mortality rates, even after adjusting for geographical accessibility to healthcare facilities. Particularly, ADI was more correlated in metropolitan areas. This suggests that implementing CVD prevention strategies that reflect regional characteristics, with a focus on reducing socioeconomic gradient in metropolitan areas and increasing healthcare accessibility in other provinces, could potentially improve population health.