Breaking the symmetry: immune enhancement increases persistence of dengue viruses in the presence of asymmetric transmission rates

Temporary Cross-Immunity as a Plausible Driver of Asynchronous Cycles of Dengue Serotypes

Many infectious diseases exist as multiple variants, with interactions between variants potentially driving epidemiological dynamics. These diseases include dengue, which infects hundreds of millions of people every year and exhibits complex multi-serotype dynamics. Antibodies produced in response to primary infection by one of the four dengue serotypes can produce a period of temporary cross-immunity (TCI) to infection by other serotypes. After this period, the remaining antibodies can facilitate the entry of heterologous serotypes into target cells, thus enhancing severity of secondary infection by a heterologous serotype. This represents antibody-dependent enhancement (ADE). In this study, we analyze an epidemiological model to provide novel insights into the importance of TCI and ADE in producing cyclic outbreaks of dengue serotypes. Our analyses reveal that without TCI, such cyclic outbreaks are synchronous across serotypes and only occur when ADE produces high transmission rates. In contrast, the presence of TCI allows asynchronous cycles of serotypes by inducing a time lag between recovery from primary infection by one serotype and secondary infection by another, with such cycles able to occur without ADE. Our results suggest that TCI is a fundamental driver of asynchronous cycles of dengue serotypes and possibly other multi-variant diseases.

Applying a multi-strain dengue model to epidemics data

Dengue disease transmission is a complex vector-borne disease, mainly due to the co-circulation of four serotypes of the virus. Mathematical models have proved to be a useful tool to understand the complexity of this disease. In this work, we extend the model studied by Esteva et al., 2003, originally proposed for two serotypes, to four circulating serotypes. Using epidemic data of dengue fever in Iquitos (Peru) and San Juan (Puerto Rico), we estimate numerically the co-circulation parameter values for selected outbreaks using a bootstrap method, and we also obtained the Basic Reproduction Number, R₀, for each serotype, using both analytical calculations and numerical simulations. Our results indicate that the impact of co-circulation of serotypes in population dynamics of dengue infection is such that there is a reduced effect from DENV-3 to DENV-4 in comparison to no-cross effect for epidemics in Iquitos. Concerning San Juan epidemics, also comparing to no-cross effect, we also observed a reduced effect from the predominant serotype DENV-3 to both DENV-2 and DENV-1 epidemics neglecting the very small number of cases of DENV-4.

Mathematical models for dengue fever epidemiology: A 10-year systematic review

Mathematical models have a long history in epidemiological research, and as the COVID-19 pandemic progressed, research on mathematical modeling became imperative and very influential to understand the epidemiological dynamics of disease spreading. Mathematical models describing dengue fever epidemiological dynamics are found back from 1970. Dengue fever is a viral mosquito-borne infection caused by four antigenically related but distinct serotypes (DENV-1 to DENV-4). With 2.5 billion people at risk of acquiring the infection, it is a major international public health concern. Although most of the cases are asymptomatic or mild, the disease immunological response is complex, with severe disease linked to the antibody-dependent enhancement (ADE) - a disease augmentation phenomenon where pre-existing antibodies to previous dengue infection do not neutralize but rather enhance the new infection. Here, we present a 10-year systematic review on mathematical models for dengue fever epidemiology. Specifically, we review multi-strain frameworks describing host-to-host and vector-host transmission models and within-host models describing viral replication and the respective immune response. Following a detailed literature search in standard scientific databases, different mathematical models in terms of their scope, analytical approach and structural form, including model validation and parameter estimation using empirical data, are described and analyzed. Aiming to identify a consensus on infectious diseases modeling aspects that can contribute to public health authorities for disease control, we revise the current understanding of epidemiological and immunological factors influencing the transmission dynamics of dengue. This review provide insights on general features to be considered to model aspects of real-world public health problems, such as the current epidemiological scenario we are living in.

The ratio of single to co-colonization is key to complexity in interacting systems with multiple strains

The high number and diversity of microbial strains circulating in host populations have motivated extensive research on the mechanisms that maintain biodiversity. However, much of this work focuses on strain-specific and cross-immunity interactions. Another less explored mode of pairwise interaction is via altered susceptibilities to co-colonization in hosts already colonized by one strain. Diversity in such interaction coefficients enables strains to create dynamically their niches for growth and persistence, and "engineer" their common environment. How such a network of interactions with others mediates collective coexistence remains puzzling analytically and computationally difficult to simulate. Furthermore, the gradients modulating stability-complexity regimes in such multi-player endemic systems remain poorly understood. In a recent study (Madec & Gjini, Bulletin of Mathematical Biology, 82), we obtained an analytic representation for N-type coexistence in an SIS epidemiological model with co-colonization. We mapped multi-strain dynamics to a replicator equation using timescale separation. Here, we examine what drives coexistence regimes in such co-colonization system. We find the ratio of single to co-colonization, µ, critically determines the type of equilibrium and number of coexisting strains, and encodes a trade-off between overall transmission intensity R 0 and mean interaction coefficient in strain space, k. Preserving a given coexistence regime, under fixed trait variation, requires balancing between higher mean competition in favorable environments, and higher cooperation in harsher environments, and is consistent with the stress gradient hypothesis. Multi-strain coexistence tends to steady-state attractors for small µ, whereas as µ increases, dynamics tend to more complex attractors. Following strain frequencies, evolutionary dynamics in the system also display contrasting patterns with µ, interpolating between multi-stable and fluctuating selection for cooperation and mean invasion fitness, in the two extremes. This co-colonization framework could be applied more generally, to study invariant principles in collective coexistence, and to quantify how critical shifts in community dynamics get potentiated by mean-field and environmental gradients.

Vaccine- and natural infection-induced mechanisms that could modulate vaccine safety

A degraded/dysfunctional immune system appears to be the main determinant of serious/fatal reaction to viral infection (for COVID-19, SARS, and influenza alike). There are four major approaches being employed or considered presently to augment or strengthen the immune system, in order to reduce adverse effects of viral exposure. The three approaches that are focused mainly on augmenting the immune system are based on the concept that pandemics/outbreaks can be controlled/prevented while maintaining the immune-degrading lifestyles followed by much of the global population. The fourth approach is based on identifying and introducing measures aimed at strengthening the immune system intrinsically in order to minimize future pandemics/outbreaks. Specifically, the four measures are: 1) restricting exposure to virus; 2) providing reactive/tactical treatments to reduce viral load; 3) developing vaccines to prevent, or at least attenuate, the infection; 4) strengthening the immune system intrinsically, by a) identifying those factors that contribute to degrading the immune system, then eliminating/reducing them as comprehensively, thoroughly, and rapidly as possible, and b) replacing the eliminated factors with immune-strengthening factors. This paper focuses on vaccine safety. A future COVID-19 vaccine appears to be the treatment of choice at the national/international level. Vaccine development has been accelerated to achieve this goal in the relatively near-term, and questions have arisen whether vaccine safety has been/is being/will be compromised in pursuit of a shortened vaccine development time. There are myriad mechanisms related to vaccine-induced, and natural infection-induced, infections that could adversely impact vaccine effectiveness and safety. This paper summarizes many of those mechanisms.

Understanding dengue virus evolution to support epidemic surveillance and counter-measure development

Dengue virus (DENV) causes a profound burden of morbidity and mortality, and its global burden is rising due to the co-circulation of four divergent DENV serotypes in the ecological context of globalization, travel, climate change, urbanization, and expansion of the geographic range of the Ae.aegypti and Ae.albopictus vectors. Understanding DENV evolution offers valuable opportunities to enhance surveillance and response to DENV epidemics via advances in RNA virus sequencing, bioinformatics, phylogenetic and other computational biology methods. Here we provide a scoping overview of the evolution and molecular epidemiology of DENV and the range of ways that evolutionary analyses can be applied as a public health tool against this arboviral pathogen.

The Role of Serotype Interactions and Seasonality in Dengue Model Selection and Control: Insights from a Pattern Matching Approach

The epidemiology of dengue fever is characterized by highly seasonal, multi-annual fluctuations, and the irregular circulation of its four serotypes. It is believed that this behaviour arises from the interplay between environmental drivers and serotype interactions. The exact mechanism, however, is uncertain. Constraining mathematical models to patterns characteristic to dengue epidemiology offers a means for detecting such mechanisms. Here, we used a pattern-oriented modelling (POM) strategy to fit and assess a range of dengue models, driven by combinations of temporary cross protective-immunity, cross-enhancement, and seasonal forcing, on their ability to capture the main characteristics of dengue dynamics. We show that all proposed models reproduce the observed dengue patterns across some part of the parameter space. Which model best supports the dengue dynamics is determined by the level of seasonal forcing. Further, when tertiary and quaternary infections are allowed, the inclusion of temporary cross-immunity alone is strongly supported, but the addition of cross-enhancement markedly reduces the parameter range at which dengue dynamics are produced, irrespective of the strength of seasonal forcing. The implication of these structural uncertainties on predicted vulnerability to control is also discussed. With ever expanding spread of dengue, greater understanding of dengue dynamics and control efforts (e.g. a near-future vaccine introduction) has become critically important. This study highlights the capacity of multi-level pattern-matching modelling approaches to offer an analytic tool for deeper insights into dengue epidemiology and control.

The fluctuations of multi-serotype infectious diseases are often highly irregular and hard to predict. Previous theoretical approaches have attempted to disentangle the drivers that may underlie this behaviour in dengue dynamics with variable success. Here, we examine the role of such drivers using a pattern-oriented modelling (POM) approach. In POM, multiple patterns observed at different scales are used to test a model’s proficiency in capturing real-world dynamics. We examined dengue models with combinations of cross-immunity, cross-enhancement, seasonal fluctuations in the transmission rate, and with sensitivity analyses of asymmetric transmission rates between serotypes as well as the possibility for four subsequent heterologous infections. We demonstrate the ability of POM to model dynamical drivers that have gone unnoticed in single pattern or synthetic likelihood approaches. Further, our results present a determining role of seasonality in the selection and operation of these processes in governing dengue dynamics, in particular when full, heterologous immunity is assumed to occur after a secondary infection. We show that this structural model uncertainty can have important practical significance, as demonstrated by the differences in control efforts required to disrupt transmission. These results highlight the importance of localised model selection and calibration using multiple data-matching, as well as taking explicit account of model uncertainty in predicting and planning control efforts for multi-serotype diseases.

How direct competition shapes coexistence and vaccine effects in multi-strain pathogen systems

We describe an integrated modeling framework for understanding strain coexistence in polymorphic pathogen systems. Previous studies have debated the utility of neutral formulations and focused on cross-immunity between strains as a major stabilizing mechanism. Here we convey that direct competition for colonization mediates stable coexistence only when competitive abilities amongst pathogen clones satisfy certain pairwise asymmetries. We illustrate our ideas with nested SIS models of single and dual colonization, applied to polymorphic pneumococcal bacteria. By fitting the models to cross-sectional prevalence data from Portugal (before and after the introduction of a seven-valent pneumococcal conjugate vaccine), we are able to not only statistically compare neutral and non-neutral epidemiological formulations, but also estimate vaccine efficacy, transmission and competition parameters simultaneously. Our study highlights that the response of polymorphic pathogen populations to interventions holds crucial information about strain interactions, which can be extracted by suitable nested modeling.

Are we modelling the correct dataset? Minimizing false predictions for dengue fever in Thailand

Models describing dengue epidemics are parametrized on disease incidence data and therefore high-quality data are essential. For Thailand, two different sources of long-term dengue data are available, the hard copy data from 1980 to 2005, where hospital admission cases were notified, and the electronic files, from 2003 to the present, where clinically classified forms of disease, i.e. dengue fever, dengue haemorrhagic fever, and dengue shock syndrome, are notified using separate files. The official dengue notification data, provided by the Bureau of Epidemiology, Ministry of Public Health in Thailand, were cross-checked with dengue data used in recent publications, where an inexact continuous time-series was observed to be consistently used since 2003, affecting considerably the model dynamics and its correct application. In this paper, numerical analysis and simulation techniques giving insights on predictability are performed to show the effects of model parametrization by using different datasets.

Characterizing the likelihood of dengue emergence and detection in naïve populations

Background

Vector-borne disease transmission is dependent on the many nuances of the contact event between infectious and susceptible hosts. Virus acquisition from a viremic human to a susceptible mosquito is often assumed to be nearly perfect and almost always uniform across the infectious period. Dengue transmission models that have previously addressed variability in human to vector transmission dynamics do not account for the variation in infectiousness of a single individual, and subsequent infection of naïve mosquitoes. Understanding the contribution of this variability in human infectiousness is especially important in the context of introduction events where an infected individual carries the virus into a population of competent vectors. Furthermore, it could affect the ability to detect an epidemic (and the timing of detection) following introduction.

Methods

We constructed a stochastic, compartmental model to describe the heterogeneity of human viremia and calculate the probability of a successful introduction, taking into account the viremia level (and thus acquisition potential) of the index case on, and after, the day of introduction into a susceptible population and varying contact rates between the human and mosquito populations. We then compared the results of this model with those generated by a simpler model that has the same average infectiousness but only a single infectious class.

Results

We found that the infectivity of the index case as well as the contact rate affected the probability of emergence, but that contact rate had the most significant effect. We also found that the interaction between contact rate and the infectiousness of the index case affected the time to detection relative to the peak of the epidemic curve. Additionally, when compared to our model that accounts for variable infectiousness, a model with a single infectious class underestimates the probability of emergence and transmission intensity.

Conclusion

Understanding the interplay between individual human heterogeneity of infectiousness and the rate of contact with the vector population will be important when predicting the likelihood, detection, and magnitude of an outbreak.

The impact of coinfections and their simultaneous transmission on antigenic diversity and epidemic cycling of infectious diseases

Epidemic cycling in human infectious diseases is common; however, its underlying mechanisms have been poorly understood. Much effort has been made to search for external mechanisms. Multiple strains of an infectious agent were usually observed and coinfections were frequent; further, empirical evidence indicates the simultaneous transmission of coinfections. To explore intrinsic mechanisms for epidemic cycling, in this study we consider a multistrain Susceptible-Infected-Recovered-Susceptible epidemic model by including coinfections and simultaneous transmission. We show that coinfections and their simultaneous transmission widen the parameter range for coexistence and coinfections become popular when strains enhance each other and the immunity wanes quickly. However, the total prevalence is nearly independent of these characteristics and approximated by that of one-strain model. With sufficient simultaneous transmission and antigenic diversity, cyclical epidemics can be generated even when strains interfere with each other by reducing infectivity. This indicates that strain interactions within coinfections and cross-immunity during subsequent infection provide a possible intrinsic mechanism for epidemic cycling.

Potential opportunities and perils of imperfect dengue vaccines

Dengue vaccine development efforts have focused on the development of tetravalent vaccines. However, a recent Phase IIb trial of a tetravalent vaccine indicates a protective effect against only 3 of the 4 serotypes. While vaccines effective against a subset of serotypes may reduce morbidity and mortality, particular profiles could result in an increased number of cases due to immune enhancement and other peculiarities of dengue epidemiology. Here, we use a compartmental transmission model to assess the impact of partially effective vaccines in a hyperendemic Thai population. Crucially, we evaluate the effects that certain serotype heterogeneities may have in the presence of mass-vaccination campaigns. In the majority of scenarios explored, partially effective vaccines lead to 50% or greater reductions in the number of cases. This is true even of vaccines that we would not expect to proceed to licensure due to poor or incomplete immune responses. Our results show that a partially effective vaccine can have significant impacts on serotype distribution and mean age of cases.