ABSTRACT
Mosquito-borne viruses pose a perpetual public health threat to countries and territories in the Carribean due to the region’s tropical climate and seasonal reception of international tourists. Outbreaks of the emerging viruses chikungunya and Zika in 2014 and 2016, respectively, demonstrated the rapidity with which these viruses can spread between islands. At the same time, the number of reported dengue fever cases, caused by the endemic dengue virus, has steadily climbed over the past decade, and a large dengue outbreak that began sweeping through this region in 2019 continues in 2020. Sustainable disease and mosquito control measures are urgently needed to quell virus transmission in the long term and prevent future outbreaks from occurring. To improve upon current surveillance methods, we analyzed temporal and spatial patterns of chikungunya, Zika, and dengue outbreaks reported in the Dominican Republic between 2012 and 2018. The viruses that caused these outbreaks are transmitted by Aedes mosquitoes, which are sensitive to seasonal climatological variability. In this study, we evaluated whether climate and the spatio-temporal dynamics of past dengue outbreaks could inform when and where future emerging disease outbreaks might occur. We found that the temporal and spatial distribution of emerging disease outbreaks did not conform to those of seasonal dengue outbreaks. Rather, the former occurred when climatological conditions were suboptimal for Aedes activity. Provincial dengue attack rates did not correspond to those of emerging diseases. Our study also provides evidence for under-reporting of dengue cases, especially following the 2016 Zika outbreak. We advocate for the implementation of a sustainable and long-term surveillance system to monitor the spread of known mosquito-borne viruses and to identify emerging threats before they cause outbreaks. Specifically, we recommend the use of febrile illness incidence, ca se fatality rates, and serosurveys during inter-outbreak periods to better understand rates of transmission and asymptomatic infection.
INTRODUCTION
Emerging and endemic mosquito-borne viruses are a constant public health concern in the Carribean (Cao-Lormeau, 2016; Mavian et al., 2018; Musso et al., 2018). This region is especially vulnerable to the spread of the former due to its tropical climate and large tourism industry, which attracts visitors from all across the globe. The recent outbreaks of chikungunya (Cassadou et al., 2014; Dorléans et al., 2018; Halstead, 2015) and Zika (Faria et al., 2016b; Grubaugh et al., 2018; Metsky et al., 2017) in 2014 and 2016, respectively, demonstrated that viruses once associated with mild illness can re-emerge and cause devastating health outcomes. Reports that Mayaro virus, which has already been detected in the Caribbean, can be transmitted by the urbanized mosquito species Aedes indicate that future outbreaks of emerging mosquito-borne diseases may be on the horizon (Ganjian & Riviere-Cinnamond, 2020; Hotez & Murray, 2017; Long et al., 2011; Lorenz et al., 2019; Weaver & Reisen, 2010). In addition to this threat, dengue virus is endemic to many Caribbean countries and territories and has caused outbreaks with increased frequency over the past decade. Large outbreaks, which began in 2019 but have continued through 2020 (Pan-American Health Organization (PAHO), 2020a), extend an alarming trend of a rising number of dengue cases reported annually in the Americas in recent years (PAHO, 2014). However, the danger of these viruses and the diseases they cause lie not only in their debilitating and sometimes life-threatening symptoms. The unpredictability of when and where a new outbreak will occur precludes preparedness. Outbreak response strategies are inherently reactionary and, in their transience, disrupt public health systems when new initiatives are introduced and subsequently phased out. This is especially problematic in resource-limited settings where the strategic allocation of resources should be prioritized to maximize the impact of disease control efforts. New approaches centered around sustainable, long-term surveillance are needed to curtail the potential for future public health crises in the Caribbean.
One such approach is the use of climate data to evaluate the risk of viral spread. The transmission of mosquito-borne viruses, both emerging and endemic, should adhere to similar temporal and spatial patterns such that the dynamics of past outbreaks can inform those of future outbreaks. Many flaviviruses, like dengue and Zika viruses, and alphaviruses, like chikungunya virus, are transmitted by Aedes mosquitoes, which are sensitive to climatological variability (Kraemer et al., 2019; Kraemer, et al., 2015; Obolski et al., 2019). Temporal concordance between seasonal weather patterns and transmission of yellow fever virus (flavivirus) has been documented on the African continent (Hamlet et al., 2018). Similar associations have been reported for dengue fever incidence in Hanoi, Vietnam and in Myanmar (Do et al., 2014; Perez-Guzman et al., 2018) and for chikungunya and Zika incidence in the Americas (Obolski et al., 2019). However, there is evidence to suggest that this relationship does not extend to the transmission dynamics of viruses during their first year of introduction into a new population (Lourenço et al., 2017; Faria et al., 2016). Therefore, we considered whether, while climate may be a useful indicator for future endemic virus outbreaks, other factors including population size, demographics, and the timing of introduction should be considered when developing strategies to prevent future emerging disease outbreaks.
To answer these questions in the context of the Caribbean, we analyzed dengue, chikungunya, and Zika cases reported daily in the Dominican Republic between 2012 and 2018. We found that emerging disease outbreaks (chikungunya and Zika) occurred earlier in the year than dengue outbreaks, and the timing and location of introductions of emerging viruses impacted when and where corresponding outbreaks occurred. Moreover, the spread of chikungunya and Zika viruses was tolerant to sub-optimal climates for transmission by Aedes mosquitoes, likely due to the large size of the susceptible human population. Predicted mosquito abundance was similarly uninformative for the spatial distribution of emerging disease attack rates and force of infection. Moreover, provincial-level dengue attack rates were consistent between dengue outbreaks, but they did not correspond to local attack rates of chikungunya and Zika. Taken together, we demonstrate that dengue virus may not always be an appropriate model to prepare for future emerging mosquito-borne disease outbreaks. Instead, a sustainable and long-term mosquito-borne disease surveillance system is needed to facilitate proactive responses to emerging threats and to track the continued spread of known diseases including dengue, chikungunya, and Zika. We specifically propose the use of two indicators, incidence of febrile illness cases and dengue case fatality rates, to monitor surveillance performance and identify potential emerging threats.
RESULTS
Multiple endemic and emerging mosquito-borne virus outbreaks in the Dominican Republic
Between 2012 and 2018, the Dominican Republic, a country that shares the Caribbean island La Hispaniola with Haiti (Fig. 1a), experienced five disease outbreaks caused by mosquito-borne viruses (Fig. 1b,c). We delineated this period into six seasons (Seasons 1-6), each beginning in April of every year between 2012 - 2018, coincident with the start of the rainy season. Three of the five outbreaks were caused by dengue virus (Seasons 1, 2, and 4; Fig. 1b,c). Sufficient serosurvey data for the Dominican Republic are not available to determine the predominant dengue virus serotype(s) during each outbreak, but the Pan-American Health Organization reports that serotypes 1, 2, and 4 were prevalent in the Caribbean between 2012 and 2014, and serotypes 2, 3, and 4 were circulating in the region in 2015 (PAHO, 2020b). The other two outbreaks were caused by emerging viruses, chikungunya and Zika (Seasons 3 and 5, respectively; Fig. 1b,c) (Bowman et al., 2018). Those outbreaks were the first and only to be reported of either disease in the country, and the number of reported cases in the Caribbean as a whole plummeted in subsequent years (Lindsey et al., 2018; Zingman et al., 2017; PAHO, 2020c). After the season dominated by Zika virus (Season 5), only a small number of dengue cases were diagnosed (Season 6).
Clinical and demographic characteristics differed between reported dengue, chikungunya, and Zika cases. The majority of dengue cases experienced fever (98.3%), were hospitalized for their condition (93.7%), and were between the ages of 0 and 15 (62.3%), suggesting immunity in the adult population (Table 1). Patients diagnosed with chikungunya and Zika were older and predominantly female. Rash and/or arthralgia in conjunction with fever are typical clinical manifestations of chikungunya and Zika infections (Duffy et al., 2009; He et al., 2017; Langsjoen et al., 2016; Martinez et al., 2019; Pineda et al., 2016; Rosario et al., 2015). Although lower than rates reported elsewhere, arthralgia rates among chikungunya cases in our dataset were consistent with previously documented rates among individuals diagnosed with chikungunya in the Dominincan Republic (Langsjoen et al., 2016).
Chikungunya and Zika outbreak dynamics did not conform to seasonal dengue patterns
Implementing effective disease prevention and control measures requires knowledge of when their etiological viruses are most likely to emerge, re-emerge, and spread (Grubaugh, Ladner, et al., 2019). Dengue, chikungunya, and Zika are caused by RNA viruses transmitted by the mosquito vectors Aedes aegypti and Aedes albopictus. The abundance and capacity of these vectors are sensitive to climatological factors including temperature, rainfall, and humidity (Kraemer et al., 2019; Kraemer, et al., 2015), such that transmission of all three viruses has been shown to fluctuate in kind (Do et al., 2014; Hamlet et al., 2018; Obolski et al., 2019; Perez-Guzman et al., 2018; Faria et al., 2016a). We therefore hypothesized that the timing of the five outbreaks should concord with seasonal weather patterns.
We found that the emerging disease outbreaks (chikungunya [Season 3] and Zika [Season 5]) occurred earlier in the year compared to endemic outbreaks (dengue). The chikungunya and Zika outbreaks peaked 15 and 26 weeks earlier, respectively, than the averaged dengue peak (epidemiological week 41; Fig. 1b). While both chikungunya and dengue cases began to rise around epidemiological week 15, a rise in Zika cases was observed during the first weeks of the year.
Because the implementation of national-level reporting of chikungunya and Zika cases influenced when these outbreaks were detected (Fig 1c, arrows), we calculated the number of weeks elapsed between the first reported case of each outbreak and the peak number of cases within provinces (Fig. 1d). We assumed that climatological factors did not vary widely between provinces during a given season. Our reasoning followed that the rate of viral transmission is limited by the extrinsic incubation period (EIP) of the virus in the mosquito vector (Liu-Helmersson et al., 2014), and the EIP for Aedes mosquitoes is mainly influenced by temperature (Chan & Johansson, 2012; Mordecai et al., 2017; Winokur et al., 2020). Therefore, because we did not expect transmission rates to vary widely between provinces, provincial outbreaks that peaked soon after the reporting system was implemented would indicate that transmission in those provinces preceded the establishment of this system.
Our analysis identified one such instance. The majority of provincial chikungunya and Zika outbreaks peaked within 12 and 25 weeks, and the mean time to peak across provinces did not significantly differ (Fig. 1d,e). This was in contrast to the substantial heterogeneity of the timing of the three dengue outbreaks. The lack of uniformity across these outbreaks may have been due to a slower rate of spread between regions because of a pre-existing and spatially heterogeneous herd immunity landscape. Five provinces, four of which are located in the south-western part of the country, reported peak numbers of Zika cases in January and February of 2016 (Fig. 1d). For the reasons stated above and because Zika virus had been circulating in the Americas for at least two years before it was reported in the Dominican Republic (Faria et al., 2016, 2017), this observation suggested that Zika virus was already circulating in those provinces before the national reporting system was implemented at the beginning of January. We hypothesized that early Zika virus transmission could have occurred in the south-western region because of the region’s proximity to Haiti. In particular, this region shares a border with Haiti, is connected to the Haitian capital Port-au-Prince via a major roadway, and is home to a large binational market (Fig. 1e, map, star). A temporal comparison of the Zika outbreaks reported in this region, the rest of the country, and Haiti revealed that the regional outbreak was more consistent with the 2016 Zika outbreak reported in Haiti (Fig. 1e, inset). These findings indicate that the Dominican Republic experienced two geographically and temporally distinct Zika outbreaks. From reported case counts alone we cannot conclude whether Zika virus was introduced into the Dominican Republic from Haiti or vice versa; however, it is clear that binational coordination is an essential component of local mosquito-borne disease control because the viruses these vectors transmit do not recognize international borders.
Initial outbreaks of emerging mosquito-borne diseases can occur during periods of sub-optimal climatological conditions
We considered two possible explanations for the asynchronicity we observed between outbreaks of the emerging and endemic viruses (Fig. 1). (1) Either climate patterns differed between seasons and the emerging outbreaks preceded those of dengue due to seasonal stochasticity, or (2) the spread of emerging viruses was less sensitive to climatological factors compared to dengue virus. By analyzing daily climatological data collected over the duration of our study period, we found support for the latter (Fig. 2).
To investigate the relationship between climate and case incidence, we used temperature and humidity time series data to estimate the mosquito-borne transmission potential throughout our study period (Fig. 2). For this analysis, we used Index P, a metric that is calculated by incorporating climate and entomological data into a Bayesian framework to estimate the transmission potential of individual female mosquitoes (Obolski et al., 2019). We calculated transmission potential for each week of our study period using data reported in six of thirty-two provinces. We selected these six provinces based on their representativeness of the country’s sub-climates and data availability. On average, transmission potential fluctuated seasonally, rising between April and November of each year, coincident with the rainy season, and fell shortly thereafter (Fig. 2). When we compared the temporal dynamics of transmission potential to reported disease incidence, we found that the number of emerging disease cases reported weekly peaked before transmission potential had reached a seasonal maximum for both outbreaks (chikungunya and Zika), while the number of dengue cases reported weekly peaked after this point for all three dengue outbreaks (Fig. 2a). Similar climatological patterns persisted on the provincial level (Fig. 2b). Some provinces experienced clear seasonal fluctuations in transmission potential, whereas seasonal peaks and troughs were less defined in Distrito Nacional and María Trinidad Sánchez, which have low to moderately humid sub-climates. Despite this variation, the timing of provincial outbreaks conformed to the trend observed on the national level: dengue outbreaks peaked after a period characterized by high transmission potential, just as transmission potential was beginning to wane; in contrast emerging disease outbreaks occurred concurrently with increasing transmission potential. Moreover, we cannot attribute climatic anomalies to the decline in reported dengue cases in Seasons 3, 5, and 6.
The sizes of the susceptible human population influenced the speed at which emerging viruses spread
Variability in seasonal weather patterns and vectorial capacity did not account for differences in the timing of emerging disease outbreaks (Fig. 2). Rather, we hypothesized that the patterns observed for dengue were likely to have been influenced by a pre-existing and spatially heterogeneous herd immunity landscape. This would mandate that transmission potential must remain high for an extended period before epidemic growth is achieved each season.
To test the hypothesis that local susceptibility influenced the timing of the outbreaks, we examined the relationship between the size of the susceptible human population and each outbreak’s doubling time (Fig. 3). We estimated that the median estimated effective reproduction number (Reff) for the three dengue outbreaks was between 1.3 and 1.4, while those of the chikungunya and Zika outbreaks ranged from 1.6 to 2.45 (Fig. 3a). These findings are consistent with the epidemiology of the three diseases and previous estimates of Reff for chikungunya and Zika (Bowman et al., 2018; Cauchemez et al., 2014; Nishiura et al., 2016). Published estimates of the basic reproduction number (R0) for dengue, equivalent to Reff in a wholly susceptible population, vary widely (0.97-65) due to the parameter’s sensitivity to the size of the local susceptible population (Liu et al., 2020). We observed a slightly longer doubling time (smaller Reff) during the third dengue outbreak, consistent with increasing levels of herd immunity from the two previous outbreaks. In our dataset, reported dengue cases were primarily in the 0-15 age group, indicating older age groups had high levels of pre-existing immunity and were not susceptible to disease (Table 1; Fig. 3b). In contrast, we speculated that the entire population was susceptible to chikungunya and Zika, allowing these viruses to spread despite sub-optimal weather conditions, facilitating an earlier outbreak peak, and affecting a much wider range of ages (Table 1; Fig. 3b).
Spatial distribution of dengue attack rates is a poor indicator for initial outbreaks of emerging disease
Given that the timing of emerging disease outbreaks (chikungunya and Zika) did not conform to that of dengue outbreaks (Fig. 1,2), we suspected that the spatial distribution of dengue cases would be equally uninformative for chikungunya and Zika outbreaks. Specifically, we hypothesized that the relative burden of dengue within provinces during individual outbreaks would correlate with the burden of dengue during subsequent outbreaks but not with the relative burden of emerging disease cases.
To address this question, we measured the age- and sex-adjusted attack rates by province for each of the outbreaks. We found that the attack rates for individual provinces across outbreaks were well correlated between dengue outbreaks and between chikungunya and Zika outbreaks, but not across the endemic and emerging viruses (Fig. 4a, b).
Next, we investigated the role of climate and land-use in perpetuating this trend, reasoning that larger mosquito populations would facilitate higher attack rates (Kraemer, et al., 2015; Rodriguez-Barraquer, Salje, et al., 2019). To this end, we compared an Aedes aegypti suitability score (AaS; Fig. 4c), a metric that incorporates ecological variables not included in our transmission potential estimates such as vegetation levels and precipitation, to population-adjusted attack rates of disease (Fig. S1, Fig. 4a). Also unlike our transmission potential estimates, AaS is an aggregated estimate of suitability per month in any given year. Interestingly, AaS was not well correlated with the adjusted provincial attack rates or with mean age of infection (Fig. S1). Moreover, provinces in which the highest levels of mosquito abundance were predicted did not consistently report the highest burden of disease (Fig. 4a,c). One reason for these inconsistencies is that AaS relies on historical climate data from 1970 to 2000 to estimate mosquito abundance. Therefore, our estimates do not capture anomalous climatological events such as the hurricane that made landfall in the Dominican Republic in October 2012. AaS also assumes a high level of stability in land use and minimal urbanization since the thirty year period from which the estimates were made. Regardless of the true underlying reason for these discrepancies, the spatial distributions between endemic (dengue) and emerging (chikungunya and Zika) mosquito-borne virus diseases are not very well correlated.
Clinical outcomes of dengue cases and trends in febrile illness cases are indicators of disease prevalence
Up to this point, our findings demonstrate that emerging mosquito-borne disease outbreaks are not confined to the temporal and spatial patterns of endemic disease outbreaks in the Dominican Republic (Figs. 1-4). While similar inconsistencies in outbreak dynamics have been observed elsewhere in the Americas (Fig. S2) (Lourenço et al., 2017; Faria et al., 2016a), it is difficult to discern if our findings are due to differences in epidemiology or in underlying reporting biases. Therefore, we used general patterns of febrile illness and clinical characteristics to show that reporting biases did not influence our previous conclusions (Fig. 5).
Because the vast majority of cases in our dataset reported fever independent of disease (Table 1), we hypothesized that febrile illness incidence should reflect disease incidence reported during our study period. When we examined the number of febrile illness cases reported in the Dominican Republic per epidemiological week per season, we found that temporal trends in febrile illness cases in Season 3 (2014-2015) and Season 5 (2016-2017; Fig. 5a) were consistent with those we observed in our chikungunya and Zika case data, respectively (Fig. 1a). In Season 4 (2015-2016, dengue), febrile illness cases peaked in week 44, coincident with the peak of the dengue outbreak (Fig. 1a). To ensure that the spatial relationships we identified in Figure 4 were not purely the product of differential reporting practices, we compared population-adjusted attack rates of mosquito-borne disease and febrile illness reported by individual provinces during each season and found no significant correlation (Fig. S3). Equally pronounced seasonal peaks of reported febrile illness were not observed in Season 1 and Season 2 despite the known dengue outbreaks that occurred within that period. We speculate that changes made in reporting protocols for febrile illness cases during the chikungunya outbreak in Season 3 prompted a sustained increase in national reporting thereafter (Ministerio de Salud Pública, 2014a).
We then investigated why reported febrile illness cases peaked seasonally in Seasons 5 and 6 despite an apparent decline in coincident dengue cases by comparing case-fatality rates (CFRs) of dengue cases across seasons. Specifically, we considered the possibility that this decline was due to under-reporting of dengue cases during the interim periods between dengue outbreaks (Seasons 3, 5, and 6). The average CFR during these seasons was 1.87% (range: 1.60%-2.40%), while the average CFR of reported cases during the three large dengue outbreaks in Seasons 1, 2, and 4 was 0.8% (range: 0.48%-1.20%) (Fig. 5b). A wide range of dengue CFRs is reported in the literature, with a mean of 1.62% (95% CI: 0.64–4.02%) (Low et al., 2018), and without serosurveillance data, we are unable to rule out the possibility that the circulation of different dengue serotypes resulted in differences in seasonal CFRs. However, if we assume that fatality among cases during the third dengue outbreak (Season 4) were well reported, dengue cases were under-reported by 33% in Season 3 and Season 5, and by 2-fold in Season 6. These are likely conservative estimates as the current surveillance system is passive and therefore does not capture asymptomatic cases. As a result, Seasons 3, 5, and 6 may represent periods characterized by a large number of mild cases that were not captured by the surveillance system. It is also likely that the number of deaths due to chikungunya was substantially under-reported; however, our calculation of 0.19% corresponds to a corrected, post-hoc estimate of this rate (0.15%) for the Dominican Republic (Freitas et al., 2018).
DISCUSSION
Our study demonstrates that, even when transmitted by the same mosquito vector, viruses are not beholden to the same temporal and spatial outbreak dynamics. Instead, when and where the new virus is introduced, the size of the susceptible human population, and the capacity of local surveillance systems determine these dynamics are detected. In short, dengue epidemiology cannot be used to anticipate the location and timing of future emerging mosquito-borne disease outbreaks in the Dominican Republic, and likely in other Caribbean countries and territories. Instead, consistent and sustainable surveillance methods should be implemented to limit disease and prevent future outbreaks. These methods could include serosurveil lance of the population during periods between outbreaks (Imai et al., 2015), testing local mosquito vectors for viral infections (Fauver et al., 2017, 2018; Grubaugh et al., 2015), and monitoring health outcomes of travelers who visit the country (Grubaugh, et al., 2019).
Maintaining a sustainable surveillance system is critical for preventing the silent transmission of viruses that can fuel large outbreaks. Other countries in the Americas reported subsequent outbreaks of chikungunya and Zika after their initial outbreaks. We cannot conclusively determine whether the Dominican Republic experienced a similar pattern because surveillance data for chikungunya and Zika are not available for seasons following the initial outbreaks of these diseases. Elucidating whether dengue, chikungunya, and Zika are co-circulating in the country will be critical for triaging and providing appropriate clinical care to patients who present with febrile illness (Vogels et al., 2019), especially if chikungunya and Zika virus transmission is now in sync with dengue transmission (Bisanzio et al., 2018; L. P. Freitas et al., 2019).
Equally important is the identification of viruses that could spread across international borders. In our study, we observed an epidemiological link between Zika outbreaks in Haiti and the south-west region of the Dominincan Republic. This connection likely extends to other mosquito-borne viruses such as Mayaro virus and West Nile virus, which have been reported in Haiti (Blohm et al., 2019; Lednicky et al., 2016; Weppelmann et al., 2017). The former can be transmitted by Aedes mosquitoes (Kantor et al., 2019; Long et al., 2011; Pereira et al., 2020; Wiggins et al., 2018), and the latter has been isolated from zoonotic reservoirs in the Dominican Republic (Komar et al., 2005). Similarly, reports of international travelers carrying mosquito-borne viruses from the Dominican Republic back to their home countries demonstrate that gaps in surveillance have global implications (Barzon et al., 2016; Díaz-Menéndez et al., 2018; Duijster et al., 2016; Goncé et al., 2018; Millman et al., 2016; Requena-Méndez et al., 2014). Therefore, understanding the transmission patterns of viruses and developing a unified, international plan to combat them before they cause an outbreak will help mitigate the potential for such an event.
Our findings suggest that dengue cases were under-reported following the Zika outbreak in 2016; however, there are a number of possible explanations for the ostensible decline in dengue cases that should be explored. Wide-spread mosquito control measures motivated by the Zika outbreak could have limited the spread of dengue later that year. Although plausible, this line of reasoning does not explain why a similar post-Zika decline was observed in other countries in the Americas (Perez et al., 2019), nor why there was a resurgence of dengue cases in 2019. A second explanation is that Zika infections confer some level of temporary immunity to subsequent dengue infections (Gordon et al., 2019; Ribeiro et al., 2018, 2020; Rodriguez-Barraquer, Costa, et al., 2019). This theory cannot account for the small number of dengue cases reported in 2014 following the chikungunya outbreak, as the etiological agent of that disease is an alphavirus, and it assumes very high attack rates and extensive under-reporting of Zika in 2016 to have achieved sufficient levels of herd immunity. Our findings instead suggest that prior Zika infection protects against symptomatic dengue infections because we observed significant positive correlation between dengue attack rates within provinces across outbreaks. Given that cross-reacting immunity between dengue serotypes is well documented (Katzelnick et al., 2016; Tsang et al., 2019), our data suggest that a similar relationship between dengue and Zika would not result in a wide-spread decline of cases. Rather, dengue transmission could have reasonably persisted undetected if most of those infected were not hospitalized. If true, this hypothesis would explain why seasonal peaks in reported febrile illness cases persisted in 2017 and 2018, and why CFRs among reported dengue cases appeared to be elevated in Seasons 5 and 6 if the true number of cases was under-reported. To better understand these complex interactions, the collection of serotype information should be incorporated into current dengue surveillance efforts.
There are a few important limitations to our study. First, our dataset included chikungunya and Zika case data from the initial wave of each disease, and we cannot therefore compare temporal and spatial dynamics of these diseases across seasons. After these initial outbreaks, diagnostic testing for these diseases has largely ceased. While the number of cases of these diseases reported in the Dominican Republic has declined to zero the true burden of disease is unknown. Future studies should investigate whether these viruses have continued to circulate undetected in the country and whether their spatiotemporal dynamics have since synchronized with that of dengue virus. Second, the reporting system for suspected chikungunya cases differed from that used for suspected Zika cases. During the chikungunya outbreak, most febrile illnesses cases without apparent cause were initially classified as suspected chikungunya cases. For this reason, the number of cases reported by the Pan-American Health Organization (PAHO) and the Ministry of Health was significantly larger than those which we have reported here (Ministerio de Salud Pública, 2014b). Our chikungunya case data contains a disproportionate number of children in the <1 year age group, indicating that the dissemination of diagnostic testing may have been skewed towards high-risk groups. Third, our findings demonstrate that an epidemiological relationship existed between the Dominican Republic and Haiti during the Zika epidemic in 2016, but we cannot determine the directionality of cross - border virus movement without virus genomic data. However, given that this relationship exists and that mosquitoes do not recognize political boundaries, it can be assumed that bi-directional spillover of mosquito-borne diseases will occur in the future unless appropriate bi-national surveillance and control measures are implemented. Finally, our analysis primarily focused on virus transmission by Aedes aegypti mosquitoes, but it is possible that other mosquito vectors contributed to the propagation of the outbreaks we investigated. Specifically, Aedes albopictus may have played a key role in chikungunya transmission (Lamballerie et al., 2008; Vazeille et al., 2007). Limited data are available on the distribution of relevant mosquito vectors in the country (González et al., 2019, 2020), and broader entomological surveillance is needed to better address this question.
Taken together, our study demonstrates that surveillance for mosquito-borne diseases should be sustained during periods when transmission appears to be low because patterns in reported dengue cases are poor indicators of future emerging mosquito-borne virus outbreak dynamics. Reported symptoms and case demographics may be useful for identifying shifts in disease prevalence, but many clinical features, especially fever, that we have noted are likely a function of the reporting and diagnostic algorithms used during an outbreak. Active reporting of new dengue, chikungunya, and Zika cases and the broader deployment of diagnostics for newly emerged diseases are needed to ascertain more accurate case profiles. Outbreaks of emerging tropical diseases are a threat to the public health of the Carribean, and endemic diseases such as dengue precipitate health crises with increasing frequency. Given the pervasiveness of mosquito-borne diseases in tropical climates, sustainable surveillance systems rather than reactionary disease control measures should be implemented to prevent future crises.
MATERIALS AND METHODS
Description of data
Data for suspected and confirmed cases of dengue, chikungunya, and Zika reported between 2012 -2018 were extracted from the Dominican Republic Ministry of Health Weekly Reports (digepi.gob.do). Depersonalized demographic and clinical characteristics of cases were solicited from the National Statistics Directorate of the Ministry of Health (MoH) (saip.gob.do). Data were organized by date of onset of symptom report.
The National Epidemiology Directorate of the MOH collects dengue case reports by passive surveillance in a weekly manner in its Epi 1 form within their Sistema Nacional de Vigilancia Epidemiológica (SINAVE) Digital Platform for every healthcare setting across the country. These data were collected in weekly reports submitted to their website (digepi.gob.do) and organized by province and week of reported cases. The data collected by the MoH includes age, sex, province and municipality of residence, date of symptom onset, clinical outcome, and symptomatology. The MoH included mandatory reports from suspected and confirmed chikungunya infection through 2014 and Zika virus throughout 2016, without continuing to do so thereafter.
Province- and municipality-level weekly acute febrile illness data from 2012-2018 were solicited from the National Epidemiology Directorate of the Ministry of Health (saip.gob.do), which collects the data as part of a passive surveillance system. In 2014, the Ministry of Health increased its efforts to identify febrile illness cases (Ministerio de Salud Pública, 2014a).
Daily climatological variables for mean temperature and relative humidity for 5 cities (Santo Domingo, Distrito Nacional; Barahona, Barahona; La Union, Puerto Plata; Cabrera, María Trinidad Sánchez; Santiago, Santiago) was collected from the National Meteorology Office (ONAMET) database (saip.gob.do) from January 2012 through December 2018. Hourly climatological data for Bonao, Monseñor Nouel for the same timeframe were obtained from openweathermap.org.
Population data were extracted from the 9th National Population and Household Census. This census was conducted in 2010 by the National Statistics Office (National Statistics Office, 2012).
Reported cases of dengue, chikungunya, and Zika for Haiti, Bolivia, Jamaica, and Venezuela were extracted from the Pan-American Health Organization website and PLISA database (PAHO, 2020d). Dengue and Zika cases, which were reported by epidemiological week, were aggregated by month to allow for a direct comparison to chikungunya cases, which were reported by month.
Time to outbreak peak
We calculated the time to peak for provincial outbreaks by first identifying the epidemiological week in which the first case of the national outbreak was reported (chikungunya: EW 6, Zika: EW 1, dengue (3): EW 14) and then counting the number of weeks elapsed until each of the 32 provinces reported the maximum number of cases for the corresponding outbreak. The mean time to peak for each outbreak was compared using an unpaired T test implemented in Prism v8.4.2.
Transmission potential
We calculated weekly transmission potential (Index P) with the Bayesian approach developed by Obolski et al. (Obolski et al., 2019). Briefly, we extracted daily average temperature and relative humidity for 6 cities as described above. Rarely, temperature or humidity were not available for a given day. In these cases, we averaged the respective variable from the same date across the remaining 6 years. We did this for 14 days for Barahona and Puerto Plata, 3 days for Santiago, and 1 day for María Trinidad Sánchez.
We then used the R package MVSE and entomological and epidemiological priors documented in the literature (Table S1) to calculate daily transmission potential, which we then aggregated by week (Lourenço, 2019). We found that the model was reasonably robust to a range of priors for each of the parameters and therefore elected to use short human incubation and infectious period estimates to inform the model. We adopted our prior estimate of human life expectancy to the Dominican Republic based on estimates by the United Nations Development Program (UNDP, 2019), but otherwise used the same priors as reported by Oboloski et al. (Table S1) (Obolski et al., 2019).
Effective reproduction number (Reff) estimates
We estimated Reff for the five outbreaks using the method developed by Lipsitch et al. (Lipsitch et al., 2003) to fit a linear regression model to our case data (Eq. 1).
We used a range of values for latent and infectious periods that were well reported in the literature to calculate f, the proportion of the serial interval in the latent period, and v, the serial interval (Table S2). We obtained r, the epidemic growth rate, by extracting the slope of our model after fitting it to our case data. All calculations were done using R v4.0.0.
Adjusted attack rates and linear regression
Province-level attack rates by outbreak were age- and sex-adjusted using the direct standardization method, with the national population as the reference population.
We calculated Pearson’s R correlation coefficient for province-level attack rates between various pairs of outbreaks (Fig. 4a,b). An outlier analysis showed that the removal of outliers did not substantially affect the size or significance of the correlation coefficients, so we included all data points.
Aedes aegypti suitability score (AaS)
AaS was calculated for each month using both climate and land-use data on a 5km x 5 km scale collected from 1970 and 2000 and compiled by WorldClim as described by Kraemer et al. (Kraemer, et al., 2015). Monthly suitability scores were extracted and averaged by province in Rv4.0.0.
Data Availability
Case and febrile illness data were solicited from the Dominican Republic Ministry of Health (see link below). National population, age, and sex data are available through the Oficina Nacional de Estadistica (link below).
https://www.one.gob.do/demograficas/proyecciones-de-poblacion
COMPETING INTERESTS
The authors declare no competing interests.
ACKNOWLEDGEMENTS
We thank C. Vogels, A. Brito, J. Fauver, C. Kalinich, I. Ott, S. Lapidus, K. Gangavarapu, J. Pack, and S. Taylor for feedback and/or assistance with the methodology. The study was funded in part by the Hecht Global Health Faculty Network Award provided to N.D.G. from the Yale Institute of Global Health.