Immunological serotype interactions and their effect on the epidemiological pattern of dengue

1Dengue is caused by infection with any one of the four closely related virus serotypes DEN-1, DEN-2, DEN-3 and DEN-4, and transmitted to humans primarily by the Aedes aegypti mosquito vector. Its emergence and geographical spread are of major global concern and it has seen a 30-fold increase in incidence during the last 50 years. The virus currently affects over 50 million people annually in the tropical and subtropical regions of the world and causes tens of thousands of deaths. Clinical manifestations of dengue can include dengue fever (DF) and the more severe dengue shock syndrome and dengue haemorrhagic fever (DHF), a life-threatening illness with case fatality rates up to 20 per cent in the absence of treatment (which can be reduced to less than 1 per cent with proper medical management).

2Despite recent advances, there are currently no licensed vaccines or effective antiviral drugs available against the disease. One of the key problems in developing a vaccine is that pre-existing antibodies to a particular serotype can exacerbate an infection of a heterologous serotype. This effect is hypothesized to be due to the action of antibody-dependent enhancement, or ADE, whereby cross-reactive antibodies elicited by a previously encountered serotype bind to the newly infecting heterologous serotype, but fail to neutralize it. The virus–antibody complex engages Fc receptors on the surface of cells such as macrophages and dendritic cells, enabling the virus to enter the cell more easily for replication. Animal studies indicate that the ADE process results in higher viral burdens in vivo, which is likely to explain the epidemiological association of severe disease with secondary infection by a heterologous serotype and also the association of severe disease in infants with passively acquired antibody.

3The temporal epidemiological pattern of dengue is characterized by multi-annual fluctuations in DF or DHF incidence and also cyclical oscillations of the individual serotypes whereby the dominant type is sequentially replaced over time; it is unclear, however, whether and how these two features relate to one another. Various mathematical models have been proposed to illuminate these intriguing dynamics and to uncover the driving forces behind them. The consequences of the effect of ADE on the transmission dynamics of dengue were studied by, and it was shown that ADE could promote the coexistence of strains and lead to nonlinear dynamics. further demonstrated that ADE-driven dynamics can provide a good qualitative approximation to the chaotic desynchronization of strains as commonly observed within longitudinal datasets of numbers of clinical cases, although the selective advantage of ADE may be capped by threats to persistence arising from increasingly high amplitude oscillations. However, it has since been argued that either the period of oscillations or the desynchronization between serotypes generated by ADE alone is not in good agreement with those observed among reported time series of DHF cases, and that some form of cross-protective immunity or seasonal forcing (e.g. ) is necessary to explain the observed epidemic periods and disease patterns.

4Here, we present the results of an epidemiological model for dengue transmission, which employs two parameters to represent the effects of ADE. In addition to a parameter describing the enhancement of transmission during secondary infection through increased viraemia, we include the possibility of enhanced susceptibility to heterologous infections among individuals who have experienced primary infections due to the existence of non-neutralizing cross-reactive antibodies. We show that by including both facets explicitly, both the observed temporal patterns of DF outbreaks and serotype dynamics can robustly be reproduced without the need for extrinsic factors such as seasonal variation or stochasticity. The model further allows us to investigate the relationship between the epidemic pattern of dengue and the dynamics of the co-circulating serotypes. Analyses of time series produced by our model suggest that overall incidence rates and relative serotype prevalence, and especially their respective rates of change, are intrinsically linked.

5We use a simple mathematical model to investigate the transmission dynamics of four co-circulating dengue serotypes within the population, similar to the one analysed by. We assume that recovery from an infection with any one serotype provides permanent immunity against this particular serotype, but can lead to an enhancement of other serotypes upon secondary infection after which individuals acquire complete immunity against all four serotypes. (The last assumption is based on the rarity of third or fourth dengue infections; e.g..) With these assumptions, we divide the population into the following classes: s denotes the fraction of the population that has not yet been infected with any of the serotypes and is thus totally susceptible; yi is the proportion infectious with a primary infection with serotype i, ri is the proportion recovered from primary infection with serotype i; yij is the proportion infectious with serotype j, having already recovered from infection with serotype i; and, finally, r is the proportion of completely immune (those who have recovered after being exposed to two serotypes). The model equations are given as follows:

6The host population size is assumed to be constant with 1/μ(=70 years) being the average host life expectancy and 1/σ(=3.65 days) is the duration of infectiousness. As dengue infections lasts less than a week, the incidence of multiple simultaneous infections is assumed to be negligible. Within this framework, the force of infection of serotype i, λi, is given bywhere βi(=400 per year) is the transmission coefficient of serotype i. The effects of ADE are represented by two parameters that describe (i) the enhancement of susceptibility to secondary infection, γij≥1, and (ii) the increase in transmissibility during secondary infection, ϕij≥1. The model behaviour can then be analysed by varying the respective degrees of enhancement while keeping all other parameters constant.

7To investigate a possible effect of the serotype dynamics on disease prevalence within the population, we used a time-dependent Shannon index of biodiversitywhere pi(t) is the relative frequency of serotype i at time t, and calculated the cross-correlation between the two time series H(t) and total disease prevalence P(t) as follows:where μP and μH are the means, and σP and σH are the standard deviations, of time-series P and H, respectively; n is the number of data points.

8Data on yearly patterns of disease incidence, monthly patterns of dengue hospitalization and relative serotype abundance in southern Viet Nam were kindly provided by the dengue surveillance programme at the Institute Pasteur, Ho Chi Minh City.

9As in most countries in South East Asia, dengue is highly endemic in Viet Nam, particularly in the tropical southern regions of the country. Figure 1 shows the epidemiological data of total disease incidence and relative serotype prevalence for South Viet Nam from 1994 to 2007 (serological data were not available for the years 1994 and 1995), together with a more detailed temporal record of children hospitalized with dengue at monthly time intervals for the period 2004–2007. Figure 1a clearly demonstrates the two key epidemiological features of dengue: first, the total incidence rates of DF/DHF show marked variations over the years, with no indication of a trend towards either reduced or increased disease burden within this region. Second, the serotype dynamics display the characteristic fluctuations of single-type dominance over time. Interestingly, despite the period covered by the data in figure 1a being too short for a meaningful analysis of the inter-epidemic period, it is obvious that it is likely to be longer than the 3-year cycles often reported to underlie dengue incidence rates in Thailand.

10Figure 1b offers a more detailed picture of dengue epidemiology within this region over a 4-year period, revealing the strong seasonal component in dengue transmission, which peaks between July and October before dropping to very low levels over the drier winter months. Notably, though, there is no pronounced component of seasonal forcing in the longitudinal behaviour of the individual serotypes. We therefore do not take into account seasonal or stochastic variation in dengue transmission, but investigate whether immune interaction through ADE and cross-immunity can cause the cyclical behaviour in disease incidence and serotype prevalence.

11To compare the relative impact of ADE on the transmission dynamics of dengue and its four serotypes, we used a simple compartmental model and performed a series of simulations varying the two key model parameters: the enhancement of susceptibility to secondary infection by a heterologous serotype, γ, and the increased transmissibility during secondary infection due to elevated viral load, ϕ (see §2 for full model description). Figure 2 shows the various time series tracing the level of prevalence of the four serotypes within the population over time in response to changes in γ and ϕ. The system readily exhibits an abundance of different dynamical behaviours. For small degrees of enhancement (ϕ,γ<1.2), the system approaches a steady-state equilibrium (figure 2a) that is replaced by synchronized oscillations as the level of enhancement increases. Further increase in γ or ϕ desynchronizes the serotypes, resulting in out-of-phase oscillations of increasing amplitude and periods. As the effects of ADE amplify, we observe prolonged episodes of very low disease prevalence with intermittent epidemic outbreaks (figure 2i).

12Of particular interest here is the region of chaotic oscillations that most closely resemble the observed temporal pattern of DF/DHF. Within this region, we find the characteristic multi-annual cycles (π∼3–5 year) in disease prevalence (see the electronic supplementary material for the according periodograms) and, importantly, the sequential replacement of prevailing serotypes. This behaviour is exemplified in figure 3, which contrasts a typical output of our model with a previously published time record of dengue cases in Thailand. The figure shows a good qualitative agreement between the data and model output in both overall serotype dynamics and disease incidence.

13We then analysed the model dynamics within the parameter space of enhancement with respect to synchrony between individual serotypes (figure 4a), serotype replacement (figure 4b), inter-epidemic period (figure 4c) and serotype persistence (figure 4d). From the time series shown in figure 2, it is clear that the desynchronization of serotypes leads to chaotic fluctuations and the characteristic replacements of the dominant types over time. Studying the synchronization patterns within the (ϕ,γ)-parameter plane by measuring the divergence between various pairs of serotypes over time, we found that desynchronization and hence chaotic oscillations typify a significant portion of the parameter region (figure 4a). Using a measure for single-type dominance similar to the one introduced by, we also found that sequential serotype replacement where individual serotypes are temporarily dominant characterize most of the (ϕ,γ)-parameter space (figure 4b). This confirms that desynchronized fluctuations in serotype prevalence are not an artefact of a particular combination of parameter values, but rather the principal model behaviour under the action of ADE alone. Importantly, these oscillations arise even at relatively small values of ϕ and γ, and for realistic values of the basic reproductive number, R0; that is, we do not require biologically unreasonable high levels of enhancement (either in transmissibility or susceptibility) to produce a realistic epidemiological and serological pattern of dengue.

14The reduction in the degree of enhancement necessary for desynchronized oscillations achieved by the decomposition of the effects of ADE has two important implications for the system's dynamical behaviour. First, it allows for greater flexibility in terms of inter-epidemic periods with a significant proportion of the (ϕ,γ)-parameter space yielding realistic 3–5 year epidemic cycles in disease incidence (figure 4c). Second, the reduction in enhancement also lowers the threat of stochastic extinction due to large amplitude oscillations (as shown in figure 4d) as the proportion of time a particular serotype persists above a given threshold level (we chose a threshold of 1E-8, in line with ).

15We next investigated a possible association between disease prevalence and serotype dynamics by measuring diversity H(t) and total prevalence P(t) over a 50-year period. Figure 5a shows a typical time course of relative serotype frequency generated by our model, with the according time series for disease prevalence, P, shown in figure 5b. From figure 5b, we can identify the presence of a time-lagged correlation between infection prevalence and diversity, in that the fluctuations in prevalence appear to trail the fluctuations in diversity. This time-delayed correlation is even more pronounced between the time courses of the rates of change in P and H (dP/dt and dH/dt, respectively), as shown in figure 5c. We confirmed these findings by computing the cross-correlations P⊗H and dP/dt⊗dH/dt as a function of a lag between the two series (see §2). The graph in figure 5d clearly demonstrates the three principal findings: (i) there is a positive correlation between the level of serotype diversity and infection prevalence that is maximized at time-lag ΔP, (ii) the rate of change in diversity correlates inversely with the rate of change in disease prevalence, and (iii) the rate of change in diversity is positively correlated with prevalence at lag ΔdP. (This behaviour is also robust to changes in ϕ and γ, and can be found across a vast region of the enhancement spectrum; see the electronic supplementary material.) Taken together, these results strongly suggest an intricate interplay between the dynamics of individual serotypes and overall disease pattern of dengue where an increase in serotype diversity causes a subsequent increase in disease prevalence, which itself can cause a reduction in the number of circulating types.

16The limited amount of data available on the dynamics of co-circulating dengue serotypes suggest that they are complex, containing periods of synchronization between certain serotypes and desynchronization among others. Simple models of dengue transmission incorporating ADE have demonstrated that this behaviour is a natural feature of such a system, provided that heterologous immunity is established upon secondary infection. However, it has recently been posited that ADE alone cannot account for the observed patterns, and that some form of cross-immunity after primary infection, in conjunction with other extrinsic factors, must be invoked to generate the appropriate periods and synchronization profiles. Here, we have shown that this limitation can be overcome by decomposing the effects of ADE on the population dynamics of dengue into two independent forms of enhancement: susceptibility and transmission. A wide range of dynamical regimes results from the disjunction between these two parameters, which easily accommodate the asynchrony of different serotypes and inter-epidemic periods of 3–5 years that the data evince. Furthermore, decomposing the effects of ADE significantly lowers the levels of enhancement necessary for obtaining realistic disease patterns. This has the added benefit of reducing the system's propensity for exhibiting the high amplitude oscillations, which in turn lowers the risk of stochastic extinction during the low-prevalence, inter-epidemic periods.

17While our results suggest that ADE alone can give rise to observed patterns of serotype abundance through time, they do not rule out immunological cross-protection after primary infection as a contributing phenomenon. Indeed, immunological cross-protection through shared epitopes among four serotypes can give rise to periods of asynchrony, but not as reliably as through the action of ADE. Furthermore, the inclusion of a period of serotype-transcending immunity following primary infection can also cause desynchronized oscillations, even in the absence of ADE. However, in this case, the resulting periods in serotype prevalence and disease incidence are strongly coupled to the period of cross-immunity, a limitation that can easily be overcome by the inclusion of ADE. Another putative mechanism by which the dynamics of dengue can be explained has recently been suggested by; they conclude that a transient period of clinical cross-protection during which the hosts may gain immunity through heterologous infections without developing disease (together with seasonal variation in transmission) can account for observed epidemiological patterns, particularly the apparent shift towards longer periods in DHF incidence and in its age distribution. There is, however, no strong evidence as yet in favour of this mechanism.

18There are a number of methods of discriminating between these various hypotheses. Certainly, further experimental studies on the effects of ADE are necessary to confirm its role of enhancing susceptibility and transmissibility within individuals suffering from severe disease during secondary infections. Its role in DHF has also yet to be clarified, although it is clear that excess mortality alone cannot account for observed dengue patterns, but may play a role in the separation of serotypes. The extent of cross-protection conferred by primary infection can also be determined by experimental studies, and its role in the population dynamics of dengue may possibly be assessed by comparing patterns of dengue with other strain-structured pathogens where immunological cross-protection has been demonstrated. Similarly, the contribution of environmental factors will be better established as data emerge on patterns in areas of differing intensities and seasonalities of transmission. Spatial models and phylodynamic tools (as used in ) will play a critical role in connecting these processes to the patterns observed.

19One avenue we have not yet explored is the possibility of asymmetries in the levels of enhancement (or cross-protection) between the four serotypes. For example, there is evidence that DENV-2 is more often associated with severe secondary infection, which suggests that it requires heterologous priming to cause severe disease. These asymmetries could also play an important role in terms of the sequence of infection, as shown by. However, the inherent difficulty in determining infection histories from patient sera, due to high immunological cross-reactivity, prevents the construction of a coherent picture of serotype–serotype (immune) interaction. Future longitudinal cohort studies should help to illuminate the antigenic and genetic relationships between the circulating dengue serotypes. These will also help to elucidate the possible correlation between incidence of disease and serotype diversity identified by our model.

20Finally, the enhancement of transmission as opposed to an increase in susceptibility may well be attributed to a different mechanism to ADE. A recent study shows DHF patients exhibit both a higher viral load and a slower rate of viral clearance in comparison with DF patients. For the latter, viral loads drop to undetectable levels by the first day of defervescence, but they remain at high levels in DHF patients. This difference in the rate of loss of viraemia is better explained by interference in the development of the appropriate immune response either by original antigenic sin or altered peptide ligand antagonism. If so, the two parameters in our model could represent two different positive interactions between serotypes, rather than two aspects of ADE alone. Whichever the biological basis, the decomposition into effects on susceptibility and transmission affords a parsimonious explanation for the observed variety in the epidemiological behaviour of dengue.