INTRODUCTION
Dengue, a vector-borne viral disease caused by four serotypes (dengue virus [DENV]-1-4) and transmitted by Aedes aegypti mosquitoes1, is endemic in most Latin American countries, and an important cause of hospitalization and death among children and adults in the region2. All four virus serotypes are capable of causing illness and are known to circulate in the region, with cocirculation of two or more serotypes at any given time in some countries. Such cocirculation of several serotypes has the potential to cause severe disease outbreaks3.
In the Americas, 10.6 million dengue cases and 6542 deaths caused by dengue were reported to the Pan-American Health Organization (PAHO)/World Health Organization (WHO) between 2008 and 20144. Of these, 1,025,825 dengue cases and 653 dengue-associated deaths were reported in Mexico. Annually, Mexico has accounted for 2.8-21.7% of reported dengue cases and 1.7-21.1% of dengue-associated deaths in the Americas. Case fatality rates (CFR) have ranged in Mexico from 0.04% to 0.10%4. Between 2008 and 2014, the number of dengue cases reported in Mexico annually had ranged from 25,040 to 249,886. Although the PAHO data are comprehensive, there is usually significant underreporting inherent in surveillance systems, which hinders assessment of the true burden of the disease. The level of under-reporting of dengue cases in Mexico at the national level may be as high as 46-fold5. Indeed, dengue seroprevalence was reported to be about 75% in typical endemic communities in Mexico6,7, suggestive of high dengue transmission rates in the country even in areas considered of low dengue transmission based on historical epidemiological surveillance data.
Although there are numerous reports on the medical and economic burden of dengue, the impact of the disease during outbreaks in terms of health-care management outcome of other morbid conditions unrelated to dengue remains to be adequately explored. Since dengue outbreaks are usually clustered in time or location - cases have been shown to be clustered mainly in two coastal areas of Mexico, the Pacific and Gulf8 - they may place a considerable burden on health-care systems and compromise treatment quality and decisions across all treatment wards, including degrading performance of clinical laboratories in the affected regions9. The information (clinical, financial, and operational) collated in hospital databases can provide insights into resource utilization10 and could be used to help determine the impact of dengue on health-care demand and other morbid conditions at a given time and location.
The Mexican public health-care system is divided into two groups. One group comprises social security institutions such as the Mexican Social Security Institute (IMSS), Institute of Social Security and Services of State Workers (ISSSTE), Mexican Oil Company (Petróleos Mexicanos, PEMEX), and the Armed Forces (SEDENA and SEMAR), and covers workers in the formal sector of the economy. The other group comprises institutions that cover the population without social security/insurance plans such as the Ministry of Health (SSa), the State Health Services (SESA), and the IMSS-Opportunities Program11. Here, we used hospital admission records in the database of hospital services (Automated Subsystem of Hospital Discharges; SAEH) from the Ministry of Health in Mexico, which typically cover public institutions for those without social security/insurance plans11, to explore the extent of the burden imposed by dengue on hospital services in the country. Specifically, we compared dengue-associated hospital admissions with other potentially preventable hospitalizations (PPHs)12, assessed the impact of dengue activity on in-hospital mortality from other morbid conditions, and explored the link between dengue activity and dengue-related mortality.
METHODS
Data sources
SAEH contains anonymized information on cases regarding reasons for hospitalization, diagnostic and treatment procedures, patient outcome, reasons for discharge, and duration of hospital stay for all causes, grouped according to the International Statistical Classification of Diseases and Related Health Problems version 10 (ICD-10)13. For the period assessed (2008-2014), the SAEH database included data on 19.2 million hospital admission records from 3370 departments in 817 hospitals. The Ministry of Health hospitals provided 38.3% of the hospital services in the country14.
Dengue-related hospitalizations were identified using ICD codes for dengue cases: A90 (dengue) or A91 (dengue hemorrhagic fever). The selection of other PPHs for comparison to dengue was based on those proposed by the Australian Institute of Health and Welfare12. ICD codes of conditions selected are shown in supplementary information Table S1. The Australian codes were chosen as these are the most referenced in the literature, are updated annually to ensure reliability, and include vaccine PPHs (VPPHs) (Table S1)15.
The duration of hospital stay documented for each admission was used to build daily tables, per hospital departments, with calculated metrics as follows: (1) the daily dengue patient load in the hospital department was defined as the number of beds attributed to dengue according to the ICD10 classification on admission; (2) the bed capacity on any given day, calculated as the maximum number of daily beds occupied over the past 180 days in a given hospital department; and (3) the relative daily dengue patient load calculated as the ratio of the former two metrics, i.e., the percentage of beds attributed to dengue on any given day with respect to the estimated bed capacity of the hospital department.
There were 613 small departments that were unclassified (and excluded from further analysis) as they recorded very low (< 200) total admissions over the period 2009-2014. Hospital departments were categorized into high dengue activity (HDA), low dengue activity (LDA), or zero dengue activity (ZDA) based on dengue-related hospital admission data from the years 2009 to 2014. HDA departments were defined as those with at least 200 dengue admissions and a maximum relative daily dengue patient load exceeding 20% of the estimated bed capacity of the department. ZDA departments were those with no recorded dengue admissions. The LDA departments were by definition the remaining departments (neither HDA nor ZDA) and included those hospital departments with < 200 total dengue admissions (but with at least one dengue admission) or a maximum relative daily dengue patient load <20% of the estimated bed capacity of the department over 2009-2014.
Outcome measures
The main outcome measures were dengue hospital admissions (where dengue was considered the principal or secondary diagnosis) in the SAEH hospital database from 2008 to 2014, by age group, severity (ICD codes for admission were A90 [dengue] or A91 [dengue hemorrhagic fever]), whether a death occurred, and by category of hospital department (HDA, LDA). Dengue-related hospital admissions (i.e., where dengue was considered the principal diagnosis) were compared to admissions due to other PPHs or chronic conditions. The impact of dengue on general (i.e., non-dengue) in-hospital mortality of a hospitalized episode was explored using the dengue patient load (see definition below).
Statistical analysis
Analytics platforms used were MySQL Workbench 6.3.6 for storage16, and KNIME 3.4.2 and R 3.3.217 for analysis of the data. An ordinary logistic regression model18 was generated based on all documented admissions to HDA departments from 2009 to 2014 (553,000 episodes) and used to quantify if department dengue patient load was a risk factor for in-hospital mortality. The model was developed to take into account possible confounding effects: age group, principal diagnosis defined according to the ICD1013, department characteristics (size, specialty, and intensive care unit), and seasonality (admission year and quarter). The dengue patient load and department size factors assigned to each hospitalized episode were the average values during the episode determined from the admission and discharge dates.
RESULTS
Hospital department data
Overall, 68,000 hospitalized dengue admissions, where the disease was considered the principal or secondary diagnosis (A90 and A91), were included in the current analysis. There were 79 HDA departments, managing 55% of dengue admissions, 923 LDA departments managing 45% of dengue admissions, and 1755 ZDA departments from 2009 to 2014 (Fig. 1).
Dengue activity compared to other PPHs
Weekly dengue admissions where the disease was considered the principal diagnosis (A90 and A91) in the SAEH hospital database are shown in figure 2. There was a clear peak in dengue activity in the last quarter of each calendar year. The last quarters of 2009, 2012, and 2013 had the largest peak in weekly dengue admissions in the SAEH hospital database. Results from the SAEH database from 2008 to 2014 suggested that dengue consumed more hospital resources than other VPPHs (Fig. 3a). In 2013, dengue represented 87% of VPPHs (Fig. 3b-i; ranked 1 vs. other VPPHs), and dengue patients utilized more hospital beds than prevalent chronic conditions (ranked 2 vs. chronic conditions) such as chronic obstructive pulmonary disease, hypertension, or asthma, although not as many as diabetes (Fig. 3b-ii). Dengue patients occupied as many beds as prevalent acute conditions such as urinary tract infections (Fig. 3b-iii; ranked joint first with urinary tract infections), and dengue patients occupied more hospital resources than any other infectious disease (ranked 2 vs. other infectious diseases), except for intestinal diseases (Fig. 3b-iv).
In-hospital mortality related to dengue patient load
Dengue daily load was shown to be a significant risk factor for overall in-hospital mortality (i.e., non-dengue-related mortality) in HDA departments. Multivariate logistic regression analysis indicated that high dengue patient load, older patient age, and larger department size increased the odds of in-hospital mortality in HDA departments (Table S2). After controlling for possible confounding factors, the odds of in-hospital death increased significantly and monotonically with increasing dengue patient load (Table 1). Mortality odds increased by 24.5% (20.0-29.1%, p < 0.001) for dengue patient load > 3.
Dengue CFR related to hospital department category
Figure 4 shows the dengue CFR rates stratified by age, HDA and LDA departments, and according to the dengue severity (A90, A91). Higher dengue CFRs were observed in LDA than HDA departments irrespective of dengue severity. For the study period, 2009-2014, the average A90 CFRs were 0.35% for LDA and 0.18% for HDA hospital departments, and the average A91 CFRs were 1.14% and 0.55%, respectively. For all ages, there was an approximately two-fold higher dengue CFR in LDA than HDA hospital departments; this difference was greater in the ≥ 60 years of age population irrespective of dengue severity.
DISCUSSION
This study examined a large data set derived from a total of 3370 hospital departments (in 817 hospitals) across Mexico. A total of 68,000 hospitalized dengue cases were included in the analysis as well as the remaining 19.2 million admissions for comparison. Weekly data of dengue hospital admissions analyzed from the Mexican SAEH database showed peaks in dengue activity in the second half of each calendar year between 2008 and 2014. This seasonal pattern of dengue disease, with most cases occurring in the second half of the year, has previously been reported in Mexico19. The yearly variation in dengue incidence may be related to seasonal fluctuations in climate parameters influencing dengue transmission cycles; for example, an increase in sea-surface temperature, minimum temperature, and rainfall was shown to be associated with increased reported cases of dengue in the coastal municipalities of the Gulf of Mexico20. However, despite the seasonal pattern, our study showed that a number of dengue cases continued to be reported throughout the year.
As dengue admission comprised a small proportion of all hospitalizations, we chose to initially compare dengue with other preventable diseases. We found that dengue accounted for the majority of VPPHs in our study; 72% of VPPHs during our study period were due to dengue, which increased to 82% in 2013, the year with the largest peak in weekly dengue admissions. Moreover, dengue admissions were usually ranked either first or second most frequent cause of PPHs as a proportion of acute conditions, chronic conditions, and infectious diseases. These results show that dengue imposes a significant burden on hospital resources in Mexico and reflects the positive impact that the national immunization strategy has had on the public health of Mexico21.
In dengue epidemic years, the burden may be severe enough to have an adverse impact on the health-care systems performance, as physicians dedicate less time to their patients, and demand increases for laboratory resources22, which would also indirectly affect the quality of care for other illnesses/conditions and their outcomes. Indeed, a recent study suggested that dengue outbreaks had the potential to undermine medical, social, economic, and political sectors in the wider society23. Such outbreaks are generally characterized by a rapid patient influx into the health-care system which may quickly become overwhelmed, as well as placing a financial burden on governments and households that could lead to wider political and social disruption and disengagement on outbreak control measures. Other potential negative results of dengue outbreaks, as seen in other epidemics, are the reductions in admissions related with non-urgent procedures and the reduction in the availability of general hospital resources24,25. The seasonal nature of dengue may also have an impact on the management of dengue hospitalization; it has previously been observed that a larger proportion of cases was managed in LDA departments during dengue epidemic periods26. In our study, since the CFR for dengue patients was higher in LDA departments, the fact that these departments receive more patients during outbreaks implies a higher mortality risk for those dengue patients they treat. We also confirm the disruption due to dengue on the quality of care as a whole, with dengue patient load being a significant risk factor for overall non-dengue-related in-hospital mortality in HDA departments; the odds of in-hospital death increased by 24.5% for dengue patient loads > 3. We hypothesize that the worse prognosis for patients in LDA departments is related to the lack of experience in the management of severe or evolving dengue disease. To improve the prognosis, a training program for the management of dengue should be mandatory as well as compliance with management guidelines.
We also observed dengue-related CFRs for severe cases (A91) of 0.39% and 1.09% in HDA and LDA settings, respectively, which are higher than those reported for Mexico (0.02-0.32% CFR) published by PAHO (official figures were based on severe cases reported by all the health units)27. Although there have been important variations in the dengue mortality rates in Mexico over the last couple of decades, the mortality by age group has previously been shown to be mainly toward the extremes of life, with high mortality rates in that aged < 1 year and ≥ 60 years28. Moreover, dengue mortality has consistently been highest in those aged ≥ 60 years in the three decades to 2009. Our data confirm the high mortality rate in those aged ≥ 60 years and was irrespective of disease severity reported. The higher mortality rate in those aged ≥ 60 years may be related to the more atypical disease presentation which can delay diagnosis and lifesaving interventions, combined with the higher likelihood of underlying comorbidities29. The high prevalence of comorbidities in Mexico must be highlighted. Recent models have proposed that, for 2020, 26.3% of the population 50 years and older will have been medically diagnosed with diabetes, for instance30.
We observed for all age groups an approximately two-fold higher dengue CFR in LDA than HDA hospital departments; this difference was more pronounced among those aged ≥ 60 years irrespective of dengue severity (A90 and A91 dengue admissions). Dengue-related CFRs were also higher in LDA than in HDA departments in another similar study undertaken in Brazil26. There are a number of reasons why dengue-related CFRs may be higher in LDA departments. These hospital departments may have less expertise/experience in the care of patients with dengue. There is also evidence to suggest that hospitalized dengue patients treated by dedicated medical teams experienced in the clinical management of dengue would substantially improve dengue-related quality of care and clinical outcome31. Other reasons may include that LDA departments might be less prepared or equipped to provide adequate diagnostic and supportive care to dengue patients32, or lack clinical experience and adequately trained staff in the management of dengue symptoms31,33.
Overall, these results support the case for preventative measures (e.g., vector control, disease awareness/protection, and dengue vaccine) in limiting the extent of any potential dengue outbreak. However, the use of the dengue vaccine appears to confer protection among dengue-seropositive individuals but predisposes those dengue-seronegative to severe dengue34. As such, the WHO Strategic Advisory Group of Experts on immunization preferred approach to dengue vaccination in endemic settings is to screen for prior dengue infection, to ensure that only those with previous dengue infection are vaccinated35.
To the best of our knowledge, this is the largest retrospective analysis conducted in Mexico exploring dengue disease and its potential influence on hospital admissions and non-dengue patient mortality. We used the SAEH database over 7 years of continuous hospital activity, which represents a robust data set over a substantial time period. In addition, the quality of records examined was good; 99% of the admissions analyzed had complete records. The use of real-world evidence has enabled us to further characterize the true burden of dengue on hospital practice as well as on direct and indirect dengue-related mortality. However, there were some limitations to this study; the databases used captured only mortality in public hospitals run by the Mexican Ministry of Health, and so these findings may not necessarily be applicable across private hospitals or to the broader Latin American population. We used dengue admissions and mortality as indicators of disease burden, although most of the burden is associated with non-hospitalized and nonfatal cases9,36. In addition, dengue diagnosis may be challenging as a number of disease-causing pathogens (such as Zika or Chikungunya) may cause symptoms similar to dengue that are difficult to differentiate in the absence of laboratory confirmation37. The logistic regression did not account for clustering effects at the hospital or regional levels, which could lead to more accurate comparisons of the potential drivers of overall mortality at a fine geographical level. There also remain gaps in the epidemiological knowledge with regard to local serotype distribution in Mexico38.
In summary, our study confirms the important burden due to dengue in Mexican hospitals when compared to other preventable chronic and acute hospitalizations, as well as preventable hospitalizations due to other infectious diseases. The study also highlights the impact of dengue activity not only on the dengue case-fatality rate but also on the overall in-hospital mortality. It is hoped these data will help inform policy regarding resource allocation to improve the quality of care and clinical outcome of all hospitalized patients during periods of high dengue activity.