A Probabilistic Synthesis of Malaria Epidemiology: Exposure, Infection, Parasite Densities, and Detection

A Probabilistic Synthesis of Malaria Epidemiology: Exposure, Infection, Parasite Densities, and Detection | medRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-P4HH5NV'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search A Probabilistic Synthesis of Malaria Epidemiology: Exposure, Infection, Parasite Densities, and Detection John M. Henry , Austin R. Carter , View ORCID Profile Sean L. Wu , View ORCID Profile David L. Smith doi: https://doi.org/10.1101/2025.03.24.25324561 John M. Henry a College of the Environment, University of Washington , Seattle, WA, USA b Institute for Health Metrics and Evaluation, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Austin R. Carter b Institute for Health Metrics and Evaluation, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sean L. Wu b Institute for Health Metrics and Evaluation, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Sean L. Wu David L. Smith c Department of Health Metrics Sciences, School of Medicine, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for David L. Smith For correspondence: smitdave{at}uw.edu Abstract Full Text Info/History Metrics Data/Code Preview PDF Abstract The epidemiology of Plasmodium falciparum malaria presents a unique set of challenges due to the complicated dynamics of infection, immunity, disease, and detection. Studies of malaria epidemiology commonly measure malaria parasite densities or prevalence, but since malaria is so complex with so many factors to consider, a complete mathematical synthesis of malaria epidemiology has been elusive. Here, we take a new approach. From a simple model of malaria exposure and infection in human cohorts as they age, we develop random variables describing the multiplicity of infection (MoI) and the age of infection (AoI). Next, using the MoI and AoI distributions, we develop random variables describing parasite densities, parasite counts, and detection. We also derived a random variable describing the age of the youngest infection (AoY), which can be used to compute approximate parasite densities in complex infections. Finally, we derive a simple system of differential equations with hybrid variables that track the mean MoI, AoI and AoY, and we show it matches the complex probabilistic system with reasonable accuracy. We can thus compute the state of any individual chosen at random from the population in two ways. The same approach – pairing random variables and hybrid models – can be extended to model other features of malaria epidemiology, including disease, malaria immunity, treatment and chemoprotection, and infectiousness. The computational simplicity of hybrid models has some advantages over compartmental models and stochastic individual-based models, and with the supporting probabilistic framework, provide a sound basis for a synthesis of observational malaria epidemiology. Significance The probabilistic framework presented here represents a step forward in the theory for observational malaria epidemiology. Using these methods, one can connect the heterogeneous dynamics of individual infections to both observational and process uncertainty in estimates of population-level metrics such as cross-sectional surveys of prevalence by light microscopy. This serves as a basis for more complex simulation studies using the present models and extensions that include factors such as fever, treatment, and transmission. The epidemiology of human malaria caused by infection with Plasmodium falciparum is complex. Exposure to P. falciparum , as measured by the annual entomological inoculation rate (aEIR), varies by several orders of magnitude seasonally, across years, and among locations ( 1 – 3 ). The duration of detectable parasitemia for simple malaria infections varies from a few weeks to more than a year, and parasite densities are highly variable among individuals in populations depending on their age, location, and the time of year ( 4 – 6 ). In malaria therapy – a peculiar set of observations in which exposure was fully known – detection, disease and infectiousness all varied with the age of infection (AoI) ( 7 , 8 ). While these studies provide a rich source of data describing malaria, controlled infections in those institutionalized patient populations differed from natural settings in many ways. In natural settings, neither immunity nor infection prevents reinfection (also called superinfection), so malaria infections can be complex – the number of distinct parasite lineages has been called the multiplicity of infection (MoI) ( 9 – 11 ), and a person could experience dozens of bouts of moderate to severe malaria over a lifetime ( 3 , 12 ). Acquired and innate immunity modify infection and the risk and severity of disease, but protection is weak, has a poor memory, and accumulates slowly with age and exposure ( 13 – 15 ). The consensus is that malaria immunity varies in some way with age, cumulative exposure, and recent exposure ( 3 , 12 , 16 ), but there are some important data gaps and lingering uncertainty. A synthesis is needed. An important feature of malaria is the long, variable time course of infection ( 5 , 6 , 17 ). After the infective bite and a short period in the liver, parasites begin a cycle of asexual replication in red blood cells, all while their numbers are increasing geometrically ( 18 ). Immunity to malaria and other factors stop the geometric growth phase, and then the infection begins a long and variable chronic phase, during which parasite densities tend to decline and are detected less frequently until eventually the infection ends ( 8 ). Untreated malaria infections can, nevertheless, persist for a few weeks to several years. The time course of malaria infections is highly variable over the time course of an infection, and the features of those infections vary over a lifetime. Substantial attention has focused on the var genes, and var gene switching as a mechanism that allows parasite infections to persist, but parasites can also modify gene transcription in other ways ( 17 , 19 ), and there are other genes that play an important role in development of immunity. The net effect of the interplay between parasites and the human immune response is that parasite densities are highly variable over the time course of an infection, and predictable only in a statistical sense ( 8 , 17 ). Throughout an infection, the probability of detecting parasites varies, so accurate measurement of the distribution of basic features of malaria, such as the duration of infections, has been difficult. Parasite densities in blood fluctuate over several orders of magnitude, and only a fraction of blood is examined. The sensitivity of parasite counts or other estimates of parasite density and malaria prevalence varies across diagnostic methods ( 20 , 21 ). Available data thus measure the duration of patentcy, but evidence for very long carriage times comes from infections that have recrudesced after persisting for decades. A second feature of malaria epidemiology is re-exposure and superinfection. Acquired immunity to malaria does not prevent reinfection, so after subsequent exposure, so multiple, distinct parasite lineages can coinfect a single host. It is also possible for a single infective bite to transmit multiple distinct parasite clones. Several studies have attempted to measure the MoI, but the metric is difficult to interpret given the probability of detection ( 9 ). Understanding parasite density distributions in complex infections ( i . e ., those with MoI > 1) presents a daunting challenge, in part, because parasite densities are so variable. With highly variable infections, highdensity infections in one sample can mask other low-density infections. The complexity of infections may be partially revealed by sampling the same person repeatedly over time ( 22 , 23 ). The drive to understand malaria has led to a proliferation of models that address various aspects of malaria, including simple models exploring various aspects of transmission, and comprehensive and synthetic studies using individual-based and intra-host simulation models of malaria infections ( 24 – 27 ). Because the number of factors affecting malaria is large, the state space describing malaria is high-dimensional. Given the uncertainty about malaria epidemiology and its complexity, a conceptual and quantitative synthesis of malaria epidemiology must be developed around some theory and mathematical models. A practical challenge, from a mathematical or computational perspective, has been how to represent the multi-dimensional state space for infection and immunity – as models consider more factors, the size of the state space expands combinatorically and becomes difficult to manage. One way of handling the complex state space has been through development of individual-based simulation models, which replicate the basic patterns in models that simultaneously deal with the complex time course of infections, superinfection, and waxing and waning immunity ( 28 – 36 ). While IBMs are a useful computational tool for many purposes, they are often presented as black boxes due to the complexity of their design and the probabilistic nature of their output. What is needed is a study of relevant detail to help identify the critical dimensions without losing information and prioritize studies that could fill critical gaps ( ? ). Here, we take a probabilistic approach to understanding malaria epidemiology, including parasite density distributions, by formulating random variables that track the AoI, the MoI, and parasite density distributions in human birth cohorts as they age. In the simple infections observed in malaria therapy patients, most aspects of infections were well described statistically by the AoI, the time elapsed since the infective bite for each parasite clone infecting a host: the fraction persisting, the fraction still detectable, the mean and variance of logged parasite densities, the risk of fever and disease, and infectiousness ( 8 ). Using the AoI a unifying concept, we developed models of understand the density distributions of parasite clones in cohorts of humans as they age; these distribution link exposure, the MoI, the AoI, parasite densities in complex infections, and detection. Using this approach, we show that malaria epidemiology, including parasite densities, can be understood in terms of a new random variable describing the age of the youngest infection (AoY). This approach does not rely on stochastic simulations, and we have found it more amenable to analysis. Finally, for analysis and synthesis, we developed hybrid variables describing the mean AoI and mean AoY, and we show how a simple systems of differential equations holds most of the information about malaria infection dynamics, so computationally simple models of malaria epidemiology can be used to understand complex patterns minimal information loss. Altogether, the methods present a basis for a conceptual synthesis of observational malaria epidemiology. Dynamical Models for Cohorts In the following, we analyze a simple model for malaria infection in a host cohort as it ages. The analysis extends a queuing model for super-infection called M/M/ ∞ ( Fig 1 ), which has played an important role in malaria epidemiology ( 37 – 40 ). A core assumption of the model is that infections clear independently, and in the spirit of that model, we also assume that parasite densities over the time courses of two or more parasite infections are also independent. In this study, we have focused on understanding how malaria epidemiology changes in young children as they age, depending on exposure. Download figure Open in new tab Fig. 1. A diagram of the queuing model M/M/ ∞ (top), and the master equations (bottom), in which ζ i denotes the proportion of the population with MoI = i . MoI increments by one if an infection occurs at rate h and decrements by one if an infection clears at rate ri . In the following, we let a denote the age of a host birth cohort as it ages over time. A parasite clone is defined to include all the progeny of the parasites arising from a single infective mosquito bite. We are thus assuming that each bite would transmit a single clone, an assumption that could be relaxed in subsequent analyses. We let α denote the AoI: the age of a parasite clone, or the time elapsed since the infective bite that initiated the clonal infection. The random variables are all computed conditioned on exposure, which we define as the force of infection (FoI), the hazard rate for infection experienced by a cohort of susceptible humans as it ages. Starting from a model for the infection dynamics, we developed random variables describing the MoI, the AoI, and parasite densities as a function of the FoI (see Fig 2 ). We then model the observational process, which translates parasite densities into parasite counts and measures of prevalence. A synthesis was achieved through two extensions: first, we defined a new random variable describing the age of the youngest infection (AoY); next, we derived hybrid variables describing the mean AoI and the mean AoY. Download figure Open in new tab Fig. 2. An overview of the probabilistic model and development of the random variables. Model assumptions are drawn as hexagons. The random variables are developed around a model for exposure and infection dynamics. We pass the daily FoI in the population, h ( t ), and an assumption about exposure by age, ω ( a ), so that the model input is the daily FoI for a particular cohort as it ages, h τ ( a ). We then solve the equation describing the density of parasite clones of age α , distributed among a cohort of hosts as it ages, z ( α, a | h ). We compute random variables describing the MoI, M τ ( a | h ), and the AoI A τ ( a | h ). We assume mean log 10 parasite densities, ξ , are a function of the AoI ( F µ ( α )) with a probability density function given by Ω( µ ). Assuming independence, we then derive random variables describing the clonal density of parasites in blood, P τ ( a | h ), and total parasite densities, summed over all clonal infections, B τ ( a | h ). Under a model for sampling parasites from a fraction q of the blood, S ( B, q ,…), We compute the observed PR, ,and the distribution of counts, The primary computational methods for the paper are found in open source R package on github, called ramp.falciparum. Derivations for most of the formulas and some proofs are found in the Supporting Information. Vignettes in ramp.falciparum walk through the derivations, provide pseudo-code to compute most random variables, and provide numerical validation for most of the claims ( https://github.com/dd-harp/ramp.falciparum ). The software was designed to be both modular and extensible, and there are plans to extend this probabilistic approach to model immunity, disease, treatment and chemoprotection, and other factors. Exposure and Infection Differences in the rate of exposure by age and season are modelled through the FoI, a time- and age-dependent average hazard rate for new infections in a susceptible population. To make the analysis as general and applicable as possible, we have developed methods to generate a broad class of functions to describe the FoI in a population over time, h ( t ). In doing so, we are explicitly reducing the FoI to a model input, or a trace function . In many studies with mechanistic models, the FoI computed as the result of a dynamical process reflecting parasite transmission from humans to mosquitoes, and blood feeding by mosquitoes. These trace-driven simulations are one way to decouple different parts of a simulation when one aspect of a simulation can be replaced by an external input (trace) from other parts of the simulation on which it depends ( 41 ? ). Using this approach, we can investigate malaria under almost any pattern of exposure without the need to specify exactly how that pattern of exposure might have arisen. The Force of Infection The GitHub R package, ramp.falciparum , includes a flexible yet “realistic” set of functions to model exposure ( Fig 3 and Supporting Information). We let be a parameter that sets the overall level of exposure. We assume exposure has a seasonal pattern, S ( t ), and a temporal trend, T ( t ): Download figure Open in new tab Fig. 3. FoI trace functions developed to simulate exposure combine four elements: a parameter that sets the overall intensity, ) functions that modify exposure by age, ω ( a ); 3) functions that describe the seasonal pattern, S ( t ); and 4) and a long-term trend, T ( t ). The left panel illustrates representative functions for agemodification, and the right panel shows representative seasonal patterns. Motivated by the evidence for age-dependent exposure, the rate of exposure in the cohort is assumed to change with age. We use the age-weighting function, ω ( a ), consistent with the available evidence ( 42 , 43 ). Each cohort is indexed by its birth day, τ . The FoI at age a for the cohort born on day τ , which is at time t = a + τ , is: We can thus develop arbitrary exposure patterns and compare outcomes for cohorts of different ages in the same population, of the same age at different times. Infection Dynamics Let z ( α, a ) denote the density of parasite clonal infections of age α in a host cohort of age a . We assume infections clear at the constant rate r . Since infections in M/M/ ∞ are independent, we can track the dynamics for the AoI of all parasite clones with the equation, with the boundary condition z τ (0, a ) = h τ ( a ). We note that the age of the host birth cohort sets an upper limit for the age of the parasite clones, so 0 < α < a . The solution to Eq. 2 , which describes density of infections of age α in a host cohort of age a , is given by the formula: In the Supporting Information, we prove that Eq. 3 is a solution to Eq. 2 . Age and Multiplicity of Infection Under the assumptions of M/M/ ∞, all the parasite clones are distributed independently among hosts. Importantly, the mean MoI in the cohort at age a is the integral of this density function: Elsewhere, it has been shown that the distribution of the MoI converges to a Poisson ( 44 ). Since the MoI is zero at birth (a degenerate Poisson distribution), the distribution of the MoI remains Poisson as the cohort ages (in this model), and mean MoI, m τ ( a ) completely describes the probability density function for the MoI: We let F M denote the CDF, and the random number function. In ramp.falciparum, we verify that, for the same function h τ ( a ), that Eq. 4 has the same mean as M/M/ ∞. Since f M is Poisson, true prevalence ( p ) is the proportion with MoI greater than 0, which can be computed directly from m τ ( a ): While z τ ( α, a ) describes the density of clonal infections in the population, a random variable describing the density of the AoI, A , within hosts in particular cohort at a given age a is: We let F A denote the CDF, and  the random number function. From a trace function, h τ ( a ), we can directly compute and visualize the full probability distribution function for the AoI and the MoI in a cohort as it ages ( Fig 4 ). In ramp.falciparum , we developed functions to compute the probability density functions, cumulative distribution functions, and random numbers for a cohort of age a born on day τ directly from h ( t ). Download figure Open in new tab Fig. 4. The distribution of the a) multiplicity of infection (MoI), b) the age of infection (AoI), and c) age of the youngest clonal infection (AoY) in a cohort as it ages from birth up through 5 years old. The FoI was sinusoidal, with a population mean aFoI of 5 infections, per year. An important feature of this model is that the distribution of the AoI is affected by the seasonal pattern and age, but not by its overall intensity. Later, we will show that the patterns tend to be driven by youngest infections, which is affected by overall intensity ( Fig 4b,c ). Parasite Densities in Infections Parasite densities in the blood are often measured clinically and reported as parasite counts, and counts or infection prevalence are commonly used as covariates in observational studies. Early studies of malaria parasite counts observed enormous variability in a single host over time as well as among hosts ( 4 , 5 ). Evidence describing parasite densities improved in quality, granularity, and quantity during malaria therapy, when patients with neuro-syphilis were deliberately infected with malaria and carefully monitored ( 45 , 46 ). More recently, a vast amount of detail about the early stages of infection has been learned through controlled human malaria infections (CHMI), where patients are deliberately exposed, often as part of clinical trials for developmental phase drugs or vaccines ( 18 , 47 ). Here, we present a model for parasite densities and counts using the concept of the AoI, and we fit it to the malaria therapy data using maximum likelihood estimation. In the malaria therapy data, the variance in parasite counts was a power-function of the mean, suggesting that the variance in logged parasite would be a simple linear function of the logged mean. We thus developed a realistic probability distribution function that explicitly modeled the variance as a function of the mean. Under a model of detection, we fitted the model to the malariotherapy data ( Fig. 5 ). Download figure Open in new tab Fig. 5. The surface plots the model fit to the malaria-therapy data (Supplement-S5). To model parasite density distributions, we proceed in several steps: first, we formulate models for mean log 10 parasite densities in a parasite clonal infection (denoted µ ) as a function of its AoI ( α ) or F µ ( α ); second, we develop a family of probability distributions describing the number of infected red blood cells (log iRBCs, ξ ) infected by each clone (denoted Ω) as a function of the mean ( µ ), or ξ ~ Ω( µ ); third, we combine these two models, and we develop a random variable describing the density of parasite clones infecting a host cohort of age a born on day τ , denoted P τ ( a | h ) ( Fig 5 ); and fourth, we develop a new random variable that sums parasite densities over the MoI to describe total iRBCs in an infection, B τ ( a | h ) ( Fig. 6 ). Download figure Open in new tab Fig. 6. Top) the expected value for parasite densities, F µ ( α ) (black, Eq. 8 ) and the distribution of AoI, f A ( α ) (black line with colored points, Eq. 7 ), and the AoY, f Y ( α ) (grey line) for one particular cohort at 5 years old. Bottom) The total distribution of parasite clonal densities across all parasite infections, f P ( α ), integrates parasite density distributions, Ω( µ ) with expected value µ = F µ ( α ) (top panel), over the distribution of the AoI. For each colored point in the top, we plot the density distribution of of parasites clones, Ω, of each age from a host cohort. c) Parasite density distribution for a clone ( f P ) for MoI= N ( f N , N=2,…,10, colored purple to green), for the full distribution ( f B , orange). We also plot (grey). Time Course of Infections Consistent with most of observational epidemiology, and based on an analysis of malariotherapy data ( 8 ), we assume that the densities of each parasite clone vary probabilistically with its age, α . Generically, we let µ = F µ ( α ) denote some function describing the expected mean log 10 parasite clonal densities in a simple infection as a function of the AoI. The functions could be modified to include variables describing cumulative exposure, recent exposure, or maternal protection. We developed and analyzed one particular model for parasite densities that matches the features of malaria documented across multiple previous studies ( Fig 5 ). The basic features of the function describe the time course of an infection: after an infectious bite, parasites spend around 6 days in the liver ( 18 ? ); parasites enter the blood on or after day 7 of the infection and the asexual blood-stage parasites begin a period of geometric growth; symptoms can appear around day 9-10, and parasite populations become patent ( i . e . detectable) around day 14; the geometric growth phase of an infection ends around day δ ≈ 20 days when parasite densities peak, and the infection begins a chronic phase during which parasite densities fluctuate overall but tend to decline as the infection ages. Infections can persist for several months, or even years; the duration of infections is highly variable. The model is described by the equation: A key feature of the model is thus that infections have two distinct phases: an acute growth phase when parasite densities are increasing geometrically (for α < δ ); and a chronic phase when parasite densities decline ( Fig. 5 ). The parameter values here are provided as approximate values that are roughly consistent with observations, a starting point for further investigation about about which substantial uncertainty remains. This framework provides a rigorous way of formulating and testing ideas in silico . Clonal Infections Total parasite densities in a clonal infection must be bounded between 1 and the total number of RBCs, which is a function of host size and age. The log 10 number of RBCs is denoted b ( a ). We model Ω( µ ), where µ is the expected value of log 10 iRBCs (denoted ξ ) as a function of the AoI; µ = F µ ( α ). To compute Ω, we modified a beta distribution to: 1) return a value between 0 and b ( a ); with mean log 10 densities equal to µ ; and 3) the variance, which computed as a function of the mean (Supporting Information, and ramp.falciparum ). A random variable describing log 10 parasite densities in each parasite clonal infection cohort is denoted P τ ( a ), with a probability density function: Fig 6 . We let F P denote the CDF, and the random number function. The software, ramp.falciparum , was developed around a flexible design to accommodate alternative model formulations, ( i . e ., for both F µ and Ω). Ω could thus be developed for any reasonable double bounded continuous probability distribution function family. Parasite Densities To model parasite densities in complex infections, some assumption about interactions among parasites within a host is required. Do they compete for resources or immune-free space, interfere directly, facilitate one another or fluctuate independently? While there are many compelling candidate models that could be developed and tested, we have analyzed parasite densities under the simple and parsimonious assumption, consistent with the other assumptions of M/M/ ∞, that parasite densities fluctuate independently. Alternative models based on other assumptions can be explored in future studies. Computing total densities for potentially complex infections thus entails adding up the parasite densities of all parasite clones, the number of which is given by the distribution of the MoI. Let B denote a random variable describing log 10 parasite densities summed over all clonal infections in a host cohort of age a , denoted B τ ( a ) with a probability density function: We let F B denote the distribution function, and the random number function. We were unable to derive a simple formula for B τ ( a ), but we have developed algorithms to compute it and to verify the computation ( ramp.falciparum ). This equation can be understood through the pseudo-code that describes how to draw random values of parasite densities from individual hosts in a cohort of age a born on day τ , by drawing a random number of clonal infections with random values for clonal parasite densities (Algorithm 1). [1] m τ ( a | h ) Compute Mean MoI Draw Random MoI in 1 to Draw Random Draw Random Clone Sum of Clone Densities is Undefined Host is Uninfected We developed functions to compute the density and distribution functions for complex infections, as a convolution over a distribution of MoI ( Fig 7b , Fig 6c , ramp.falciparum , and Supporting Information). As an intermediate step, we developed random variables for parasite densities in infections with MoI= N , with distribution function f N ( ξ ; a, τ | h ), and functions to draw random number N -tuples from P τ ( a ). We verified the convolution-based algorithms against these density and distribution functions in complex infections with fixed MoI ( Fig 6c , Supporting Information, ramp.falciparum ). Download figure Open in new tab Fig. 7. Parasite densities and counts for the same cohort illustrated in Fig. 4. a ) The distribution of log 10 parasite densities across all cohorts of parasites, log 10 P ; b) The distribution of log 10 parasite densities in the host cohort as it ages, log 10 B ; and c) The distribution of log 10 parasite counts, conditioned on counting at least one parasite, log 10 C . Compare the patterns in panels a) and b) here to the patterns in panels b) and c) in Fig 4 . Mean log 10 parasite densities in complex infections can be understood in terms of the AoY. Since B is a log 10 transformed value of a sum of values of 10 P , it follows that B is close to the maximum value of P , and we can bound the sum: The total density would be near the upper bound only if all clonal parasite densities were approximately equal. If N clones had approximately the same density, then the total density would be N times the maximum value of P, but after log transforming parasite densities, the effect adds log 10 N . For N = 2, parasite densities are approximately 0.3 higher in log 10 parasite densities, and for N = 5, approximately 0.7 higher. Compared to the daily variance in parasite densities observed in malariotherapy patients ( Fig 5 ), these numbers are small ( 8 ). Superinfection does have a big effect on parasite density distributions in this model, but mainly because parasite densities a sum of multiple infections, which is bounded below by the one with the highest density (See Fig 6c and compare Fig 7a and Fig 7b ). Since parasite density distributions fluctuate over several powers of ten, the sums of parasite densities can be guided by an understanding of extreme values. It follows that B is likely to follow the youngest infection, not the AoI, which motivated development of a new random variable describing the youngest infection ( Fig 6 ). A potential problem with this model formulation is that while the parasite densities for individual clones are bounded well below 10 b ( a ) , total parasite densities summed over multiple clones need not be. To deal with extreme chance events in the upper tail, the mean of the distribution of parasite densities was bounded well below the total number of iRBCs, such that computed parasite densities exceeding 1% of all RBCs are extremely rare. Although this does break independence of infections, this distortion arises rarely and does not have a large impact on calculations on detection as this saturates quickly at high densities. Practically if parasite densities in a patient were close to the absolute maximum, catastrophic health risks would be likely. Hyperparasitemia, a form of severe malaria defined by having more that 5% infected RBCs in immunologically-naïve individuals, is associated with severe, acute clinical symptoms. Age of the Youngest Infection (AoY) We can take advantage of the extreme value arithmetic that govern parasite density distributions, as suggested by Eq. 11 , by developing a new random variable, Y τ ( a ), that describes the age of the youngest infection (AoY) with the distribution function: Logically, since the youngest infections tend to have the highest densities (unless a parasite is in the liver or early acute growth phase), expected parasite densities in infections would be dominated either by the oldest clone in the acute phase or the youngest parasite clone in the chronic phase. In this derivation of the AoY, we ignore the inaccuracies introduced by the early acute phase infections. To compute Y , we derive a formula for the AoY in terms of the AoI density and distribution functions, and the MoI: We let F Y denote the CDF, and Ŷ the random number function. Based on the arguments made above, we can use Y to compute a lower bound on B τ ( a ), a distribution that would be observed if the entire distribution was determined by the youngest clone, with a distribution function: In effect, is computed by substituting Y for A (also see Fig. 6a,b ), and We note that the distribution of the AoY has a very different shape than the AoI (contrast Fig 4b,c , and see Fig. 5 ). While the approximation is not exact, it provides a useful heuristic; the shape of Y , not A , explains the shapes of the distribution of parasite densities for simple and complex infections (compare Fig 7b,c ). Hybrid Variables Using random variables, we can compute parasite density distributions in a cohort of any age given a function that describes the history of exposure (the FoI, h τ ( a )), an assumption about parasite clearance, r , and a function that describes parasite densities as infections age, Ω( F µ ( α )). A disadvantage of the probabilistic approach is that it can be computationally intensive, and it is difficult integrate these equations into models for the transmission dynamics and control of malaria. To address these limitations, we developed new hybrid variables that are computationally much simpler ( 44 ). It has been shown that the dynamic changes in mean MoI, summarizing much of the behavior of the system of equations in Fig 1 , can be modeled as a simple ordinary differential equation ( 44 ). Noting that m τ ( a ) is the mean MoI, we also define the mean AoI, x τ ( a ) and the mean AoY, y τ ( a ): Next, we derived equations describing the dynamics of the variables x and y (Supporting Information), which when combined with the hybrid model for the MoI ( 44 ), allow us to compute changes over time using a simple system of differential equations: The function ϕ ( r, m ) captures the dynamics of changes in the AoY in complex infections, in which the loss of the youngest infection doesn’t result in clearance but an increase in the age of the youngest clone (see the Supporting Information). Further, the mean MoI, AoI, and AoY can be used compute approximate mean logged parasite densities: While the approximations are not exact, they are reasonably close ( Fig. 8 ). These approximations, and the underlying heuristic, thus serves as a useful conceptual bridge connecting exposure to mean parasite densities. Download figure Open in new tab Fig. 8. Mean log10 parasite densities (red hues) and clone densities (blue hues) computed exactly and using the hybrid variable approximations Eqs 18 . In both cases, the hybrid variable approximation is the slightly lighter shade. The approximation performs comparatively poorly for the first 200 days because of the strongly non-linear relationship between the AoI and parasite densities around the peak. Thereafter, the approximations differ slightly from the exact values. Measuring Infection The random variables describing parasite infections and parasite densities are unobservable (latent) states, so to understand malaria infections, we need to model the observation process. What is the relationship between our latent random variables and observable quantities? Studies in malaria epidemiology often rely on parasite prevalence or parasite counts. We thus developed random variables describing parasite counts from examination of a blood sample, from which we derived a related variable describing the probability of counting at least one parasite. We developed random variables for parasite counts on light microscopy. Sources of sampling error include variation in parasite densities in the sampled blood, and errors introduced through counting, which involves examination of a set of fields on the blood slide ( 48 ). These differences can be affected by differences among light microscopists. We note, as an aside, that we have not considered a set of issues associated with false positives. Similar kinds of concerns would apply to modeling PCR or RDTs, but new random variables would need to be developed around a sampling model that was adapted to the diagnostic method. It is reasonable to expect that the fraction of samples that were positive would be different based on the sampling model as well as the sensitivity of the method ( e . g ., by light microscopy, qPCR or rapid diagnostic tests which detect a biomarker that is believed to be strongly correlated with parasite densities). Sampling To diagnose malaria by light microscopy, a small fraction of total blood is sampled and examined for parasites. Parasites are thus often not detected, even when they are present. The relationship between the true parasite densities and parasite counts depends on the total fraction of blood examined, q . For light microscopy, where approximately 1 ≈ µ l, of blood is drawn and a fraction of it examined. If s 6 were the log number of RBCs in a blood sample, then the fraction of RBCs examined would be approximately: The log 10 parasite counts in a sample, ,is a random variable that depends on the log 10 parasite densities in blood, ξ . We let S denote a probability density function family to model raw counts as a function of parasite densities, conditioned on being infected. We let Q denote the probability of counting at least one parasite, conditioned on sampling an infected host: Next, we define a function describing the distribution of log 10 positive counts ( i . e ., conditioned on detecting at least one parasite) In developing models for counts in complex infections, we identified some issues that we could not resolve from a priori considerations: What probability density function family is appropriate for sampling parasites and parasite counts data? and should counts be modeled for independent detection of individual parasite clones, or for aggregate parasite densities in an infection? On the choice of sampling distributions, if blood were drawn and parasites counted repeatedly, would the counts data follow a Poisson distribution, negative binomial distribution, or some other probability distribution function family? On the question of whether we should treat parasite densities in the aggregate, or whether to treat clones independently, we could find no strong a priori reason to prefer one or the other, so we developed both methods, which very similar, but slightly different predictions. In the associated software, ramp.falciparum, functions were developed to model detection based on a Poisson or negative binomial family of distributions, and the software design is extensible, making it possible to add other suitable probability density function families. In the following, we use the negative binomial distribution. We fit a model for parasite density distributions and detection to the malaria-therapy data using maximum likelihood estimation (Supplement-S5, Fig 9 ). The analysis gave different answers depending on the value of the size parameter of the negative binomial parameter, and parameters describing the variance in parasite population densities among individuals as a function of the mean. If the size parameter was larger than one, then the sampling was reasonably accurate, andwparasite densities among individuals were more variable. If the size parameter was smaller (≈ 0.2), then the sampling was extremely noisy, and parasite population densities were less variable. A few studies that resampled infected individuals at short intervals suggested the sampling is reasonably accurate ( 22 , 23 ). Download figure Open in new tab Fig. 9. Sampling error for the detection process. a) For values of log 10 parasite densities ( x -axis), vertical slices show the proportion of parasite counts in a range defined by powers of 10 ( y -axis). To make this plot, we assume the proportion of blood examined is q = 10 −6 , and the size parameter was fitted to the malariatherapy data. b) For each integer value of B , we show a histogram of a vertical slice of the figure above it. Each one is a histogram, binned by powers of 10 (colors and categories match the top figure) for log 10 ( B ) for integer values from 3 to 11. Download figure Open in new tab Fig. 10. The prevalence over time computed using random variables ( Eq. 23 ), and an approximation using hybrid variables ( Eq. 28 ). The errors in this approximation might be too large to be useful, but the comparison is a useful heuristic for understanding detection, and it can serve as a basis for translating hybrid variables into observations using the probabilistic framework. Counts For detection based on aggregate parasite densities, we can simply use C to compute parasite counts over the distribution of parasite densities, B τ ( a ) ( e . g . see Fig 9 ): Since B describes the distribution of zero counts, and positive counts for slides from infected hosts, the number of zero counts over all slides examined is thus inflated by p τ ( a ). Detection Using B τ ( a ), apparent prevalence is Alternatively, we developed a model around detection of independent clones, which aligns in spirit with the model M/M/ ∞. We define a function, F Q , that computes the probability of detection as a function of the age of infection α , using the density distribution function, Ω and F µ : If present, a clone would be detected with probability: so the apparent MoI, ,if individual clones could be counted, would be Poisson distributed with mean and alternative formula for detection of parasites in a complex infection is thus: Hybrid Detection Dynamics We can also show that mean logged parasite counts for parasite clones scale in a reasonably simple way, with interesting counts distributions ( Fig 7c ). We can also show that we can compute a good approximation to the observed prevalence using hybrid variables: Once again, these approximations are useful as a heuristic to explore understanding the relationship between true and observed prevalence. Discussion Using a probabilistic framework, we showed how exposure over time translates into patterns of infection and malaria parasite densities and prevalence in cohort of children as they age. The approach emphasizes that malaria infection and all the metrics we use to measure it hold a memory of recent exposure. Since parasite densities are highly stochastic, the memory of exposure is probabilistic, observable only by sampling from cohorts in a population. We also showed that this memory is well-approximated by a simple system with hybrid variables that describe the mean MoI, AoI and AoY. The quantitative logic of these two approaches, relating probabilistic models to hybrid variables, makes it possible to draw a clear line from exposure to parasite density distributions. In the statistical sense, parasite densities are wellpredicted by the AoY. While the quantitative logic connecting the two is long and sometimes mathematically cumbersome, it illustrates how a detailed and highly realistic probabilistic description of a system with variables captures much of the complexity. While the models we have presented lack some of the core features of malaria – immunity, disease, treatment and chemoprotection, and infectiousness – these challenges can also be handled by extending the probabilistic approach. The approaches are planned for future studies. Altogether, this synthesis represents a new way of understanding malaria epidemiology represents a new way of understanding malaria epidemiology. Although the medical approach to understanding malaria epidemiology relies heavily on parasite densities, it is difficult to understand patterns without grappling with the large uncertainty associated with this data and the underlying process. Without doing so, we are likely to overlook existing but messy relationships such as that between malaria incidence and febrile infections in the growth phase. In developing the functions, µ ( α ), our attention was drawn to a minor, but perhaps critically important feature of infections: some features of malaria epidemiology could occur early in the growth phase, while parasite densities are still low ( 8 ). There is a narrow window of time, probably less a week, in which the growth phase of a parasite was observed to cause disease, but when parasite densities are comparatively low ( 8 ). Severe disease may be therefore be associated with these growth-phase acute infections. The observation begs a serious question about using parasite densities as a covariate in studies of malaria. Symptoms that arise from parasite infections in the acute phase may not undetectable or they present with low densities, while parasite densities from old infections may be a source of noise. This study takes a semi-mechanistic approach that starts from the premise that parasite densities can be computed, in a statistical sense, using information about infection age. The probabilistic nature of the model acknowledges the nature of infection age as predicting - not determining - the parasitemia in an individual chosen at random. Because we wished to model observable quantities such as measured prevalence which depend on the latent state of the parasitemia in an infected individual chosen at random, it was not enough to know the expectation alone. As we model each individual as being equally likely to be measured here for simplicity, the predicted rate of detection - a function of parasitemia - must be weighted according to its likelihood of being observed. In order to model variables which can be predicted through other functions of parasitemia, such as the probability of fever, more information about the parasite density distribution may be required. The approach does not address the question of what drives the shape of curves describing parasite densities as they age, but by finding functions that emulate the behavior of other models, including individual-based within-host models, it could provide a rigorous basis for comparison. While there are many parameters in the full model, each sub-model is parametrically simple. To describe parasite densities with a high degree of realism, we have we developed several linked sub-models, each one of which is highly realistic. The dynamics of infection under the model M/M/ ∞ have one fixed parameter ( r ) describing the loss of infections, a decay rate that describes how prevalence reflects the memory of past exposure. The model for parasite densities and detection requires the most: five parameters describing the trajectory of the expected value (two parameters are bounds, one describes time in the liver, one describes parasite clonal density at peak, and one describes the decline in parasite densities with the parasite clone’s age. The distribution function required three parameters to express the variance as a function of the mean. Two parameters were needed to describing detection, of which one is context dependent ( i . e . fraction of blood examined). The methods and the software can be extended to other mathematical forms. In presenting this framework, we emphasize that the claims are not specific inferential claims but general ones. The general claim is that models of this sort naturally replicate the sorts of quantitative patterns that have been observed in observational malaria epidemiological studies. The claims are not about any specific mechanism, but the method provides a useful way of understanding and modeling malaria epidemiology. Since the semi-Markov models are extending an observed pattern, it is not surprising that they should give seemingly reasonable results. The approach taken here also emphasizes an a posteriori approach with a priori tools ( 49 ). That is, after observing patterns in the data for parasitemia which is difficult to model (a posteriori), we created an a priori model of infection age and used the statistical model linking it to parasitemia to then produce a simpler model of parasitemia. To emphasize the structure and flow of the model, we largely implemented computational approaches here. Approximations suggest some robustness in the structure of the models as output is weakened only slightly by using the age of the youngest infection. Many of the complex features of malaria epidemiology are made more coherent when seen as emergent properties of the biological needs of the malaria parasite. Sexual reproduction occurs in the mosquito, and outcrossing occurs only if the multiplicity of infection (MoI) is greater than one. Frequent outcrossing, in turn, supports a genetically diverse parasite population, which is partly responsible for the weakness of the immune responses and other features of malaria epidemiology. Most outcomes of being exposed to malaria are expected to vary with an individual’s age and previous history of exposure and infection. The core challenges for applied malaria epidemiology are to understand or predict likely outcomes and to identify modifiable factors affecting the burden of malaria and the prospects for malaria elimination. Framed broadly, outcomes are affected by host infection status ( e . g . parasite densities), acquired immunity resulting from a lifetime of exposure, other host factors (anemia, genetics, nutrition, etc .), parasite factors (genotype and phenotype), and case management with anti-malarial drugs. The difficulty of understanding and quantifying outcomes is increased by the challenge of simultaneously measuring recent exposure, infection, disease, immunity acquired over a lifetime of exposure ( e . g . through serology), and other important covariates in enough settings to quantify their relationships. The contemporary view of malaria epidemiology has been pieced to-gether from a heterogeneous set of studies conducted over roughly 140 years, but the number of studies that has been done is small compared to the number that would be needed. Important gaps in our understanding of malaria epidemiology remain. Conclusion Parasite infection dynamics and detection are but a starting point for understanding malaria epidemiology. We have introduced a new modeling framework which bridges the gap between stochastic individual-level behavior and population-level observable quantities. This is a novel simulation method for infectious disease that takes a non-compartmental approach to managing the complexity and uncertainty inherent in the mechanisms that control parasitemia. However as with all models, this complexity and stochasticity can hinder our ability to isolate mechanism almost as much as the system it is meant to model. Future directions include extending the framework to model immunity, disease, treatment, and transmission. We also plan to explore the effects of disease-prompted treatment on age-prevalence curves; and comparing the evolutionary dynamics of the parasite with strain-specific immunity. Data Availability All data and code are available online from: https://dd-harp.github.io/ramp.falciparum/ or https://github.com/dd-harp/memory.core https://dd-harp.github.io/ramp.falciparum/ https://github.com/dd-harp/memory.core Author Contributions JMH and DLS designed the study, carried out all the analyses. JMH, DLS, ARC, and SLW wrote the first draft. All authors contributed to the final draft. ACKNOWLEDGEMENTS This research was supported by a grant from the National Institute of Allergies and Infectious Diseases (R01 AI163398), which supported DLS, JMH, and ARC. Additional support for DLS from the Bill and Melinda Gates Foundation (INV 030600). References 1. ↵ Simon I. Hay , D. J. Rogers , J. F. Toomer , and Robert W. Snow . Annual Plasmodium fal-ciparum entomological inoculation rates (EIR) across Africa: Literature survey, internet access and review . Trans. R. Soc. Trop. Med. Hyg ., 94 ( 2 ): 113 – 127 , March 2000 . ISSN 0035-9203 . doi: 10.1016/s0035-9203(00)90246-3 . OpenUrl CrossRef PubMed 2. Louise A. Kelly-Hope and F. Ellis McKenzie . The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa . Malar. J ., 8 : 19 , January 2009 . ISSN 1475-2875 . doi: 10.1186/1475-2875-8-19 . OpenUrl CrossRef PubMed 3. ↵ Walther H. Wernsdorfer and Ian A. McGregor Louis Molineaux , D. A. Muir , H. C. Spencer , and W. H. Wernsdorfer . The epidemiology of malaria and its measurement . In Walther H. Wernsdorfer and Ian A. McGregor , editors, Malaria: Principles and Practice of Malariology , pages 999 – 1089 . Churchill Livingstone, New York , March 1988 . ISBN 978-0-443-02417-7 . 4. ↵ Ronald Ross and David Thomson . Some enumerative studies on malaria fever . Proceedings of the Royal Society of London B: Biological Sciences , 83 ( 562 ): 159 – 173 , December 1910 . OpenUrl CrossRef 5. ↵ Don E. Eyles and Martin D. Young . The duration of untreated or inadequately treated Plasmodium falciparum infections in the human host . Bulletin of the World Health Organization , 10 ( 4 ): 557 – 619 , December 1951 . ISSN 0042-9686 . OpenUrl 6. ↵ Elizabeth A. Ashley and Nicholas J. White . The duration of Plasmodium falciparum infections . Malar. J ., 13 ( 1 ): 500 , December 2014 . ISSN 1475-2875 . doi: 10.1186/1475-2875-13-500 . OpenUrl CrossRef PubMed 7. ↵ Ingrid M. Daey Ouwens , C. Elisabeth Lens , Aernoud T. L. Fiolet , Alewijn Ott , Peter J. Koehler , Piet A. Kager , and Willem M. A. Verhoeven . Malaria Fever Therapy for General Paralysis of the Insane: A Historical Cohort Study . European Neurology , 78 ( 1-2 ): 56 – 62 , June 2017 . ISSN 0014-3022 . doi: 10.1159/000477900 . OpenUrl CrossRef PubMed 8. ↵ John M. Henry , Austin Carter , and David L. Smith . Infection age as a predictor of epidemiological metrics for malaria . Malar. J ., 21 ( 1 ): 117 , April 2022 . ISSN 1475-2875, 1475-2875 . doi: 10.1186/s12936-022-04134-5 . OpenUrl CrossRef PubMed 9. ↵ Daibin Zhong , Cristian Koepfli , Liwang Cui , and Guiyun Yan . Molecular approaches to determine the multiplicity of Plasmodium infections . Malar. J ., 17 ( 1 ): 172 , April 2018 . ISSN 1475-2875 . doi: 10.1186/s12936-018-2322-5 . OpenUrl CrossRef PubMed 10. Martina Eldh , Ulf Hammar , David Arnot , Hans-Peter Beck , André Garcia Anne Liljander , Odile Mercereau-Puijalon , Florence Migot-Nabias , Ivo Mueller , Francine Ntoumi , Amanda Ross , Thomas Smith , Klara Sondén , Manijeh Vafa Homann , Victor Yman , Ingrid Felger , and Anna Färnert . Multiplicity of Asymptomatic Plasmodium falciparum Infections and Risk of Clinical Malaria: A Systematic Review and Pooled Analysis of Individual Participant Data . The Journal of Infectious Diseases , 221 ( 5 ): 775 – 785 , February 2020 . ISSN 0022-1899 . doi: 10.1093/infdis/jiz510 . OpenUrl CrossRef PubMed 11. ↵ Kristan Alexander Schneider , Henri Christian Junior Tsoungui Obama , George Kamanga , Loyce Kayanula , and Nessma Adil Mahmoud Yousif . The many definitions of multiplicity of infection . Frontiers in Epidemiology , 2 , October 2022 . ISSN 2674-1199 . doi: 10.3389/fepid.2022.961593 . OpenUrl CrossRef 12. ↵ Walther H. Wernsdorfer and Ian A. McGregor Louis Molineaux . The epidemiology of human malaria as an explanation of its distribution, including some implications for its control . In Walther H. Wernsdorfer and Ian A. McGregor , editors, Malaria: Principles and Practice of Malariology , pages 913 – 998 . Churchill Livingstone , New York , March 1988 . ISBN 978-0-443-02417-7 . 13. ↵ Siske S. Struik and Eleanor M. Riley . Does malaria suffer from lack of memory? Immunol. Rev ., 201 ( 1 ): 268 – 290 , October 2004 . ISSN 0105-2896 . doi: 10.1111/j.0105-2896.2004.00181.x . OpenUrl CrossRef PubMed Web of Science 14. Denise L. Doolan , Carlota Dobaño, and J. Kevin Baird . Acquired immunity to malaria . Clin. Microbiol. Rev ., 22 ( 1 ): 13 – 36 , Table of Contents, January 2009 . ISSN 0893-8512, 1098-6618 . doi: 10.1128/CMR.00025-08 . OpenUrl Abstract / FREE Full Text 15. ↵ Caterina Guinovart , Carlota Dobaño, Quique Bassat , Augusto Nhabomba , Llorenç Quintó , Maria Nélia Manaca , Ruth Aguilar , Mauricio H. Rodríguez , Arnoldo Barbosa , John J. Aponte , Alfredo G. Mayor , Montse Renom , Cinta Moraleda , David J. Roberts , Evelin Schwarzer , Peter N. L. Souëf , Louis Schofield , Chetan E. Chitnis , Denise L. Doolan , and Pedro L. Alonso . The role of age and exposure to Plasmodium falciparum in the rate of acquisition of naturally acquired immunity: a randomized controlled trial . PLoS One , 7 ( 3 ): e32362 , March 2012 . ISSN 1932-6203 . doi: 10.1371/journal.pone.0032362 . OpenUrl CrossRef PubMed 16. ↵ Brian M. Greenwood , K. Marsh , and R. Snow . Why do some African children develop severe malaria? Parasitol. Today , 7 ( 10 ): 277 – 281 , October 1991 . ISSN 0169-4758 . doi: 10.1016/0169-4758(91)90096-7 . OpenUrl CrossRef PubMed Web of Science 17. ↵ Lauren M. Childs and Caroline O. Buckee . Dissecting the determinants of malaria chronicity: why within-host models struggle to reproduce infection dynamics . J. R. Soc. Interface , 12 ( 104 ): 20141379 – 20141379 , March 2015 . ISSN 1742-5689 . doi: 10.1098/rsif.2014.1379 . OpenUrl CrossRef PubMed 18. ↵ Luc E. Coffeng , Cornelus C. Hermsen , Robert W. Sauerwein , and Sake J. de Vlas . The Power of Malaria Vaccine Trials Using Controlled Human Malaria Infection . PLOS Computational Biology , 13 ( 1 ): e1005255 , January 2017 . ISSN 1553-7358 . doi: 10.1371/journal.pcbi.1005255 . OpenUrl CrossRef PubMed 19. ↵ Carolina M. Andrade , Hannah Fleckenstein , Richard Thomson-Luque , Safiatou Doumbo , Nathalia F. Lima , Carrie Anderson , Julia Hibbert , Christine S. Hopp , Tuan M. Tran , Shanping Li , Moussa Niangaly , Hamidou Cisse , Didier Doumtabe , Jeff Skinner , Dan Sturdevant , Stacy Ricklefs , Kimmo Virtaneva , Muhammad Asghar , Manijeh Vafa Homann , Louise Turner , Joana Martins , Erik L. Allman , Marie-Esther N’Dri , Volker Winkler , Manuel Llinás , Catherine Lavazec , Craig Martens , Anna Färnert , Kassoum Kayentao , Aissata Ongoiba , Thomas Lavstsen , Nuno S. Osório , Thomas D. Otto , Mario Recker , Boubacar Traore , Peter D. Crompton , and Silvia Portugal . Increased circulation time of Plasmodium falciparum underlies persistent asymptomatic infection in the dry season . Nat. Med ., October 2020 . ISSN 1078-8956, 1546-170X . doi: 10.1038/s41591-020-1084-0 . OpenUrl CrossRef PubMed 20. ↵ Lucy C. Okell , Azra C. Ghani , Emily Lyons , and Christopher J. Drakeley . Submicroscopic infection in Plasmodium falciparum-endemic populations: a systematic review and meta-analysis . J. Infect. Dis ., 200 ( 10 ): 1509 – 1517 , November 2009 . ISSN 0022-1899, 1537-6613 . doi: 10.1086/644781 . OpenUrl CrossRef PubMed Web of Science 21. ↵ Bonnie Mappin , Ewan Cameron , Ursula Dalrymple , Daniel J. Weiss , Donal Bisanzio , Samir Bhatt , and Peter W. Gething . Standardizing Plasmodium falciparum infection prevalence measured via microscopy versus rapid diagnostic test . Malar. J ., 14 ( 1 ): 460 , November 2015 . ISSN 1475-2875 . doi: 10.1186/s12936-015-0984-9 . OpenUrl CrossRef PubMed 22. ↵ Anna Farnert , Georges Snounou , Ingegerd Rooth , and Anders Bjorkman . Daily Dynamics of Plasmodium falciparum Subpopulations in Asymptomatic Children in a Holoendemic Area . The American Journal of Tropical Medicine and Hygiene , 56 ( 5 ): 538 – 547 , May 1997 . ISSN 0002-9637 . doi: 10.4269/ajtmh.1997.56.538 . OpenUrl Abstract / FREE Full Text 23. ↵ Véronique Delley , Paul Bouvier , Norman Breslow , Ogobara Doumbo , Issaka Sagara , Mahamadou Diakite , Anne Mauris , Amagana Dolo , and André Rougemont . What does a single determination of malaria parasite density mean? A longitudinal survey in Mali . Tropical Medicine & International Health , 5 ( 6 ): 404 – 412 , 2000 . ISSN 1365-3156 . doi: 10.1046/j.1365-3156.2000.00566.x . _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1046/j.1365-3156.2000.00566.x . OpenUrl CrossRef PubMed Web of Science 24. ↵ Robert C. Reiner , T. Alex Perkins , Christopher M. Barker , Tianchan Niu , Luis Fernando Chaves , Alicia M. Ellis , Dylan B. George , Arnaud Le Menach , Juliet R. C. Pulliam , Donal Bisanzio , Caroline Buckee , Christinah Chiyaka , Derek A. T. Cummings , Andres J. Garcia , Michelle L. Gatton , Peter W. Gething , David M. Hartley , Geoffrey Johnston , Eili Y. Klein , Edwin Michael , Steven W. Lindsay , Alun L. Lloyd , David M. Pigott , William K. Reisen , Nick Ruktanonchai , Brajendra K. Singh , Andrew J. Tatem , Uriel Kitron , Simon I. Hay , Thomas W. Scott , and David L. Smith . A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970-2010 . J. R. Soc. Interface , 10 ( 81 ): 20120921 , April 2013 . OpenUrl CrossRef PubMed 25. Neal R. Smith , James M. Trauer , Manoj Gambhir , Jack S. Richards , Richard J. Maude , Jonathan M. Keith , and Jennifer A. Flegg . Agent-based models of malaria transmission: A systematic review . Malar. J ., 17 ( 1 ): 299 , August 2018 . ISSN 1475-2875 . doi: 10.1186/s12936-018-2442-y . OpenUrl CrossRef PubMed 26. Louis Molineaux and Klaus Dietz . Review of intra-host models of malaria . Parassitologia , 41 ( 1-3 ): 221 – 231 , September 1999 . ISSN 0048-2951 . OpenUrl PubMed 27. ↵ Flavia Camponovo , Tamsin E. Lee , Jonathan R. Russell , Lydia Burgert , Jaline Gerardin , and Melissa A. Penny . Mechanistic within-host models of the asexual Plasmodium falciparum infection: a review and analytical assessment . Malaria Journal , 20 ( 1 ): 309 , July 2021 . ISSN 1475-2875 . doi: 10.1186/s12936-021-03813-z . OpenUrl CrossRef PubMed 28. ↵ Jamie T. Griffin , T. Deirdre Hollingsworth , Lucy C. Okell , Thomas S. Churcher , Michael White , Wes Hinsley , Teun Bousema , Chris J. Drakeley , Neil M. Ferguson , María-Gloria Basáñez , and Azra C. Ghani . Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies . PLoS Med ., 7 ( 8 ): e1000324 , August 2010 . ISSN 1549-1277, 1549-1676 . doi: 10.1371/journal.pmed.1000324.s010 . OpenUrl CrossRef PubMed 29. Jamie T. Griffin , Neil M. Ferguson , and Azra C. Ghani . Estimates of the changing ageburden of Plasmodium falciparum malaria disease in sub-Saharan Africa . Nat. Commun ., 5 : 3136 , 2014 . ISSN 2041-1723 . doi: 10.1038/ncomms4136 . OpenUrl CrossRef PubMed 30. Philip A. Eckhoff . A malaria transmission-directed model of mosquito life cycle and ecology . Malar. J ., 10 ( 1 ): 303 , October 2011 . ISSN 1475-2875 . doi: 10.1186/1475-2875-10-303 . OpenUrl CrossRef PubMed 31. Philip A. Eckhoff . Malaria parasite diversity and transmission intensity affect development of parasitological immunity in a mathematical model . Malar. J ., 11 ( 1 ): 419 , December 2012 . doi: 10.1186/1475-2875-11-419 . OpenUrl CrossRef PubMed 32. Philip A. Eckhoff . Mathematical models of within-host and transmission dynamics to determine effects of malaria interventions in a variety of transmission settings . Am. J. Trop. Med. Hyg ., 88 ( 5 ): 817 – 827 , May 2013 . ISSN 0002-9637 . doi: 10.4269/ajtmh.12-0007 . OpenUrl Abstract / FREE Full Text 33. Edward A. Wenger and Philip A. Eckhoff . A mathematical model of the impact of present and future malaria vaccines . Malar. J ., 12 ( 1 ): 126 , April 2013 . ISSN 1475-2875 . doi: 10.1186/1475-2875-12-126 . OpenUrl CrossRef PubMed 34. Thomas A. Smith , Gerry F. Killeen , Nicolas Maire , Amanda Ross , Louis Molineaux , Fabrizio Tediosi , Guy Hutton , Juerg Utzinger , Klaus Dietz , and Marcel Tanner . Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: Overview . Am. J. Trop. Med. Hyg ., 75 ( 2 Suppl ): 1 – 10 , August 2006 . ISSN 0002-9637 . doi: 10.4269/ajtmh.2006.75.2_suppl.0750001 . OpenUrl FREE Full Text 35. Thomas A. Smith , N. Maire , A. Ross , M. Penny , Nakul Chitnis , A. Schapira , A. Studer , B. Genton , C. Lengeler , F. Tediosi , D. de Savigny , and Marcel Tanner . Towards a comprehensive simulation model of malaria epidemiology and control . Parasitology , 135 ( 13 ): 1507 – 1516 , November 2008 . ISSN 0031-1820, 1469-8161 . doi: 10.1017/S0031182008000371 . OpenUrl CrossRef PubMed Web of Science 36. ↵ Thomas A. Smith , Amanda Ross , Nicolas Maire , Nakul Chitnis , Alain Studer , Diggory Hardy , Alan Brooks , Melissa Penny , and Marcel Tanner . Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine . PLoS Med ., 9 ( 1 ): e1001157 , January 2012 . ISSN 1549-1277, 1549-1676 . doi: 10.1371/journal.pmed.1001157 . OpenUrl CrossRef PubMed 37. ↵ K. Dietz , Louis Molineaux , and A. Thomas . A malaria model tested in the African savannah . Bull World Health Organ , 50 ( 3-4 ): 347 – 357 , 1974 . OpenUrl PubMed Web of Science 38. Paul E. M. Fine . Superinfection a problem in formulating a problem . Tropical Diseases Bulletin , 75 ( 6 ): 475 – 488 , 1975 . OpenUrl 39. Norman T. J. Bailey . The Biomathematics of Malaria. Charles Griffin & Company Ltd., Oxford , 1982 . ISBN 978-0-85264-266-5 . 40. ↵ W. Wernsdorfer and Ian A. McGregor Klaus Dietz . Mathematical models for transmission and control of malaria . In W. Wernsdorfer and Ian A. McGregor , editors, Malaria: Principles and Practice of Malariology , pages 1091 – 1133 . Churchill Livingstone , New York , March 1988 . ISBN 978-0-443-02417-7 . 41. ↵ Lawrence M. Leemis and Stephen Keith Park . Discrete-event simulation: A first course . Pearson Prentice Hall Upper Saddle River, NJ , 2006 . 42. ↵ G. R. Port , P. F. L. Boreham , and Joan H. Bryan . The relationship of host size to feeding by mosquitoes of the Anopheles gambiae Giles complex (Diptera: Culicidae) . Bulletin of Entomological Research , 70 ( 01 ): 133 – 144 , March 1980 . ISSN 0007-4853, 1475-2670 . doi: 10.1017/S0007485300009834 . OpenUrl CrossRef 43. ↵ P. Carnevale , J. L. Frézil , M. F. Bosseno , F. Le Pont , and J. Lancien . [The aggressiveness of Anopheles gambiae A in relation to the age and sex of the human subjects] . Bull. World Health Organ ., 56 ( 1 ): 147 – 154 , 1978 . ISSN 0042-9686 . OpenUrl PubMed Web of Science 44. ↵ Ingemar Nåsell . Hybrid Models of Tropical Infections, volume 59 of Lecture Notes in Biomathematics. Springer-Verlag, Berlin, 1st edition edition , 1985 . ISBN 978-0-387-15978-2 . 45. ↵ Georges Snounou and Jean-Louis Pérignon . Malariotherapy–insanity at the service of malariology . Adv. Parasitol ., 81 : 223 – 255 , 2013 . ISSN 0065-308X . doi: 10.1016/B978-0-12-407826-0.00006-0 . OpenUrl CrossRef PubMed 46. ↵ William E. Collins and Geoffrey M. Jeffery . A retrospective examination of mosquito infection on humans infected with Plasmodium falciparum . Am. J. Trop. Med. Hyg ., 68 ( 3 ): 366 – 371 , March 2003 . ISSN 0002-9637 . OpenUrl Abstract / FREE Full Text 47. ↵ Xi Zen Yap , Matthew B. B. McCall , and Robert W. Sauerwein . Fast and fierce versus slow and smooth: Heterogeneity in immune responses to Plasmodium in the controlled human malaria infection model . Immunol. Rev ., 293 ( 1 ): 253 – 269 , January 2020 . ISSN 0105-2896, 1600-065X . doi: 10.1111/imr.12811 . OpenUrl CrossRef PubMed 48. ↵ Wendy Prudhomme O’Meara , F. Ellis McKenzie , Alan J. Magill , J. Russ Forney , Barnyen Permpanich , Carmen Lucas , Robert A. Gasser , and Chansuda Wongsrichanalai . Sources of variability in determining malaria parasite density by microscopy . Am. J. Trop. Med. Hyg ., 73 ( 3 ): 593 – 598 , September 2005 . ISSN 0002-9637 . OpenUrl Abstract / FREE Full Text 49. ↵ Paul E. M. Fine . Ross’s a priori pathometry - a perspective . Proc. R. Soc. Med ., 68 ( 9 ): 547 – 551 , September 1975 . ISSN 0035-9157 . OpenUrl PubMed View the discussion thread. Back to top Previous Next Posted March 25, 2025. Download PDF Data/Code Email Thank you for your interest in spreading the word about medRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following A Probabilistic Synthesis of Malaria Epidemiology: Exposure, Infection, Parasite Densities, and Detection Message Subject (Your Name) has forwarded a page to you from medRxiv Message Body (Your Name) thought you would like to see this page from the medRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share A Probabilistic Synthesis of Malaria Epidemiology: Exposure, Infection, Parasite Densities, and Detection John M. Henry , Austin R. Carter , Sean L. Wu , David L. Smith medRxiv 2025.03.24.25324561; doi: https://doi.org/10.1101/2025.03.24.25324561 Share This Article: Copy Citation Tools A Probabilistic Synthesis of Malaria Epidemiology: Exposure, Infection, Parasite Densities, and Detection John M. Henry , Austin R. Carter , Sean L. Wu , David L. Smith medRxiv 2025.03.24.25324561; doi: https://doi.org/10.1101/2025.03.24.25324561 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Infectious Diseases (except HIV/AIDS) Subject Areas All Articles Addiction Medicine (568) Allergy and Immunology (863) Anesthesia (300) Cardiovascular Medicine (4436) Dentistry and Oral Medicine (444) Dermatology (382) Emergency Medicine (608) Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1509) Epidemiology (15229) Forensic Medicine (30) Gastroenterology (1124) Genetic and Genomic Medicine (6600) Geriatric Medicine (668) Health Economics (997) Health Informatics (4538) Health Policy (1368) Health Systems and Quality Improvement (1613) Hematology (542) HIV/AIDS (1264) Infectious Diseases (except HIV/AIDS) (15916) Intensive Care and Critical Care Medicine (1103) Medical Education (623) Medical Ethics (146) Nephrology (667) Neurology (6599) Nursing (346) Nutrition (998) Obstetrics and Gynecology (1144) Occupational and Environmental Health (957) Oncology (3333) Ophthalmology (974) Orthopedics (369) Otolaryngology (420) Pain Medicine (436) Palliative Medicine (130) Pathology (663) Pediatrics (1693) Pharmacology and Therapeutics (691) Primary Care Research (711) Psychiatry and Clinical Psychology (5447) Public and Global Health (9232) Radiology and Imaging (2198) Rehabilitation Medicine and Physical Therapy (1370) Respiratory Medicine (1196) Rheumatology (593) Sexual and Reproductive Health (712) Sports Medicine (530) Surgery (712) Toxicology (99) Transplantation (289) Urology (265) (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'a00f96e8ef3b8650',t:'MTc3OTY2MDE5Nw=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();