key: cord-0884134-mwjk2dkb
authors: Hanifi, Syed Manzoor Ahmed; Alam, Sayed Saidul; Shuma, Sanjida Siddiqua; Reidpath, Daniel D.
title: Insights Into Excess Mortality During the First Months of the COVID-19 Pandemic From a Rural, Demographic Surveillance Site in Bangladesh
date: 2021-07-29
journal: Front Public Health
DOI: 10.3389/fpubh.2021.622379
sha: 64f35b3dfb4e27e6b2b432113a0ebdb49bd8b31a
doc_id: 884134
cord_uid: mwjk2dkb

Background: Coronavirus disease 2019 (COVID-19) has spread globally, and the government of each affected country is publishing the number of deaths every day. This official figure is an underestimate as it excludes anybody who did not die in a hospital, who did not test positive, who had a false result, or those who recovered on their own without a test. Objective: This study aimed to measure the community level excess mortality using health and demographic surveillance in a rural area of Bangladesh. Method: The study was conducted in Matlab, in a rural area of Bangladesh, with a Health and Demographic Surveillance System (HDSS) covering a population of 239,030 individuals living in 54,823 households in 142 villages. We examined the mortality in January-April from 2015 to 2020 and compared the mortality in 2020 with the historical trend of 2015–2019. Between 2015 and 2020, we followed 276,868 people until migration or death, whichever occurred first. We analyzed mortality using crude mortality rate ratio (MRR) and adjusted MRR (aMRR) from a Cox proportional hazard model. Mortality was analyzed according to age, sex, and period. Results: During follow-up, 3,197 people died. The mortality rate per 1,000 person-years increased from 10 in 2019 to 12 in 2020. Excess mortality was observed among the elderly population (aged 65 years and above). The elderly mortality rate per 1,000 person-years increased from 80 in 2019 to 110 in 2020, and the aMRR was 1.40 (95% CI: 1.19–1.64). Although an increasing tendency in mortality was observed between 2015 and 2019, it was statistically insignificant. Conclusions: The study reported a 28% increase in excess deaths among the elderly population during the first months of the pandemic. This all-cause mortality estimation at the community level will urge policymakers, public health professionals, and researchers to further investigate the causes of death and the underlying reasons for excess deaths in the older age-group.

Coronavirus disease 2019 , caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (1), first emerged in Wuhan city of Hubei province in China on December 31, 2020, when Chinese health officials informed the WHO about a cluster of 41 patients with mysterious pneumonia, supposedly connected to Huanan Seafood Wholesale Market (2) . The detection time of the first cases varied from country to country; some counties detected early in the pandemic and some later. For instance, Bangladesh announced the first cases on March 8, 2020 (Table 1) , and 2 days later on March 11, 2020, WHO declared COVID-19 as a pandemic (2) .

The number of people migrating overseas from Bangladesh for employment annually is more than 400,000 (15) . During the first months of 2020, more countries were placing lockdowns after the detection of the first cases of COVID-19 surfaced in each country. Bangladeshi migrant workers residing in these countries were compelled to return home to rural Bangladesh due to lack of income (16, 17) . In addition, when Bangladesh announced its own lockdown on March 22, 2020 (18) , after detecting the first three cases, an increasing number of people were leaving the capital city (19) . Such movement of people opens the window for the spreading of the virus to more places in the country.

Official numbers of COVID-19 cases and deaths in Bangladesh are likely to be an underestimate of the real scenario because this mainly accounts for cases that have tested positive of coronavirus through laboratory confirmation (20) and deaths that are recorded in hospitals. What the numbers have missed could be cases of deaths and infections before testing started, and the first cases were confirmed when only six PCR laboratories were available throughout the country (21) . After the tests were made available, false test results could have eliminated any real cases and unaccounted for people who did not step forward for tests fearing isolation and stigma (22) or those that recovered from taking treatment at home (23) . In addition, it is also difficult to determine the cause of death in many instances to find out whether the person was COVID-19 positive if they died prior to testing (24) . Weak civil and vital registration statistics (CVRS) for the great majority of low-and middle-income countries (LMIC) (25) means that we have remarkably little insight into the magnitude of the excess mortality associated with the COVID-19 pandemic. This is in marked contrast to some high-income countries; for instance, during the first months of pandemic the United Kingdom published preliminary all-cause mortality data and COVID-19-related mortality data with a lag of only a few weeks unlike Bangladesh (26) .

Although most countries are submitting daily data on the number of COVID-19 deaths to the WHO (27) , without comparable all-cause mortality data, and all-cause mortality data for the equivalent time period over preceding years, it is difficult to estimate the excess mortality attributable to the pandemic. The excess mortality, of course, is not restricted to COVID-19 deaths alone; it also includes the non-COVID-19 deaths that arise from the loss of adequate care, as health systems become overstretched coping with COVID-19 cases.

Health and Demographic Surveillance System (HDSS) conducts surveillance of geographically prescribed populations for extended periods of time (28) . They have, historically, provided estimates for many LMICs about the underlying birth, death, and fertility rates in the absence of an effective CVRS system (28) . The Matlab HDSS is the longest running HDSS which has provided Bangladesh with some of its earliest data on rates of births and deaths (29) . Utilizing these ongoing surveillance data, we are able to estimate the age and sex of mortality rates over the past 6 years and thereby to estimate the excess mortality in 2020 associated with the COVID-19 pandemic.

The Matlab HDSS covers a population of 239,030 individuals living in 53,823 households in 142 villages in a rural district of Bangladesh that is situated 60 km south of Dhaka (30) . All households are visited every 3 months to enquire about marriages, pregnancies, births, migrations, and deaths. Annually in the Matlab HDSS area, 5,298 babies are born, 1,687 people died, 2,671 people moved out of Bangladesh, and 1,371 people returned from abroad (30) .

A web-based software application has been designed and developed for Matlab HDSS. Thirty-one tablets (smartphones) are connected to the mobile internet through the network of mobile operators. Traditionally, community health research workers (CHRWs) visit households every 3 months and record Frontiers in Public Health | www.frontiersin.org health and demographic events using these devices, and data are stored in the central database server. In the COVID-19 pandemic, to follow the precautionary guidelines in Bangladesh, CHRWs continued registering birth and death through mobile phones instead of household visits between March 25, 2020, and November 9, 2020. During data collection through mobile phones, CHRWs reached 85% of households in the first contact. Information of absent households (in the first contact) was updated in the next round.

We analyzed mortality data in the study area between January 2020 and April 2020 to take into account the transmission of coronavirus, not only related to the movement of population from the capital city to rural areas after declaring the first lockdown on March 22, 2020, in Bangladesh but also the spread of the virus exacerbated by the migrants returning from abroad during the first months of the pandemic (January 2020-April 2020) (24, 31) . The mortality rate is compared with the mortality rates for the period of January 1 and April 30 in the years from 2015 to 2019. Between 2015 and 2020, we followed 278,290 people until migration or death, whichever occurred first (Figure 1) . We calculated the mortality rate between January and April as seasonality in mortality was observed in the study area (30) (Figure 1) . Deaths due to COVID-19 vary by age (32) and sex (33) . We presented age-specific and sex-specific mortality and population size as the denominator. We used Cox proportional hazard models with age as the underlying timescale to estimate the aMRR. We included sex (33) , period (30) , and village (30) as potential confounders in the adjusted models.

In the Matlab HDSS area between January 2020 and April 2020, 1,008 people returned from abroad, which is 2.39 times (95% CI: 2.29-2.55) higher than in 2019. Moreover, in 2020, in-migration was very high, which is 3.87 times higher (95% CI: 3.10-4.14) in January and February compared to March and April (Figure 2 ).

During the follow-up period from January to April from 2015 to 2020, a total of 3,197 people died. In 2020, the number of people migrating back to Matlab from abroad was 406 in January, 375 in February, 207 in March, and 17 in April. From 2015 to 2019, the crude mortality rates (CMRs) per 1,000 person-years were 7.37 (95% CI: 6.75-8.05), 7.63 (95% CI: 7.01-8.31), 7.81 (95% CI: 7.18-8.50), 8.98 (95% CI: 8.29-9.73), and 9.70 (95% CI: 8.99-10.48), respectively. In 2020, the CMR was 11.61 (95% CI: 10.51-12.82) ( Table 2 ). Figure 3 shows the age by year mortality rates. Between 2015 and 2019, there is little evidence of substantial yearto-year variation in the mortality rate in any age-group. In 2020, the mortality rate in the 65+ age-group (older adults) is markedly higher than in the previous 5 years.

The mortality rate in older adults per 1,000 person-years increased from 80 in 2019 to 110 in 2020, and the aMRR was 1.39 (95% CI: 1.17-1.66) among older adults. The aMRR was 

Compared to the five previous years, there is a clear, statistically significant difference in the mortality rate in the 65+ age-group during the COVID-19 period in Matlab. This result is strong indication of a pandemic effect. No such difference was observed in any of the younger age-groups. The age by mortality effect is consistent with the global data on excess mortality that shows the highest mortality risk in those who are aged 65 years or older (32, 34, 35) .

The COVID-19 pandemic has had a significant negative impact on mortality rates in older people living in rural Bangladesh. It is not possible at this stage to determine the cause of deathalthough a verbal autopsy process, to be conducted in a few months, may give some insight into this. However, given the age effect in the excess mortality, the results strongly indicate either a direct or an indirect COVID-19 effect. The indirect effect may be attributable to a decreased likelihood of seeking life-saving care or a decreased capacity of the health system to manage non-COVID-19-related healthcare needs (36) . For instance, the prevalence of hypertension is 53% among older age-groups (aged 65 years and above) in Bangladesh (37) and these hypertensive patients may not have been able to avail regular checkups or acquire medicines during the pandemic. Moreover, about 50% of deaths occur in the study area per year (30) due to conditions of chronic disease like heart disease and stroke. The pandemic situation may be responsible for deaths among such patients by preventing them from traveling the distance to a hospital that could provide them the immediate intensive care services that they required.

We analyzed mortality from January to April, 2015-2020, which reduced any seasonality bias (30) . Mortality rates were calculated using person-time techniques (38) that remain the basic epidemiological approach to estimating mortality, yet one of them is frequently missed in the calculation of mortality (39) . The data were only collected from one rural area, which may not reflect the situation in all rural Bangladesh. Indeed, Matlab appears to have better health outcomes than other rural areas of Bangladesh, and this may indicate that excess mortality rates would be worse elsewhere. The relatively short period of observation was to ensure the timely reporting of data.

We cited a number of news media-released reports as there was no detailed timeline information on these events during the first months of the pandemic.

Globally, the COVID-19 pandemic impacted the mortality of the overall population. COVID-19 pandemic attributed 30 deaths per 1,000 among the older age-group in the study area. We did not determine the deaths related to COVID-19 and non-COVID-19 causes. A further cause of death analysis will provide an estimate of excess deaths associated with COVID-19 and non-COVID-19 causes. It is important to examine whether the excess deaths as a result of access to healthcare and how the national COVID-19 policy operates on the decisions and actions of the people. Bangladesh needs to strengthen the CVRS system and national health statistics to monitor timely morbidity and mortality, especially in an epidemic or pandemic situation. Bangladesh should strive to strengthen its health systems by providing additional resources to make healthcare services more accessible to its residents irrespective of geographical locations.

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

The studies involving human participants were reviewed and approved by Ethical Review Committee of icddr,b. The patients/participants provided their written informed consent to participate in this study.

SH and DR conceived and designed the study and are the guarantors of the study. SA prepared the data file. SH analyzed the data and wrote the first draft of the manuscript. SS contributed to the literature review. All authors contributed to the final version of the manuscript. All authors had full access to all of the data (including statistical reports and tables) in the study and can take responsibility for the integrity of the data and the accuracy of the data analysis.

Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study

World Health Organisation. Coronavirus Disease (COVID-19) Outbreak

A Comprehensive Timeline of the Coronavirus Pandemic at 1 Year, From China's First Case to the Present Online

World Health Organisation. Novel Coronavirus (2019-nCoV). World Health Organization

Timeline of the UK's First Recorded Covid Cases Last Year Online: The Guardian

Germany's First Coronavirus Case Is Human-to-Human Transmission Online: Al-Jazeera

First confirmed case of COVID-19 infection in India: a case report

First Middle East Cases of Coronavirus Confirmed in the UAE Online: CNBC. (2020)

Timeline -How the Coronavirus Spread in Spain

A Timeline of the Coronavirus Pandemic

COVID-19)

Saudi Arabia Announces First Case of Coronavirus

South Africa's War on COVID-19 Online: Think Global Health

Bangladesh Confirms Its First Three Cases of Coronavirus

htm#:~:text=Each%20year%2C%20more%20than %20400%2C000,the%20Bangladesh%20for%20overseas%20employment

Bangladeshs Sees a Huge Surge in Last Three Weeks Online (2020)

International Organization for Migration. IOM Reports that 70 per cent of Returning Migrants to Bangladesh Struggle to Find Employment

Available online at

Lockdowns in Asia Have Sparked a Stampede Home Online: The Economist (2020) Available online at

Coronavirus Pandemic (COVID-19). Our World in Data

Covid-19 Testing and Health Sector Resource Mobilisation

Bangladesh: Low COVID-19 Testing Rate Raises Concerns Online

COVID-19 Patients at Home Missing From DGHS Record in Bangladesh

Tracking Covid-19 Excess Deaths Across Countries

Civil Registration and Vital Statistics Benefit Health, Child Protection, and Governance: A Case Study on the Return on Investment for CRVS Systems

Available online at

World Health Organization. WHO Coronavirus Disease (COVID-19) Dashboard. (2020)

Health and demographic surveillance systems: a step towards full civil registration and vital statistics system in sub-Sahara Africa?

Health and demographic surveillance system (HDSS) in Matlab, Bangladesh

Available online at

Demographic science aids in understanding the spread and fatality rates of COVID-19

Gender differences in patients with COVID-19: focus on severity and mortality

Estimating excess 1-year mortality associated with the COVID-19 pandemic according to underlying conditions and age: a population-based cohort study

Excess mortality estimation during the COVID-19 pandemic: preliminary data from Portugal

Healthcare crisis in Bangladesh during the COVID-19 pandemic

National Institute of population research and training (NIPORT), ICF. Bangladesh Demographic and Health Survey 2017-18

Stata Survival Analysis and Epidemiological Tables Reference Manual

Covid-19 worldwide: we need precise data by age group and sex urgently

Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest

This research study was funded by core donors who provide unrestricted support to the icddr,b for its operations and research. icddr,b is grateful to the Government of Bangladesh, Canada, Sweden, and the UK for providing core/unrestricted support.

The authors are grateful to icddr,b for funding this study and the villagers for their cooperation in providing invaluable information. The untiring efforts of the team members of the Matlab HDSS in maintaining the surveillance system are gratefully acknowledged.