New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Region-wide synchrony and traveling waves of dengue across eight countries in Southeast Asia
Edited by Rino Rappuoli, GSK Vaccines, Siena, Italy, and approved August 25, 2015 (received for review January 21, 2015)

Significance
Persons living in the tropics and subtropics are at risk for dengue fever and dengue hemorrhagic fever, and large epidemics occur unexpectedly that can overburden healthcare systems. The spatial and temporal dynamics of dengue transmission are poorly understood, limiting disease control efforts. We compiled a large-scale dataset and analyzed continental-scale patterns of dengue in Southeast Asia. Our analysis shows that periods of elevated temperatures can drive the occurrence of synchronous dengue epidemics across the region. This multicountry collaborative study improved insight that may lead to improved prediction of dengue transmission patterns and more effective disease surveillance and control efforts.
Abstract
Dengue is a mosquito-transmitted virus infection that causes epidemics of febrile illness and hemorrhagic fever across the tropics and subtropics worldwide. Annual epidemics are commonly observed, but there is substantial spatiotemporal heterogeneity in intensity. A better understanding of this heterogeneity in dengue transmission could lead to improved epidemic prediction and disease control. Time series decomposition methods enable the isolation and study of temporal epidemic dynamics with a specific periodicity (e.g., annual cycles related to climatic drivers and multiannual cycles caused by dynamics in population immunity). We collected and analyzed up to 18 y of monthly dengue surveillance reports on a total of 3.5 million reported dengue cases from 273 provinces in eight countries in Southeast Asia, covering ∼107 km2. We detected strong patterns of synchronous dengue transmission across the entire region, most markedly during a period of high incidence in 1997–1998, which was followed by a period of extremely low incidence in 2001–2002. This synchrony in dengue incidence coincided with elevated temperatures throughout the region in 1997–1998 and the strongest El Niño episode of the century. Multiannual dengue cycles (2–5 y) were highly coherent with the Oceanic Niño Index, and synchrony of these cycles increased with temperature. We also detected localized traveling waves of multiannual dengue epidemic cycles in Thailand, Laos, and the Philippines that were dependent on temperature. This study reveals forcing mechanisms that drive synchronization of dengue epidemics on a continental scale across Southeast Asia.
Dengue virus (DENV) is an arbovirus transmitted by Aedes mosquitos in the tropics and subtropics of the world. The virus causes an estimated 390 million infections per year, resulting in 96 million clinically symptomatic cases (1). DENV has four serotypes (DENV-1, DENV-2, DENV-3, and DENV-4) that each circulate worldwide. The spatial propagation of dengue transmission at short distances by the mosquito vector is well-understood, but the mechanism of long-distance spread has remained unclear. Disease transmission over large geographical distances is difficult to measure directly, but epidemiological coupling of locations revealed by synchrony in population-level disease patterns has been used successfully in the past to infer mechanisms of spread (2⇓–4). For example, synchrony and its spatial hierarchies indicated that measles in the United Kingdom spread from urban centers to rural areas through a mechanism of fadeout and reintroduction (2). Other studies have suggested that influenza in the United States spreads through workforce commuting (3) and that dengue spreads along a major road in Cambodia (5). Studying epidemic synchrony requires data at high spatiotemporal resolution for a large sample of locations. Data limitations have restricted previous studies on disease spread and synchrony to small geographical areas within country boundaries. Given the increased (cross-border) mobility of populations, strong evidence of global warming, and potential for rapid, global spread of highly pathogenic infectious diseases, a better understanding of the mechanisms of long-distance disease spread and spatial synchrony is becoming essential for global health security.
Some infectious diseases in endemic settings, such as dengue or influenza, occur annually in well-defined cycles that depend on climate factors, such as precipitation. In addition to an annual cycle, significant variability at multiannual periodicities has been observed for dengue (6, 7). This multiannual periodicity is thought to be driven by cycling of immunity in the host population and has been observed to vary over time (6). To better understand the spatiotemporal dynamics of large epidemics, previous studies used various time series decomposition methods to isolate multiannual oscillations from background annual cycles and higher frequency noise. For example, previous work on disease spread focused on the 1.5–3 y cycle for measles (2), the 3.5–4.5 and 5–6 y cycles for pertussis (8), and the 2–3 y (9, 10) and 3–4 y cycles (6) for dengue.
Whereas mechanisms that cause spatial patterns of multiannual cycles of diseases, such as measles and pertussis, are known (2, 4, 11), these mechanisms remain unclear for dengue. Previous studies have suggested that immunity-driven extinction–reintroduction dynamics of DENV serotypes can play a role, particularly around urban centers (6, 12). The role of multiyear climate variation has also been studied but without consistent results (13).
We studied the synchrony of multiannual dengue cycles across a large geographical area of eight countries in Southeast Asia that span 3,500 km east to west by 2,500 km north to south, with a combined population of 320 million in 2010. We used monthly dengue surveillance data that represent ∼3.5 million reported cases at the provincial level. High dengue transmission rates across all countries combined with extensive diversity in population density, climate, and geology make this region ideal to investigate the long-distance spread of major dengue epidemics that occurred in this region during the past decades.
Results
Strong Region-Wide Synchrony of Dengue Transmission.
We found strong synchrony across the entire region of multiannual dengue cycles and also, for annual cycles and unfiltered incidence rates (IRs) (Figs. 1A and 2). The average power (amplitude) of statistically significant multiannual cycles changed over time and was the highest in 1993–2004 (Fig. S1A). In comparison, the average power of annual cycles was more constant over time but reduced in 1997–2001 (Fig. S1B). We found few statistically significant multiannual cycles in the northern provinces and few annual cycles in the north and the south. Temperatures were too low in the north and too constant in the equatorial south to support cyclical dengue transmission. Synchrony of multiannual cycles changed over time, with the strongest median synchrony of 0.59 (interquartile range = 0.48–0.76) in the 1996–2000 time window (Fig. 3). In comparison, the regional median for annual cycles was consistent over time, ranging between 0.50 and 0.64. The median synchrony of the unfiltered IRs fell between that for the annual and multiannual cycles. Using the average wavelet coherency of dengue cycles as an alternative metric for synchrony over time, we found equivalent patterns (Fig. S1 C and D). Using the entire time series, the average synchrony decreased as the distance between province pairs increased up to ∼1,000 km. (Fig. S2). Synchrony of multiannual cycles peaked at 0.55 [95% confidence interval (95% CI) = 0.45–0.67] compared with annual cycles at 0.70 (95% CI = 0.64–0.75) at short distances.
Monthly dengue IRs (per 100,000 people) and longitudinal climate indicators. Monthly dengue IRs for each province ranked by latitude and monthly climate indicators for corresponding latitudes and time periods. Upper shows median values across provinces or latitudes. NA, not available. (A) Monthly dengue IRs per 100,000 people that have been centered and reduced into z scores, log10-transformed, detrended, and imputed. We imputed missing data by random draws from values of the same months but for different years (Fig. S3). (B) Map of the study provinces by latitude. (C) Average monthly temperature in degrees Celsius from gridded data covering the entire region averaged by latitude and centered and reduced into z scores. (D) The same as C but for total monthly precipitation (millimeters).
Wavelet transforms. Reconstructed periodic cycles of monthly dengue IRs for provinces ranked by latitude. Upper shows monthly distributions across provinces. NA, not available. (A) Reconstructed multiannual cycles. (B) Reconstructed annual cycles.
Synchrony of dengue cycles over time. We computed the average synchrony for moving, overlapping 5-y windows to detect changes over time. (A) Distributions of average synchrony per province per time window plotted at the midyear of each window for multiannual and annual cycles and unfiltered IRs. (B) The average synchrony of multiannual dengue cycles per province for four time windows. (C) The same as B but for annual cycles.
Changes in periodicity over time in months as shown by wavelet transforms and wavelet coherency. For each province, we computed the average power of statistically significant wavelet transforms per month in the multiannual or annual periodicity band. We also computed for each province the percentage of other provinces that had statistically significant wavelet coherency with this province. Upper shows monthly distributions. NA, not available. (A) Average power of statistically significant wavelet transforms in the multiannual periodicity band per month for each province ranked by latitude. (B) The same as A but for the annual periodicity band. (C) For each province ranked by latitude, the percentage of other provinces that had statistically significant wavelet coherency with this province for periodicities within the multiannual band. (D) The same as C but for the annual periodicity band.
Cross-correlation function of geographical distance vs. pairwise Pearson correlation. (A) Average cross-correlation functions (solid lines) and regional averages (dashed lines) of annual and multiannual cycles and log10 IRs centered and reduced to z scores. (B) Average cross-correlations and 95% CIs (solid lines) and the regional averages (dashed lines) for log10 IRs centered and reduced to z scores. (C) The same as B but for annual dengue cycles (6–18 mo). (D) The same as B but for multiannual dengue cycles (19–60 mo).
Dengue IRs and transformations. Monthly values for each province ranked by latitude in color coding. NA, not available. The distributions across provinces per month are shown in Upper. (A) Reported dengue IRs per 100,000 people. (B) Log10-transformed IRs. (C) Log10-transformed IRs that were detrended by subtracting fitted values of a linear model. (D) The same as in C but with missing data imputed.
Climate Forcing of Synchrony.
The strong synchronization of multiannual dengue cycles in 1997–2001 coincided with high temperatures across most latitudes (Fig. 1C) but not with an anomaly in precipitation (Fig. 1D). High temperatures in these years were related to the strongest El Niño episode of the past century (14, 15). We measured a strong wavelet coherency (>0.8) between multiannual dengue cycles and the Oceanic Niño Index (ONI) during 1993–2002 and during 2009–2010 for almost all provinces (Fig. 4). Wavelet coherency for the latter period should be interpreted with caution, because it is at the end of the time series and subject to edge effects. Synchrony of multiannual dengue cycles in the 1997–2001 time window was reduced in a corridor running from Laos to eastern Cambodia (Fig. 3B and Fig. S4A). This area is characterized by high altitude and low temperatures and includes the Annamite Mountain Range (Fig. S4). When studying synchrony between each of 12 major cities in the region and all other provinces, we found two clusters of cities: one cluster consisting of Bangkok, Singapore, Zamboanga (the Philippines), Davao, Cebu, and Taipei, that was synchronous with the Annamite corridor, and a second cluster consisting of Phnom Penh, Vientiane, Hanoi, and Kuala Lumpur, that was synchronous with the Annamite corridor (Fig. S4). These clusters of different synchrony suggest at least two separate networks of epidemiologically connected areas in this region, possibly determined by temperature. Indeed, using a linear regression, we found that synchrony of multiannual cycles was stronger at higher temperatures: synchrony increased with 0.029 for each 1 °C increase in temperature (95% CI = 0.026−0.031) (Table S1). For annual cycles, this coefficient was much lower at 0.004 (95% CI = 0.002−0.005). These results suggest that temperature plays a significant role in the spread of major dengue epidemics. Interestingly, we also found a negative association between population density and synchrony for both annual (β = −0.088/log10 population density per km2) and multiannual (β = −0.037/log10 population density per km2) cycles. This association suggests that more densely populated areas may be able to determine their own nonsynchronous dynamics as independent “pacemakers” instead of phase-locking dynamics with other areas.
Wavelet coherency between the ONI and multiannual DENV cycles. The monthly average statistically significant wavelet coherency between ONI and multiannual DENV cycles across (Lower) the multiannual periodicity band for each province ranked by latitude. Upper shows the distributions (medians and interquartile ranges) of province average wavelet coherency per month. NA, not available.
Average synchrony of multiannual dengue cycles between each of 10 major cities and all other provinces across the entire study period. We found two clusters of cities with different synchrony dependent on temperature. (A) Average synchrony of multiannual dengue cycles per province for cities that were not synchronous with the Annamite region. Two of the largest metropolitan areas, Manila and Ho Chi Minh City, had relatively low synchrony with most of the other provinces in the region, including the Annamite region, and are not shown. (B) Average annual temperature per province in degrees Celsius. (C) Average synchrony of multiannual dengue cycles per province for cities synchronous with the Annamite region.
Linear models of annual and multiannual synchrony vs. climate and population covariates
Traveling Waves of Multiannual Cycles.
We detected traveling waves of multiannual dengue cycles in various parts of the region. We used the phase difference θ to determine the difference (in months) in epidemic timing between provinces. A province could have either a positive- or negative θ compared with another province. A positive θ indicated that a province was timed earlier (leading ahead) vs. the other, and a negative θ indicated that a province was timed later (lagging behind). Outgoing traveling waves can emerge from a province with epidemic dynamics timed ahead of others. In contrast, a province with epidemics lagging behind others could experience an incoming traveling wave. For each province separately, we tested for the presence of local incoming or outgoing traveling waves. For provinces that were lagging behind (θ < 0), we defined an incoming traveling wave as a decreasing lag time with decreasing distance. For provinces that were leading ahead (θ > 0), we defined an outgoing traveling wave as an increasing lag time with increasing distance; 28 provinces had statistically significant incoming traveling waves of multiannual dengue cycles, and 33 had outgoing traveling waves (Fig. 5). Provinces with outgoing traveling waves were located in west Thailand and the Bangkok area, central Laos (Savannakhét and Khammouan), and southern Philippines (Bohol). We found fewer incoming traveling waves for annual cycles concentrated in the northern Philippines (n = 8) and central-eastern Thailand (n = 21). The presence of multiannual waves but not annual waves was statistically significantly associated with temperature and precipitation. Provinces with outgoing multiannual waves had an average of 1.5 °C higher temperature (95% CI = 1.0–2.1 °C) and 53.1 mm (95% CI = 39.7–66.6 mm) lower precipitation compared to provinces without waves. Provinces with incoming multiannual waves had an average of 39.5 mm (95% CI = 22.4–56.5 mm) lower precipitation, but no significant temperature difference compared to other provinces.
Traveling waves of synchrony across provinces. For each province, we fitted a linear model of the phase difference θ of multiannual and annual dengue cycles vs. geographical distance (kilometers). A negative θ indicated that a province epidemic cycle was timed later than another province, possibly experiencing an incoming traveling wave (decreasing θ with decreasing distance). A positive θ indicated that a province was timed earlier than another province, possibly experiencing a positive traveling wave (increasing θ with increasing distance). For θ < 0, we inversed the distance for more intuitive displays. (A) Fitted values of linear models of θ for multiannual cycles vs. distance for each province. We fitted models separately for incoming and outgoing waves. Fitted values are only shown for provinces with a statistically significant model coefficient. We used a Bonferroni-corrected significance level (P < 2e−4) for each province but also, showed fitted values for models with significant coefficients at the 0.05 level (gray lines). The fitted values for the regional average model are shown as black lines. (B) The same as A but for annual cycles. (C) Provinces with statistically significant incoming (red) and outgoing (blue) waves of multiannual cycles. (D) The same as C but for annual cycles.
Discussion
This analysis of large-scale surveillance data revealed strong region-wide synchrony in multiannual dengue cycles. We used a “synoptic epidemiology” approach that spans a large geographical scale but includes granular detail, providing an instantaneous picture of region-wide and local disease dynamics. Synchrony of multiannual cycles changed over time, with a maximal region-wide synchronization occurring during 1997–2001, whereas synchrony of the annual cycles was consistent over time. Synchrony of multiannual dengue cycles during this period coincided with the highest temperatures of the study period across most latitudes and the strongest El Niño event of the century (14, 15). We measured strong wavelet coherency between multiannual dengue cycles and the ONI across most provinces during 1993–2002, but this coherency decreased afterward. A previous study found identical nonstationary wavelet coherency during this period between multiannual dengue cycles and El Niño Southern Oscillation indices for one province in Vietnam (Binh Thuan) (9). The transient nature of this association over time suggests a threshold effect, where the spread of major DENV epidemics may be facilitated by abnormally high temperatures or high temperatures for an abnormally long period. Indeed, we found that synchrony of multiannual cycles increased as temperature increased. Throughout 1997 and 1998, high temperatures across the region could have sustained high levels of dengue transmission, leading to a depletion of susceptibles and low transmission in the following years (Fig. 1A). This hypothesis is consistent with the biology of the Aedes aegypti vector, which reproduces faster and transmits DENV more efficiently at higher temperatures (16, 17). Also, in 1998, a new DENV-2 strain (Cosmopolitan genotype) emerged in Asia (18). High temperatures across a sufficiently large geographical area combined with a large pool of susceptibles could enable the spread of major synchronous dengue epidemics when new DENV types emerge. Indeed, we found a 1.5 °C higher temperature in provinces with outgoing multiannual traveling waves compared to provinces without these waves.
Synchrony of multiannual as well as annual cycles was inversely associated with population density. This association could suggest an extinction–reinvasion mechanism, where synchronous cycles emerge in areas of low population density after a period of low transmission. Densely populated areas are less prone to such “fade-outs,” because they supply a constant pool of susceptibles that can sustain ongoing transmission of all four DENV serotypes. These urban centers could act as independent pacemakers of epidemic dengue cycles into the surrounding areas (6, 12). This mechanism is consistent with the two clusters of synchrony among major cities in the region: one synchronous with the Annamite Mountain Range and the other synchronous with the west-central Thailand area. These areas are also the two main areas for which we detected outgoing traveling waves for multiannual cycles, suggesting that forcing mechanisms, such as temperature, act independently in each of these areas.
The role of human movement in the spread of disease epidemics has been a strong focus of recent research, greatly facilitated by the emergence of novel data sources, such as mobile phone and flight data (19, 20). Strong synchrony and traveling waves of dengue epidemics across Southeast Asia could be driven by (cross-border) population movement. We were unable to formally test this because of limited data on human movement. However, a high degree of population connectivity and spread of DENV across country borders in Southeast Asia has been shown by phylogenetic studies, showing that DENV genotypes circulating in countries, such as Thailand and the Philippines, were isolated across the entire region (18, 21⇓–23). We found spatial structuring in time lags of dengue cycles across provinces that was consistent with traveling waves of dengue incidence. Outgoing traveling waves were concentrated in central Thailand, the east Mekong, and the southern Philippines. It was also in these countries that dengue was first recognized in the 1950s (24). Provinces with outgoing traveling waves had higher temperatures and lower precipitation compared to other provinces. Local climate, virus, and population conditions in these areas may have ignited the emergence and spread of new DENV types that resulted in region-wide synchronous dynamics through widespread high temperatures during a strong El Niño episode and population movement.
Powerful forcing mechanisms in Southeast Asia, particularly sustained high temperatures, can drive the synchronized spread of major dengue epidemics on a continental scale. This analysis improves opportunities for future studies on the causal mechanisms and for predictive modeling of large-scale dynamics of dengue as well as other infectious diseases. This study also demonstrates the advantages of multicountry collaboration to advance infectious disease surveillance, analysis, and control.
Materials and Methods
Data.
Monthly dengue surveillance data and corresponding population (25) and climate data (26⇓–28) at the provincial level (29) were available for 273 provinces in Thailand, Cambodia, Laos, Vietnam, Malaysia, Singapore, the Philippines, and Taiwan (more details on data sources and inclusion are in SI Materials and Methods). We computed monthly dengue IRs per 100,000 people for 1993–2010 for provinces in Thailand, Malaysia, and Singapore; 1994–2010 for the Philippines and Vietnam; and 1998–2010 for the other countries (Figs. S3 and S5 and Movie S1). All data are publicly available through Project Tycho (www.tycho.pitt.edu).
Monthly data available by province ranked by country. For each province, the months between 1993 and 2010 were included in the analysis for a total of 273 provinces. Upper shows the total number of provinces included per month. The length of the time series determined the maximum periodic cycle that could be studied by wavelet analysis. To reduce edge effects, we defined the maximum periodicity for each province as the number of observations divided by 2.5 (i.e., the time series should be able to contain at least 2.5 repeats of a periodic cycle).
Wavelet Transforms.
We used wavelet methods to decompose reported dengue IRs into multiannual and annual cycles as described previously (2, 8, 9, 30, 31). Wavelet transforms are appropriate to characterize nonstationary signals with multiple periodicities, such as dengue IRs (6). We used a Morlet wavelet with a nondimensional frequency ω0 = 6 as used previously (8, 30). This wavelet is complex, enabling the extraction of phase angles to study epidemic timing. We explored the influence of the selected value of ω0 on synchrony in a sensitivity analysis and found that only extreme values influenced study results (Fig. S6). We found that most provinces had statistically significant annual cycles with a periodicity of 6–18 mo and statistically significant multiannual cycles with a periodicity of 19–60 mo (Fig. S7). We analyzed reconstructions of these cycles for provinces with statistically significant cycles in the annual or multiannual periodicity range.
Sensitivity analysis of average synchrony per province for different values of
Variance explained by different periodicities for each province. We used wavelet analysis for periodicities ranging from 2 mo to the maximum periodicity for a province, with a 1-mo step size along a linear scale. (A) The global wavelet power spectrum using only statistically significant wavelet transforms per province ranked by latitude. Upper shows the distributions across provinces. NA, not available. (B) The average multiannual periodicity across periodicities ranging from 19 to the maximum supported by a province time series weighted by the average power for each periodicity.
Synchrony.
We used pairwise Pearson correlation coefficients of multiannual and annual dengue cycles and unfiltered IRs between provinces to measure synchrony ρ. Pearson correlation indicates similarity in both timing and amplitude of epidemic cycles. We computed the average ρ for a province as the average across all province pairs that included that province weighted by the number of pairs with non-missing data. We computed this average using the entire time series of dengue cycles but also, for moving, overlapping 5-y time windows. We also computed the average wavelet coherency between province pairs in the annual and multiannual periodicity bands as described previously (7, 32) using parameter values for the wavelet transforms as described. Wavelet coherency (ranging from zero to one) describes the phase relationship between two time series localized in a time–periodicity spectrum. For strong wavelet coherency, statistically significant cycles of a specific periodicity need to be phase-locked (positively or negatively). We also used wavelet coherency to assess the association between multiannual dengue cycles and the ONI. The ONI identifies El Niño and La Niña events in the tropical Pacific based on sea surface temperature in the Niño 3.4 Region. To measure the dependency of synchrony on the average level of population density, temperature, and precipitation of province pairs, we used a multivariate linear regression with synchrony as the dependent variable and these covariates as independent variables.
Phase Angles.
We used phase angle transforms of multiannual and annual cycles to study epidemic timing as described previously (2, 7, 31, 32). We expressed the pairwise phase difference θ between province pairs in months by assuming a 12-mo periodicity for annual cycles and 39-mo periodicity for multiannual cycles. We defined a traveling wave for a province as a statistically significantly linear association between θ and geographical distances for that province vs. all others (Fig. S8). A negative θ indicated that a province epidemic cycle was timed later vs. others and that this province could have an incoming traveling wave. A positive θ indicated that a province was timed earlier and could have an outgoing wave. For each province, we tested the presence of “local” incoming traveling waves (decreasing θ with decreasing distance for θ < 0) and local outgoing traveling waves (increasing θ with increasing distance for θ > 0) using a linear regression model:
Local traveling waves for Bangkok and four hypothetical example provinces. (A) Location of Bangkok and four hypothetical provinces at varying distances. (B) IRs (per 100,000 people) for each location. (C) Wavelet reconstruction of multiannual (19–60 m) dengue cycles. (D) Phase angles of multiannual dengue cycles. (E) Linear model fit of the association between lag time in months and geographical distance from Bangkok. For negative phase angle differences, we inversed the geographical distance for a more intuitive display. Incoming traveling waves into Bangkok were defined as decreasing lag times with decreasing distance (other provinces timed ahead of Bangkok). Outgoing waves emerging from Bangkok were defined as increasing lag times with increasing geographical distance (other provinces timed after Bangkok). We measured the presence of these incoming or outgoing traveling waves for each province in our study (Fig. 5).
where θp,q is the lag time between provinces p and q for annual or multiannual cycles, and dp,q is the distance in kilometers between provinces. We inversed the sign of distance for negative lag times for more intuitive displays of incoming waves. We defined local as distances ≤1,000 km (multiannual cycles) or ≤1,500 km (annual cycles). Using a linear model of θ vs. distance including all province pairs, we found that, after these distances, θ did not continue to statistically significantly change with distance. We defined a statistically significant traveling wave as a positive
The entire analysis was conducted in the R Statistical Package, version 3.2.1.
SI Materials and Methods
1.1. Data Collection.
We collected monthly dengue surveillance data at the provincial level for Thailand, Cambodia, Laos, Vietnam, Malaysia, Singapore, the Philippines, and Taiwan. Data for Thailand were provided by the Thai Ministry of Public Health Bureau of Epidemiology (1968–2010). Data for Cambodia were provided by the Cambodia Ministry of Health (MOH) National Dengue Control Program (1998–2009) and supplemented with data from the WHO DengueNet Database (1998–2010) and the WHO Western Pacific Regional Office (1998–2001) where needed. Data for Laos were provided by the Lao People's Democratic Republic MOH National Center for Laboratory and Epidemiology (1998–2010) and supplemented with data from the WHO DengueNet Database (1998–2007) and the WHO Western Pacific Regional Office (1998–2001) where needed. Data for Vietnam were provided by the National Institute of Hygiene and Epidemiology (1994–2010). Data for Malaysia were provided by the Malaysia MOH Disease Control Division (1988–2010). Data for Singapore were provided by the Singapore National Environment Agency (1993–2010). Data for the Philippines were provided by the Philippines MOH National Epidemiology Center (1993–2010), and data for Taiwan were provided by the Taiwan Centers for Disease Control (1998–2010). Dengue surveillance systems were similar across countries and mostly passive, except for the Philippines and Cambodia (sentinel system). In most countries, dengue is predominantly reported among children. Cambodia was the only country for which dengue surveillance was explicitly restricted to children <15 y old. Despite these slight differences, population-level disease patterns were comparable across all countries. For each country, data were collected at the provincial level. We defined a province as the smallest geographical area for which data were available (the first or second administrative level). Where needed, surveillance reports for dengue fever and dengue hemorrhagic fever were combined into the total number of dengue cases. Population counts were obtained for provinces by year for each country. Population counts for Thailand were provided by the Thai Ministry of Public Health Bureau of Epidemiology (1968–2010); population counts for Cambodia were provided by the National Institute of Statistics (1998) and the MOH National Dengue Control Program (2002–2010). Population counts for Laos were obtained from the Laos Statistics Bureau (1995 and 2000–2010). Population counts for Vietnam were provided by the General Statistics Office of Vietnam (1995–2010). Population counts for Malaysia were provided by the Department of Statistics (2000 and 2003–2010). Population counts for Singapore were obtained from the US Census Bureau International Database (25). Population counts for the Philippines were obtained from the Philippines Census Bureau (2000, 2005, and 2010), and population counts for Taiwan were provided by the Taiwan Centers for Disease Control (1998–2010). Administrative boundary files were obtained from the Global Administrative Areas (GADM) Database (29). Monthly gridded (
1.2. Data Exclusions.
Province names were standardized based on the administrative level one or two divisions in the GADM Database as of 2013 (29). Provinces that could not be located were excluded (e.g., because of a split after the last update of the boundary definition files in the GADM Database). The study period was determined separately for each province. No data before 1993 were included because of lack of data for most countries. For each province, the first and last years of a study period were defined by the first and last years with at least six monthly observations. Provinces were excluded if 12 or more consecutive monthly observations were missing, two-thirds or more of all observations were missing, or if the study period was less than 4 y (2.5 cycles of the smallest multiannual period of 19 mo). These in- and exclusions resulted in a total of 273 of 289 provinces included (94%).
1.3. Data Management.
For multiple countries, the numbers and names of provinces changed over time because of merging or splitting events. We used province boundaries as in the GADM Database as of 2013; we merged data for provinces that merged and only used data after a splitting event for provinces that had split. Geographical distances between all province pairs were computed using latitude/longitude of province centroids according to the World Geodetic System 1984 Revision Coordinate System. We computed the average air temperature and precipitation per latitude for each month. We centered and reduced the time series for each latitude to z scores around the latitude mean. We also extracted the average value of air temperature, precipitation, and population density for province polygons from the gridded data. Population estimates were not available for the entire study period of each province. We applied country-specific, nationwide annual population growth rates obtained from the US Census Bureau International Database (25) to estimate the provincial-level population for missing years. The growth rate
We used dengue case data and population estimates to compute monthly IRs per 100,000 people for each province.
1.4. Data Transformation.
IRs were log-transformed as
1.5. Reconstructing Annual and Multiannual Cycles.
We used wavelet analysis to isolate cycles with multiannual and annual periodicities from the reported IRs. Wavelet analysis is well-suited to characterize epidemiological time series containing cyclical variability with periodicities that change over time (nonstationarity). Wavelet analysis has been used previously to study a wide range of infectious diseases, such as measles, pertussis, and dengue (2, 5, 8, 9, 11, 13). Detailed methods for wavelet analysis have been described previously (30, 31). We computed Morlet wave transforms for each province using a nondimensional frequency
For each scale s and time interval
For each province, we computed the average of wavelet transforms for each periodic intervals over all time points n (global wavelet spectrum) according to the work by Torrence and Campo (30):
Similarly, for each province, we also computed the average of statistically significant wavelet transforms per time interval
We reconstructed annual epidemic cycles defined as the average of 6–18 mo periodic cycles and multiannual cycles defined as the average of 19–60 mo periodic cycles. We reconstructed these cycles using the filter described in the work by Torrence and Campo (30):
1.6. Estimating Synchrony.
We defined synchrony between provinces as the pairwise Pearson correlation coefficient, taking into account both timing and amplitude of the signals. We computed synchrony for both annual and multiannual cycles. To assess the association between synchrony and geographical distance between provinces, we used a previously described nonparametric covariance function that estimates the underlying correlation without assuming any particular shape (2, 33). It uses a smoothing spline (bandwidth of 300 data points) to estimate the curve and bootstrapping (n = 1,000) to estimate confidence limits:
We measured changes in synchrony over time by computing synchrony separately for parts of the time series within moving and overlapping 5-y time windows. For each time window, we computed the average synchrony per province weighted by the number of province pairs included. We did this for annual and multiannual cycles and the unfiltered IRs (Fig. 3). We used a linear model to estimate the association between synchrony and the pairwise average log10 population density, precipitation and temperature for province pairs, adjusted for geographical distance between pairs:
1.7. Wavelet Coherency.
As an alternative metric for synchrony over time and for comparison with the wave correlation analysis, we measured wavelet coherency between province pairs. As previously described, wavelet coherency uses wave transforms of two time series to indicate their localized phase relationship in a time–frequency spectrum (7, 32). Wavelet coherency ranges from zero to one, and high wavelet coherency requires that statistically significant cycles of a particular periodicity are detected in both time series and that these cycles are phase-dependent (positively or negatively):
We studied changes over time in wavelet coherency for each province by computing the percentages of other provinces that were statistically significantly coherent with this province at each time point (Fig. S1 C and D). In addition, we used wavelet coherency to estimate the association between the ONI and the multiannual dengue cycle for each province. For each province, we averaged all statistically significant values of this wavelet coherency matrix within the multiannual periodicity band (19–60 mo).
1.8. Phase Differences.
As previously described, we computed phase angles for annual and multiannual cycles as indicators of epidemic timing (2, 30, 31). For each province and for annual and multiannual periodicities separately,
We computed the average phase angle-
We used these phase angles to compute the phase angle difference θ between provinces for annual and multiannual cycles. The phase angle difference θ was constrained between
1.9. Traveling Waves.
We defined a local traveling wave for a province as a statistically significant positive linear association between θ and geographical distances between that province and the other provinces. A province could have a negative or positive θ with other provinces. A negative θ (lag time) indicated that a province was timed later compared to the other province (lagging behind), and a positive θ (lead time) indicated that the province was timed earlier compared to the other province (leading ahead). Outgoing traveling waves can emerge from provinces that are timed ahead of others. Incoming traveling waves can occur for provinces that are timed behind others. For provinces that were lagging behind (negative θ), we defined an incoming traveling wave as a decreasing θ with decreasing distance. For provinces that were leading ahead (positive θ), we defined an outgoing traveling wave as an increasing θ with increasing distance. For each province and separately for annual and multiannual cycles, we used a linear model to detect traveling waves:
1.10. Sensitivity Analysis.
We conducted a sensitivity analysis on the effect of the value for the wavelet parameters
1.11. Computing Environment.
All analyses were conducted using the R System, version 3.2.1. We used the package dplR for the wavelet analysis with a modified Morlet function with a linear periodicity scale. We used the package biwavelet for wavelet coherency with modified wtc and wt functions with a linear periodicity scale. We used the function Sncf for the nonparametric analysis of the correlation between synchrony and geographical distance.
Acknowledgments
The authors thank surveillance data managers in each country for support in compiling dengue surveillance data. We also thank Dr. Nyphonh (Lao People's Democratic Republic), Ms. Hannah Lewis (WHO), Dr. Torréton (Institut de Recherche pour le Développement, Vietnam), and Dr. Sok Touch (CDC Cambodia) for support in facilitating collaborations and Mr. Guido Camargo for proofreading of this manuscript. W.G.v.P., D.A.T.C., X.X., N.S.C., M.C., and D.S.B. would like to thank the Bill and Melinda Gates Foundation (Grant 49276, Evaluation of Candidate Vaccine Technologies Using Computational Models) and the US National Institute of General Medical Sciences (Grant 5U54GM088491, Computational Models of Infectious Disease Threats) for their support. S.K.L. would like to thank the University of Malaya for their support (UM.C/HIR/MOHE/CHAN/15). M.C. would like to thank the Groupe de Recherche International (GDRI) “Biodiversity and infectious diseases in Southeast Asia.” D.A.T.C. was funded by a Career Award at the Scientific Interface from the Burroughs Wellcome Fund and the US NIH R01AI102939 and R01AI114703.
Footnotes
- ↵1To whom correspondence should be addressed. Email: wav10{at}pitt.edu.
Author contributions: W.G.v.P., M.C., D.S.B., and D.A.T.C. designed research; W.G.v.P., M.C., X.X., N.S.C., P.A., S.I., S.K.L., C.K.C., F.C.L., B.P., P.V., K.B., H.R., N.T.H., P.Q.T., T.N.D., J.-H.C., Y.-L.L., L.-C.N., Y.S., E.A.T., V.G.R., L.L.L.S., R.G.J., R.V.G., J.M.S.V., I.-K.Y., D.S.B., and D.A.T.C. performed research; W.G.v.P., M.C., X.X., N.S.C., and D.A.T.C. analyzed data; and W.G.v.P., M.C., S.I., S.K.L., F.C.L., K.B., H.R., P.Q.T., Y.-L.L., R.G.J., R.V.G., D.S.B., and D.A.T.C. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1501375112/-/DCSupplemental.
Freely available online through the PNAS open access option.
References
- ↵
- ↵
- ↵.
- Viboud C, et al.
- ↵.
- Rohani P,
- Earn DJ,
- Grenfell BT
- ↵.
- Teurlai M, et al.
- ↵
- ↵.
- Cazelles B,
- Chavez M,
- Magny GC,
- Guégan J-F,
- Hales S
- ↵.
- Choisy M,
- Rohani P
- ↵
- ↵
- ↵.
- Magpantay FMG,
- Rohani P
- ↵
- ↵
- ↵
- ↵
- ↵.
- Watts DM,
- Burke DS,
- Harrison BA,
- Whitmire RE,
- Nisalak A
- ↵
- ↵.
- Salda LTD, et al.
- ↵
- ↵.
- Colizza V,
- Barrat A,
- Barthélemy M,
- Vespignani A
- ↵
- ↵
- ↵
- ↵.
- Hammon WM,
- Rudnick A,
- Sather GE
- ↵US Census Bureau (2014) International Database. Available at www.census.gov/population/international/data/idb/informationGateway.php. Accessed July 14, 2014.
- ↵.
- Fan Y,
- van den Dool H
- ↵
- ↵Center for International Earth Science Information Network (CIESIN), Columbia University; Centro Internacional de Agricultura Tropical (CIAT) (2005) Gridded Population of the World Version 3 (GPWv3) [NASA Socioeconomic Data and Applications Center (SEDAC), Palisades, NY].
- ↵GADM (2013) Database of Global Administrative Areas. Available at gadm.org/home. Accessed July 21, 2015.
- ↵
- ↵
- ↵
- ↵
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Biological Sciences
- Medical Sciences