Abstract
Introduction Ensuring health equity, especially for vulnerable populations in less developed settings with poor health system is essential for the current and future global health threats. This study examined the heterogeneity of COVID-19 mortality and its association with population health characteristics, health care capacity in responding pandemic, and socio-economic characteristics across 514 districts in Indonesia.
Methods This nationwide ecological study included aggregated COVID-19 cases data from all 514 districts in Indonesia, recorded in the National COVID-19 Task Force database, during the first two years of the epidemic, from 1 March 2020 to 27 February 2022. We calculated incidence and mortality rate by time, sex, and age. We extracted district-level socio-demographics, population health, and health care capacity data from government official sources. We used multivariable linear regression to examine factors associated with higher mortality rate.
Results Of total 5,539,333 reported cases, 148,034 (2·7%) died, and 5,391,299 (97.4%) were recovered. The national mortality rate was 55 per 100,000 population, ranged from 13 per 100,000 population in Papua to 156 per 100,000 population in East Kalimantan province. At district-level, higher mortality rate was associated with higher COVID-19 incidence (p<0.0001), higher proportion of ≥60 years old population (p<0.0001), higher prevalence of diabetes mellitus (p<0.0001), lower prevalence of obesity (p=0.019), lower number of doctors per population (p=0.001), higher life expectancy at birth (p=0.035), and lower formal education (p=0.021). There was no association between COVID-19 mortality rate with expenditure, prevalence of hypertension and pneumonia, vaccine coverage for ≥60 years old population, number of nurses, midwives, and hospitals per population (p>0.05 each).
Conclusion COVID-19 mortality rate in Indonesia was highly heterogeneous and associated with different prevalence of pre-existing comorbidity, healthcare capacity in responding the pandemic, and socio-economic characteristics. This study revealed the need of health capacity strengthening and better resource allocation to ensure optimal health outcomes for vulnerable population.
What is already known on this topic
The severity of COVID-19 illness and clinical outcomes can be affected by the concentration of comorbidities in susceptible groups in communities, and through disparities of access to health care for preventive measures or prompt diagnosis and treatment.
However, evidence on the heterogeneity of COVID-19 impact from low- and middle-income country (LMIC) where differences in age distribution, comorbidities, access to quality health services, and other factors, may greatly influence mortality risk, are limited.
What this study adds
This study affirmed that COVID-19 disproportionately affected areas with high proportion of elder population, high prevalence of diabetes mellitus, lower doctor to population ratio, higher life expectancy at birth, and lower level of formal education.
These findings indicate that vulnerability to death associated with COVID-19 in LMIC includes not only elder and comorbid, but also males and communities living in area with lower health care capacity and with lower level of education.
How this study might affect research, practice and/or policy
These findings may inform decisions on health resource allocation against COVID-19 delivering the greatest possible health dividends by prioritising interventions, including even distribution of essential health care need such as doctors, and a tailored risk communication and community engagement for the most vulnerable communities in LMIC, especially with decentralised health systems like in Indonesia.
Background
Globally, coronavirus disease 2019 (COVID-19) pandemic has caused millions of deaths [1], public health systems collapse [2], and socio-economic disruptions [3,4]. As of 17 June 2022, there have been 535,863,950 total confirmed cases and 6,314,972 deaths reported to the World Health Organization (WHO) [1]. It was reported that employment loss rates have been high since the beginning of the pandemic, and significantly higher among women than men (26% vs 20% in September 2021), with a steadily decreasing trend worldwide [3]. Income loss was reported at 58·4%, with higher rates seen in less developed settings such as Sub-Saharan Africa, Southeast Asia, east Asia, and Oceania [3]. As per 15 June 2022, 11,902,271,619 vaccine doses have been administered [1], yet the pandemic continues to affect the global population, especially the most vulnerable communities.
Ensuring health equity, especially for vulnerable populations in less developed settings with poor health systems is essential for the current COVID-19 pandemic and future global health threats [5-10]. At community-level, it is known that the severity of illness and clinical outcomes can be affected by the concentration of comorbidities in susceptible groups in communities [11-14], and through disparities in access to health care for preventive measures or prompt diagnosis and treatment [14,15]. Recent findings in US, Chile and Brazil suggested that COVID-19 poor outcome was concentrated in groups with higher socio-demographics and health system vulnerability [14, 16-20]. However, the heterogeneity of COVID-19 impacts from low- and middle-income countries (LMIC) where differences in age distribution, comorbidities, access to quality health services, and other factors, may significantly influence mortality risk, are limited.
Indonesia has suffered the highest number of COVID-19 confirmed cases and deaths in Southeast Asia, second only to India in all of Asia [21], at 6,068,075 cases and 156,687 deaths (2·6% case fatality rate (CFR)) up to 19 June 2022 [22]. The first SARS-CoV-2 epidemic wave occurred from 2 March 2020 to 30 April 2021, and a more intense second wave dominated by Delta variant peaked in July 2021 [23], followed by the third wave peaked in February 2022 [22]. The majority of cases and deaths in Indonesia were reported in Java Island, a more developed setting populated by 152 million individuals (56% of the total Indonesia’s population). Recent studies from Indonesia’s capital city of Jakarta suggested that COVID-19 disproportionately affected individuals with older age and pre-existing chronic comorbidities, as well as those areas with lower vaccine coverage, and higher poverty and population density [24,25]. However, data on the impact of COVID-19 across 514 districts of Indonesia that have different pre-existing burdens of major infectious diseases such as malaria, tuberculosis, HIV and other tropical infections [26], as well as non-communicable diseases like cardiovascular diseases, cancers, chronic pulmonary diseases, diabetes, and others [27], are scarce.
Indonesia is the fourth most populous country (population 270 million) and the LMIC featuring great geographic, cultural and socio-economic diversity across the archipelago. The 2020 human development index (HDI) ranged from 0·32 in Nduga District, Papua to 0·87 in Yogyakarta city, Yogyakarta [28]. In addition, substantial proportions of the Indonesian population face barriers in accessing quality health care services due to under-resourced and fragile health systems [26]. A heavily decentralised health systems [30] has resulted in distinct public health capacity across 514 districts of Indonesia. For example, the 2018 Public Health Development Index (PHDI) ranges from 35% in Paniai district, Papua province to 75% in Gianyar district, Bali province [31]. That heterogeneity and the large number of COVID-19 cases and deaths provides insights directly relevant to the national public health response to the COVID-19 crisis, and other LMIC settings. In this study, we assessed the COVID-19 heterogeneity, district-level socio-economics, population health-related conditions, and health care capacity among 514 districts of Indonesia and how those factors were associated with COVID-19 mortality rate during the first 24 months of the epidemic in Indonesia (March 2020 through February 2022).
Methods
Study design and participants
This was a nationwide ecological study to assess COVID-19 burden, heterogeneity, and factors associated with mortality rate in all 514 districts in Indonesia. The study analysed aggregated data of individuals diagnosed with COVID-19 based on either rapid antigen diagnostic test (Ag RDT) or polymerase chain reaction (PCR) recorded by the National COVID-19 Task Force from 1 March 2020 to 27 February 2022. In accordance with Indonesia’s national COVID-19 guideline, individuals are categorized as COVID-19 confirmed cases if tested positive by Ag RDT or PCR.
The main outcome variable in this study was district-level COVID-19 mortality rate per 100,000 population. The determinants assessed were gender, proportion of ≥60 years old individuals, prevalence of hypertension, prevalence of diabetes mellitus, prevalence of obesity, prevalence of pneumonia, COVID-19 vaccine coverage for ≥60 years old population, number of doctors, nurses, midwives, and hospitals per 10,000 population, expenditure, life expectancy at birth, and length of formal education.
Patient and public involvement
This study did not include patients and public in the design, or conduct, or reporting, or dissemination plans. This study was a secondary analysis of aggregated routine surveillance data conducted under a formal collaboration between the National COVID-19 Task Force (Komite Pengendalian COVID-19 dan Pemulihan Ekonomi Nasional) and United Nation Development Program, Indonesia.
Data collection
The aggregated data of weekly number of COVID-19 cases, incidence rate per 100,000 population, number of COVID-19 deaths among confirmed cases, mortality rate per 100,000 population, by district, sex, and age group from 1 March 2020 to 27 February 2022 were collected from the government official COVID-19 database managed by the National COVID-19 Task Force.
District-level data were collected from various government sources. Data on the number of populations were collected from Statistics Bureau Database available in each province. Data on the number of doctors, nurses, midwives, and hospitals per October 2021 were collected from the Indonesia Ministry of Health records. Data on prevalence of hypertension, diabetes mellitus, obesity, and pneumonia were collected from the National Public Health Development Index 2018 Report [31]. Data on the COVID-19 two-dose vaccine coverage for ≥60 years old population by district per 27 February 2022 were collected from the National COVID-19 Vaccination Database. Data on expenditure, life expectancy at birth, and length of formal education were collected from the Human Development Index 2020 Report [28].
Statistical analysis
Incidence and mortality rate per 100,000 population was calculated by time (weekly), sex, and province. Male to female incidence and mortality rate ratio (RR) was calculated to compare risk between males and females, and by province. Maps showing spatial distribution of male to female incidence and mortality rate ratio were generated using QGIS 3.20 software.
District-level incidence and mortality rate per 100,000 population, proportion of ≥60 years old population, prevalence of hypertension (%), prevalence of diabetes mellitus (%), prevalence of obesity (%), prevalence of pneumonia (%), COVID-19 vaccine coverage for ≥60 years old population (%), number of doctors, nurses, midwives, and hospitals per 100,000 population, expenditure (millions IDR), life expectancy at birth (years), and length of formal education (years) were calculated then categorised into quartiles. Descriptive statistics included proportions and the chi-squared test to compare district-level characteristics between different quartiles. Spearman’s correlation tests were done to assess correlation between each district-level variable. We used bivariable and multivariable linear regression models to determine factors associated with higher mortality rate at district level, expressed as beta, with 95% confidence intervals. All independent variables with p-value <0.10 in bivariable analysis were included in the multivariable models. Final model selection was informed by likelihood ratio tests. We set statistical significance at 0·05, and all tests were two-sided. All analyses were done in Stata/IC 15.1 (StataCorp, College Station, TX, USA). This study is reported as per Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines [31].
Ethics
This study was a secondary analysis of aggregated routine surveillance data collected by COVID-19 control program, with district as the analysis unit. Use of the data, which was aggregated, and at the population level is permitted by the Indonesian Ministry of Health under Regulation Number 45 (2014), Article 3, paragraphs 1 and 2. There was no ethical approval sought for this study.
Role of the funding source
The funder of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all data and the final responsibility to submit for publication.
Results
Cumulative COVID-19 incidence and mortality rate
Table 1 presents the reported cumulative COVID-19 incidence and mortality rate in Indonesia between 1 March 2020 and 27 February 2022. A total of 5,539,333 COVID-19 cases were recorded on the Indonesia National COVID-19 Database. Of those, 148,034 (2·7%) were deceased, and 5,391,299 (97.4%) were recovered. The overall incidence and mortality rates were 2,050 and 55 per 100,000 population, respectively. The incidence and mortality rate varied widely by province. The highest incidence rate was reported in DKI Jakarta (10,764 cases per 100,000 population), and the lowest was reported in West Nusa Tenggara (667 cases per 100,000 population). The highest mortality rate was reported in East Kalimantan (156 per 100,000 population), while the lowest was reported in Papua (13 per 100,000 population).
Heterogeneity of COVID-19 incidence and mortality rate at province-level
Figure 1 shows the weekly incidence and mortality at the national and province-level. The data suggest that Indonesia had experienced three pandemic waves with the highest mortality recorded in week 73rd (12 to 18 July 2021) where cases were predominantly infected by Delta variant (Figure 1). The province-level data show that DKI Jakarta had consistently experienced the highest incidence rate over time (Figure 2A). Overall, the mortality rate increased over time, with the highest rate seen in East Kalimantan, followed by Yogyakarta, DKI Jakarta, and North Kalimantan (Figure 2B).
The mortality rate was higher for males than females across 34 provinces in Indonesia (Table 2), with the highest male to female mortality rate ratio (RR) seen in East Nusa Tenggara (RR 3.00), followed by Yogyakarta (RR 2.61), Southeast Sulawesi (RR 2.52), Banten (RR 2.51), and West Papua (RR 2.50). The lowest male to female mortality was reported in Aceh province (RR 1.15). The incidence rate was also higher for males than females in the majority of provinces (94.12%), except for in Aceh and Bali (RR 0.92 and 0.99, respectively) (Supplementary Table 1). The highest male to female mortality rate ratio was seen in East Kalimantan (RR 2.42), followed by North Kalimantan (RR 2.17), Riau islands, Papua, and West Papua (RR 2.12 each).
Heterogeneity and factors associated with higher mortality rate at district-level
Table 3 presents the summary of district-level COVID-19 burden, prevalence of health-related conditions, and health care capacity associated with COVID-19 mortality rate across 514 districts in Indonesia. The incidence and mortality rate were highly heterogeneous, ranged from 8.13 to 10,626.09 per 100,000 population, and 0.28 to 283.82 per 100,000 population, respectively. Compared to eastern Indonesia, the mortality rate tends to be higher in the western areas, especially in districts and cities located in Java and Kalimantan Island (Figure 3A). Higher incidence rate was seen in Kalimantan, Sulawesi, and Papua Island (Figure 3B). Based on quartiles, the mortality rates ranged from 0.28 to 16.43 per 100,000 population (quartile 1), 16.43 to 31.47 per 100,000 population (quartile 2), 31.47 to 64.68 per 100,000 population (quartile 3), and 64.68 to 283.82 per 100,000 population (quartile 4). Of 129 districts in the highest mortality rate quartile, 71.32% had the highest incidence rate.
Regarding population health-related conditions, 44.19% of districts with the highest mortality rate had the highest proportion of ≥60 years old population (10.43-20.38%), 34.11% had the highest prevalence of hypertension (34.90-50.00%), 48.06% had the highest prevalence of diabetes mellitus (2.50-50.00%), 37.98% had the highest prevalence of obesity (35.10-50.00%), and 13.95% had the highest prevalence of pneumonia (2.40-13.00%). Regarding health care capacity, 47.29% of districts with the highest mortality rate had the highest COVID-19 vaccine coverage for ≥60 years old population (45.79-86.20%), 13.28% had the lowest number of doctors per 100,000 population (0.15-5.43), 24.22% had the lowest number of nurses per 100,000 population (0.51-10.59), 12.50% had the lowest number of midwives per 100,000 population (0.28-11.46), and 43.31% had the lowest number of hospitals per 100,000 population (0.29-3.90). Regarding socio-economic, only 3.10% of districts with the highest mortality rate had the lowest expenditure (IDR 3.98-8.51 million), 3.88% had the lowest life expectancy at birth (55.27-67.24 years), and 15.50% had the lowest mean of length of formal education (1.13-7.43 years) (See Table 3 for details).
In bivariable linear regression analysis (Table 4), district-level COVID-19 mortality rate was associated with COVID-19 incidence rate, proportion of ≥60 years old individuals, prevalence of hypertension, prevalence of diabetes mellitus, prevalence of obesity, prevalence of pneumonia, vaccine coverage for ≥60 years old population, doctors per population, hospital per population, expenditure, life expectancy, and length of formal education. We found no association between mortality rate with number of nurses and midwives per population. Correlation matrix of COVID-19 burden, prevalence of health-related conditions, vaccine coverage for >60 years old population, health care capacity, and socio-economic characteristics can be seen in Figure 4.
Significance level of the correlation can be seen in Supplementary Table 2.
In the multivariable model, higher mortality rate was associated with higher COVID-19 incidence (p<0.0001), higher proportion of ≥60 years old population (p<0.0001), higher prevalence of diabetes mellitus (p<0.0001), lower prevalence of obesity (p=0.019), lower number of doctors per population (p=0.001), higher life expectancy at birth (p=0.035), and shorter length of formal education (p=0.021). There was no association between mortality rate with prevalence of hypertension, prevalence of pneumonia, vaccine coverage for ≥60 years ole population, number of nurses, midwives, and hospitals per population, and expenditure (p>0.05 each).
Discussion
This nationwide ecological study described the epidemiological surveillance data of confirmed COVID-19 cases reported from 514 districts, during three epidemic waves in the first 24 months of the SARS-CoV-2 transmission in Indonesia. The overall case fatality rate was 2·7% (148,034/5,539,333), which equate to a mortality rate of 55 per 100,000 population. The mortality rate was highly heterogeneous and increased over time, with the highest rate seen in East Kalimantan, followed by Yogyakarta, DKI Jakarta, and North Kalimantan. The mortality rate tends to be higher for males than female in the majority of provinces in Indonesia, with the highest male-to-female mortality rate ratio seen in eastern Indonesia provinces. At district-level, higher mortality rate was associated with higher COVID-19 incidence, higher proportion of ≥60 years old population, higher prevalence of diabetes mellitus, lower prevalence of obesity, lower number of doctors per population, higher life expectancy at birth, and shorter length of formal education.
Findings from this study revealed the extent of health inequity and its association with COVID-19 mortality across 514 districts in Indonesia. In contrast to a previous province-level analysis (n=34) found the absence of correlation between mortality rate with doctor to population ratio during the first six months of the epidemic [33], our district-level analysis found that a higher COVID-19 mortality rate was strongly associated with a lower number of doctors per 100,000 population at district-level. This finding contrasts a previous finding from a higher-income country, reporting an association between higher COVID-19 mortality rate with higher physicians density in France [34]. In Indonesia, the lack of doctors, combined with the high number of infections [35] and deaths among health care workers and doctors [36], especially during the early phase of the epidemic, had weakened the health system resilience. In the context of a heavily decentralized health system such as in Indonesia [30], coordination and prioritisation of available resources and public health intervention will be critical to ensure optimal health outcomes for vulnerable communities, especially for those areas with lower number of doctors per population. Short-term solutions such as reallocating doctors to districts with the lowest doctor to population ratio may assist local government in mitigating this inequity in health care access. Improving access to medical doctor training for local residents, and providing incentives to doctors in rural areas can provide long-term solution in ensuring future health equity.
Consistent with individual-level risk factors at from diverse settings [18, 25, 37-38], our district-level analysis affirms that higher proportion of elder population (also reflected by higher life expectancy at birth), males gender, and higher prevalence of diabetes mellitus were significantly associated with higher COVID-19 mortality rate in the general population of Indonesia. This was also consistent with findings from a previous study suggesting that individual-level risk factors such as older age, male sex, and pre-existing hypertension and diabetes were associated with an elevated risk of COVID-19 death during the early epidemic phase in DKI Jakarta, Indonesia [25]. These findings indicate the need of controlling both SARS-Cov-2 transmission and non-communicable diseases, especially in areas with high prevalence of diabetes mellitus. Reducing COVID-19 mortality in such settings may necessitate comprehensive and specific interventions such as improving diagnosis and case management of those known non-communicable diseases, community awareness, as well as a sustainable and accessible social security network that may reduce vulnerability of these communities.
Our analysis reveals the heterogeneity of COVID-19 burden across 34 provinces and 514 districts in Indonesia. The province-level findings suggesting that mortality and incidence rate was higher for males than females, may reflect the cultural phenomenon in Indonesia where males were more likely to get more exposure to the SARS-Cov-2 due to their domestic responsibility to work even during the social restrictions were implemented. Our following findings suggesting higher mortality rate among districts with lower-level of formal education (reflected by shorter mean length of formal education) was consistent with findings from more developed settings such as in the US [11, 39] and Sweden [40]. A lower educational level might be associated with a lower health literacy that can lead to lower access and understanding of public health information. A tailored risk communication and community engagement strategy targeting population living in districts with lower education level is warranted in addressing health inequity in the future.
A previous cohort study in DKI Jakarta reported that higher COVID-9-related mortality risk was significantly associated with lower COVID-19 vaccine coverage at subdistrict-level [24]. By contrast, our current study found that there was no association between district-level COVID-9 mortality rate with COVID-19 vaccine coverage among >60 years old population in Indonesia. This absence of association could possibly be explained by the very low vaccine coverage across the 514 districts. The vaccine coverage was highly varied from 0.03% to 86.20%. The majority of districts (82.3%) still had a vaccine coverage of less than 50% after one year of vaccine roll-out (27 February 2022), thus further highlights the extent of inequity across the country. A previous study from Brazil reported that rapid scaling up of vaccination coverage among elderly Brazilians was associated with significant declines in relative mortality compared with younger individuals, in a setting where the gamma variant predominated [41]. Those findings from DKI Jakarta and Brazil, highlight that rapid vaccination roll outs targeting the most vulnerable is crucial in reducing COVID-19-related deaths. As per 19 June 2022, the coverage for two-dose and three-doses COVID-19 vaccination in Indonesia was 80.14% (166,911,457/208,265,720 targeted population) and 21.26% (44,273,456/208,265,720), respectively [42].
This study had several limitations. Firstly, this study was based on routine surveillance data derived from district-level health office case reports. The imperfect contact tracing, testing, and reporting activities could result in underreporting of cases, especially those asymptomatic and mild cases, which could cause an overestimation of case fatality rate and an underestimation of mortality rate per population in this study. Secondly, as this study was based on aggregated data available at district level, our analysis was unable to capture and adjust the effect of important individual-level risk factors. Therefore, interpretation of the study findings should be restricted to district-level risk factors. Finally, previous studies in Indonesia have suggested the association between COVID-19 mortality rate with population density [24] and ratio of hospitals bed [33]. However, these variables were not evaluated in our study, due to limited access to the relevant data.
In conclusion, our analysis suggested that COVID-19 disproportionately affected districts with high proportion of elder population, high prevalence of diabetes mellitus, lower doctor to population ratio, higher life expectancy at birth, and lower level of formal education. These findings indicate that vulnerability to death associated with COVID-19 includes not only elders and comorbid, but also males and communities living in area with lower health care capacity and with lower level of education. These findings may inform decisions on health resource allocation against COVID-19 delivering the greatest possible health dividends by prioritising interventions, including even distribution of essential health care need such as doctors, and a tailored risk communication and community engagement for the most vulnerable communities in LMIC, especially with decentralised health systems like in Indonesia. Future nationwide studies incorporating individual and district-level data to assess vulnerability associated with COVID-19-related morbidity and mortality are needed to better comprehend the COVID-19 impact and to better prioritise interventions for the most vulnerable communities.
Data Availability
After publication, the datasets used for this study will be made available to others on reasonable requests to the corresponding author, including a detailed research proposal, study objectives and statistical analysis plan.
Contributors
HS was the principal investigator of this study. HS designed the study, did the analysis, and had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. HS, DP, NNA, AIP, AAS, and EP contributed to data collection and verification. LR, DG, and MG supervised the data collection and verification. HS and IRFE drafted the paper. All authors critically revised the manuscript for important intellectual content and all authors gave final approval for the version to be published.
Funding
There was no specific funding for this study. HS, DP, and NNA were supported by funding from UNDP Indonesia. HS, EP, and IRFE were funded by the Wellcome (UK) Africa Asia Programme Vietnam (106680/Z/14/Z).
Conflict of interests
We declare no competing interests.
Data sharing statement
After publication, the datasets used for this study will be made available to others on reasonable requests to the corresponding author, including a detailed research proposal, study objectives and statistical analysis plan.
Acknowledgments
We acknowledge all health care workers involved in the care for the COVID-19 patients, as well as those involved in the field data collection.
Acknowledgment
We acknowledge the Ministry of Health of Republic of Indonesia, the National COVID-19 Task Force (Komite Pengendalian COVID-19 dan Pemulihan Ekonomi Nasional), and all health care workers involved in the care for the COVID-19 patients, as well as those involved in the field data collection.
References
