High Plasmodium infection intensity in naturally infected malaria vectors in Africa

Revisiting the Plasmodium sporozoite inoculum and elucidating the efficiency with which malaria parasites progress through the mosquito

Malaria is initiated when infected anopheline mosquitoes inoculate sporozoites as they probe for blood. It is thought that all infected mosquitoes are equivalent in terms of their infectious potential, with parasite burden having no role in transmission success. In this study, using mosquitoes harboring the entire range of salivary gland sporozoite loads observed in the field, we demonstrate a strong and highly significant correlation between mosquito parasite burden and inoculum size. We then link the inoculum data to oocyst counts, the most commonly-used metric to assess mosquito infection in the field, and determine the efficiency with which oocyst sporozoites enter mosquito salivary glands. Taken together our data support the conclusion that mosquitoes with higher parasite burdens are more likely to initiate infection and contribute to onward transmission. Overall these data may account for some of the unexplained heterogeneity in transmission and enable more precise benchmarks for transmission-blocking interventions.

Dual effector population modification gene-drive strains of the African malaria mosquitoes, Anopheles gambiae and Anopheles coluzzii

Proposed genetic approaches for reducing human malaria include population modification, which introduces genes into vector mosquitoes to reduce or prevent parasite transmission. We demonstrate the potential of Cas9/guide RNA (gRNA)-based gene-drive systems linked to dual antiparasite effector genes to spread rapidly through mosquito populations. Two strains have an autonomous gene-drive system coupled to dual anti-Plasmodium falciparum effector genes comprising single-chain variable fragment monoclonal antibodies targeting parasite ookinetes and sporozoites in the African malaria mosquitoes Anopheles gambiae (AgTP13) and Anopheles coluzzii (AcTP13). The gene-drive systems achieved full introduction within 3 to 6 mo after release in small cage trials. Life-table analyses revealed no fitness loads affecting AcTP13 gene-drive dynamics but AgTP13 males were less competitive than wild types. The effector molecules reduced significantly both parasite prevalence and infection intensities. These data supported transmission modeling of conceptual field releases in an island setting that shows meaningful epidemiological impacts at different sporozoite threshold levels (2.5 to 10 k) for human infection by reducing malaria incidence in optimal simulations by 50 to 90% within as few as 1 to 2 mo after a series of releases, and by ≥90% within 3 mo. Modeling outcomes for low sporozoite thresholds are sensitive to gene-drive system fitness loads, gametocytemia infection intensities during parasite challenges, and the formation of potentially drive-resistant genome target sites, extending the predicted times to achieve reduced incidence. TP13-based strains could be effective for malaria control strategies following validation of sporozoite transmission threshold numbers and testing field-derived parasite strains. These or similar strains are viable candidates for future field trials in a malaria-endemic region.

Intervention reducing malaria parasite load in vector mosquitoes: No impact on Plasmodium falciparum extrinsic incubation period and the survival of Anopheles gambiae

In the fight against malaria, transmission blocking interventions (TBIs) such as transmission blocking vaccines or drugs, are promising approaches to complement conventional tools. They aim to prevent the infection of vectors and thereby reduce the subsequent exposure of a human population to infectious mosquitoes. The effectiveness of these approaches has been shown to depend on the initial intensity of infection in mosquitoes, often measured as the mean number of oocysts resulting from an infectious blood meal in absence of intervention. In mosquitoes exposed to a high intensity of infection, current TBI candidates are expected to be ineffective at completely blocking infection but will decrease parasite load and therefore, potentially also affect key parameters of vector transmission. The present study investigated the consequences of changes in oocyst intensity on subsequent parasite development and mosquito survival. To address this, we experimentally produced different intensities of infection for Anopheles gambiae females from Burkina Faso by diluting gametocytes from three natural Plasmodium falciparum local isolates and used a newly developed non-destructive method based on the exploitation of mosquito sugar feeding to track parasite and mosquito life history traits throughout sporogonic development. Our results indicate the extrinsic incubation period (EIP) of P. falciparum and mosquito survival did not vary with parasite density but differed significantly between parasite isolates with estimated EIP50 of 16 (95% CI: 15-18), 14 (95% CI: 12-16) and 12 (95% CI: 12-13) days and median longevity of 25 (95% CI: 22-29), 15 (95% CI: 13-15) and 18 (95% CI: 17-19) days for the three isolates respectively. Our results here do not identify unintended consequences of the decrease of parasite loads in mosquitoes on the parasite incubation period or on mosquito survival, two key parameters of vectorial capacity, and hence support the use of transmission blocking strategies to control malaria.

Evolutionary consequences of vector-borne transmission: how using vectors shapes host, vector and pathogen evolution

Transmission mode is a key factor that influences host–parasite coevolution. Vector-borne pathogens are among the most important disease agents for humans and wildlife due to their broad distribution, high diversity, prevalence and lethality. They comprise some of the most important and widespread human pathogens, such as yellow fever, leishmania and malaria. Vector-borne parasites (in this review, those transmitted by blood-feeding Diptera) follow unique transmission routes towards their vertebrate hosts. Consequently, each part of this tri-partite (i.e. parasite, vector and host) interaction can influence co- and counter-evolutionary pressures among antagonists. This mode of transmission may favour the evolution of greater virulence to the vertebrate host; however, pathogen–vector interactions can also have a broad spectrum of fitness costs to the insect vector. To complete their life cycle, vector-borne pathogens must overcome immune responses from 2 unrelated organisms, since they can activate responses in both vertebrate and invertebrate hosts, possibly creating a trade-off between investments against both types of immunity. Here, we assess how dipteran vector-borne transmission shapes the evolution of hosts, vectors and the pathogens themselves. Hosts, vectors and pathogens co-evolve together in a constant antagonistic arms race with each participant's primary goal being to maximize its performance and fitness.

Omitting age-dependent mosquito mortality in malaria models underestimates the effectiveness of insecticide-treated nets

Mathematical models of vector-borne infections, including malaria, often assume age-independent mortality rates of vectors, despite evidence that many insects senesce. In this study we present survival data on insecticide-resistant Anopheles gambiae s.l. from experiments in Côte d’Ivoire. We fit a constant mortality function and two age-dependent functions (logistic and Gompertz) to the data from mosquitoes exposed (treated) and not exposed (control) to insecticide-treated nets (ITNs), to establish biologically realistic survival functions. This enables us to explore the effects of insecticide exposure on mosquito mortality rates, and the extent to which insecticide resistance might impact the effectiveness of ITNs. We investigate this by calculating the expected number of infectious bites a mosquito will take in its lifetime, and by extension the vectorial capacity. Our results show that the predicted vectorial capacity is substantially lower in mosquitoes exposed to ITNs, despite the mosquitoes in the experiment being highly insecticide-resistant. The more realistic age-dependent functions provide a better fit to the experimental data compared to a constant mortality function and, hence, influence the predicted impact of ITNs on malaria transmission potential. In models with age-independent mortality, there is a great reduction for the vectorial capacity under exposure compared to no exposure. However, the two age-dependent functions predicted an even larger reduction due to exposure, highlighting the impact of incorporating age in the mortality rates. These results further show that multiple exposures to ITNs had a considerable effect on the vectorial capacity. Overall, the study highlights the importance of including age dependency in mathematical models of vector-borne disease transmission and in fully understanding the impact of interventions.

Interventions against malaria are most commonly targeted on the adult mosquitoes, which transmit the infection from person to person. One of the most important interventions are bed-nets, treated with insecticides. Unfortunately, extensive exposure of mosquitoes to insecticide has led to widespread evolution of insecticide resistance, which might threaten control strategies. Piecing together the overall impact of resistance on the efficacy of insecticide-treated nets is complex, but can be informed by the use of mathematical models. However, there are some assumptions that the models frequently use which are not realistic in terms of the mosquito biology. In this paper, we formulate a model that includes age-dependent mortality rates, an important parameter in vector control since control strategies most commonly aim to reduce the lifespan of the mosquitoes. By using novel data collected using field-derived insecticide-resistant mosquitoes, we explore the effects that the presence of insecticides on nets have on the mortality rates, as well as the difference incorporating age dependency in the model has on the results. We find that including age-dependent mortality greatly alters the anticipated effects of insecticide-treated nets on mosquito transmission potential, and that ignoring this realism potentially overestimates the negative impact of insecticide resistance.

Streamlining sporozoite isolation from mosquitoes by leveraging the dynamics of migration to the salivary glands

Background

Sporozoites isolated from the salivary glands of Plasmodium-infected mosquitoes are a prerequisite for several basic and pre-clinical applications. Although salivary glands are pooled to maximize sporozoite recovery, insufficient yields pose logistical and analytical hurdles; thus, predicting yields prior to isolation would be valuable. Preceding oocyst densities in the midgut is an obvious candidate. However, it is unclear whether current understanding of its relationship with sporozoite densities can be used to maximize yields, or whether it can capture the potential density-dependence in rates of sporozoite invasion of the salivary glands.

Methods

This study presents a retrospective analysis of Anopheles stephensi mosquitoes infected with two strains of the rodent-specific Plasmodium berghei. Mean oocyst densities were estimated in the midguts earlier in the infection (11–15 days post-blood meal), with sporozoites pooled from the salivary glands later in the infection (17–29 days). Generalized linear mixed effects models were used to determine if (1) mean oocyst densities can predict sporozoite yields from pooled salivary glands, (2) whether these densities can capture differences in rates of sporozoite invasion of salivary glands, and (3), if the interaction between oocyst densities and time could be leveraged to boost overall yields.

Results

The non-linear effect of mean oocyst densities confirmed the role of density-dependent constraints in limiting yields beyond certain oocyst densities. Irrespective of oocyst densities however, the continued invasion of salivary glands by the sporozoites boosted recoveries over time (17–29 days post-blood meal) for either parasite strain.

Conclusions

Sporozoite invasion of the salivary glands over time can be leveraged to maximize yields for P. berghei. In general, however, invasion of the salivary glands over time is a critical fitness determinant for all Plasmodium species (extrinsic incubation period, EIP). Thus, delaying sporozoite collection could, in principle, substantially reduce dissection effort for any parasite within the genus, with the results also alluding to the potential for changes in sporozoites densities over time to modify infectivity for the next host.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04270-y.

R H: a genetic metric for measuring intrahost Plasmodium falciparum relatedness and distinguishing cotransmission from superinfection

Multiple-strain (polygenomic) infections are a ubiquitous feature of Plasmodium falciparum parasite population genetics. Under simple assumptions of superinfection, polygenomic infections are hypothesized to be the result of multiple infectious bites. As a result, polygenomic infections have been used as evidence of repeat exposure and used to derive genetic metrics associated with high transmission intensity. However, not all polygenomic infections are the result of multiple infectious bites. Some result from the transmission of multiple, genetically related strains during a single infectious bite (cotransmission). Superinfection and cotransmission represent two distinct transmission processes, and distinguishing between the two could improve inferences regarding parasite transmission intensity. Here, we describe a new metric, R H, that utilizes the correlation in allelic state (heterozygosity) within polygenomic infections to estimate the likelihood that the observed complexity resulted from either superinfection or cotransmission. R H is flexible and can be applied to any type of genetic data. As a proof of concept, we used R H to quantify polygenomic relatedness and estimate cotransmission and superinfection rates from a set of 1,758 malaria infections genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode. Contrary to expectation, we found that cotransmission was responsible for a significant fraction of 43% to 53% of the polygenomic infections collected in three distinct epidemiological regions in Senegal. The prediction that polygenomic infections frequently result from cotransmission stresses the need to incorporate estimates of relatedness within polygenomic infections to ensure the accuracy of genomic epidemiology surveillance data for informing public health activities.

Xeno-monitoring of molecular drivers of artemisinin and partner drug resistance in P. falciparum populations in malaria vectors across Cameroon

•

High Plasmodium infection rate in the major Anopheles vectors across Cameroon.

•

Emerging signal of the R575I polymorphism in the k13 propeller domain backbone.

•

Dominance of the N⁸⁶F¹⁸⁴mdr1 variants in natural P. falciparum populations.

•

Low k13 and mdr1 genetic diversity in P. falciparum-infected mosquitoes.

Background

Monitoring of drug resistance in Plasmodium populations is crucial for malaria control. This has primarily been performed in humans and rarely in mosquitoes where parasites genetic recombination occurs. Here, we characterized the Plasmodium spp populations in wild Anopheles vectors by analyzing the genetic diversity of the P. falciparum kelch13 and mdr1 gene fragments implicated in artemisinin and partner drug resistance across Cameroon in three major malaria vectors.

Methods

Anopheles mosquitoes were collected across nine localities in Cameroon and dissected into the head/thorax (H/T) and abdomen (Abd) after species identification. A TaqMan assay was performed to detect Plasmodium infection. Fragments of the Kelch 13 and mdr1 genes were amplified in P. falciparum positive samples and directly sequenced to assess their drug resistance polymorphisms and genetic diversity profile.

Results

The study revealed a high Plasmodium infection rate in the major Anopheles vectors across Cameroon. Notably, An. funestus vector recorded the highest sporozoite (8.0%) and oocyst (14.4%) infection rates. A high P. falciparum sporozoite rate (80.08%) alongside epidemiological signatures of significant P. malariae (15.9%) circulation were recorded in these vectors. Low genetic diversity with six (A578S, R575I, G450R, L663L, G453D, N458D) and eight (H53H, V62L, V77E, N86Y, G102G, L132I, H143H, Y184F) point mutations were observed in the k13 and mdr1 backbones respectively. Remarkably, the R575I (4.4%) k13 and Y184F (64.2%) mdr1 mutations were the predominant variants in the P. falciparum populations.

Conclusion

The emerging signal of the R575I polymorphism in the Pfk13 propeller backbone entails the regular surveillance of molecular markers to inform evidence-based policy decisions. Moreover, the high frequency of the ⁸⁶N¹⁸⁴F allele highlights concerns on the plausible decline in efficacy of artemisinin-combination therapies (ACTs); further implying that parasite genotyping from mosquitoes can provide a more relevant scale for quantifying resistance epidemiology in the field.

The potential impact of Anopheles stephensi establishment on the transmission of Plasmodium falciparum in Ethiopia and prospective control measures

BACKGROUND

Sub-Saharan Africa has seen substantial reductions in cases and deaths due to malaria over the past two decades. While this reduction is primarily due to an increasing expansion of interventions, urbanisation has played its part as urban areas typically experience substantially less malaria transmission than rural areas. However, this may be partially lost with the invasion and establishment of Anopheles stephensi. A. stephensi, the primary urban malaria vector in Asia, was first detected in Africa in 2012 in Djibouti and was subsequently identified in Ethiopia in 2016, and later in Sudan and Somalia. In Djibouti, malaria cases have increased 30-fold from 2012 to 2019 though the impact in the wider region remains unclear.

METHODS

Here, we have adapted an existing model of mechanistic malaria transmission to estimate the increase in vector density required to explain the trends in malaria cases seen in Djibouti. To account for the observed plasticity in An. stephensi behaviour, and the unknowns of how it will establish in a novel environment, we sample behavioural parameters in order to account for a wide range of uncertainty. This quantification is then applied to Ethiopia, considering temperature-dependent extrinsic incubation periods, pre-existing vector-control interventions and Plasmodium falciparum prevalence in order to assess the potential impact of An. stephensi establishment on P. falciparum transmission. Following this, we estimate the potential impact of scaling up ITN (insecticide-treated nets)/IRS (indoor residual spraying) and implementing piperonyl butoxide (PBO) ITNs and larval source management, as well as their economic costs.

RESULTS

We estimate that annual P. falciparum malaria cases could increase by 50% (95% CI 14-90) if no additional interventions are implemented. The implementation of sufficient control measures to reduce malaria transmission to pre-stephensi levels will cost hundreds of millions of USD.

CONCLUSIONS

Substantial heterogeneity across the country is predicted and large increases in vector control interventions could be needed to prevent a major public health emergency.

A Survey of Potentially Pathogenic-Incriminated Arthropod Vectors of Health Concern in Botswana

Arthropod vectors play a crucial role in the transmission of many debilitating infections, causing significant morbidity and mortality globally. Despite the economic significance of arthropods to public health, public knowledge on vector biology, ecology and taxonomic status remains anecdotal and largely unexplored. The present study surveyed knowledge gaps regarding the biology and ecology of arthropod vectors in communities of Botswana, across all districts. Results showed that communities are largely aware of individual arthropod vectors; however, their 'potential contribution' in disease transmission in humans, livestock and wildlife could not be fully attested. As such, their knowledge was largely limited with regards to some aspects of vector biology, ecology and control. Communities were strongly concerned about the burden of mosquitoes, cockroaches, flies and ticks, with the least concerns about fleas, bedbugs and lice, although the same communities did not know of specific diseases potentially vectored by these arthropods. Knowledge on arthropod vector control was mainly limited to synthetic chemical pesticides for most respondents, regardless of their location. The limited knowledge on potentially pathogen-incriminated arthropod vectors reported here has large implications for bridging knowledge gaps on the bio-ecology of these vectors countrywide. This is potentially useful in reducing the local burden of associated diseases and preventing the risk of emerging and re-emerging infectious diseases under global change.

Detection of Plasmodium falciparum in laboratory-reared and naturally infected wild mosquitoes using near-infrared spectroscopy

There is an urgent need for high throughput, affordable methods of detecting pathogens inside insect vectors to facilitate surveillance. Near-infrared spectroscopy (NIRS) has shown promise to detect arbovirus and malaria in the laboratory but has not been evaluated in field conditions. Here we investigate the ability of NIRS to identify Plasmodium falciparum in Anopheles coluzzii mosquitoes. NIRS models trained on laboratory-reared mosquitoes infected with wild malaria parasites can detect the parasite in comparable mosquitoes with moderate accuracy though fails to detect oocysts or sporozoites in naturally infected field caught mosquitoes. Models trained on field mosquitoes were unable to predict the infection status of other field mosquitoes. Restricting analyses to mosquitoes of uninfectious and highly-infectious status did improve predictions suggesting sensitivity and specificity may be better in mosquitoes with higher numbers of parasites. Detection of infection appears restricted to homogenous groups of mosquitoes diminishing NIRS utility for detecting malaria within mosquitoes.

Immunofocusing humoral immunity potentiates the functional efficacy of the AnAPN1 malaria transmission-blocking vaccine antigen

Malaria transmission-blocking vaccines (TBVs) prevent the completion of the developmental lifecycle of malarial parasites within the mosquito vector, effectively blocking subsequent infections. The mosquito midgut protein Anopheline alanyl aminopeptidase N (AnAPN1) is the leading, mosquito-based TBV antigen. Structure-function studies identified two Class II epitopes that can induce potent transmission-blocking (T-B) antibodies, informing the design of the next-generation AnAPN1. Here, we functionally screened new immunogens and down-selected to the UF6b construct that has two glycine-linked copies of the T-B epitopes. We then established a process for manufacturing UF6b and evaluated in outbred female CD1 mice the immunogenicity of the preclinical product with the human-safe adjuvant Glucopyranosyl Lipid Adjuvant in a liposomal formulation with saponin QS21 (GLA-LSQ). UF6b:GLA-LSQ effectively immunofocused the humoral response to one of the key T-B epitopes resulting in potent T-B activity, underscoring UF6b as a prime TBV candidate to aid in malaria elimination and eradication efforts.

Predicting the public health impact of a malaria transmission-blocking vaccine

Transmission-blocking vaccines that interrupt malaria transmission from humans to mosquitoes are being tested in early clinical trials. The activity of such a vaccine is commonly evaluated using membrane-feeding assays. Understanding the field efficacy of such a vaccine requires knowledge of how heavily infected wild, naturally blood-fed mosquitoes are, as this indicates how difficult it will be to block transmission. Here we use data on naturally infected mosquitoes collected in Burkina Faso to translate the laboratory-estimated activity into an estimated activity in the field. A transmission dynamics model is then utilised to predict a transmission-blocking vaccine's public health impact alongside existing interventions. The model suggests that school-aged children are an attractive population to target for vaccination. Benefits of vaccination are distributed across the population, averting the greatest number of cases in younger children. Utilising a transmission-blocking vaccine alongside existing interventions could have a substantial impact against malaria.

Plasmodium oocysts respond with dormancy to crowding and nutritional stress

Malaria parasites develop as oocysts in the mosquito for several days before they are able to infect a human host. During this time, mosquitoes take bloodmeals to replenish their nutrient and energy reserves needed for flight and reproduction. We hypothesized that these bloodmeals are critical for oocyst growth and that experimental infection protocols, typically involving a single bloodmeal at the time of infection, cause nutritional stress to the developing oocysts. Therefore, enumerating oocysts disregarding their growth and differentiation state may lead to erroneous conclusions about the efficacy of transmission blocking interventions. Here, we examine this hypothesis in Anopheles coluzzii mosquitoes infected with the human and rodent parasites Plasmodium falciparum and Plasmodium berghei, respectively. We show that oocyst growth and maturation rates decrease at late developmental stages as infection intensities increase; an effect exacerbated at very high infection intensities but fully restored with post infection bloodmeals. High infection intensities and starvation conditions reduce RNA Polymerase III activity in oocysts unless supplemental bloodmeals are provided. Our results suggest that oocysts respond to crowding and nutritional stress with a dormancy-like strategy, which urges the development of alternative methods to assess the efficacy of transmission blocking interventions.

Estimating the extrinsic incubation period of malaria using a mechanistic model of sporogony

During sporogony, malaria-causing parasites infect a mosquito, reproduce and migrate to the mosquito salivary glands where they can be transmitted the next time blood feeding occurs. The time required for sporogony, known as the extrinsic incubation period (EIP), is an important determinant of malaria transmission intensity. The EIP is typically estimated as the time for a given percentile, x, of infected mosquitoes to develop salivary gland sporozoites (the infectious parasite life stage), which is denoted by EIPx. Many mechanisms, however, affect the observed sporozoite prevalence including the human-to-mosquito transmission probability and possibly differences in mosquito mortality according to infection status. To account for these various mechanisms, we present a mechanistic mathematical model, which explicitly models key processes at the parasite, mosquito and observational scales. Fitting this model to experimental data, we find greater variation in the EIP than previously thought: we estimated the range between EIP10 and EIP90 (at 27°C) as 4.5 days compared to 0.9 days using existing statistical methods. This pattern holds over the range of study temperatures included in the dataset. Increasing temperature from 21°C to 34°C decreased the EIP50 from 16.1 to 8.8 days. Our work highlights the importance of mechanistic modelling of sporogony to (1) improve estimates of malaria transmission under different environmental conditions or disease control programs and (2) evaluate novel interventions that target the mosquito life stages of the parasite.

Last-come, best served? Mosquito biting order and Plasmodium transmission

A pervasive characteristic of parasite infections is their tendency to be overdispersed. Understanding the mechanisms underlying this overdispersed distribution is of key importance as it may impact the transmission dynamics of the pathogen. Although multiple factors ranging from environmental stochasticity to inter-individual heterogeneity may explain parasite overdispersion, parasite infection is also overdispersed in an inbred host population maintained under laboratory conditions, suggesting that other mechanisms are at play. Here, we show that the aggregated distribution of malaria parasites within mosquito vectors is partially explained by a temporal heterogeneity in parasite infectivity triggered by the bites of mosquitoes. Parasite transmission tripled between the mosquito's first and last blood feed in a period of only 3 h. Surprisingly, the increase in transmission is not associated with an increase in parasite investment in production of the transmissible stage. Overall, we highlight that Plasmodium is capable of responding to the bites of mosquitoes to increase its own transmission at a much faster pace than initially thought and that this is partly responsible for overdispersed distribution of infection. We discuss the underlying mechanisms as well as the broader implications of this plastic response for the epidemiology of malaria.