Stability and bifurcations in a model of antigenic variation in malaria

Plasmodium knowlesi: a superb in vivo nonhuman primate model of antigenic variation in malaria

Antigenic variation in malaria was discovered in Plasmodium knowlesi studies involving longitudinal infections of rhesus macaques (M. mulatta). The variant proteins, known as the P. knowlesi Schizont Infected Cell Agglutination (SICA) antigens and the P. falciparum Erythrocyte Membrane Protein 1 (PfEMP1) antigens, expressed by the SICAvar and var multigene families, respectively, have been studied for over 30 years. Expression of the SICA antigens in P. knowlesi requires a splenic component, and specific antibodies are necessary for variant antigen switch events in vivo. Outstanding questions revolve around the role of the spleen and the mechanisms by which the expression of these variant antigen families are regulated. Importantly, the longitudinal dynamics and molecular mechanisms that govern variant antigen expression can be studied with P. knowlesi infection of its mammalian and vector hosts. Synchronous infections can be initiated with established clones and studied at multi-omic levels, with the benefit of computational tools from systems biology that permit the integration of datasets and the design of explanatory, predictive mathematical models. Here we provide an historical account of this topic, while highlighting the potential for maximizing the use of P. knowlesi - macaque model systems and summarizing exciting new progress in this area of research.

Chimera states in multi-strain epidemic models with temporary immunity

We investigate a time-delayed epidemic model for multi-strain diseases with temporary immunity. In the absence of cross-immunity between strains, dynamics of each individual strain exhibit emergence and annihilation of limit cycles due to a Hopf bifurcation of the endemic equilibrium, and a saddle-node bifurcation of limit cycles depending on the time delay associated with duration of temporary immunity. Effects of all-to-all and non-local coupling topologies are systematically investigated by means of numerical simulations, and they suggest that cross-immunity is able to induce a diverse range of complex dynamical behaviors and synchronization patterns, including discrete traveling waves, solitary states, and amplitude chimeras. Interestingly, chimera states are observed for narrower cross-immunity kernels, which can have profound implications for understanding the dynamics of multi-strain diseases.

Global properties of nested network model with application to multi-epitope HIV/CTL dynamics

Mathematical modeling and analysis can provide insight on the dynamics of ecosystems which maintain biodiversity in the face of competitive and prey-predator interactions. Of primary interests are the underlying structure and features which stabilize diverse ecological networks. Recently Korytowski and Smith (Theor Ecol 8(1):111-120, 2015) proved that a perfectly nested infection network, along with appropriate life history trade-offs, leads to coexistence and persistence of bacteria-phage communities in a chemostat model. In this article, we generalize their model in order to apply it to the within-host dynamics virus and immune response, in particular HIV and CTL (Cytotoxic T Lymphocyte) cells. Our model can describe sequential viral escape from dominant immune responses and rise in subdominant immune responses, consistent with observed patterns of HIV/CTL evolution. We find a Lyapunov function for the system which leads to rigorous characterization of persistent viral and immune variants, along with informing upon equilibria stability and global dynamics. Results are interpreted in the context of within-host HIV/CTL evolution and numerical simulations are provided.

Capturing the dynamics of pathogens with many strains

Pathogens that consist of multiple antigenic variants are a serious public health concern. These infections, which include dengue virus, influenza and malaria, generate substantial morbidity and mortality. However, there are considerable theoretical challenges involved in modelling such infections. As well as describing the interaction between strains that occurs as a result cross-immunity and evolution, models must balance biological realism with mathematical and computational tractability. Here we review different modelling approaches, and suggest a number of biological problems that are potential candidates for study with these methods. We provide a comprehensive outline of the benefits and disadvantages of available frameworks, and describe what biological information is preserved and lost under different modelling assumptions. We also consider the emergence of new disease strains, and discuss how models of pathogens with multiple strains could be developed further in future. This includes extending the flexibility and biological realism of current approaches, as well as interface with data.

Novel methods for quantifying individual host response to infectious pathogens for genetic analyses

We propose two novel approaches for describing and quantifying the response of individual hosts to pathogen challenge in terms of infection severity and impact on host performance. The first approach is a direct extension of the methodology for estimating group tolerance (the change in performance with respect to changes in pathogen burden in a host population) to the level of individuals. The second approach aims to capture the dynamic aspects of individual resistance and tolerance over the entire time course of infections. In contrast to the first approach, which provides a means to disentangle host resistance from tolerance, the second approach focuses on the combined effects of both characteristics. Both approaches provide new individual phenotypes for subsequent genetic analyses and come with specific data requirements. In particular, both approaches rely on the availability of repeated performance and pathogen burden measurements of individuals over the time course of one or several episodes of infection. Consideration of individual tolerance also highlights some of the assumptions hidden within the concept of group tolerance, indicating where care needs to be taken in trait definition and measurement.

Synchronous versus asynchronous oscillations for antigenically varying Plasmodium falciparum with host immune response

We consider a deterministic intra-host model for Plasmodium falciparum (Pf) malaria infection, which accounts for antigenic variation between n clonal variants of PfEMP1 and the corresponding host immune response (IR). Specifically, the model separates the IR into two components, specific and cross-reactive, respectively, in order to demonstrate that the latter can be a mechanism for the sequential appearance of variants observed in actual Pf infections. We show that a strong variant-specific IR relative to the cross-reactive IR favours the asynchronous oscillations (sequential dominance) over the synchronous oscillations in a number of ways. The decay rate of asynchronous oscillations is smaller than that for the synchronous oscillations, allowing for the parasite to survive longer. With the introduction of a delay in the stimulation of the IR, we show that only a small delay is necessary to cause persistent asynchronous oscillations and that a strong variant-specific IR increases the amplitude of the asynchronous oscillations.