Quantitative, model-based estimates of variability in the generation and serial intervals of Plasmodium falciparum malaria

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in High-Recombining Plasmodium falciparum Genomes

Genomic surveillance is crucial for identifying at-risk populations for targeted malaria control and elimination. Identity-by-descent (IBD) is increasingly being used in Plasmodium population genomics to estimate genetic relatedness, effective population size (N e ), population structure, and signals of positive selection. Despite its potential, a thorough evaluation of IBD segment detection tools for species with high recombination rates, such as P. falciparum, remains absent. Here, we perform comprehensive benchmarking of IBD callers - probabilistic (hmmIBD, isoRelate), identity-by-state-based (hap-IBD, phased IBD) and others (Refined IBD) - using population genetic simulations tailored for high recombination, and IBD quality metrics at both the IBD segment level and the IBD-based downstream inference level. Our results demonstrate that low marker density per genetic unit, related to high recombination relative to mutation, significantly compromises the accuracy of detected IBD segments. In genomes with high recombination rates resembling P. falciparum, most IBD callers exhibit high false negative rates for shorter IBD segments, which can be partially mitigated through optimization of IBD caller parameters, especially those related to marker density. Notably, IBD detected with optimized parameters allows for more accurate capture of selection signals and population structure; IBD-based N e inference is very sensitive to IBD detection errors, with IBD called from hmmIBD uniquely providing less biased estimates of N e in this context. Validation with empirical data from the MalariaGEN Pf 7 database, representing different transmission settings, corroborates these findings. We conclude that context-specific evaluation and parameter optimization are essential for accurate IBD detection in high-recombining species and recommend hmmIBD for quality-sensitive analysis, such as estimation of N e in these species. Our optimization and high-level benchmarking methods not only improve IBD segment detection in high-recombining genomes but also enhance overall genomic analysis, paving the way for more accurate genomic surveillance and targeted intervention strategies for malaria.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia

BACKGROUND

Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries.

METHODS

We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection.

RESULTS

We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure.

CONCLUSIONS

Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

No evidence of sustained nonzoonotic Plasmodium knowlesi transmission in Malaysia from modelling malaria case data

Reported incidence of the zoonotic malaria Plasmodium knowlesi has markedly increased across Southeast Asia and threatens malaria elimination. Nonzoonotic transmission of P. knowlesi has been experimentally demonstrated, but it remains unknown whether nonzoonotic transmission is contributing to increases in P. knowlesi cases. Here, we adapt model-based inference methods to estimate RC, individual case reproductive numbers, for P. knowlesi, P. falciparum and P. vivax human cases in Malaysia from 2012-2020 (n = 32,635). Best fitting models for P. knowlesi showed subcritical transmission (RC < 1) consistent with a large reservoir of unobserved infection sources, indicating P. knowlesi remains a primarily zoonotic infection. In contrast, sustained transmission (RC > 1) was estimated historically for P. falciparum and P. vivax, with declines in RC estimates observed over time consistent with local elimination. Together, this suggests sustained nonzoonotic P. knowlesi transmission is highly unlikely and that new approaches are urgently needed to control spillover risks.

Measurably recombining malaria parasites

Genomic epidemiology has guided research and policy for various viral pathogens and there has been a parallel effort towards using genomic epidemiology to combat diseases that are caused by eukaryotic pathogens, such as the malaria parasite. However, the central concept of viral genomic epidemiology, namely that of measurably mutating pathogens, does not apply easily to sexually recombining parasites. Here we introduce the related but different concept of measurably recombining malaria parasites to promote convergence around a unifying theoretical framework for malaria genomic epidemiology. Akin to viral phylodynamics, we anticipate that an inferential framework developed around recombination will help guide practical research and thus realize the full public health potential of genomic epidemiology for malaria parasites and other sexually recombining pathogens.

Life table analysis of Anopheles balabacensis, the primary vector of Plasmodium knowlesi in Sabah, Malaysia

BACKGROUND

Plasmodium knowlesi has become a major public health concern in Sabah, Malaysian Borneo, where it is now the only cause of indigenous malaria. The importance of P. knowlesi has spurred on a series of studies on this parasite, as well as on the biology and ecology of its principal vector, Anopheles balabacensis. However, there remain critical knowledge gaps on the biology of An. balabacensis, such as life history data and life table parameters. To fill these gaps, we conducted a life table study of An. balabacensis in the laboratory. Characterising vector life cycles and survival rates can inform more accurate estimations of the serial interval, the time between two linked cases, which is crucial to understanding and monitoring potentially changing transmission patterns.

METHODS

Individuals of An. balabacensis were collected in the field in Ranau district, Sabah to establish a laboratory colony. Induced mating was used, and the life history parameters of the progeny were recorded. The age-stage, two-sex life table approach was used in the analysis. The culture conditions in the laboratory were 9 h light:15 h dark, mean temperature 25.7 °C ± 0.05 and relative humidity 75.8% ± 0.31.

RESULTS

The eggs hatched within 2 days, and the larval stage lasted for 10.5 days in total, with duration of instar stages I, II, III and IV of 2.3, 3.7, 2.3, 2.2 days, respectively. The maximum total fecundity was 729 for one particular female, while the maximum female age-specific mean fecundity (m_x) was 142 at age 59 days. The gross reproductive rate or number of offspring per individual was about 102. On average, each female laid 1.81 ± 0.19 (range 1-7) batches of eggs, with 63% of the females producing only one batch; only one female laid six batches, while one other laid seven. Each batch comprised 159 ± 17.1 eggs (range 5-224) and the female ratio of offspring was 0.28 ± 0.06. The intrinsic rate of increase, finite rate of increase, net reproductive rate, mean generation time and doubling time were, respectively, 0.12 ± 0.01 day^-1, 1.12 ± 0.01 day^-1, 46.2 ± 14.97, 33.02 ± 1.85 and 5.97 days.

CONCLUSIONS

Both the net reproductive rate and intrinsic rate of increase of An. balabacensis are lower than those of other species in published studies. Our results can be used to improve models of P. knowlesi transmission and to set a baseline for assessing the impacts of environmental change on malaria dynamics. Furthermore, incorporating these population parameters of An. balabacensis into spatial and temporal models on the transmission of P. knowlesi would provide better insight and increase the accuracy of epidemiological forecasting.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

The importance of the generation interval in investigating dynamics and control of new SARS-CoV-2 variants

Inferring the relative strength (i.e. the ratio of reproduction numbers) and relative speed (i.e. the difference between growth rates) of new SARS-CoV-2 variants is critical to predicting and controlling the course of the current pandemic. Analyses of new variants have primarily focused on characterizing changes in the proportion of new variants, implicitly or explicitly assuming that the relative speed remains fixed over the course of an invasion. We use a generation-interval-based framework to challenge this assumption and illustrate how relative strength and speed change over time under two idealized interventions: a constant-strength intervention like idealized vaccination or social distancing, which reduces transmission rates by a constant proportion, and a constant-speed intervention like idealized contact tracing, which isolates infected individuals at a constant rate. In general, constant-strength interventions change the relative speed of a new variant, while constant-speed interventions change its relative strength. Differences in the generation-interval distributions between variants can exaggerate these changes and modify the effectiveness of interventions. Finally, neglecting differences in generation-interval distributions can bias estimates of relative strength.

Receptivity to malaria: meaning and measurement

"Receptivity" to malaria is a construct developed during the Global Malaria Eradication Programme (GMEP) era. It has been defined in varied ways and no consistent, quantitative definition has emerged over the intervening decades. Despite the lack of consistency in defining this construct, the idea that some areas are more likely to sustain malaria transmission than others has remained important in decision-making in malaria control, planning for malaria elimination and guiding activities during the prevention of re-establishment (POR) period. This manuscript examines current advances in methods of measurement. In the context of a decades long decline in global malaria transmission and an increasing number of countries seeking to eliminate malaria, understanding and measuring malaria receptivity has acquired new relevance.

Inferring person-to-person networks of Plasmodium falciparum transmission: are analyses of routine surveillance data up to the task?

BACKGROUND

Inference of person-to-person transmission networks using surveillance data is increasingly used to estimate spatiotemporal patterns of pathogen transmission. Several data types can be used to inform transmission network inferences, yet the sensitivity of those inferences to different data types is not routinely evaluated.

METHODS

The influence of different combinations of spatial, temporal, and travel-history data on transmission network inferences for Plasmodium falciparum malaria were evaluated.

RESULTS

The information content of these data types may be limited for inferring person-to-person transmission networks and may lead to an overestimate of transmission. Only when outbreaks were temporally focal or travel histories were accurate was the algorithm able to accurately estimate the reproduction number under control, R_c. Applying this approach to data from Eswatini indicated that inferences of R_c and spatiotemporal patterns therein depend upon the choice of data types and assumptions about travel-history data.

CONCLUSIONS

These results suggest that transmission network inferences made with routine malaria surveillance data should be interpreted with caution.

Plasmodium falciparum importation does not sustain malaria transmission in a semi-arid region of Kenya

Human movement impacts the spread and transmission of infectious diseases. Recently, a large reservoir of Plasmodium falciparum malaria was identified in a semi-arid region of northwestern Kenya historically considered unsuitable for malaria transmission. Understanding the sources and patterns of transmission attributable to human movement would aid in designing and targeting interventions to decrease the unexpectedly high malaria burden in the region. Toward this goal, polymorphic parasite genes (ama1, csp) in residents and passengers traveling to Central Turkana were genotyped by amplicon deep sequencing. Genotyping and epidemiological data were combined to assess parasite importation. The contribution of travel to malaria transmission was estimated by modelling case reproductive numbers inclusive and exclusive of travelers. P. falciparum was detected in 6.7% (127/1891) of inbound passengers, including new haplotypes which were later detected in locally-transmitted infections. Case reproductive numbers approximated 1 and did not change when travelers were removed from transmission networks, suggesting that transmission is not fueled by travel to the region but locally endemic. Thus, malaria is not only prevalent in Central Turkana but also sustained by local transmission. As such, interrupting importation is unlikely to be an effective malaria control strategy on its own, but targeting interventions locally has the potential to drive down transmission.

Inference of malaria reproduction numbers in three elimination settings by combining temporal data and distance metrics

Individual-level geographic information about malaria cases, such as the GPS coordinates of residence or health facility, is often collected as part of surveillance in near-elimination settings, but could be more effectively utilised to infer transmission dynamics, in conjunction with additional information such as symptom onset time and genetic distance. However, in the absence of data about the flow of parasites between populations, the spatial scale of malaria transmission is often not clear. As a result, it is important to understand the impact of varying assumptions about the spatial scale of transmission on key metrics of malaria transmission, such as reproduction numbers. We developed a method which allows the flexible integration of distance metrics (such as Euclidian distance, genetic distance or accessibility matrices) with temporal information into a single inference framework to infer malaria reproduction numbers. Twelve scenarios were defined, representing different assumptions about the likelihood of transmission occurring over different geographic distances and likelihood of missing infections (as well as high and low amounts of uncertainty in this estimate). These scenarios were applied to four individual level datasets from malaria eliminating contexts to estimate individual reproduction numbers and how they varied over space and time. Model comparison suggested that including spatial information improved models as measured by second order AIC (ΔAICc), compared to time only results. Across scenarios and across datasets, including spatial information tended to increase the seasonality of temporal patterns in reproduction numbers and reduced noise in the temporal distribution of reproduction numbers. The best performing parameterisations assumed long-range transmission (> 200 km) was possible. Our approach is flexible and provides the potential to incorporate other sources of information which can be converted into distance or adjacency matrices such as travel times or molecular markers.

Malaria Elimination in Costa Rica: Changes in Treatment and Mass Drug Administration

Costa Rica is a candidate to eliminate malaria by 2020. The remaining malaria transmission hotspots are located within the Huétar Norte Region (HNR), where 90% of the country's 147 malaria cases have occurred since 2016, following a 33-month period without transmission. Here, we examine changes in transmission with the implementation of a supervised seven-day chloroquine and primaquine treatment (7DCPT). We also evaluate the impact of a focal mass drug administration (MDA) in January 2019 at Boca Arenal, the town in HNR reporting the greatest local transmission. We found that the change to a seven-day treatment protocol, from the prior five-day program, was associated with a 98% reduction in malaria transmission. The MDA helped to reduce transmission, keeping the basic reproduction number, R_T, significantly below 1, for at least four months. However, following new imported cases from Nicaragua, autochthonous transmission resumed. Our results highlight the importance of appropriate treatment delivery to reduce malaria transmission, and the challenge that highly mobile populations, if their malaria is not treated, pose to regional elimination efforts in Mesoamerica and México.

Tracking progress towards malaria elimination in China: Individual-level estimates of transmission and its spatiotemporal variation using a diffusion network approach

In order to monitor progress towards malaria elimination, it is crucial to be able to measure changes in spatio-temporal transmission. However, common metrics of malaria transmission such as parasite prevalence are under powered in elimination contexts. China has achieved major reductions in malaria incidence and is on track to eliminate, having reporting zero locally-acquired malaria cases in 2017 and 2018. Understanding the spatio-temporal pattern underlying this decline, especially the relationship between locally-acquired and imported cases, can inform efforts to maintain elimination and prevent re-emergence. This is particularly pertinent in Yunnan province, where the potential for local transmission is highest. Using a geo-located individual-level dataset of cases recorded in Yunnan province between 2011 and 2016, we introduce a novel Bayesian framework to model a latent diffusion process and estimate the joint likelihood of transmission between cases and the number of cases with unobserved sources of infection. This is used to estimate the case reproduction number, Rc. We use these estimates within spatio-temporal geostatistical models to map how transmission varied over time and space, estimate the timeline to elimination and the risk of resurgence. We estimate the mean Rc between 2011 and 2016 to be 0.171 (95% CI = 0.165, 0.178) for P. vivax cases and 0.089 (95% CI = 0.076, 0.103) for P. falciparum cases. From 2014 onwards, no cases were estimated to have a Rc value above one. An unobserved source of infection was estimated to be moderately likely (p>0.5) for 19/ 611 cases and high (p>0.8) for 2 cases, suggesting very high levels of case ascertainment. Our estimates suggest that, maintaining current intervention efforts, Yunnan is unlikely to experience sustained local transmission up to 2020. However, even with a mean of 0.005 projected up to 2020, locally-acquired cases are possible due to high levels of importation.

Falciparum malaria from coastal Tanzania and Zanzibar remains highly connected despite effective control efforts on the archipelago

BACKGROUND

Tanzania's Zanzibar archipelago has made significant gains in malaria control over the last decade and is a target for malaria elimination. Despite consistent implementation of effective tools since 2002, elimination has not been achieved. Importation of parasites from outside of the archipelago is thought to be an important cause of malaria's persistence, but this paradigm has not been studied using modern genetic tools.

METHODS

Whole-genome sequencing (WGS) was used to investigate the impact of importation, employing population genetic analyses of Plasmodium falciparum isolates from both the archipelago and mainland Tanzania. Ancestry, levels of genetic diversity and differentiation, patterns of relatedness, and patterns of selection between these two populations were assessed by leveraging recent advances in deconvolution of genomes from polyclonal malaria infections.

RESULTS

Significant decreases in the effective population sizes were inferred in both populations that coincide with a period of decreasing malaria transmission in Tanzania. Identity by descent analysis showed that parasites in the two populations shared long segments of their genomes, on the order of 5 cM, suggesting shared ancestry within the last 10 generations. Even with limited sampling, two of isolates between the mainland and Zanzibar were identified that are related at the expected level of half-siblings, consistent with recent importation.

CONCLUSIONS

These findings suggest that importation plays an important role for malaria incidence on Zanzibar and demonstrate the value of genomic approaches for identifying corridors of parasite movement to the island.

Thermal biology of mosquito-borne disease

Mosquito-borne diseases cause a major burden of disease worldwide. The vital rates of these ectothermic vectors and parasites respond strongly and nonlinearly to temperature and therefore to climate change. Here, we review how trait-based approaches can synthesise and mechanistically predict the temperature dependence of transmission across vectors, pathogens, and environments. We present 11 pathogens transmitted by 15 different mosquito species - including globally important diseases like malaria, dengue, and Zika - synthesised from previously published studies. Transmission varied strongly and unimodally with temperature, peaking at 23-29ºC and declining to zero below 9-23ºC and above 32-38ºC. Different traits restricted transmission at low versus high temperatures, and temperature effects on transmission varied by both mosquito and parasite species. Temperate pathogens exhibit broader thermal ranges and cooler thermal minima and optima than tropical pathogens. Among tropical pathogens, malaria and Ross River virus had lower thermal optima (25-26ºC) while dengue and Zika viruses had the highest (29ºC) thermal optima. We expect warming to increase transmission below thermal optima but decrease transmission above optima. Key directions for future work include linking mechanistic models to field transmission, combining temperature effects with control measures, incorporating trait variation and temperature variation, and investigating climate adaptation and migration.

A practical generation-interval-based approach to inferring the strength of epidemics from their speed

Infectious disease outbreaks are often characterized by the reproduction number R and exponential rate of growth r. R provides information about outbreak control and predicted final size, but estimating R is difficult, while r can often be estimated directly from incidence data. These quantities are linked by the generation interval - the time between when an individual is infected by an infector, and when that infector was infected. It is often infeasible to obtain the exact shape of a generation-interval distribution, and to understand how this shape affects estimates of R. We show that estimating generation interval mean and variance provides insight into the relationship between R and r. We use examples based on Ebola, rabies and measles to explore approximations based on gamma-distributed generation intervals, and find that use of these simple approximations are often sufficient to capture the r-R relationship and provide robust estimates of R.

Robustness of the reproductive number estimates in vector-borne disease systems

BACKGROUND

The required efforts, feasibility and predicted success of an intervention strategy against an infectious disease are partially determined by its basic reproduction number, R0. In its simplest form R0 can be understood as the product of the infectious period, the number of infectious contacts and the per-contact transmission probability, which in the case of vector-transmitted diseases necessarily extend to the vector stages. As vectors do not usually recover from infection, they remain infectious for life, which places high significance on the vector's life expectancy. Current methods for estimating the R0 for a vector-borne disease are mostly derived from compartmental modelling frameworks assuming constant vector mortality rates. We hypothesised that some of the assumptions underlying these models can lead to unrealistic high vector life expectancies with important repercussions for R0 estimates.

METHODOLOGY AND PRINCIPAL FINDINGS

Here we used a stochastic, individual-based model which allowed us to directly measure the number of secondary infections arising from one index case under different assumptions about vector mortality. Our results confirm that formulas based on age-independent mortality rates can overestimate R0 by nearly 100% compared to our own estimate derived from first principles. We further provide a correction factor that can be used with a standard R0 formula and adjusts for the discrepancies due to erroneous vector age distributions.

CONCLUSION

Vector mortality rates play a crucial role for the success and general epidemiology of vector-transmitted diseases. Many modelling efforts intrinsically assume these to be age-independent, which, as clearly demonstrated here, can lead to severe over-estimation of the disease's reproduction number. Our results thus re-emphasise the importance of obtaining field-relevant and species-dependent vector mortality rates, which in turn would facilitate more realistic intervention impact predictions.

Emerging implications of policies on malaria treatment: genetic changes in the Pfmdr-1 gene affecting susceptibility to artemether-lumefantrine and artesunate-amodiaquine in Africa

Artemether–lumefantrine (AL) and artesunate–amodiaquine (AS-AQ) are the most commonly used artemisinin-based combination therapies (ACT) for treatment of Plasmodium falciparum in Africa. Both treatments remain efficacious, but single nucleotide polymorphisms (SNPs) in the Plasmodium falciparum multidrug resistance 1 (Pfmdr1) gene may compromise sensitivity. AL and AS-AQ exert opposing selective pressures: parasites with genotype 86Y, Y184 and 1246Y are partially resistant to AS-AQ treatment, while N86, 184 F and D1246 are favoured by AL treatment. Through a systematic review, we identified 397 surveys measuring the prevalence of Pfmdr1 polymorphisms at positions 86 184 or 1246 in 30 countries in Africa. Temporal trends in SNP frequencies after introduction of AL or AS-AQ as first-line treatment were analysed in 32 locations, and selection coefficients estimated. We examined associations between antimalarial policies, consumption, transmission intensity and rate of SNP selection. 1246Y frequency decreased on average more rapidly in locations where national policy recommended AL (median selection coefficient(s) of −0.083), compared with policies of AS-AQ or both AL and AS-AQ (median s=−0.035 and 0.021, p<0.001 respectively). 86Y frequency declined markedly after ACT policy introduction, with a borderline significant trend for a more rapid decline in countries with AL policies (p=0.055). However, these trends could also be explained by a difference in initial SNP frequencies at the time of ACT introduction. There were non-significant trends for faster selection of N86 and D1246 in areas with higher AL consumption and no trend with transmission intensity. Recorded consumption of AS-AQ was low in the locations and times Pfmdr1 data were collected. SNP trends in countries with AL policies suggest a broad increase in sensitivity of parasites to AS-AQ, by 7–10 years after AL introduction. Observed rates of selection have implications for planning strategies to cycle drugs or use multiple first-line therapies to maintain drug efficacy.

Mathematical analysis to prioritise strategies for malaria elimination

Malaria and some other tropical diseases are currently targeted for elimination and eventually eradication. Since resources are limited, prioritisation of countries or areas for elimination is often necessary. However, this prioritisation is frequently conducted in an ad hoc manner. Lower transmission areas are usually targeted for elimination first, but for some areas this necessitates long and potentially expensive surveillance programs while transmission is eliminated from neighbouring higher transmission areas. We use a mathematical model to compare the implications of prioritisation choices in reducing overall burden and costs. We show that when the duration of the elimination program is independent of the transmission potential, burden is always reduced most by targeting high transmission areas first, but to reduce costs the optimal ordering depends on the actual transmission levels. In general, when overall transmission potential is low and the surveillance cost per secondary case compared to the cost per imported case is low, targeting the higher transmission area for elimination first is favoured.

Estimating spatiotemporally varying malaria reproduction numbers in a near elimination setting

In 2016 the World Health Organization identified 21 countries that could eliminate malaria by 2020. Monitoring progress towards this goal requires tracking ongoing transmission. Here we develop methods that estimate individual reproduction numbers and their variation through time and space. Individual reproduction numbers, R_c, describe the state of transmission at a point in time and differ from mean reproduction numbers, which are averages of the number of people infected by a typical case. We assess elimination progress in El Salvador using data for confirmed cases of malaria from 2010 to 2016. Our results demonstrate that whilst the average number of secondary malaria cases was below one (0.61, 95% CI 0.55–0.65), individual reproduction numbers often exceeded one. We estimate a decline in R_c between 2010 and 2016. However we also show that if importation is maintained at the same rate, the country may not achieve malaria elimination by 2020.

Twenty one countries have been identified for malaria elimination by 2020 and their progress needs to be constantly evaluated. Here, the authors present a method that estimates individual reproduction numbers and their variation through time and space and use it to monitor elimination success in El Salvador between 2010 and 2016.

Estimating the effective reproduction number of dengue considering temperature-dependent generation intervals

The effective reproduction number, R_t, is a measure of transmission that can be calculated from standard incidence data to timely detect the beginning of epidemics. It has being increasingly used for surveillance of directly transmitted diseases. However, current methods for R_t estimation do not apply for vector borne diseases, whose transmission cycle depends on temperature. Here we propose a method that provides dengue's R_t estimates in the presence of temperature-mediated seasonality and apply this method to simulated and real data from two cities in Brazil where dengue is endemic. The method shows good precision in the simulated data. When applied to the real data, it shows differences in the transmission profile of the two cities and identifies periods of higher transmission.

Mapping multiple components of malaria risk for improved targeting of elimination interventions

There is a long history of considering the constituent components of malaria risk and the malaria transmission cycle via the use of mathematical models, yet strategic planning in endemic countries tends not to take full advantage of available disease intelligence to tailor interventions. National malaria programmes typically make operational decisions about where to implement vector control and surveillance activities based upon simple categorizations of annual parasite incidence. With technological advances, an enormous opportunity exists to better target specific malaria interventions to the places where they will have greatest impact by mapping and evaluating metrics related to a variety of risk components, each of which describes a different facet of the transmission cycle. Here, these components and their implications for operational decision-making are reviewed. For each component, related mappable malaria metrics are also described which may be measured and evaluated by malaria programmes seeking to better understand the determinants of malaria risk. Implementing tailored programmes based on knowledge of the heterogeneous distribution of the drivers of malaria transmission rather than only consideration of traditional metrics such as case incidence has the potential to result in substantial improvements in decision-making. As programmes improve their ability to prioritize their available tools to the places where evidence suggests they will be most effective, elimination aspirations may become increasingly feasible.

Temperature modulates dengue virus epidemic growth rates through its effects on reproduction numbers and generation intervals

Epidemic growth rate, r, provides a more complete description of the potential for epidemics than the more commonly studied basic reproduction number, R₀, yet the former has never been described as a function of temperature for dengue virus or other pathogens with temperature-sensitive transmission. The need to understand the drivers of epidemics of these pathogens is acute, with arthropod-borne virus epidemics becoming increasingly problematic. We addressed this need by developing temperature-dependent descriptions of the two components of r—R₀ and the generation interval—to obtain a temperature-dependent description of r. Our results show that the generation interval is highly sensitive to temperature, decreasing twofold between 25 and 35°C and suggesting that dengue virus epidemics may accelerate as temperatures increase, not only because of more infections per generation but also because of faster generations. Under the empirical temperature relationships that we considered, we found that r peaked at a temperature threshold that was robust to uncertainty in model parameters that do not depend on temperature. Although the precise value of this temperature threshold could be refined following future studies of empirical temperature relationships, the framework we present for identifying such temperature thresholds offers a new way to classify regions in which dengue virus epidemic intensity could either increase or decrease under future climate change.

Recurrent, rapidly intensifying epidemics of dengue–the world’s most prevalent mosquito-borne viral disease–pose a challenge to healthcare systems throughout the tropical and subtropical world. An acute disease that tends to respond well to proper treatment, the sometimes intense nature of dengue epidemics has been known to overwhelm healthcare systems and elevate the morbidity and mortality of patients left without adequate medical treatment under peak epidemic conditions. Here, we quantify the temperature dependence of dengue epidemic intensity by quantifying two distinct determinants of epidemic growth rate: the average number of secondary infections arising from each primary infection and the average time between successive infections in humans. Our results show that the time between successive infections in humans decreases steadily with increasing temperatures, whereas the average number of secondary infections peaks at intermediate temperatures. Altogether, this suggests a peak temperature for dengue epidemic intensity. Applying this result to global temperature projections under future climate change scenarios suggests that dengue epidemics in many regions of the world could become more intense under future temperature increases.