Introduction

Chikungunya is an emerging arbovirus that reached the Western Hemisphere at the end of 2013. On December 9, 2013, the Pan American Health Organization (PAHO) has issued an alert about the transmission of chikungunya in the Americas1. Since then, the transmission of CHIKV were confirmed in 44 countries and territories in the region, with more than 2 million reported cases and 403 deaths. In Brazil, the transmission was confirmed in September 2014, but the worst epidemic occurred in 2016 with 216,102 cases reported in 696 municipalities and 91 deaths.

Studies in Indic Ocean and India suggest that passive surveillance system cannot recognize many of deaths associated with chikungunya, which can be inferred by an increase in the overall mortality observed during chikungunya epidemics 3,4.

The lethality observed in Brazil and in all America, has been much lower than that reported in other countries1, perhaps due to difficulty in diagnosis and lack of adequate tools for surveillances1. The objective of this study was to assess the mortality associated with chikungunya epidemics in Brazil by comparing the mortality of pre-chikungunya with mortality during chikungunya epidemics2,3.

Methods

We studied the 3 Brazilian states with higher chikungunya incidence rates (IR), Pernambuco (9,410,772 inhabitants), Rio Grande do Norte (3,474,998 inhabitants), and Bahia (15,498,733 inhabitants)5. Bahia is the largest and most populated state, with a heterogeneous incidence of chikungunya, and therefore we analyzed separately each one of 9 administrative regions of Bahia. We retrospectively analyzed anonymous, publicly available data from the Brazilian Ministry of Health.

To evaluate the possible effect of chikungunya on all-cause mortality, we compared the number of deaths expected for the current estimated population with that observed during the epidemic, during 2015 and 2016. Given the lack of previous data of the CHIKV circulation in the studied regions, we arbitrarily considered as epidemic period onset on first of at least two consecutive months with incidence rates higher than 50 cases per 100 thousand population; and for the end of epidemic period we defined as incidence rate lower than 50 per 100 thousand for two consecutive months for all studied regions. We plot the incidence rate of dengue during 2015 and 2016 in all studied regions for comparing.

The expected mortality rates were calculated based on the mean monthly Age Standardized Mortality Rate (ASMR, using the Brazilian population, census 2010) of the same months between 2011 to 2013, years prior to chikungunya transmission. Values conservative 99% of the upper limit of the confidence interval were calculated for each month as a maximum threshold of the monthly-expected deaths for each studied region. We defined excess deaths as the difference between the observed and expected deaths during the epidemics months.

We also compared the expected and observed ASMR for each month during the chikungunya epidemic period and calculated change in mortality rate per 100,000 and percentage change (%) in each region. Similar methodology has been used previously2,3,4. In most affected areas, we also calculated the mortality rate for epidemic period and 99% of the upper limit of the confidence interval by age group. We used the mean monthly mortality rate by age group of the same periods of the non-epidemic years 2011-2013 as expected mortality rates.

We obtained population data (census and projections) from the official Brazilian Institute of Geography and Statistics (IBGE), and mortality data from the official Mortality Information System. Probable cases of chikungunya, dengue, and zika virus infection (laboratory and clinical epidemiological criteria) were available in the Brazilian Epidemiological Surveillance System (SINAN), based on cases reported by healthcare professionals. For the processing of data and figures and the statistical analysis, we used Excel-2013® and SPSS v. 20.

Results

Using our criteria of epidemics, it was possible to identify one wave CHIKV outbreaks in: Pernambuco (Jan-Apr 2016), Rio Grande do Norte (Feb-Jun 2016), and Noth-center Bahia (Feb-Mar 2016), Northeast Bahia (Jan-Mar 2016), North-Bahia (Sep 2015 – Apr 2016), South-Bahia (Jan-Apr 2016); and two waves outbreaks in: East-center Bahia (first Mar-Jul 2015, and second Feb-Mar 2016). Using the same criteria no epidemics was confirmed in other regions of Bahia (Table 1).

During CHIKV epidemics there was an increase in all-cause mortality in all analyzed regions, almost simultaneously with the peak of epidemics (Figure 1). Excess mortality rates were higher in the four region with the highest incidence rate of the chikungunya cases (table-1). In states of Pernambuco, South-Bahia, and Rio Grande do Norte incidence of dengue were high in 2015 and 2016. In Pernambuco, the highest incidence of dengue occurred in 2015, a year in which there were no excess deaths (Figure 1).

The excess mortality rates per 100,000 population and the excess deaths in the 4 critical regions studied were: Pernambuco (47.9 and excess death = 4,505) corresponding to the 21% increase in deaths, Rio Grande do Norte (42.5 and 1,478), South-Bahia (32.0 and 546), and North-Bahia (30.0 and 339) (Table 1). Table 2 shows excess mortality rates by age group for each region and Figure 2 presents time series incidence rate of chikungunya and mortality rate by age groups in Pernambuco.

Table 1. Observed and expected deaths, all-causes mortality rates (deaths/100.000 population), excess mortality, change in mortality rate and incidence rate of chikungunya, zika, and dengue virus in 3 states of Northeast Brazil.

Table 1: Observed and expected deaths, all-causes mortality rates (deaths/100.000 population), excess mortality, change in mortality rate and incidence rate of chikungunya, zika, and dengue virus in 3 states of Northeast Brazil.

Figure 1 2nd version

Fig. 1: Monthly observed Age Standardized Mortality Rate (ASMR, blue bar) and expected ASMR (white bar), 99% confidence interval (red line), monthly incidence rate of chikungunya (black line) and dengue (blue line) in Pernambuco, Rio Grande do Norte and Bahia (two regions).

OBK-17-0020 Table 2

Table 2: Observed, expected all-causes mortality rates (deaths/100,000 population) and excess mortality by age group during epidemics periods in Pernambuco, Rio Grande do Norte and Bahia (two most affected regions, North and South).

Figure 2 2nd version

Fig. 2: Monthly all causes mortality rate by age group and monthly incidence rate of chikungunya from 2015 to 2016 in Pernambuco.

Discussion

The understanding of the current chikungunya burden depends on a better understanding of the mortality associated with this arboviruses. Our findings show excess mortality during chikungunya epidemics in three most affected states in Brazil. The correlation between the incidence of chikungunya and increase in all-cause mortality reinforce the possibility that chikungunya virus has an impact on mortality rates. Similar results have already been found in other chikungunya epidemics in the past decades2,3,4. Incidence of dengue virus was high in all studied areas in the years 2011-2013, used as a baseline5. Therefore, we suggest that the excess deaths should be related to chikungunya infections and not another arbovirus that occurs at the same season. This is the first time that the excess mortality associated with chikungunya by age group has been calculated, and the results are very close to that found in the Reunion Islands epidemic using the data obtained in the death certificates16.

Since the chikungunya arrive in Americas, the virus has been identified in 44 countries and territories, leading to more than 2 million cases and 403 deaths6. In Brazil, the first case was confirmed in September 2014; however, the worst outbreak occurred in 2016, with 236,287 reported cases and 120 officially confirmed deaths5. The Brazilian Epidemiological Surveillance System detected far fewer deaths compared to our estimates in Pernambuco (54 deaths), Rio Grande do Norte (19 deaths), and all of Bahia (5 deaths)5. Mortality identified by the passive surveillance system was approximately 50 times lower than estimated from the excess of deaths. These differences may be due to failure to diagnose and difficulty recognizing severe forms of chikungunya1,15. Deaths caused by chikungunya occur more frequently in the elderly as a complication of preexisting diseases as reported in other hospital-based studies7,8. The present study showed that the increase in mortality was higher than the 99% confidence interval threshold in all regions in the age groups above 50 years. In fact the elderly are known to be the age group with the highest risk of severe forms and deaths due to chikungunya7. Because chikungunya is a disease not yet well known by health professionals in Brazil, the association of exacerbation of the underlying disease with chikungunya may have been underestimated leading to underreporting. In the age range between 15 and 19 years, there was an increase in mortality rate above the 99% confidence interval in 3 of the 4 regions studied. Clinical complications due to CHIKV infection have been reported in adults, not specifically in adolescents. Though, a systematic review on systemic infection due to CHIKV reports 52.2% and 37.1% articles respectively about cardiovascular and neurological manifestations11. On the other hand, severe cases and deaths of children and adolescents with postnatal zika virus infection are not commonly described in the literature9. However, Azevedo et al. (2016) draws attention to deaths in young people (16, 20 and 35 years old) with autoimmune diseases and immunodeficiency, besides the syndrome Guillain Barre10. Future investigations are necessary to understand better this occurrence.

There was also an increase, greater than 99% confidence interval threshold, in the mortality rate in children under 1 year. These latter findings had not been previously described and need to be better understood, however, it is not possible to exclude the possibility that the excess mortality in these age groups was related to the circulation of zika virus12.

The arbitrary definition of the epidemic period from 50 cases/100,000 could be less sensitive to express an epidemic context at lower incidence rates, however it was a parameter to enable comparisons between different periods and regions.

We must analyze results of time-series studies very carefully because of the method’s own difficulties, which include time-coincident confounding factors. However, the periods of excess mortality observed do not coincide exactly with the seasons of influenza in the states of the Northeast region. Influenza AH1N1 was the most frequent viral subtype in 2016 and has not been found to be associated with high elderly mortality. Zika virus was not a reportable disease in 2015 there is no official data to calculate the IR of zika virus in that year, but even so the mortality associated with zika is admittedly low13. Other factors could be associated with excess mortality in that period; but no changes in the epidemiological profile of other diseases or events were reported in the studied regions. Furthermore, the overall arbovirus incidence may be underestimated due to sub diagnosis and under-reporting by health professionals. The impact of arbovirus co-circulation is still poorly understood14. The interaction of arboviruses, could theoretically result in more intense viremia or other immunological alterations.

Considering the magnitude and rapid territorial expansion of chikungunya transmission, we believe that global health systems should be prepared to manage the severity of the disease and associated deaths, not only acute and chronic joint disorders that have been the main concern related to chikungunya historically. Excess mortality is a useful indicator to quantify the impact of a health-related event and has traditionally been used to describe the increase of deaths during influenza season and extreme climates events8,10,15. All causes mortality during chikungunya epidemics could be monitored as a strategic tool, beyond individual case reporting to the epidemiological surveillance system, to assess excess mortality and the overall burden of chikungunya.

Competing Interests

The authors have declared that no competing interests exist.

Data Availability

Data are available from the figshare repository: https://doi.org/10.6084/m9.figshare.5397160.

Corresponding Author

André Ricardo Ribas Freitas: [email protected]