Polymorphisms in Anopheles gambiae immune genes associated with natural resistance to Plasmodium falciparum

The Molecular Toll Pathway Repertoire in Anopheline Mosquitoes

Innate immunity in mosquitoes has received much attention due to its potential impact on vector competence for vector-borne disease pathogens, including malaria parasites. The nuclear factor (NF)-κB-dependent Toll pathway is a major regulator of innate immunity in insects. In mosquitoes, this pathway controls transcription of the majority of the known canonical humoral immune effectors, mediates anti-bacterial, anti-fungal and anti-viral immune responses, and contributes to malaria parasite killing. However, besides initial gene annotation of putative Toll pathway members and genetic analysis of the contribution of few key components to immunity, the molecular make-up and function of the Toll pathway in mosquitoes is largely unexplored. To facilitate functional analyses of the Toll pathway in mosquitoes, we report here manually annotated and refined gene models of Toll-like receptors and all putative components of the intracellular signal transduction cascade across 19 anopheline genomes, and in two culicine genomes. In addition, based on phylogenetic analyses, we identified differing levels of evolutionary constraint across the intracellular Toll pathway members, and identified a recent radiation of TOLL1/5 within the An. gambiae complex. Together, this study provides insight into the evolution of TLRs and the putative members of the intracellular signal transduction cascade within the genus Anopheles.

Expanding the CRISPR base editing toolbox in Drosophila melanogaster

CRISPR base editors can introduce point mutations into DNA precisely, and cytosine base editors (CBEs) catalyze C to T transitions. While CBEs have been thoroughly explored in cell culture and organisms such as mice, little is known about DNA base editing in insects. In this study, we evaluated germline editing rates of three different CBEs expressed under actin (ubiquitous) or nanos (germline) promoters utilizing Drosophila melanogaster. The original Rattus norvegicus-derived cytosine deaminase APOBEC1 (rAPO-1) displayed high base editing rates (~99%) with undetectable indel formation. Additionally, we show that base editors can be used for generating male sterility and female lethality. Overall, this study highlights the importance of promoter choice and sex-specific transmission for efficient base editing in flies while providing new insights for future genetic biocontrol designs in insects.

Homing gene drives can transfer rapidly between Anopheles gambiae strains with minimal carryover of flanking sequences

CRISPR-Cas9 homing gene drives are designed to induce a targeted double-stranded DNA break at a wild type allele ('recipient'), which, when repaired by the host cell, is converted to the drive allele from the homologous ('donor') chromosome. Germline localisation of this process leads to super-Mendelian inheritance of the drive and the rapid spread of linked traits, offering a novel strategy for population control through the deliberate release of drive individuals. During the homology-based DNA repair, additional segments of the recipient chromosome may convert to match the donor, potentially impacting carrier fitness and strategy success. Using Anopheles gambiae strains with variations around the drive target site, here we assess the extent and nature of chromosomal conversion. We show both homing and meiotic drive contribute as mechanisms of inheritance bias. Additionally, over 80% of homing events resolve within 50 bp of the chromosomal break, enabling rapid gene drive transfer into locally-adapted genetic backgrounds.

Revisiting the determinants of malaria transmission

Malaria parasites have coevolved with humans over thousands of years, mirroring their migration out of Africa. They persist to this day, despite continuous elimination efforts worldwide. These parasites can adapt to changing environments during infection of human and mosquito, and when expanding the geographical range by switching vector species. Recent studies in the human malaria parasite, Plasmodium falciparum, identified determinants governing the plasticity of sexual conversion rates, sex ratio, and vector competence. Here we summarize the latest literature revealing environmental, epigenetic, and genetic determinants of malaria transmission.

The Anopheles leucine-rich repeat protein APL1C is a pathogen binding factor recognizing Plasmodium ookinetes and sporozoites

Leucine-rich repeat (LRR) proteins are commonly involved in innate immunity of animals and plants, including for pattern recognition of pathogen-derived elicitors. The Anopheles secreted LRR proteins APL1C and LRIM1 are required for malaria ookinete killing in conjunction with the complement-like TEP1 protein. However, the mechanism of parasite immune recognition by the mosquito remains unclear, although it is known that TEP1 lacks inherent binding specificity. Here, we find that APL1C and LRIM1 bind specifically to Plasmodium berghei ookinetes, even after depletion of TEP1 transcript and protein, consistent with a role for the LRR proteins in pathogen recognition. Moreover, APL1C does not bind to ookinetes of the human malaria parasite Plasmodium falciparum, and is not required for killing of this parasite, which correlates LRR binding specificity and immune protection. Most of the live P. berghei ookinetes that migrated into the extracellular space exposed to mosquito hemolymph, and almost all dead ookinetes, are bound by APL1C, thus associating LRR protein binding with parasite killing. We also find that APL1C binds to the surface of P. berghei sporozoites released from oocysts into the mosquito hemocoel and forms a potent barrier limiting salivary gland invasion and mosquito infectivity. Pathogen binding by APL1C provides the first functional explanation for the long-known requirement of APL1C for P. berghei ookinete killing in the mosquito midgut. We propose that secreted mosquito LRR proteins are required for pathogen discrimination and orientation of immune effector activity, potentially as functional counterparts of the immunoglobulin-based receptors used by vertebrates for antigen recognition.

Evaluation of diffuse reflectance spectroscopy for predicting age, species, and cuticular resistance of Anopheles gambiae s.l under laboratory conditions

Mid-infrared spectroscopy (MIRS) combined with machine learning analysis has shown potential for quick and efficient identification of mosquito species and age groups. However, current technology to collect spectra is destructive to the sample and does not allow targeting specific tissues of the mosquito, limiting the identification of other important biological traits such as insecticide resistance. Here, we assessed the use of a non-destructive approach of MIRS for vector surveillance, micro diffuse reflectance spectroscopy (µDRIFT) using mosquito legs to identify species, age and cuticular insecticide resistance within the Anopheles gambiae s.l. complex. These mosquitoes are the major vectors of malaria in Africa and the focus on surveillance in malaria control programs. Legs required significantly less scanning time and showed more spectral consistence compared to other mosquito tissues. Machine learning models were able to identify An. gambiae and An. coluzzii with an accuracy of 0.73, two ages groups (3 and 10 days old) with 0.77 accuracy and we obtained accuracy of 0.75 when identifying cuticular insecticide resistance. Our results highlight the potential of different mosquito tissues and µDRIFT as tools for biological trait identification on mosquitoes that transmit malaria. These results can guide new ways of identifying mosquito traits which can help the creation of innovative surveillance programs by adapting new technology into mosquito surveillance and control tools.

Measuring the Impact of Genetic Heterogeneity and Chromosomal Inversions on the Efficacy of CRISPR-Cas9 Gene Drives in Different Strains of Anopheles gambiae

The human malaria vector Anopheles gambiae is becoming increasingly resistant to insecticides, spurring the development of genetic control strategies. CRISPR-Cas9 gene drives can modify a population by creating double-stranded breaks at highly specific targets, triggering copying of the gene drive into the cut site ("homing"), ensuring its inheritance. The DNA repair mechanism responsible requires homology between the donor and recipient chromosomes, presenting challenges for the invasion of laboratory-developed gene drives into wild populations of target species An. gambiae species complex, which show high levels of genome variation. Two gene drives (vas2-5958 and zpg-7280) were introduced into three An. gambiae strains collected across Africa with 5.3-6.6% variation around the target sites, and the effect of this variation on homing was measured. Gene drive homing across different karyotypes of the 2La chromosomal inversion was also assessed. No decrease in gene drive homing was seen despite target site heterology, demonstrating the applicability of gene drives to wild populations.

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.

Impacts of dual active-ingredient bed nets on the behavioural responses of pyrethroid resistant Anopheles gambiae determined by room-scale infrared video tracking

BACKGROUND

The success of insecticide treated bed nets (ITNs) for malaria vector control in Africa relies on the behaviour of various species of Anopheles. Previous research has described mosquito behavioural alterations resulting from widespread ITN coverage, which could result in a decrease in net efficacy. Here, behaviours were compared including timings of net contact, willingness to refeed and longevity post-exposure to two next-generation nets, PermaNet^® 3.0 (P3 net) and Interceptor^® G2 (IG2 net) in comparison with a standard pyrethroid-only net (Olyset Net™ (OL net)) and an untreated net.

METHODS

Susceptible and resistant Anopheles gambiae mosquitoes were exposed to the nets with a human volunteer host in a room-scale assay. Mosquito movements were tracked for 2 h using an infrared video system, collecting flight trajectory, spatial position and net contact data. Post-assay, mosquitoes were monitored for a range of sublethal insecticide effects.

RESULTS

Mosquito net contact was focused predominantly on the roof for all four bed nets. A steep decay in activity was observed for both susceptible strains when P3 net and OL net were present and with IG2 net for one of the two susceptible strains. Total mosquito activity was higher around untreated nets than ITNs. There was no difference in total activity, the number, or duration, of net contact, between any mosquito strain, with similar behaviours recorded in susceptible and resistant strains at all ITNs. OL net, P3 net and IG2 net all killed over 90% of susceptible mosquitoes 24 h after exposure, but this effect was not seen with resistant mosquitoes where mortality ranged from 16 to 72%. All treated nets reduced the willingness of resistant strains to re-feed when offered blood 1-h post-exposure, with a more pronounced effect seen with P3 net and OL net than IG2 net.

CONCLUSION

These are the first results to provide an in-depth description of the behaviour of susceptible and resistant Anopheles gambiae strains around next-generation bed nets using a room-scale tracking system to capture multiple behaviours. These results indicate that there is no major difference in behavioural responses between mosquito strains of differing pyrethroid susceptibility when exposed to these new ITNs under the experimental conditions used.

Molecular and functional properties of human Plasmodium falciparum CSP C-terminus antibodies

Human monoclonal antibodies (mAbs) against the central repeat and junction domain of Plasmodium falciparum circumsporozoite protein (PfCSP) have been studied extensively to guide malaria vaccine design compared to antibodies against the PfCSP C terminus. Here, we describe the molecular characteristics and protective potential of 73 germline and mutated human mAbs against the highly immunogenic PfCSP C-terminal domain. Two mAbs recognized linear epitopes in the C-terminal linker with sequence similarity to repeat and junction motifs, whereas all others targeted conformational epitopes in the α-thrombospondin repeat (α-TSR) domain. Specificity for the polymorphic Th2R/Th3R but not the conserved RII+/CS.T3 region in the α-TSR was associated with IGHV3-21/IGVL3-21 or IGLV3-1 gene usage. Although the C terminus specific mAbs showed signs of more efficient affinity maturation and class-switching compared to anti-repeat mAbs, live sporozoite binding and inhibitory activity was limited to a single C-linker reactive mAb with cross-reactivity to the central repeat and junction. The data provide novel insights in the human anti-C-linker and anti-α-TSR antibody response that support exclusion of the PfCSP C terminus from malaria vaccine designs.

Glycosylated nanoparticle-based PfCSP vaccine confers long-lasting antibody responses and sterile protection in mouse malaria model

The development of an effective and durable vaccine remains a central goal in the fight against malaria. Circumsporozoite protein (CSP) is the major surface protein of sporozoites and the target of the only licensed Plasmodium falciparum (Pf) malaria vaccine, RTS,S/AS01. However, vaccine efficacy is low and short-lived, highlighting the need for a second-generation vaccine with superior efficacy and durability. Here, we report a Helicobacter pylori apoferritin-based nanoparticle immunogen that elicits strong B cell responses against PfCSP epitopes that are targeted by the most potent human monoclonal antibodies. Glycan engineering of the scaffold and fusion of an exogenous T cell epitope enhanced the anti-PfCSP B cell response eliciting strong, long-lived and protective humoral immunity in mice. Our study highlights the power of rational vaccine design to generate a highly efficacious second-generation anti-infective malaria vaccine candidate and provides the basis for its further development.

Distribution of Anopheles gambiae thioester-containing protein 1 alleles along malaria transmission gradients in The Gambia

BACKGROUND

Thioester-containing protein 1 (TEP1) is a highly polymorphic gene playing an important role in mosquito immunity to parasite development and associated with Anopheles gambiae vectorial competence. Allelic variations in TEP1 could render mosquito either susceptible or resistant to parasite infection. Despite reports of TEP1 genetic variations in An. gambiae, the correlation between TEP1 allelic variants and transmission patterns in malaria endemic settings remains unclear.

METHODS

TEP1 allelic variants were characterized by PCR from archived genomic DNA of > 1000 An. gambiae mosquitoes collected at 3 time points between 2009 and 2019 from eastern Gambia, where malaria transmission remains moderately high, and western regions with low transmission.

RESULTS

Eight common TEP1 allelic variants were identified at varying frequencies in An. gambiae from both transmission settings. These comprised the wild type TEP1, homozygous susceptible genotype, TEP1s; homozygous resistance genotypes: TEP1rA and TEP1rB, and the heterozygous resistance genotypes: TEP1srA, TEP1srB, TEP1rArB and TEP1srArB. There was no significant disproportionate distribution of the TEP1 alleles by transmission setting and the temporal distribution of alleles was also consistent across the transmission settings. TEP1s was the most common in all vector species in both settings (allele frequencies: East = 21.4-68.4%. West = 23.5-67.2%). In Anopheles arabiensis, the frequency of wild type TEP1 and susceptible TEP1s was significantly higher in low transmission setting than in high transmission setting (TEP1: Z = - 4.831, P < 0.0001; TEP1s: Z = - 2.073, P = 0.038).

CONCLUSIONS

The distribution of TEP1 allele variants does not distinctly correlate with malaria endemicity pattern in The Gambia. Further studies are needed to understand the link between genetic variations in vector population and transmission pattern in the study settings. Future studies on the implication for targeting TEP1 gene for vector control strategy such as gene drive systems in this settings is also recommended.

A closer look at the WHO cone bioassay: video analysis of the hidden effects of a human host on mosquito behaviour and insecticide contact

Background

The WHO cone test is one of three tests currently used to evaluate the efficacy of insecticide-treated bed nets (ITNs). It generates two test outputs, knockdown and 24-h mortality, both indicative of immediate toxicity but that reveal little about the nature of mosquito and ITN interaction or how results translate to real-world settings.

Methods

A human arm held 5 mm behind the net surface acted as a host attractant during cone tests and a smartphone was used to capture mosquito behaviour in the cone. Post-exposure blood feeding and survival for nine days were recorded; ingested blood meal size was determined by measuring excreted haematin. Four strains of Anopheles gambiae (insecticide susceptible: Kisumu and N’gousso; insecticide resistant: Banfora and VK7) were tested with and without the host attractant using untreated, Permanet 2.0 and Olyset nets. Video recordings were scan sampled every five seconds to record mosquito positions on either the net, in flight or in contact with the cone. Generalized estimating equations were used to analyse all data except survival within nine days which was analysed using Weighted Cox Regression.

Results

Net contact was the most frequently recorded behaviour in all Anopheles spp. strains on all nets. Adding the human host as attractant triggered excitatory behaviours: in all strains, the magnitude of net contact was significantly decreased compared to tests without a host. ITN exposure altered the observed behaviour of the two susceptible strains, which exhibited a decreased response to the host during ITN tests. The resistant strains did not alter their behaviour during ITN tests. Significantly less net contact was observed during Olyset Net tests compared to Permanet 2.0. The host presence affected survival after exposure: Banfora and VK7 mosquitoes exposed to Permanet 2.0 with a host lived longer compared to tests performed without a host. However, mosquitoes that blood-fed and survived long enough to digest the blood meal did not exhibit significantly reduced longevity regardless of the presence of the host attractant.

Conclusions

Simple modifications to the WHO cone test and extension of post-test monitoring beyond the current 24 h enable detailed behavioural characterizations of individual ITNs to be compiled. The effects observed from testing with a host and including blood feeding suggest that more representative estimates of true of ITN efficacy are gained with these modifications than when using the current testing protocol.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04232-4.

Genetics and immunity of Anopheles response to the entomopathogenic fungus Metarhizium anisopliae overlap with immunity to Plasmodium

Entomopathogenic fungi have been explored as a potential biopesticide to counteract the insecticide resistance issue in mosquitoes. However, little is known about the possibility that genetic resistance to fungal biopesticides could evolve in mosquito populations. Here, we detected an important genetic component underlying Anopheles coluzzii survival after exposure to the entomopathogenic fungus Metarhizium anisopliae. A familiality study detected variation for survival among wild mosquito isofemale pedigrees, and genetic mapping identified two loci that significantly influence mosquito survival after fungus exposure. One locus overlaps with a previously reported locus for Anopheles susceptibility to the human malaria parasite Plasmodium falciparum. Candidate gene studies revealed that two LRR proteins encoded by APL1C and LRIM1 genes in this newly mapped locus are required for protection of female A. coluzzii from M. anisopliae, as is the complement-like factor Tep1. These results indicate that natural Anopheles populations already segregate frequent genetic variation for differential mosquito survival after fungal challenge and suggest a similarity in Anopheles protective responses against fungus and Plasmodium. However, this immune similarity raises the possibility that fungus-resistant mosquitoes could also display enhanced resistance to Plasmodium, suggesting an advantage of selecting for fungus resistance in vector populations to promote naturally diminished malaria vector competence.

Sympatric Populations of the Anopheles gambiae Complex in Southwest Burkina Faso Evolve Multiple Diverse Resistance Mechanisms in Response to Intense Selection Pressure with Pyrethroids

Simple Summary

Targeting mosquitoes with insecticides is one of the most effective methods to prevent malaria transmission. Although numbers of malaria cases have declined substantially this century, this pattern is not universal and Burkina Faso has one of the highest burdens of malaria; it is also a hotspot for the evolution of insecticide resistance in malaria vectors. We have established laboratory colonies from multiple species within the An. gambiae complex, the most efficient group of malaria vectors in the world, from larval collections in southwest Burkina Faso. Using bioassays with different insecticides widely used to control public health pests, we provide a profile of insecticide resistance in each of these colonies and, using molecular tools, reveal the genetic changes underpinning this resistance. We show that, whilst many resistance mechanisms are shared between species, there are some important differences which may affect resistance to current and future insecticide classes. The complexity, and diversity of resistance mechanisms highlights the importance of screening any potential new insecticide intended for use in malaria control against a wide range of populations. These stable laboratory colonies provide a valuable resource for insecticide discovery, and for further studies on the evolution and dispersal of insecticide resistance within and between species.

Abstract

Pyrethroid resistance in the Anopheles vectors of malaria is driving an urgent search for new insecticides that can be used in proven vector control tools such as insecticide treated nets (ITNs). Screening for potential new insecticides requires access to stable colonies of the predominant vector species that contain the major pyrethroid resistance mechanisms circulating in wild populations. Southwest Burkina Faso is an apparent hotspot for the emergence of pyrethroid resistance in species of the Anopheles gambiae complex. We established stable colonies from larval collections across this region and characterised the resistance phenotype and underpinning genetic mechanisms. Three additional colonies were successfully established (1 An. coluzzii, 1 An. gambiae and 1 An. arabiensis) to add to the 2 An. coluzzii colonies already established from this region; all 5 strains are highly resistant to pyrethroids. Synergism assays found that piperonyl butoxide (PBO) exposure was unable to fully restore susceptibility although exposure to a commercial ITN containing PBO resulted in 100% mortality. All colonies contained resistant alleles of the voltage gated sodium channel but with differing proportions of alternative resistant haplotypes. RNAseq data confirmed the role of P450s, with CYP6P3 and CYP6Z2 elevated in all 5 strains, and identified many other resistance mechanisms, some found across strains, others unique to a particular species. These strains represent an important resource for insecticide discovery and provide further insights into the complex genetic changes driving pyrethroid resistance.

Multi-insecticide resistant malaria vectors in the field remain susceptible to malathion, despite the presence of Ace1 point mutations

Insecticide resistance in Anopheles mosquitoes is seriously threatening the success of insecticide-based malaria vector control. Surveillance of insecticide resistance in mosquito populations and identifying the underlying mechanisms enables optimisation of vector control strategies. Here, we investigated the molecular mechanisms of insecticide resistance in three Anopheles coluzzii field populations from southern Côte d’Ivoire, including Agboville, Dabou and Tiassalé. All three populations were resistant to bendiocarb, deltamethrin and DDT, but not or only very weakly resistant to malathion. The absence of malathion resistance is an unexpected result because we found the acetylcholinesterase mutation Ace1-G280S at high frequencies, which would typically confer cross-resistance to carbamates and organophosphates, including malathion. Notably, Tiassalé was the most susceptible population to malathion while being the most resistant one to the pyrethroid deltamethrin. The resistance ratio to deltamethrin between Tiassalé and the laboratory reference colony was 1,800 fold. By sequencing the transcriptome of individual mosquitoes, we found numerous cytochrome P450-dependent monooxygenases – including CYP6M2, CYP6P2, CYP6P3, CYP6P4 and CYP6P5 – overexpressed in all three field populations. This could be an indication for negative cross-resistance caused by overexpression of pyrethroid-detoxifying cytochrome P450s that may activate pro-insecticides, thereby increasing malathion susceptibility. In addition to the P450s, we found several overexpressed carboxylesterases, glutathione S-transferases and other candidates putatively involved in insecticide resistance.

Insecticide-based mosquito control has saved millions of lives from malaria and other vector-borne diseases. However, the emergence and increase of insecticide resistant Anopheles populations seriously threaten to derail malaria control programmes. Surveillance of insecticide resistance and understanding the underlying molecular mechanisms are key for choosing effective vector control strategies. Here, we characterised the degree and mechanisms of resistance in three malaria vector populations from Côte d’Ivoire. Our key finding was that these multi-insecticide resistant malaria vectors largely remained susceptible to malathion, despite the presence of a mutation in the target enzyme of this organophosphate insecticide that would typically confer resistance. Intriguingly, we found overexpression of metabolic P450 enzymes that are known to detoxify insecticides and activate pro-insecticides such as malathion. It is highly probable that, here, we observed P450-mediated negative cross-resistance for the first time in Anopheles field populations. Negative cross-resistance merits further investigation as advantage could be taken of this phenomenon in the fight against multi-resistant malaria vectors.

Assessment of Plasmodium falciparum Infection and Fitness of Genetically Modified Anopheles gambiae Aimed at Mosquito Population Replacement

Genetically modified (GM) mosquitoes expressing anti-plasmodial effectors propagating through wild mosquito populations by means of gene drive is a promising tool to support current malaria control strategies. The process of generating GM mosquitoes involves genetic transformation of mosquitoes from a laboratory colony and, often, interbreeding with other GM lines to cross in auxiliary traits. These mosquito colonies and GM lines thus often have different genetic backgrounds and GM lines are invariably highly inbred, which in conjunction with their independent rearing in the laboratory may translate to differences in their susceptibility to malaria parasite infection and life history traits. Here, we show that laboratory Anopheles gambiae colonies and GM lines expressing Cas9 and Cre recombinase vary greatly in their susceptibility to Plasmodium falciparum NF54 infection. Therefore, the choice of mosquitoes to be used as a reference when conducting infection or life history trait assays requires careful consideration. To address these issues, we established an experimental pipeline involving genetic crosses and genotyping of mosquitoes reared in shared containers throughout their lifecycle. We used this protocol to examine whether GM lines expressing the antimicrobial peptide (AMP) Scorpine in the mosquito midgut interfere with parasite infection and mosquito survival. We demonstrate that Scorpine expression in the Peritrophin 1 (Aper1) genomic locus reduces both P. falciparum sporozoite prevalence and mosquito lifespan; both these phenotypes are likely to be associated with the disturbance of the midgut microbiota homeostasis. These data lead us to conclude that the Aper1-Sco GM line could be used in proof-of-concept experiments aimed at mosquito population replacement, although the impact of its reduced fitness on the spread of the transgene through wild populations requires further investigation.

Susceptibility of Field-Collected Nyssorhynchus darlingi to Plasmodium spp. in Western Amazonian Brazil

Mosquito susceptibility to Plasmodium&nbsp;spp. infection is of paramount importance for malaria occurrence and sustainable transmission. Therefore, understanding the genetic features underlying the mechanisms of susceptibility traits is pivotal to assessing malaria transmission dynamics in endemic areas. The aim of this study was to investigate the susceptibility of Nyssorhynchus darlingi-the dominant malaria vector in Brazil-to Plasmodium spp. using a reduced representation genome-sequencing protocol. The investigation was performed using a genome-wide association study (GWAS) to identify mosquito genes that are predicted to modulate the susceptibility of natural populations of the mosquito to Plasmodium infection. After applying the sequence alignment protocol, we generated the variant panel and filtered variants; leading to the detection of 202,837 SNPs in all specimens analyzed. The resulting panel was used to perform GWAS by comparing the pool of SNP variants present in Ny.&nbsp;darlingi infected with Plasmodium spp. with the pool obtained in field-collected mosquitoes with no evidence of infection by the parasite (all mosquitoes were tested separately using RT-PCR). The GWAS results for infection status showed two statistically significant variants adjacent to important genes that can be associated with susceptibility to Plasmodium infection: Cytochrome P450 (cyp450) and chitinase. This study provides relevant knowledge on malaria transmission dynamics by using a genomic approach to identify mosquito genes associated with susceptibility to Plasmodium infection in Ny. darlingi in western Amazonian Brazil.

Malaria-Resistant Mosquitoes (Diptera: Culicidae); The Principle is Proven, But Will the Effectors Be Effective?

Over the last few decades, a substantial number of anti-malarial effector genes have been evaluated for their ability to block parasite infection in the mosquito vector. While many of these approaches have yielded significant effects on either parasite intensity or prevalence of infection, just a few have been able to completely block transmission. Additionally, many approaches, while effective against the parasite, also disrupt or alter important aspects of mosquito physiology, leading to corresponding changes in lifespan, reproduction, and immunity. As the most promising approaches move towards field-based evaluation, questions of effector gene robustness and durability move to the forefront. In this forum piece, we critically evaluate past effector gene approaches with an eye towards developing a deeper pipeline to augment the current best candidates.

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.

Impact of field-realistic doses of glyphosate and nutritional stress on mosquito life history traits and susceptibility to malaria parasite infection

Glyphosate is the world's most widely used herbicide. The commercial success of this molecule is due to its nonselectivity and its action, which would supposedly target specific biosynthetic pathways found mainly in plants. Multiple studies have however provided evidence for high sensitivity of many nontarget species to glyphosate and/or to formulations (glyphosate mixed with surfactants). This herbicide, found at significant levels in aquatic systems through surface runoffs, impacts life history traits and immune parameters of several aquatic invertebrates' species, including disease-vector mosquitoes. Mosquitoes, from hatching to emergence, are exposed to aquatic chemical contaminants. In this study, we first compared the toxicity of pure glyphosate to the toxicity of glyphosate-based formulations for the main vector of avian malaria in Europe, Culex pipiens mosquito. Then we evaluated, for the first time, how field-realistic dose of glyphosate interacts with larval nutritional stress to alter mosquito life history traits and susceptibility to avian malaria parasite infection. Our results show that exposure of larvae to field-realistic doses of glyphosate, pure or in formulation, did not affect larval survival rate, adult size, and female fecundity. One of our two experimental blocks showed, however, that exposure to glyphosate decreased development time and reduced mosquito infection probability by malaria parasite. Interestingly, the effect on malaria infection was lost when the larvae were also subjected to a nutritional stress, probably due to a lower ingestion of glyphosate.

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

Comparison of infectivity of Plasmodium vivax to wild-caught and laboratory-adapted (colonized) Anopheles arabiensis mosquitoes in Ethiopia

BACKGROUND

Mosquito-feeding assays that assess transmission of Plasmodium from man-to-mosquito typically use laboratory mosquito colonies. The microbiome and genetic background of local mosquitoes may be different and influence Plasmodium transmission efficiency. In order to interpret transmission studies to the local epidemiology, it is therefore crucial to understand the relationship between infectivity in laboratory-adapted and local mosquitoes.

METHODS

We assessed infectivity of Plasmodium vivax-infected patients from Adama, Ethiopia, using laboratory-adapted (colony) and wild-caught (wild) mosquitoes raised from larval collections in paired feeding experiments. Feeding assays used 4-6 day-old female Anopheles arabiensis mosquitoes after starvation for 12 h (colony) and 18 h (wild). Oocyst development was assessed microscopically 7 days post-feeding. Wild mosquitoes were identified morphologically and confirmed by genotyping. Asexual parasites and gametocytes were quantified in donor blood by microscopy.

RESULTS

In 36 paired experiments (25 P. vivax infections and 11 co-infections with P. falciparum), feeding efficiency was higher in colony (median: 62.5%; interquartile range, IQR: 47.0-79.0%) compared to wild mosquitoes (median: 27.8%; IQR: 17.0-38.0%; Z = 5.02; P < 0.001). Plasmodium vivax from infectious individuals (51.6%, 16/31) infected a median of 55.0% (IQR: 6.7-85.7%; range: 5.5-96.7%; n = 14) of the colony and 52.7% (IQR: 20.0-80.0%; range: 3.2-95.0%; n = 14) of the wild mosquitoes. A strong association (ρ(16) = 0.819; P < 0.001) was observed between the proportion of infected wild and colony mosquitoes. A positive association was detected between microscopically detected gametocytes and the proportion of infected colony (ρ(31) = 0.452; P = 0.011) and wild (ρ(31) = 0.386; P = 0.032) mosquitoes.

CONCLUSIONS

Infectivity assessments with colony and wild mosquitoes yielded similar infection results. This finding supports the use of colony mosquitoes for assessments of the infectious reservoir for malaria in this setting whilst acknowledging the importance of mosquito factors influencing sporogonic development of Plasmodium parasites.

Gene copy number and function of the APL1 immune factor changed during Anopheles evolution

Background

The recent reference genome assembly and annotation of the Asian malaria vector Anopheles stephensi detected only one gene encoding the leucine-rich repeat immune factor APL1, while in the Anopheles gambiae and sibling Anopheles coluzzii, APL1 factors are encoded by a family of three paralogs. The phylogeny and biological function of the unique APL1 gene in An. stephensi have not yet been specifically examined.

Methods

The APL1 locus was manually annotated to confirm the computationally predicted single APL1 gene in An. stephensi. APL1 evolution within Anopheles was explored by phylogenomic analysis. The single or paralogous APL1 genes were silenced in An. stephensi and An. coluzzii, respectively, followed by mosquito survival analysis, experimental infection with Plasmodium and expression analysis.

Results

APL1 is present as a single ancestral gene in most Anopheles including An. stephensi but has expanded to three paralogs in an African lineage that includes only the Anopheles gambiae species complex and Anopheles christyi. Silencing of the unique APL1 copy in An. stephensi results in significant mosquito mortality. Elevated mortality of APL1-depleted An. stephensi is rescued by antibiotic treatment, suggesting that pathology due to bacteria is the cause of mortality, and indicating that the unique APL1 gene is essential for host survival. Successful Plasmodium development in An. stephensi depends upon APL1 activity for protection from high host mortality due to bacteria. In contrast, silencing of all three APL1 paralogs in An. coluzzii does not result in elevated mortality, either with or without Plasmodium infection. Expression of the single An. stephensi APL1 gene is regulated by both the Imd and Toll immune pathways, while the two signaling pathways regulate different APL1 paralogs in the expanded APL1 locus.

Conclusions

APL1 underwent loss and gain of functions concomitant with expansion from a single ancestral gene to three paralogs in one lineage of African Anopheles. We infer that activity of the unique APL1 gene promotes longevity in An. stephensi by conferring protection from or tolerance to an effect of bacterial pathology. The evolution of an expanded APL1 gene family could be a factor contributing to the exceptional levels of malaria transmission mediated by human-feeding members of the An. gambiae species complex in Africa.

A sensory appendage protein protects malaria vectors from pyrethroids

Pyrethroid-impregnated bed nets have driven considerable reductions in malaria-associated morbidity and mortality in Africa since the beginning of the century1. The intense selection pressure exerted by bed nets has precipitated widespread and escalating resistance to pyrethroids in African Anopheles populations, threatening to reverse the gains that been made by malaria control2. Here we show that expression of a sensory appendage protein (SAP2), which is enriched in the legs, confers pyrethroid resistance to Anopheles gambiae. Expression of SAP2 is increased in insecticide-resistant populations and is further induced after the mosquito comes into contact with pyrethroids. SAP2 silencing fully restores mortality of the mosquitoes, whereas SAP2 overexpression results in increased resistance, probably owing to high-affinity binding of SAP2 to pyrethroid insecticides. Mining of genome sequence data reveals a selective sweep near the SAP2 locus in the mosquito populations of three West African countries (Cameroon, Guinea and Burkina Faso) with the observed increase in haplotype-associated single-nucleotide polymorphisms mirroring the increasing resistance of mosquitoes to pyrethroids reported in Burkina Faso. Our study identifies a previously undescribed mechanism of insecticide resistance that is likely to be highly relevant to malaria control efforts.

Influence of genetic polymorphism on transcriptional enhancer activity in the malaria vector Anopheles coluzzii

Enhancers are cis-regulatory elements that control most of the developmental and spatial gene expression in eukaryotes. Genetic variation of enhancer sequences is known to influence phenotypes, but the effect of enhancer variation upon enhancer functional activity and downstream phenotypes has barely been examined in any species. In the African malaria vector, Anopheles coluzzii, we identified candidate enhancers in the proximity of genes relevant for immunity, insecticide resistance, and development. The candidate enhancers were functionally validated using luciferase reporter assays, and their activity was found to be essentially independent of their physical orientation, a typical property of enhancers. All of the enhancers segregated genetically polymorphic alleles, which displayed significantly different levels of functional activity. Deletion mutagenesis and functional testing revealed a fine structure of positive and negative regulatory elements that modulate activity of the enhancer core. Enhancer polymorphisms segregate in wild A. coluzzii populations in West Africa. Thus, enhancer variants that modify target gene expression leading to likely phenotypic consequences are frequent in nature. These results demonstrate the existence of naturally polymorphic A. coluzzii enhancers, which may help explain important differences between individuals or populations for malaria transmission efficiency and vector adaptation to the environment.

Interaction of RNA viruses of the natural virome with the African malaria vector, Anopheles coluzzii

Mosquitoes are colonized by a little-studied natural virome. Like the bacterial microbiome, the virome also probably influences the biology and immunity of mosquito vector populations, but tractable experimental models are lacking. We recently discovered two novel viruses in the virome of wild Anopheles and in colonies of the malaria vector Anopheles coluzzii: Anopheles C virus and Anopheles cypovirus. Here, we describe biological interactions between these two viruses and An. coluzzii mosquitoes. Viral abundance varies reproducibly during mosquito development. DNA forms of these viruses were not detected, and thus viral persistence is likely based on vertical transmission of RNA genomes. At least Anopheles C virus is vertically transmitted by an intraembryonic route. Relative abundance of the two viruses is inversely correlated in individual mosquitoes. One possible mechanism for this could be interactions with host immunity, and functional genomic analysis indicated differential influence of at least the Toll and JAK/STAT immune signaling pathways upon the viruses. The nonrandom distributions and interactions with host immunity suggest that these and other members of the natural virome may constitute a source of unrecognized heterogeneity in mosquito vector populations.

Transcriptomic meta-signatures identified in Anopheles gambiae populations reveal previously undetected insecticide resistance mechanisms

Increasing insecticide resistance in malaria-transmitting vectors represents a public health threat, but underlying mechanisms are poorly understood. Here, a data integration approach is used to analyse transcriptomic data from comparisons of insecticide resistant and susceptible Anopheles populations from disparate geographical regions across the African continent. An unbiased, integrated analysis of this data confirms previously described resistance candidates but also identifies multiple novel genes involving alternative resistance mechanisms, including sequestration, and transcription factors regulating multiple downstream effector genes, which are validated by gene silencing. The integrated datasets can be interrogated with a bespoke Shiny R script, deployed as an interactive web-based application, that maps the expression of resistance candidates and identifies co-regulated transcripts that may give clues to the function of novel resistance-associated genes.

Increasing insecticide resistance of mosquitoes represents a public health threat, and underlying mechanisms are poorly understood. Here, Ingham et al. identify putative insecticide resistance genes in Anopheles gambiae populations across Africa and develop a web-based application that maps their expression.

Highly focused transcriptional response of Anopheles coluzzii to O'nyong nyong arbovirus during the primary midgut infection

Background

Anopheles mosquitoes are efficient vectors of human malaria, but it is unknown why they do not transmit viruses as well as Aedes and Culex mosquitoes. The only arbovirus known to be consistently transmitted by Anopheles mosquitoes is O’nyong nyong virus (ONNV, genus Alphavirus, family Togaviridae). The interaction of Anopheles mosquitoes with RNA viruses has been relatively unexamined.

Results

We transcriptionally profiled the African malaria vector, Anopheles coluzzii, infected with ONNV. Mosquitoes were fed on an infectious bloodmeal and were analyzed by Illumina RNAseq at 3 days post-bloodmeal during the primary virus infection of the midgut epithelium, before systemic dissemination. Virus infection triggers transcriptional regulation of just 30 host candidate genes. Most of the regulated candidate genes are novel, without known function. Of the known genes, a significant cluster includes candidates with predicted involvement in carbohydrate metabolism. Two candidate genes encoding leucine-rich repeat immune (LRIM) factors point to possible involvement of immune protein complexes in the mosquito antiviral response. The primary ONNV infection by bloodmeal shares little transcriptional response in common with ONNV infection by intrathoracic injection, nor with midgut infection by the malaria parasites, Plasmodium falciparum or P. berghei. Profiling of A. coluzzii microRNA (miRNA) identified 118 known miRNAs and 182 potential novel miRNA candidates, with just one miRNA regulated by ONNV infection. This miRNA was not regulated by other previously reported treatments, and may be virus specific. Coexpression analysis of miRNA abundance and messenger RNA expression revealed discrete clusters of genes regulated by Imd and JAK/STAT, immune signaling pathways that are protective against ONNV in the primary infection.

Conclusions

ONNV infection of the A. coluzzii midgut triggers a remarkably limited gene regulation program of mostly novel candidate genes, which likely includes host genes deployed for antiviral defense, as well as genes manipulated by the virus to facilitate infection. Functional dissection of the ONNV-response candidate genes is expected to generate novel insight into the mechanisms of virus-vector interaction.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4918-0) contains supplementary material, which is available to authorized users.

The Anopheles gambiae ATP-binding cassette transporter family: phylogenetic analysis and tissue localization provide clues on function and role in insecticide resistance

The role of ATP-binding cassette (ABC) transporters in conferring insecticide resistance has received much attention recently. Here we identify ABC transporters differentially expressed in insecticide-resistant populations of the malaria vector, Anopheles gambiae. Although we found little evidence that the orthologues of the multidrug resistance proteins described in other species are associated with resistance in An. gambiae we did identify a subset of ABC proteins consistently differentially expressed in pyrethroid-resistant populations from across Africa. We present information on the phylogenetic relationship, primary sites of expression and potential role of ABC transporters in mediating the mosquito's response to insecticides. Furthermore we demonstrate that a paralogous group of eight ABCG transporters, clustered on chromosome 3R, are highly enriched in the legs of An. gambiae mosquitoes, consistent with a proposed role for this ABC subfamily in transport of lipids to the outer surface of the cuticle. Finally, antibodies raised against one of the most highly expressed ABC transporters in adult females, ABCG7 (AGAP009850), localized this transporter to the pericardial cells. These data will help prioritize members of this gene family for further localization and functional validation studies to identify the in vivo function of these transporters in the mosquito and determine whether elevated expression of members of this family contribute to insecticide resistance.

Transmission traits of malaria parasites within the mosquito: Genetic variation, phenotypic plasticity, and consequences for control

Evaluating the risk of emergence and transmission of vector‐borne diseases requires knowledge of the genetic and environmental contributions to pathogen transmission traits. Compared to the significant effort devoted to understanding the biology of malaria transmission from vertebrate hosts to mosquito vectors, the strategies that malaria parasites have evolved to maximize transmission from vectors to vertebrate hosts have been largely overlooked. While determinants of infection success within the mosquito host have recently received attention, the causes of variability for other key transmission traits of malaria, namely the duration of parasite development and its virulence within the vector, as well as its ability to alter mosquito behavior, remain largely unknown. This important gap in our knowledge needs to be bridged in order to obtain an integrative view of the ecology and evolution of malaria transmission strategies. Associations between transmission traits also need to be characterized, as they trade‐offs and constraints could have important implications for understanding the evolution of parasite transmission. Finally, theoretical studies are required to evaluate how genetic and environmental influences on parasite transmission traits can shape malaria dynamics and evolution in response to disease control.

Microbial Pre-exposure and Vectorial Competence of Anopheles Mosquitoes

Anopheles female mosquitoes can transmit Plasmodium, the malaria parasite. During their aquatic life, wild Anopheles mosquito larvae are exposed to a huge diversity of microbes present in their breeding sites. Later, adult females often take successive blood meals that might also carry different micro-organisms, including parasites, bacteria, and viruses. Therefore, prior to Plasmodium ingestion, the mosquito biology could be modulated at different life stages by a suite of microbes present in larval breeding sites, as well as in the adult environment. In this article, we highlight several naturally relevant scenarios of Anopheles microbial pre-exposure that we assume might impact mosquito vectorial competence for the malaria parasite: (i) larval microbial exposures; (ii) protist co-infections; (iii) virus co-infections; and (iv) pathogenic bacteria co-infections. In addition, significant behavioral changes in African Anopheles vectors have been associated with increasing insecticide resistance. We discuss how these ethological modifications may also increase the repertoire of microbes to which mosquitoes could be exposed, and that might also influence their vectorial competence. Studying Plasmodium–Anopheles interactions in natural microbial environments would efficiently contribute to refining the transmission risks.

Natural Parasite Exposure Induces Protective Human Anti-Malarial Antibodies

Antibodies against the NANP repeat of circumsporozoite protein (CSP), the major surface antigen of Plasmodium falciparum (Pf) sporozoites, can protect from malaria in animal models but protective humoral immunity is difficult to induce in humans. Here we cloned and characterized rare affinity-matured human NANP-reactive memory B cell antibodies elicited by natural Pf exposure that potently inhibited parasite transmission and development in vivo. We unveiled the molecular details of antibody binding to two distinct protective epitopes within the NANP repeat. NANP repeat recognition was largely mediated by germline encoded and immunoglobulin (Ig) heavy-chain complementarity determining region 3 (HCDR3) residues, whereas affinity maturation contributed predominantly to stabilizing the antigen-binding site conformation. Combined, our findings illustrate the power of exploring human anti-CSP antibody responses to develop tools for malaria control in the mammalian and the mosquito vector and provide a molecular basis for the structure-based design of next-generation CSP malaria vaccines.

The Anopheles gambiae 2La chromosome inversion is associated with susceptibility to Plasmodium falciparum in Africa

Chromosome inversions suppress genetic recombination and establish co-adapted gene complexes, or supergenes. The 2La inversion is a widespread polymorphism in the Anopheles gambiae species complex, the major African mosquito vectors of human malaria. Here we show that alleles of the 2La inversion are associated with natural malaria infection levels in wild-captured vectors from West and East Africa. Mosquitoes carrying the more-susceptible allele (2L+^a) are also behaviorally less likely to be found inside houses. Vector control tools that target indoor-resting mosquitoes, such as bednets and insecticides, are currently the cornerstone of malaria control in Africa. Populations with high levels of the 2L+^a allele may form reservoirs of persistent outdoor malaria transmission requiring novel measures for surveillance and control. The 2La inversion is a major and previously unappreciated component of the natural malaria transmission system in Africa, influencing both malaria susceptibility and vector behavior.

Interactive cost of Plasmodium infection and insecticide resistance in the malaria vector Anopheles gambiae

Insecticide resistance raises concerns for the control of vector-borne diseases. However, its impact on parasite transmission could be diverse when considering the ecological interactions between vector and parasite. Thus we investigated the fitness cost associated with insecticide resistance and Plasmodium falciparum infection as well as their interactive cost on Anopheles gambiae survival and fecundity. In absence of infection, we observed a cost on fecundity associated with insecticide resistance. However, survival was higher for mosquito bearing the kdr mutation and equal for those with the ace-1(R) mutation compared to their insecticide susceptible counterparts. Interestingly, Plasmodium infection reduced survival only in the insecticide resistant strains but not in the susceptible one and infection was associated with an increase in fecundity independently of the strain considered. This study provides evidence for a survival cost associated with infection by Plasmodium parasite only in mosquito selected for insecticide resistance. This suggests that the selection of insecticide resistance mutation may have disturbed the interaction between parasites and vectors, resulting in increased cost of infection. Considering the fitness cost as well as other ecological aspects of this natural mosquito-parasite combination is important to predict the epidemiological impact of insecticide resistance.

Population genetics of the Asian tiger mosquito Aedes albopictus, an invasive vector of human diseases

The Asian tiger mosquito Aedes albopictus is currently one of the most threatening invasive species in the world. Native to Southeast Asia, the species has spread throughout the world in the past 30 years and is now present in every continent but Antarctica. Because it was the main vector of recent Dengue and Chikungunya outbreaks, and because of its competency for numerous other viruses and pathogens such as the Zika virus, A. albopictus stands out as a model species for invasive diseases vector studies. A synthesis of the current knowledge about the genetic diversity of A. albopictus is needed, knowing the interplays between the vector, the pathogens, the environment and their epidemiological consequences. Such resources are also valuable for assessing the role of genetic diversity in the invasive success. We review here the large but sometimes dispersed literature about the population genetics of A. albopictus. We first debate about the experimental design of these studies and present an up-to-date assessment of the available molecular markers. We then summarize the main genetic characteristics of natural populations and synthesize the available data regarding the worldwide structuring of the vector. Finally, we pinpoint the gaps that remain to be addressed and suggest possible research directions.

The potential for adaptive maintenance of diversity in insect antimicrobial peptides

Genes involved in immune defence are among the fastest evolving in the genomes of many species. Interestingly, however, genes encoding antimicrobial peptides (AMPs) have shown little evidence for adaptive divergence in arthropods, despite the centrality of these peptides in direct killing of microbial pathogens. This observation, coupled with a failure to detect phenotypic consequence of genetic variation in AMPs, has led to the hypothesis that individual AMPs make minor contributions to overall immune defence and that AMPs instead act as a collective cocktail. Recent data, however, have suggested an alternative explanation for the apparent lack of adaptive divergence in AMP genes. Molecular evolutionary and phenotypic data have begun to suggest that variant AMP alleles may be maintained through balancing selection in invertebrates, a pattern similar to that observed in several vertebrate AMPs. Signatures of balancing selection include high rates of non-synonymous polymorphism, trans-species amino acid polymorphisms, and convergence of amino acid states across the phylogeny. In this review, we revisit published literature on insect AMP genes and analyse newly available population genomic datasets in Drosophila, finding enrichment for patterns consistent with adaptive maintenance of polymorphism.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Identification and Characterization of Two Novel RNA Viruses from Anopheles gambiae Species Complex Mosquitoes

Mosquitoes of the Anopheles gambiae complex display strong preference for human bloodmeals and are major malaria vectors in Africa. However, their interaction with viruses or role in arbovirus transmission during epidemics has been little examined, with the exception of O'nyong-nyong virus, closely related to Chikungunya virus. Deep-sequencing has revealed different RNA viruses in natural insect viromes, but none have been previously described in the Anopheles gambiae species complex. Here, we describe two novel insect RNA viruses, a Dicistrovirus and a Cypovirus, found in laboratory colonies of An. gambiae taxa using small-RNA deep sequencing. Sequence analysis was done with Metavisitor, an open-source bioinformatic pipeline for virus discovery and de novo genome assembly. Wild-collected Anopheles from Senegal and Cambodia were positive for the Dicistrovirus and Cypovirus, displaying high sequence identity to the laboratory-derived virus. Thus, the Dicistrovirus (Anopheles C virus, AnCV) and Cypovirus (Anopheles Cypovirus, AnCPV) are components of the natural virome of at least some anopheline species. Their possible influence on mosquito immunity or transmission of other pathogens is unknown. These natural viruses could be developed as models for the study of Anopheles-RNA virus interactions in low security laboratory settings, in an analogous manner to the use of rodent malaria parasites for studies of mosquito anti-parasite immunity.

Exposure to deltamethrin affects development of Plasmodium falciparum inside wild pyrethroid resistant Anopheles gambiae s.s. mosquitoes in Uganda

Background

Pyrethroid resistance in African vector mosquitoes is a threat to malaria control. Resistant mosquitoes can survive insecticide doses that would normally be lethal. We studied effects of such doses on Plasmodium falciparum development inside kdr-resistant Anopheles gambiae s.s. in Uganda.

Methods

We collected An. gambiae s.s. homozygous for kdr-L1014S mutation, fed them on blood samples from 42 P. falciparum-infected local patients, then exposed them either to nets treated with sub-lethal doses of deltamethrin or to untreated nets. After seven days, we dissected 692 mosquitoes and examined their midguts for oocysts. Prevalence (proportion infected) and intensity of infection (number of oocysts per infected mosquito) were recorded for each group.

Results

Both prevalence and intensity of infection were significantly reduced in deltamethrin-exposed mosquitoes, compared to those exposed to untreated nets. With low doses (2.5–5.0 mg/m²), prevalence was reduced by 59 % (95 % CI = 22 %-78 %) and intensity by 41 % (95 % CI = 25 %-54 %). With high doses (10–16.7 mg/m²), prevalence was reduced by 80 % (95 % CI = 67 %-88 %) and intensity by 34 % (95 % CI = 20 %-46 %).

Conclusions

We showed that, with locally-sampled parasites and mosquitoes, doses of pyrethroids that are sub-lethal for resistant mosquitoes can interfere with parasite development inside mosquitoes. This mechanism could enable pyrethroid-treated nets to prevent malaria transmission despite increasing vector resistance.

An Evolution-Based Screen for Genetic Differentiation between Anopheles Sister Taxa Enriches for Detection of Functional Immune Factors

Nucleotide variation patterns across species are shaped by the processes of natural selection, including exposure to environmental pathogens. We examined patterns of genetic variation in two sister species, Anopheles gambiae and Anopheles coluzzii, both efficient natural vectors of human malaria in West Africa. We used the differentiation signature displayed by a known coordinate selective sweep of immune genes APL1 and TEP1 in A. coluzzii to design a population genetic screen trained on the sweep, classified a panel of 26 potential immune genes for concordance with the signature, and functionally tested their immune phenotypes. The screen results were strongly predictive for genes with protective immune phenotypes: genes meeting the screen criteria were significantly more likely to display a functional phenotype against malaria infection than genes not meeting the criteria (p = 0.0005). Thus, an evolution-based screen can efficiently prioritize candidate genes for labor-intensive downstream functional testing, and safely allow the elimination of genes not meeting the screen criteria. The suite of immune genes with characteristics similar to the APL1-TEP1 selective sweep appears to be more widespread in the A. coluzzii genome than previously recognized. The immune gene differentiation may be a consequence of adaptation of A. coluzzii to new pathogens encountered in its niche expansion during the separation from A. gambiae, although the role, if any of natural selection by Plasmodium is unknown. Application of the screen allowed identification of new functional immune factors, and assignment of new functions to known factors. We describe biochemical binding interactions between immune proteins that underlie functional activity for malaria infection, which highlights the interplay between pathogen specificity and the structure of immune complexes. We also find that most malaria-protective immune factors display phenotypes for either human or rodent malaria, with broad specificity a rarity.

Anopheles gambiae and Anopheles coluzzii are the primary mosquito vectors of human malaria in West Africa. Both of these closely related species efficiently transmit the disease, although they display ecological differences. Previous work showed that A. coluzzii displays distinct genetic patterns in genes important for mosquito immunity. Here, we use this genetic pattern as a filter to examine a panel of potential immune genes, and show that the genetic pattern is strongly predictive for genes that play a functional role in immunity when tested with malaria parasites.

Association mapping by pooled sequencing identifies TOLL 11 as a protective factor against Plasmodium falciparum in Anopheles gambiae

Background

The genome-wide association study (GWAS) techniques that have been used for genetic mapping in other organisms have not been successfully applied to mosquitoes, which have genetic characteristics of high nucleotide diversity, low linkage disequilibrium, and complex population stratification that render population-based GWAS essentially unfeasible at realistic sample size and marker density.

Methods

We designed a novel mapping strategy for the mosquito system that combines the power of linkage mapping with the resolution afforded by genetic association. We established founder colonies from West Africa, controlled for diversity, linkage disequilibrium and population stratification. Colonies were challenged by feeding on the infectious stage of the human malaria parasite, Plasmodium falciparum, mosquitoes were phenotyped for parasite load, and DNA pools for phenotypically similar mosquitoes were Illumina sequenced. Phenotype-genotype mapping was carried out in two stages, coarse and fine.

Results

In the first mapping stage, pooled sequences were analysed genome-wide for intervals displaying relativereduction in diversity between phenotype pools, and candidate genomic loci were identified for influence upon parasite infection levels. In the second mapping stage, focused genotyping of SNPs from the first mapping stage was carried out in unpooled individual mosquitoes and replicates. The second stage confirmed significant SNPs in a locus encoding two Toll-family proteins. RNAi-mediated gene silencing and infection challenge revealed that TOLL 11 protects mosquitoes against P. falciparum infection.

Conclusions

We present an efficient and cost-effective method for genetic mapping using natural variation segregating in defined recent Anopheles founder colonies, and demonstrate its applicability for mapping in a complex non-model genome. This approach is a practical and preferred alternative to population-based GWAS for first-pass mapping of phenotypes in Anopheles. This design should facilitate mapping of other traits involved in physiology, epidemiology, and behaviour.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2009-z) contains supplementary material, which is available to authorized users.

The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing

Background

Members of the Anophelesgambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anophelescoluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known.

Methods

Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility.

Results

A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers.

Conclusions

Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0924-8) contains supplementary material, which is available to authorized users.

Connecting genotypes to medically relevant phenotypes in major vector mosquitoes

Transmission of mosquito-borne human disease relies on vectors maintaining strong human host preference and continued susceptibility to disease-causing pathogens or parasites. These traits are affected by the genetics and the environments of all involved organisms, and genotypic interactions are common between parasite and vector, and between virus and vector. A recent study on Aedes host preference has exploited natural genetic variation to make great progress. Here I review our current understanding of the genetic basis of transmission-relevant traits in Anopheles and Aedes, highlighting additional research areas that would benefit from the integration of natural genetic variation.

Three-way interactions between mosquito population, viral strain and temperature underlying chikungunya virus transmission potential

Interactions between pathogens and their insect vectors in nature are under the control of both genetic and non-genetic factors, yet most studies on mosquito vector competence for human pathogens are conducted in laboratory systems that do not consider genetic and/or environmental variability. Evaluating the risk of emergence of arthropod-borne viruses (arboviruses) of public health importance such as chikungunya virus (CHIKV) requires a more realistic appraisal of genetic and environmental contributions to vector competence. In particular, sources of variation do not necessarily act independently and may combine in the form of interactions. Here, we measured CHIKV transmission potential by the mosquito Aedes albopictus in all combinations of six worldwide vector populations, two virus strains and two ambient temperatures (20°C and 28°C). Overall, CHIKV transmission potential by Ae. albopictus strongly depended on the three-way combination of mosquito population, virus strain and temperature. Such genotype-by-genotype-by-environment (G × G × E) interactions question the relevance of vector competence studies conducted with a simpler set of conditions. Our results highlight the need to account for the complex interplay between vectors, pathogens and environmental factors to accurately assess the potential of vector-borne diseases to emerge.

Biophysical analysis of anopheles gambiae leucine-rich repeat proteins APL1A1, APL1B [corrected] and APL1C and their interaction with LRIM1

Natural infection of Anopheles gambiae by malaria-causing Plasmodium parasites is significantly influenced by the APL1 genetic locus. The locus contains three closely related leucine-rich repeat (LRR) genes, APL1A, APL1B and APL1C. Multiple studies have reported the participation of APL1A-C in the immune response of A. gambiae to invasion by both rodent and human Plasmodium isolates. APL1C forms a heterodimer with the related LRR protein LRIM1 via a C-terminal coiled-coil domain that is also present in APL1A and APL1B. The LRIM1/APL1C heterodimer protects A. gambiae from infection by binding the complement-like protein TEP1 to form a stable and active immune complex. Here we report solution x-ray scatting data for the LRIM1/APL1C heterodimer, the oligomeric state of LRIM1/APL1 LRR domains in solution and the crystal structure of the APL1B LRR domain. The LRIM1/APL1C heterodimeric complex has a flexible and extended structure in solution. In contrast to the APL1A, APL1C and LRIM1 LRR domains, the APL1B LRR domain is a homodimer. The crystal structure of APL1B-LRR shows that the homodimer is formed by an N-terminal helix that complements for the absence of an N-terminal capping motif in APL1B, which is a unique distinction within the LRIM1/APL1 protein family. Full-length APL1A1 and APL1B form a stable complex with LRIM1. These results support a model in which APL1A1, APL1B and APL1C can all form an extended, flexible heterodimer with LRIM1, providing a repertoire of functional innate immune complexes to protect A. gambiae from a diverse array of pathogens.

Seasonal variation in wing size and shape between geographic populations of the malaria vector, Anopheles coluzzii in Burkina Faso (West Africa)

The mosquito, Anopheles coluzzii is a major vector of human malaria in Africa with widespread distribution throughout the continent. The species hence populates a wide range of environments in contrasted ecological settings often exposed to strong seasonal fluctuations. In the dry savannahs of West Africa, this mosquito population dynamics closely follows the pace of surface water availability: the species pullulates during the rainy season and is able to reproduce throughout the dry season in areas where permanent water bodies are available for breeding. The impact of such environmental fluctuation on mosquito development and the phenotypic quality of emerging adults has however not been addressed in details. Here, we examined and compared phenotypic changes in the duration of pre-imaginal development, body dry mass at emergence and wing size, shape and surface area in young adult females An. coluzzii originated from five distinct geographic locations when they are reared in two contrasting conditions mimicking those experienced by mosquitoes during the rainy season (RS) and at the onset of the dry season (ODS) in Burkina Faso (West Africa). Our results demonstrated strong phenotypic plasticity in all traits, with differences in the magnitude and direction of changes between RS and ODS depending upon the geographic origin, hence the genetic background of the mosquito populations. Highest heterogeneity within population was observed in Bama, where large irrigation schemes allow year-round mosquito breeding. Further studies are needed to explore the adaptive value of such phenotypic plasticity and its relevance for local adaptation in An. coluzzii.

An overview of malaria transmission from the perspective of Amazon Anopheles vectors

In the Americas, areas with a high risk of malaria transmission are mainly located in the Amazon Forest, which extends across nine countries. One keystone step to understanding the Plasmodium life cycle in Anopheles species from the Amazon Region is to obtain experimentally infected mosquito vectors. Several attempts to colonise Anopheles species have been conducted, but with only short-lived success or no success at all. In this review, we review the literature on malaria transmission from the perspective of its Amazon vectors. Currently, it is possible to develop experimental Plasmodium vivax infection of the colonised and field-captured vectors in laboratories located close to Amazonian endemic areas. We are also reviewing studies related to the immune response to P. vivax infection of Anopheles aquasalis, a coastal mosquito species. Finally, we discuss the importance of the modulation of Plasmodium infection by the vector microbiota and also consider the anopheline genomes. The establishment of experimental mosquito infections with Plasmodium falciparum, Plasmodium yoelii and Plasmodium berghei parasites that could provide interesting models for studying malaria in the Amazonian scenario is important. Understanding the molecular mechanisms involved in the development of the parasites in New World vectors is crucial in order to better determine the interaction process and vectorial competence.

Antiviral immunity of Anopheles gambiae is highly compartmentalized, with distinct roles for RNA interference and gut microbiota

Arboviruses are transmitted by mosquitoes and other arthropods to humans and animals. The risk associated with these viruses is increasing worldwide, including new emergence in Europe and the Americas. Anopheline mosquitoes are vectors of human malaria but are believed to transmit one known arbovirus, o'nyong-nyong virus, whereas Aedes mosquitoes transmit many. Anopheles interactions with viruses have been little studied, and the initial antiviral response in the midgut has not been examined. Here, we determine the antiviral immune pathways of the Anopheles gambiae midgut, the initial site of viral infection after an infective blood meal. We compare them with the responses of the post-midgut systemic compartment, which is the site of the subsequent disseminated viral infection. Normal viral infection of the midgut requires bacterial flora and is inhibited by the activities of immune deficiency (Imd), JAK/STAT, and Leu-rich repeat immune factors. We show that the exogenous siRNA pathway, thought of as the canonical mosquito antiviral pathway, plays no detectable role in antiviral defense in the midgut but only protects later in the systemic compartment. These results alter the prevailing antiviral paradigm by describing distinct protective mechanisms in different body compartments and infection stages. Importantly, the presence of the midgut bacterial flora is required for full viral infectivity to Anopheles, in contrast to malaria infection, where the presence of the midgut bacterial flora is required for protection against infection. Thus, the enteric flora controls a reciprocal protection tradeoff in the vector for resistance to different human pathogens.

Experimental study of the relationship between Plasmodium gametocyte density and infection success in mosquitoes; implications for the evaluation of malaria transmission-reducing interventions

The evaluation of transmission reducing interventions (TRI) to control malaria widely uses membrane feeding assays. In such assays, the intensity of Plasmodium infection in the vector might affect the measured efficacy of the candidates to block transmission. Gametocyte density in the host blood is a determinant of the infection success in the mosquito, however, uncertain estimates of parasite densities and intrinsic characteristics of the infected blood can induce variability. To reduce this variation, a feasible method is to dilute infectious blood samples. We describe the effect of diluting samples of Plasmodium-containing blood samples to allow accurate relative measures of gametocyte densities and their impact on mosquito infectivity and TRI efficacy. Natural Plasmodium falciparum samples were diluted to generate a wide range of parasite densities, and fed to Anopheles coluzzii mosquitoes. This was compared with parallel dilutions conducted on Plasmodium berghei infections. We examined how blood dilution influences the observed blocking activity of anti-Pbs28 monoclonal antibody using the P. berghei/Anopheles stephensi system. In the natural species combination P. falciparum/An. coluzzii, blood dilution using heat-inactivated, infected blood as diluents, revealed positive near linear relationships, between gametocyte densities and oocyst loads in the range tested. A similar relationship was observed in the P. berghei/An. stephensi system when using a similar dilution method. In contrast, diluting infected mice blood with fresh uninfected blood dramatically increases the infectiousness. This suggests that highly infected mice blood contains inhibitory factors or reduced blood moieties, which impede infection and may in turn, lead to misinterpretation when comparing individual TRI evaluation assays. In the lab system, the transmission blocking activity of an antibody specific for Pbs28 was confirmed to be density-dependent. This highlights the need to carefully interpret evaluations of TRI candidates, regarding gametocyte densities in the P. berghei/An. stephensi system.