Climatic-driven seasonality of emerging dengue fever in Hanoi, Vietnam

Background Dengue fever (DF) has been emerging in Hanoi over the last decade. Both DF epidemiology and climate in Hanoi are strongly seasonal. This study aims at characterizing the seasonality of DF in Hanoi and its links to climatic variables as DF incidence increases from year to year. Methods Clinical suspected cases of DF from the 14 central districts of Hanoi were obtained from the Ministry of Health over a 8-year period (2002–2009). Wavelet decompositions were used to characterize the main periodic cycles of DF and climatic variables as well as the mean phase angles of these cycles. Cross-wavelet spectra between DF and each climatic variables were also computed. DF reproductive ratio was calculated from Soper’s formula and smoothed to highlight both its long-term trend and seasonality. Results Temperature, rainfall, and vapor pressure show strong seasonality. DF and relative humidity show both strong seasonality and a sub-annual periodicity. DF reproductive ratio is increasing through time and displays two clear peaks per year, reflecting the sub-annual periodicity of DF incidence. Temperature, rainfall and vapor pressure lead DF incidence by a lag of 8–10 weeks, constant through time. Relative humidity leads DF by a constant lag of 18 weeks for the annual cycle and a lag decreasing from 14 to 5 weeks for the sub-annual cycle. Conclusion Results are interpreted in terms of mosquito population dynamics and immunological interactions between the different dengue serotypes in the human compartment. Given its important population size, its strong seasonality and its dengue emergence, Hanoi offers an ideal natural experiment to test hypotheses on dengue serotypes interactions, knowledge of prime importance for vaccine development.


Background
Dengue Fever (DF) has recently been recognized by WHO as the fastest spreading tropical disease across all continents. Bhatt et al. [1] estimated that the global number of new infections per year (390 millions, 95% confidence interval: 284-528) is largely underestimated: only 96 millions (95% confidence interval: 67-136) being declared yearly. DF is one of the many symptoms (ranging from mild fever to hemorrhagic fever and shock syndrome) caused by one of the four serotypes of dengue virus (Flaviridae family). Even though recovery from dengue confers life-long immunity against the infecting serotype, immunological interactions between the different serotypes are not fully understood. The virus is transmitted by bites of female Aedes aegypti or albopictus mosquitoes in the inter-tropical regions of the world. In absence of vaccine (under development), mosquito control is the only available method of prophylaxy.
In Vietnam, dengue is recognized as a major cause of mortality and morbidity and ranks amongst the top ten communicable diseases in terms of overall health burden [2]. All four dengue virus serotypes have been found circulating in Vietnam with the dominant one varying over time. Reports from the National Institute of Hygiene and Epidemiology show that DENV-1 and DENV-2 have been the predominant circulating viruses almost every year. DENV-3 emerged in the late 1990s and was responsible for the large outbreak of 1998, whereas DENV-4 was also detected between 1999 and 2003 [3]. Dengue transmission occurs throughout the year in Vietnam, with peaks in the numbers of cases (72% of total cases) reported between June and November [4]. There are regional variations in the seasonality of dengue epidemiology in Vietnam. In the Northern and Central Highland regions, dengue notifications are low during the winter time from December to March, while in southern regions, dengue transmission occurs throughout the year, even if it sharply increases during the rainy season from July to September. Given that dengue is a vector-born disease and that vector population dynamics are strongly dependent on climatic factors, the diversity of climates in Vietnam may explain the observed diversity of dengue epidemiological dynamics.
Hanoi, the capital of Vietnam, is located in the North of the country and known as a low transmission setting of DF [5]. Hanoi experiences annual seasonal dengue outbreaks with the pinnacle of epidemics usually falling in September/October and ending in November/December. Over the last decade, the number of DF cases has been increasing from year to year, reaching a peak in 2009. According to the Ministry of Health's statistics, the outbreak in Hanoi in 2009 is the most important outbreak of the last decade, with 384 notified cases per 100,000 individuals. Interestingly, 2009 was also the year El Niño increased actively [6,7]. There are only few studies published on dengue epidemiology in Hanoi that are based on the public health surveillance data routinely collected through the Ministry of Health's notifiable diseases surveillance program. Toan et al. [8] show that there are spatio-temporal clusters of DF limited to a radius of 1,000 m and a duration of 29 days. This study also demonstrates that most of the DF cases occur between June and November, during which the rainfall and temperatures are highest. Cuong et al. [5] use wavelet analysis to relate dengue incidence to climatic variables and suggest that all the tested local climatic variables (total rainfall; mean wind velocity; mean, maximum and minimum temperatures; relative humidity) are significantly associated with dengue incidence around the annual periodicity: on average, dengue incidence follows the seasonal peak of rainfall and mean temperature with a lag of 1 to 2 months.
Other studies have been carried out on the correlation between climate and DF in other parts of Vietnam as well as in other parts of the world, using a wide spectrum of mathematical and statistical modeling methods [6,[9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25]. In Vietnam, most of the studies have been carried out in the south and the center of the country, and showed significant associations between climatic variables and dengue incidence. A wavelet analysis of monthly dengue cases from the province of Binh Thuan has shown a non-stationary relationship between El Niño Southern Oscillation indice and dengue incidence in the 2-3 year periodic band [6]. Meanwhile, a correlation study carried out on monthly dengue cases from the province of Daklak has found the risk of dengue to be associated with high temperature, high relative humidity and rainfall, but inversely associated with duration of sunshine [22]. Findings from other studies in many other parts of the world also show climatic variables to have an effect on dengue transmission. Studies in Thailand [11], Barbados [12], Taiwan [16], Guangzhou, China [18], the French West Indies [21] and in Colombia [26] showed a positive correlation between dengue incidence and minimum and maximum temperatures, precipitation and relative and absolute humidities. However, depending on the approach of analysis and the areas, these correlations were more or less strong. In Barbados, the strongest correlation was found at a lag of 6, 12 and 16 weeks for vapour pressure and minimum and maximum temperatures respectively, whereas in Taiwan the highest correlations were found with maximum temperature at a lag of 5 weeks and with total precipitation at a lag of 7 weeks. In Colombia, Eastin et al. [26]'s results suggest that DF cases increase 2 to 5 weeks after the daily temperature range remains for an extended period within the temperature range optimal for vector survival and disease transmission. Nagao et al. [11] in Thailand and Yi et al. [27] in Guangdong, China, demonstrated that the distributions of Aedes species and dengue cases were positively associated with high absolute humidity, which itself increases with high temperature and rainfall. In San Juan, Puerto Rico, Schreiber [9] used a water budgeting technique and showed that high levels of dengue are associated with reduced actual evapotranspiration, minimum temperature and with high levels of precipitation. In Taiwan, using autoregressive integrated moving average models, Wu et al. [20], found a negative association of dengue incidence with temperature and relative humidity. Finally, in the city of Noumea (New Caledonia), Descloux et al. [24] recently documented a high seasonality of dengue incidence, with an epidemic peak (March-April) lagging the warmest temperature by 1 to 2 months and in phase with maximum precipitations, relative humidity and entomological indices.
In the present study, we consider vapor pressure and relative humidities, temperature and rainfall in order to identify which of these variables are most critical for the onset of dengue epidemics. Compared to Cuong et al. [5], who also investigated the links between dengue and climatic variables in Hanoi, Vietnam, we here consider vapor pressure in addition to relative humidity. Vapor pressure is a measure of absolute humidity and this climatic variable is often neglected in the studies investigating the links between climate and disease transmission, even though it has been proved to play a role more important than relative humidity for the transmission of some diseases such as influenza (e.g. [28]). A second difference with Cuong et al. [5] is that we here work on weekly incidence aggregates instead of monthly aggregates. With this finer temporal resolution we aim at investigating intra-annual patterns of seasonality. We first quantify the synchrony between weekly dengue incidence and the four climatic variables between 2002 and 2009. We then estimate the values of reproductive ratio of dengue fever through time and characterize the trend and seasonality of reproductive ratio. We finally discuss the results in light of immunological and entomological factors specific to dengue epidemiology in Hanoi, Vietnam.

Methods
Setting: Hanoi (1,760 km 2 ) is located in the Red River Delta, in the centre of North Vietnam (21 o 2 ' N, 105°51 ' E). The city experiences the typical climate of northern Vietnam, where summers are hot and humid, and winters are, by national standards, relatively cold and dry.

Data collection
The data used in this study have been previously published by Cuong et al. [5] and Toan et al. [8]. Clinical suspected cases of DF in old Hanoi (14 districts of Hanoi before merging with Ha Tay in 2008) are reported to the surveillance system of Hanoi Center for Preventive Health. The criteria for notification of DF disease are based on the guidelines of the Ministry of Health (2006) on surveillance, diagnosis and treatment of dengue, in which suspected dengue cases are based on acute febrile illness (≥38°C) of 2-7 days duration with at least two of the following non-specific manifestations of dengue fever: headache, retro-orbital pain, myalgia, arthralgia, rash, hemorrhagic manifestations, and leucopenia [29]. The data analyzed here include all reported cases from January 2002 to December 2009, aggregated by week.
Daily weather data for Hanoi from 2002 to 2009 were provided by the National Centre for Hydrometeorological Forecasting. They include the records of mean, maximum, and minimum temperatures (T, in°C), rainfall (in mm) and relative humidity (in %).

Vapor pressure
We used vapor pressure (VP, in mb) as a proxy of absolute humidity. VP was calculated from relative humidity (RH, in %) and temperature T using the Clausius-Clapeyron formula [28,30]: where L = 2,257 J/g is the latent heat of evaporation for water, R v = 416.5 J/(kg K) is the gas constant for water vapor, T 0 = 273.15 K and VP 0 = 6.11 mb is the vapor pressure at which water would change phase between vapor and liquid if the temperature was T 0 . Note that in the above formula temperature T is expressed in°C instead of K as in Shaman and Kohn [28].
Climatic variables were aggregated by week using sums for rainfall and mean values for all the other variables.

Wavelet analysis
Epidemiological data can be substantially non-stationary [31][32][33] as is the case for dengue in Hanoi where it is emerging (see in particular the increase in mean and amplitude over time on Figure 1B). Here we performed wavelet decomposition, a time-series statistical analysis allowing to efficiently deal with non-stationary data. Specifically, we used the Morlet wavelet [34], classically used in ecology with a non-dimensional frequency ω 0 = 6. An advantage of using the Morlet wavelet is that it is a complex wavelet, allowing to quantify the phase and thus calculate time lags between different time series.
Coherence based on wavelets allows to perform similar analysis as cross-correlation but for potentially nonstationary signals. Wavelet coherences were calculated to examine the association between two time series, both in time and frequency. Coherence spectra allow to investigate whether different periodic modes of two time series tend to oscillate simultaneously and, if yes, to identify the periodicity around which this association takes place.
Significance levels were calculated by a Chi-square test assuming that the wavelet coefficients are normally distributed as described in Torrence and Compo [35]. The detailed theory for wavelet analysis has been described elsewhere [4]. Before wavelet decomposition, time series were square-root transformed in order to mitigate the weights of high values. They were also zero-padded to the next power of 2 of their length (i.e. 512), in order to minimize edge effects [35].

Reproductive ratio
The reproductive ratio R is the expected number of infections caused by one infected individual. It is maximal at the start of an epidemic when the host population is fully susceptible and then decreases over the course of the epidemics. It reaches the equilibrium value of 1 at epidemic peak and decreases below 1 after it. The initial and maximum value of the reproductive ratio is called the basic reproductive ratio R 0 and is the classical epidemiological statistics as its value relative to 1 informs about the potential for an epidemic to occur. We followed Soper [36] as reported in Keeling and Rohani [37] and approximated the reproductive ratio R by: where C t+1 and C t are the numbers of cases at times t + 1 and t respectively and α is a parameter that expresses the infection generation length (i.e. the sum of the infectious and latent periods) in the same units as the data time steps (here 1 week). The infection generation length for dengue being approximately 2 weeks [38], we set this parameter α to the value of 2. The other part of Soper [36]'s formula expresses the reproductive ratio as a function of the number of susceptibles in the population: where X* and X t+1 are the numbers of susceptibles before the epidemics and at time t + 1, respectively and k is a parameter reflecting some potential external forcings (such as climatic ones). This latter part of the Soper [36]'s equation shows that variations in the reproductive ratio are due either to variations in the number of susceptibles in the population, or to some external forcings directly affecting the transmissibility of the disease (in our case climatic factors acting on the vector population dynamics and density). Given that dengue is emerging in Hanoi and that the basic reproductive ratio of dengue is generally low [38], we expect the depletion of susceptible in the population to be very slow. We expect it to be even slower given the fact that dengue can actually be caused by four different serotypes with no permanent cross-immunity between them. Hence, among the two above-cited factors that can affect the seasonality of the reproductive ratio (susceptible depletion and external forcing), we assume that susceptible depletion is negligible before any external forcing such as climatic drivers on the mosquito population dynamics.
Time series of the calculated seasonal factor k were smoothed by lowess regression with smoothing factors equal to 0.05 and 0.90 in order to reveal its seasonality and its long-term trend respectively. Confidence intervals were calculated by assuming normal distribution of errors.
Given the uncertainty on the infection generation length and the assumption made on the number of susceptibles, estimates of the reproductive ratio will be treated with caution: only their trend and seasonality will be interpreted, not their absolute values, on which we will have limited confidence.
All analyses were conducted in R (R Core Team [39]) and wavelet analyses were performed using the "biwavelet" R package [40].

Dengue incidence and its reproductive ratio in Hanoi
From 2002 to 2009, 23,195 DF cases were reported in Hanoi with the average annual incidence rate of 69.22/ 100,000. Overall, the incidence of DF increased over the 8 years of the study, with a sharp increase during the period 2005-2009. This period of increasing incidence is visible on the wavelet spectrum of Figure 1C. The highest peak of 3,697 cases was recorded in September 2009. Over the 5 years from 2005 to 2009, annual (1 year) and subannual (6 months) periodicities were significant. These two periodicities correspond to a slow increase of DF incidence from the beginning of the year to weeks 22-24 (June), followed by a rapid increase of incidence until weeks 44-46 (November) which ends by a sharp decrease in incidence at the end of the rainy season. Figure 1A shows a long-term increase in the reproductive ratio (orange curve) as well as a non-stationary sub-annual periodicity with two peaks of this reproductive ratio per year (blue curve). These peaks are of roughly equal magnitude from 2005 to 2008, the second peak is substantially higher than the first one in 2004 and 2009, and the number of cases are too low before 2004 for any clear pattern to be visible (large confidence intervals on Figure 1B).

Meteorological variables in Hanoi
The wavelet power spectra of temperature, precipitation, and humidities (relative humidity and vapor pressure) in Hanoi during the study period are shown in Figure 2. Temperature, precipitation and vapor pressure in Hanoi show significant annual periodicities that are constant through time, whereas relative humidity shows both annual and sub-annual periodicities, as observed on the DF incidence time series ( Figure 1B).
The cross-correlation coefficients among four climatic variables in Hanoi are presented in Figure 3. It shows that mean temperature, rainfall and vapor pressure are much correlated and in phase, whereas there is a significant lag between these three variables and relative humidity.
Over the 8 years of the study, the average temperature in Hanoi was found to be lowest in January and February (16.07 ± 1.65°C standard deviation) and highest in June and July (29.80 ± 0.47°C). Over the studied period, the coolest temperature was observed in week 5 (February) of 2008 and the warmest was observed in week 27 (July) of 2009 (10.11°C and 32.84°C respectively). July to September was the time of the year receiving the majority of the annual rainfall (271.9 ± 111.9 mm in average over one week).
Week 44 (November) in 2008 recorded a rainfall extreme of 576.6 mm brought by the typhoon Maysak (Center for Excellence in Disaster Management and Humanitarian Assistance, [41]). The relative humidity in Hanoi is quite high (78.62 ± 4.28%), and is usually higher in February and March (cool but very rainy: 84.12 ± 2.56), and August and September (drier but very hot: 79.81 ± 3.54), than the rest of the year (76.80 ± 3.99).

Coherences between meteorological variables and DF incidence in Hanoi
Results of wavelet coherences between DF incidence and climate variables are shown in Figure 4. Significant coherences were observed between DF incidence and temperature, precipitation and humidity for the annual periodicity from approximately 2005 to 2009 ( Figure 4A-D). Moreover a weaker, but still significant association between relative humidity and DF incidence was also seen for the sub-annual periodicity from 2006 to 2009 ( Figure 4D).
In analyzing the phase difference at the annual cycle between DF incidence and the climatic variables, we found that dengue incidence was consistently trailing temperature, rainfall, vapor pressure and relative humidities with a delay of 9.37 ± 0.02 (standard error), 8.71 ± 0.02, 10.29 ± 0.03, and 18.05 ± 0.24 weeks respectively (see Figure 4). Interestingly, when looking at the statistical association between dengue incidence and relative humidity for the sub-annual cycle, it appears that the time delay of dengue incidence compared to relative humidity decreases from 14.30 to 5.27 weeks over the 8 years of the study.

Discussion and conclusion
Using monthly aggregated data, Cuong et al. [5] showed a clear annual cycle for dengue transmission in Hanoi from 1998 to 2009. In the present study, using weekly data allowed us to further characterize a sub-annual periodicity, in addition to the annual one. This sub-annual periodicity is reflected in the DF incidence time series, by a slow increase of incidence from the beginning of the year to the weeks 22-24 (June), followed by a faster increase of incidence until weeks 44-46 (November), which ends by a sharp decrease of incidence at the end of the rainy season. A potential drawback of working on weekly instead of monthly data is that it decreases the incidence values and thus increases the noise. However, given that we can still detect clear periodicities in our wavelet spectra, this does not seem to affect our analysis too much.
When characterizing the reproductive ratio throughout the studied period, it displays, most of the years, two peaks per year, which is in accordance with the sub-annual periodicity of DF incidence. Among the climatic variables that we investigated (temperature, rainfall, relative humidities and vapor pressure), all of them expectedly displayed strong annual periodicities with temperature, rainfall and vapor pressure leading DF incidence by a constant delay of 8 to 10 weeks. In addition to this strong annual periodicity, relative humidity displays a sub-annual periodicity, as observed on DF incidence. The annual periodicity of relative humidity leads the annual periodicity of DF incidence by a constant delay of 18 weeks whereas the sub-annual periodicity of relative humidity leads the sub-annual periodicity of DF by a delay that decreases from 14.30 weeks in 2002 to 5.27 weeks in 2009 at an almost constant rate of 1.13 week per year. These results are in general agreement with the findings of other studies that climatic factors play a role in the transmission cycles of DF. Interestingly, these two incidence peaks per year that we observed in Hanoi with periods of low incidence occurring in January and February (the coldest months in Hanoi) and in June and July (the warmest months in Hanoi) are in accordance with Eastin et al. [26]'s observation in Columbia where they noted a significant decreases of DF cases soon after extended periods of either very cool or very hot temperatures. Likewise, a study in Taiwan found three turning points of DF that occurred around early August, late August/early September, and late October/early November. The first two turning points were shown to relate with two typhoons around early to mid August in Taiwan that witnessed a sharp drop in temperature and substantial rainfall after it [16]. Similarly, other studies in Thailand and Singapore also revealed significant associations between climatic variables and dengue incidence ( [13,14,42]; Tipayamongkholgul [43,44]). For example, Tipayamongkholgul [43] conducted a study in the Gulf of Thailand and showed that the monthly average local relative humidity in the previous 3-6 months was negatively associated with epidemics of dengue and incidence of dengue cases. Woongkon et al. [44] in Chiang Rai, Thailand, showed that all climatic factors including minimum, maximum temperature, minimum and average relative humidity, evaporation, wind speed and rainfall lead increasing DF incidence by 0-2 months. In our study, DF incidence is characterized by 2 quasicycles with a periodicity of 6 months with the first one showing a slow and constant increase and the second one showing a marked epidemic. The sharp rupture between these phases can be explained by the fact that the reproductive ratio is not constant throughout the year but actually exhibits two peaks per year, with the second peak at least as high as the first one. Given that dengue is vector-born, the factor limiting its transmission is either due to the mosquito population (mostly its population size), or the human population (mostly its proportion of susceptibles). Winter climatic conditions in Hanoi are not favorable to adult mosquitoes and most of the mosquito population survive the winter either as larvae or eggs [45]. In the spring, when weather conditions become favorable again, eggs hatch and adults emerge, probably causing the first peak on the reproductive ratio and the consequent DF incidence increase. The second peak on the reproductive ratio could be due to the second mosquito generation of the year (issued from the first one), hence its potential to be higher than the first one and even partially conceal it. This second peak of higher magnitude would be the cause of the epidemic peak observed on DF incidence during the second half of the year. This epidemic peak would thus be due more to an increase of the number of infected people than to an increase in the mosquito population size and the dengue reproductive ratio. Indeed, dengue epidemic peak appears even when the second peak on the reproductive ratio is not higher than the first one. The right column shows the phase angles of the climatic variables (blue, left y-axis) and DF (red, left y-axis), as well as their difference (black, right y-axis). These phase angles are calculated on signals that have been filtered around the period of maximal power in the spectra of the left column, i.e. annual periodicity for all the climatic variables, as well as also the semi-annual periodicity for the relative humidity. In spectra of the left column, the black contour lines show the regions of power significant at the alpha-risk of 0.05, the paled region of the spectrum delineates the cone of influence due to the zero-padding of the time series, and the power increases from dark blue to dark red. hypothesis to explain the mechanism of dengue epidemiology in Hanoi can be tested by collecting entomological data (larvae and adult densities estimates) in Hanoi all year-round and translating it into mathematical equations. This would allow to check whether a model based on this hypothesis can generate epidemiological patterns that are in accordance with the ones observed on DF incidence data.
The 2 above-mentioned 6-month cycles observed on time series of DF incidence translate into the sub-annual periodicity that we have characterized in addition to the annual periodicity, with the annual periodicity of DF incidence mostly accounting for the high epidemic peak of the second half of the year, and the sub-annual periodicity mostly accounting for the slow and constant increase of DF incidence of the first half of the year. An interesting result of our analysis is that relative humidity also shows these two annual and sub-annual periodicities and that the sub-annual periodicity of relative humidity leads the sub-annual periodicity of DF incidence by a lag that decreases from 14.30 weeks in 2002 to 5.27 weeks in 2009. As interpreted above, this sub-annual periodicity reflects the first peak of the reproductive ratio that we interpreted in the paragraph above as the first mosquito hatching of the year. Explaining the observed shift in the timing of this first peak by a shift in mosquito hatching is biologically unrealistic. Alternatively, we propose that this shift is due to (i) the building-up of the human population immunity from year to year and (ii) the interactions between dengue serotypes (antibody-dependent enhancement, ADE) as explained below.
Most of primary dengue infections are asymptomatic [1]. Before the emergence of dengue in Hanoi (in 2002), most of the human population may have been susceptible to the 4 dengue serotypes and hence most of the dengue cases may have been primary infections, most likely asymptomatic and thus unnoticed by the surveillance system. As the disease progressively emerges in Hanoi, population immunity to different dengue serotypes increases, thus increasing the number of secondary infections relative to primary ones, and thus increasing the number of symptomatic detected cases. Expected consequences of this mechanism is not only an increase in the number of detected cases from year to year (as visible through the upward trend of DF incidence), but also an earlier detection of the epidemics. The latter would explain this observed shift in the sub-annual periodicity of the DF incidence. This mechanism can potentially be reinforced by some ADE-related mechanisms. Indeed, potential epidemiological consequences of the ADE hypothesis that have been proposed in the literature are that it increases the susceptibility to secondary infections and/or the transmissibility from individuals suffering from secondary infections (see for example [46]). Thus, such a mechanism could also explain the number of detected cases from year to year, the earlier detection of the epidemic, and thus the shift in the sub-annual periodicity of the DF incidence mentioned above. Such a hypothesis could be tested by collecting immunological data from the human population of Hanoi (by an aged-stratified seroprevalence survey for example) and investigating whether a mathematical model built on this hypothesis does generate the trend in DF incidence mean and timing that we observe on the data.
Both DF incidence and relative humidity exhibit conspicuous annual and sub-annual periodicities and these periodicities happen to be strongly correlated. However, we warn against over-interpretation of such correlations in term of biological causation. One reason for such a caution in particular is that relative humidity is a variable that depends on both absolute humidity and temperature (the former being naturally strongly influenced by rainfalls). In case where absolute humidity (or rainfalls) and temperature are not perfectly correlated (which is most likely the case), we do expect that relative humidity exhibits annual and biannual periodicities, as the resultant of two periodic signals that are not perfectly in phase. Thus, instead of looking for a mechanistic link between DF incidence and relative humidity, it may be more relevant to look for two links: (i) one between DF incidence and absolute humidity and (ii) one between DF incidence and temperature, possibly accounting for a possible interaction between the two climatic variables. This particular point will be the topic of a subsequent study.
In conclusion, our analysis on the links between climatic variables and DF incidence in Hanoi raises a number of questions of general interest on the relationships between climate and infectious diseases epidemiology. Because of its highly seasonal climate (and thus potentially highly seasonal dengue transmission too), its important population size and density, and its dengue epidemiological transition (current emergence), Hanoi appears as the ideal set-up to test hypotheses about interaction between serotypes. This is an issue both under-understood and potentially of high relevance for vaccine development. Further investigations on dengue in Hanoi call for additional entomological and immunological data, as well as for theoretical developments.