Transcriptome profiles of Anopheles gambiae harboring natural low-level Plasmodium infection reveal adaptive advantages for the mosquito

Anopheline mosquitoes are the sole vectors for the Plasmodium pathogens responsible for malaria, which is among the oldest and most devastating of human diseases. The continuing global impact of malaria reflects the evolutionary success of a complex vector-pathogen relationship that accordingly has been the long-term focus of both debate and study. An open question in the biology of malaria transmission is the impact of naturally occurring low-level Plasmodium infections of the vector on the mosquito’s health and longevity as well as critical behaviors such as host-preference/seeking. To begin to answer this, we have completed a comparative RNAseq-based transcriptome profile study examining the effect of biologically salient, salivary gland transmission-stage Plasmodium infection on the molecular physiology of Anopheles gambiae s.s. head, sensory appendages, and salivary glands. When compared with their uninfected counterparts, Plasmodium infected mosquitoes exhibit increased transcript abundance of genes associated with olfactory acuity as well as a range of synergistic processes that align with increased fitness based on both anti-aging and reproductive advantages. Taken together, these data argue against the long-held paradigm that malaria infection is pathogenic for anophelines and, instead suggests there are biological and evolutionary advantages for the mosquito that drive the preservation of its high vectorial capacity.

The vast majority of anopheline mosquitoes that ingest Plasmodium gametocyte-infected blood meals never become infectious 15 . In controlled laboratory studies, only 5-10% of ingested Plasmodium gametocytes successfully differentiate into ookinetes; of those, less than 0.3% form oocysts, giving rise to approximately 1-2 oocysts per mosquito midgut 15,16 . This closely mirrors West African field studies indicating that the majority of malaria-infected mosquitoes are within that range and only rarely are trapped mosquitoes observed with more than 5 oocysts per midgut 17 . In that light, it is reasonable to suggest that when investigating salient biological relationships between Anopheles and Plasmodium in the laboratory, the relative intensity of mosquito infections should be low to better align with these natural levels. Paradoxically, the extremely high malaria infection intensities that are often employed in laboratory studies to provide robust levels of infection prevalence are likely to significantly alter the Anopheles-Plasmodium paradigm potentially giving rise to biologically irrelevant shifts in mosquito gene expression, physiology, and behavior.
In order to assess the impact of malaria infection on chemosensory processes of An. gambiae we have deliberately employed low gametocytemia and, thereby, biologically relevant P. falciparum blood meals to examine how salivary-gland-stage sporozoite infection influences the transcriptome of female An. gambiae heads (including sensory appendages) and salivary glands, henceforth referred to as 'head' . These data suggest that natural levels of malaria infection align with enhanced olfactory sensitivity that likely impacts chemosensory driven host preference. Moreover, and surprisingly, natural levels of Plasmodium infection also appear to provide a broad and synergistically beneficial modulation of mosquito fitness that underlies and may indeed promote the evolutionary stability of the Anopheles-Plasmodium paradigm that continues to drive global malaria transmission.

Results and discussion
Experimental design and establishment of optimal gametocytemia. Our study was designed to rigorously compare the impact of biologically salient P. falciparum infections on the head (encompassing brain, sensory appendages, and salivary glands) transcriptome of a laboratory strain of the Afro-tropical malaria vector mosquito An. gambiae s.s. While every effort has been made to rigorously examine the impact of naturally occurring P. falciparum infections on vector mosquitoes, it must be noted there are likely to be inherent differences between established laboratory colonies and diverse field populations of both pathogen and vectors. We nevertheless believe the broad implications of this study remain salient to understanding the natural biology of malaria infected anophelines.
An experimental paradigm incorporating triplicate biological replicates (21 An. gambiae heads per replicate) was designed wherein 5-day-old females were blood fed on either P. falciparum-infected (treatment) or naïve (control) human blood meals and thereafter reared for 18 days prior to collection and tissue dissection for downstream molecular analyses. This allows for completion of the mosquito sporogonic cycle as well as sporozoite invasion of the salivary glands in treated (infected) mosquitoes (Fig. 1A), thus ensuring that mosquitoes harboring sporozoites in their salivary glands are examined. Despite of the experimental challenges resulting from significant reductions in prevalence we prioritized utilizing low-level intensities representative of natural biological infections 17 . To achieve this goal, approximately 14% of the treatment (infected) mosquitoes were sampled 8 days after blood meals, and their midguts examined for the presence of P. falciparum oocysts following established protocols 18 . These studies indicated that mosquitoes that had blood fed on human blood supplemented with P.   RNAseq/Bioinformatics. Transcriptome profiling was conducted using total RNA extracted from triplicate An. gambiae treatment (infected; T1-3) and control (uninfected; C1-3) heads. Following sequence trimming and quality control bioinformatic pipelines, individual data sets were aligned to the An. gambiae PEST strain genome. Sequencing of the T1, T2, and T3 P. falciparum-infected (treatment) An. gambiae head libraries generated a total of over 254 million reads mapping to 83%, 83%, and 85% respectively of the An. gambiae genome (VectorBase-53_AgambiaePEST; Table 2). Likewise, Illumina sequencing of the C1, C2, and C3 control (uninfected) An. gambiae head libraries generated a total of 426 million reads mapping to 81%, 80%, and 81% respectively of the An. gambiae genome (VectorBase-53_AgambiaePEST; Table 2). Uniquely mapped reads were counted using the program HTSeq. DESeq2 statistical analysis identified 5765 annotated transcripts present in significantly different abundances between treatment and control replicates (p < 0.05; Fig. 3A). These annotated transcripts represent a transcriptional shift of 36% of the An. gambiae head transcriptome. Of these, 2776 transcripts were present in significantly higher abundances in treatment replicates and the remaining 2989 transcripts present in significantly higher abundance in control replicates (p < 0.05). With such a large transcriptional shift observed between treatment and control An. gambiae heads, both principal component analysis (PCA) and Pearson's correlation coefficient (PCC) calculations were performed to verify validity. PCC calculations for the three treatment and three control An. gambiae head libraries returned r values of 0.9429 and 0.9770, respectively, indicating a strong association among biological replicates (Table 2). Additionally, PCA determined that 83% of observed variance between the six library datasets was attributed to the presence/absence of P. falciparum infection, with negligible variance attributed to biological replication (Fig. 3B). GO annotation using OmicsBox assigned 4976 GO-terms (Level II) to 2832 transcripts of the 5765 differentially abundant transcripts (control and treatment libraries). Of the total annotated transcripts, 1291 were elevated in treatment libraries and 1541 transcripts elevated in control libraries, with approximately 10% of the    4). Transcripts with no clear molecular function via GO were further BLASTed using Blast2GO and the GeneBank non-redundant protein database to assist in further functional analysis. The resultant assigned GO annotations are representative of a diverse array of biological processes and molecular functions that make interpretation difficult. Taken together, these analyses reveal a broad transcriptome profile shift that spans a wide range of cellular processes, some of which are likely to impact the behavior, physiology and ultimately the vectorial capacity of An. gambiae carrying biologically relevant P. falciparum sporozoite infections.
Plasmodium manipulation of mosquito chemosensory mechanisms. To examine the potential for behavioral modification of infected mosquitoes, our initial focus was to examine transcriptome modulations directly relevant to the chemosensory systems of An. gambiae. These systems collectively play an essential role in establishing and maintaining critically important behaviors, including the blood-meal host preference of female mosquitoes 21,22 .
Odorant receptors. The transcript abundance of several An. gambiae odorant receptors (AgORs) were highly differentiated between P. falciparum sporozoite-infected (sporozoite; treatment) and uninfected (control) heads. Indeed, 60 of the 71 AgORs and the OR co-receptor (Orco) that were detected in sequencing libraries were consistently more abundant in heads of infected mosquitoes, with 19 displaying statistically significant induction (Fig. 5A). Transcripts encoding the remaining 11 AgORs (AgOR1, AgOR3, AgOR4, AgOR26, AgOR36, AgOR45, AgOR46, AgOR52, AgOR60, AgOR62 and AgOR72) showed an insignificant trend of less abundance in samples from infected mosquitoes (Fig. 5A). With the near complete representation of AgORs (95% of annotated An. gambiae ORs), it is evident that sporozoite infection significantly impacts the levels of these ORs within mosquito chemosensory structures. The majority of anopheline ORs are functionally characterized, and we have previously integrated OR functional RNAseq data to model the receptivity profile for the antennae of An. gambiae following their initial bloodmeal 23 . Transcriptome profiling studies together with behavioral analyses support the hypothesis that the sensitivity of the anopheline olfactory system is directly associated with the relative abundance of olfactory transcripts [23][24][25] . Applying this approach here suggests that the large number of significantly abundant AgOR transcripts associated with infected An. gambiae indicates substantial increases in odorant receptivity (Fig. 6). In fact, AgOR receptivity modelling identified 70 odorants, many of which are components of human sweat, with significantly increased receptivity in Plasmodium-infected mosquitoes than in uninfected controls. Of these, several [including L-lactic-acid, 1-octen-3-ol and 4-methylphenol] are human sweat components that have been shown to be attractive kairomones for host-seeking female anophelines 26,27 . Notably, lactic acid showed a 217% increase in AgOR receptivity in treatments. Furthermore, the significant increases in receptivity to 1-butanol (247%), 3-methyl-1-butanol (260%), and 2,3-butanedione (252%) align with laboratory olfactometry and semi-field studies, indicating that these compounds significantly increase the number of anopheline collected (vs. catches) when added to MMX traps baited with CO 2 and a basic blend of amine/carboxylic acid attractants 28 . These putative increases of olfactory acuity in Plasmodium-infected anopheline mosquitoes are consistent with behavioral bioassay and surveillance studies that determined Plasmodium infections to increase mosquito bite frequency while maintaining host-preference plasticity 29,30 . That said, our receptivity modelling does not align with recent studies that examined peripheral olfactory responses in P. falciparum-infected An. gambiae 10 . It is likely that the significantly higher infection intensities used in that study over-challenged An. gambiae infection-responsive systems, such as the innate immune system, resulting in unintended fitness consequences that  www.nature.com/scientificreports/ altered host-seeking and blood-feeding behaviors 31 . Because our studies avoid impacts arising from non-natural infection intensities, they can more confidently reveal a clear effect of biologically relevant malaria infections on the anopheline olfactory system. The effect of increasing host seeking and host feeding synergistically support both parasite and mosquito survivorship as well as maximizing vectorial capacity. Interestingly, the presence of polycistronic OR genes in Anopheles 32 raises the hypothesis that Plasmodium parasites may have evolved to specifically activate the transcription factors regulating polycistronic ORs to efficiently induce large-scale olfactory modulation. Indeed, all of the AgORs identified in polycistronic cluster 1B 32 are present in our sequencing libraries, and one of these AgORs17 was significantly over-represented in malarial mosquitoes; the increased abundance of the remaining five AgORs (13,15,16,47,55) are just below statistical significance (Fig. 5A).
Odorant binding proteins/chemosensory proteins. Odorant binding proteins (OBPs), chemosensory proteins (CSPs) and sensory appendage proteins (SAPs) have been postulated to mediate the solubility and transport of odorant molecules within the sensilla lymph 33 . A total of 57 AgOBPs/AgCSPs/AgSAPs were identified in sequencing datasets, with 47 AgOBPs/AgCSPs more abundant in heads of infected mosquitoes (17 significantly) while two AgOBPs (9 and 25) and a single AgCSP5 were significantly under-represented relative to uninfected controls (Fig. 5B). These findings are consistent with microarray-based transcriptome studies in Ae. aegypti that determined dengue viral infection to increase abundance of a subset of OBPs that was correlated to increased host feeding and vectorial capacity 34 . It is reasonable to postulate that Plasmodium infection similarly impacts anopheline mosquitoes to potentially increase host feeding, synergistically supporting parasite and mosquito survivorship as well as maximizing vectorial capacity. AgCSP5 (− 0.72-fold; p-adj = 0.014; AGAP029127), AgOBP9 (− 0.54-fold; p-adj = 0.030; AGAP000278), and AgOBP25 (− 0.40-fold; p-adj = 0.018; AGAP012320) are the only significantly under-represented OBP/CSP transcripts (Fig. 5B). In any case, the precise role of AgCSP5, AgOBP9, and AgOBP25 along with those of other OBP/CSP genes in An. gambiae olfaction remains enigmatic.   23,37,38 . In contrast, AgGR22 (AGAP00999), AgGR23 (AGAP003098), and AgGR24 (AGAP001915), which make up a complex that is required for volatile carbon dioxide (CO 2 ) detection critical for blood-feeding behaviors, have been well characterized 39 . The differential transcriptome profiles for AgGR22 (0.28-fold), AgGR23 (0.31-fold), and AgGR24 (0.34-fold) are only slightly more abundant in sporozoite-infected head transcriptomes than in controls (Fig. 7A, red bar) and, not surprisingly, these changes fail to meet the rigorous significance thresholds employed here. This suggests that Plasmodium infection does not significantly impact CO 2 sensitivity and is consistent with previous electrophysiology-based studies that showed P. berghei infections had no effect on An. stephensi maxillary palp neuron sensitivity to CO 2 40 . Taken together, it appears likely that P. falciparum parasites do not manipulate early, long-range anopheline host-seeking behaviors, such as CO 2 detection 41 .
Beyond the CO 2 receptors, we observed a considerable shift in the AgGR transcriptome profile between heads of uninfected and those of sporozoite-infected mosquitoes. A total of 44 AgGRs were differentially detected; transcripts for 40 AgGRs showed a higher abundance and the remaining 4 were less abundant in heads of sporozoite-infected mosquitoes than in controls (Fig. 7A). Of these differentially abundant AgGRs, 13 were present at significantly higher levels and none at significantly lower levels in treatments. While the most significantly  www.nature.com/scientificreports/ of non-CO 2 responsive AgGRs, our data are consistent with previous studies highlighting the overabundance of sugar transporters in the salivary glands of Plasmodium-infected An. gambiae and their importance for sporozoite survival 13 . The remaining transcript, AgGR47 (2.60-fold; p-adj = 0.05; AGAP005514), found significantly abundant in treatments is a DmGR66a ortholog and putatively functions in bitter/caffeine detection 35 . It is reasonable to speculate that Plasmodium sporozoites in Anopheles mosquitoes upregulate the deterrent effects of such AgGR neuronal responses to prevent ingestion of detrimental bitter nutrient sources that have been shown to temporally inhibit mosquito sugar and water neurons and reduce nectar feeding of An. quadrimaculatus 36 . Shifting mosquito gustatory profiles to increase sensitivity to bitter and behaviorally aversive semiochemicals therefore increases the likelihood of blood-meal (and or nectar) feeding and consequently pathogen transmission. It is likely that this effect contributes to increased sensitivity and the marked decrease in mortality rate of An. gambiae harboring P. falciparum sporozoites when exposed to DEET (N, N-diéthyl-3-méthylbenzamide), a widely used repellent, in comparison with uninfected An. gambiae and other pyrethroid test substrates 44,45 . Lastly, the increased sensitivity of AgGR47 observed in sporozoite-infected An. gambiae relative to naïve controls may reflect age-related nutrient deficiencies in controls. In aging Drosophila, nutrient deprivation dramatically alters feeding behaviors, including the depotentiation of responses to bitter 46 .
Ionotropic receptors. Ionotropic receptors (IRs) are an ancient family of chemosensory receptors evolutionarily derived from ionotropic glutamate receptors 47 that are associated with olfactory responses to amines and carboxylic acids that have been identified as important semiochemicals closely associated with host-preference and seeking 21 . As is true in other insects, AgIR8a, AgIR25a and AgIR76b likely function as co-receptors along with odorant-recognizing AgIRs that 'tune' responses to host cues 48,49 . In An. gambiae, AgIR25a together with AgIR76b and AgIR8a function in the detection of amines and carboxylic acids, respectively 49 both of these chemical classes have been directly implicated in anopheline host seeking 22,27 . Interestingly, the abundance of AgIR76b transcripts was significantly enhanced (0.55-fold; p-adj = 0.05; AGAP011968) in infected mosquito heads while the levels of AgIR25a and AgIR8a were not significantly different (Fig. 7B). Transcripts for 38 of the 49 odorant-tuning AgIRs detected here were found to be more abundant in sporozoite-infected head transcriptomes; only 11 were downregulated compared with uninfected controls (Fig. 7B). Of those upregulated, 18 were significantly more abundant but none were significantly less abundant in treatment groups than controls. Inasmuch as many AgIRs are uncharacterized, including those with the most pronounced and significant differential abundance (AgIR149, AgIR168, and AgIR203), the implications of these differentially up-/downregulated IRs are difficult to assess. Nevertheless, functional information for a small subset of AgIRs is available and others can be derived from analyses of clear orthologs in Ae. aegypti, An. gambiae, An. sinensis and Drosophila melanogaster. This allows us to identify several AgIRs significantly more abundant in treatment groups (Fig. 7B) to be associated with both gustation (AgIR7x (3.11-fold; p-adj = 0.03; AGAP013520), AgIR93a (0.52-fold; p-adj = 0.05; AGAP000256) 50 , AgIR60a (2.48-fold; p-adj = 0.05; AGAP011943) 51 ) and host-seeking behaviors (AgIR75g (1.96fold; p-adj = 0.003; AGAP013085) 50 ). In addition, we observed three receptors localized to the antenna (AgIR31a (1.43-fold; p-adj = 0.03; AGAP009014), AgIR68a (2.00-fold; p-adj = 0.003; AGAP007951), and AgIR41t1 (1.44fold; p-adj = 0.03; AGAP004432) 52 ). While persistence of host seeking and blood feeding might well be advantageous for Plasmodium, it would only be beneficial for their anopheline hosts if these behavioral shifts promoted the acquisition of a successful blood meal. Transient receptor potential channels. As is true for all insects, An. gambiae transient receptor potential (AgTRP) family members are associated with a variety of sensory modalities such as chemosensation, gustation, thermosensation, hygrosensation, mechanosensation, vision, and intracellular signaling 21 . A total of nine AgTRPs were identified in sequencing libraries, five of which are more abundant (three significantly) in heads of sporozoite-infected mosquitoes while four AgTRPs (three significantly) are under-represented relative to uninfected controls (Fig. 7C). Of these, AgTRPC (canonical; 0.50-fold; p-adj = 0.0002; AGAP000349), AgTRPC (canonical; 0.67-fold; p-adj = 0.0004; AGAP008435), and AgTRPM (melastatin; 0.75-fold; p-adj = 0.01; AGAP006825) are all active in thermosensation and feeding. TRPC has been identified as having several physiological roles in insects, including signal transduction, phototransduction, and gustatory sensitivity to CO 2 53 . TRPM mediates Drosophila nociception behaviors and putatively critical behaviors for An. gambiae viability and survival 54 . Given the critical roles that AgTRPC and AgTRPM are likely to play in mosquito host seeking and blood feeding, manipulation of these transcripts would likely be advantageous for Plasmodium to ensure maximum host contact rate and transmission. TRP-Painless and TRP-Pyrexia function in thermosensation and have been documented to mediate avoidance of noxious high temperatures in D. melanogaster 55 . It is possible the significantly reduced abundance of AgTRP-Painless (− 0.78-fold; p-adj = 0.01; AGAP001243) and AgTRP-Pyrexia (− 0.75-fold; p-adj = 0.0003; AGAP010269) in treatment groups compared with naïve controls is related to a loss of nociception.
The remaining An. gambiae TRP identified as significantly differentially abundant between treatments and controls functions as a mechanoreceptor with roles in Drosophila hygrosensation. AgTRP-Waterwitch which is significantly less abundant in treatment tissues (− 0.70-fold; p-adj = 0.001; AGAP000361) is homologous to the DmTRP-Waterwitch that functions in dry-air detection 56 . Moisture detection, or importantly the lack thereof, mediates flight and host-seeking behavior in anopheline mosquitoes such that augmenting mosquito traps with moist air significantly increases flight activity and attraction of wild Aedes, Anopheles and Culex mosquitoes 57 . It is likely that upregulation of chemo/thermosensory TRPs as well as the inverse regulation of Ag-Waterwitch optimizes persistent host seeking and blood feeding and is consistent with studies showing that Anopheles mosquitoes harboring Plasmodium sporozoites exhibit more persistent host-seeking behavior than uninfected controls 12 . www.nature.com/scientificreports/ Pickpocket channel proteins/sensory neuron membrane proteins. We identified 18 An. gambiae pickpocket channel proteins (AgPPKs) of which transcripts were predominately more abundant in heads of sporozoiteinfected mosquitoes (Fig. 7D). Two of these AgPPKs, AgPPK13-like (0.97-fold; p-adj = 0.004; AGAP007945) and AgPPK31-like (1.47-fold; p-adj = 0.05; AGAP000657), showed significant enrichment. In D. melanogaster, PPKs have been implicated to function in salt and water taste, although the behavioral implications of these characterized PPKs is unclear 58 . Similarly, while aedine PPKs are hypothesized to function in a variety of sensory modalities, such as chemosensation, hygrosensation and mechanosensation 37 , anopheline PPKs remain largely uncharacterized and the implications of these differentially abundant PPKs remain unclear. Sensory neuron membrane proteins (SNMPs) are both colocalized with ORs and broadly expressed throughout D. melanogaster and are postulated to have both chemosensory and non-chemosensory functions 59 with particular roles in dipteran and lepidopteran pheromone detection 60 . Both AgSNMP1 and AgSNMP were identified in sequencing libraries, although without significance (Fig. 7E). These and other anopheline SNMPs remain uncharacterized, and the implications of differentially abundant AgSNMPs are therefore unclear.
Similarly, transcripts for NMDA and non-NMDA types of anopheline glutamate receptors were identified in sequencing libraries and predominately under-represented in heads of infected An. gambiae. Of the eight www.nature.com/scientificreports/ glutamate receptor transcripts identified, none was significantly abundant in heads of sporozoite-infected mosquitoes; in contrast, six NMDA and non-NMDA glutamate receptors-AgGluRIIa (− 1.27-fold; p-adj = < 0.0001; AGAP000803), AgGluRIId (− 0.54-fold; p-adj = 0.003; AGAP002797) AgGluRIIe (− 0.81-fold; p-adj = < 0.0001; AGAP012447), AgNMDAR1 (− 1.23-fold; p-adj = < 0.0001; AGAP001478), AgNMDAR2 (− 0.32-fold; p-adj = 0.005; AGAP012429) and AgNMDAR3 (− 0.98-fold; p-adj = 0.002; AGAP005527)-were significantly diminished in sporozoite-infected mosquitoes compared with naïve controls (Fig. 8C). It is possible that in light of the ability of these receptors to function as immunomodulators of inflammatory responses 64 , the overabundance of NMDA and non-NMDA glutamate receptors in naïve mosquitoes reflects either an aging-related immune system dysfunction or neuronal pathology. Transcripts encoding the receptors for several peptide biogenic amines, including dopamine, histamine, octopamine, serotonin and tyramine, which function as anopheline neuromodulators and neurotransmitters were identified in our libraries. With the exception of two pairs of histamine and octopamine receptor transcripts that were similarly up-/downregulated, most biogenic amine receptors showed no significant differential abundance in heads of sporozoite-infected An. gambiae relative to naïve controls (Fig. 8D-G). Transcripts encoding several biogenic amine processors and transporters were identified in similarly up-/ downregulated pairs with no clear distinction between treatment and naïve controls (Fig. 9A, D). Insect neuropeptides are essential neuromodulators that regulate a variety of physiological and behavioral functions. Several neuropeptide transcripts identified in sequencing libraries were shown to be significantly overabundant in the heads of sporozoite-infected An. gambiae (Fig. 9B). Interestingly, there were also four neuropeptide transcripts and their correlating receptors that were significantly diminished in treatment heads compared with naïve controls, allatostatin 1 (AgAs1; − 0.50-fold; p-adj = 0.001; AGAP003712), allatostatin 2 (AgAs2; − 0.41-fold; p-adj = 0.01; AGAP010157), FMRFamide (AgFMRF, − 1.34-fold; p-adj = < 0.0001; AGAP005518) and sickie (AgSickie; − 1.66-fold; p-adj = < 0.0001; AGAP009424), that are directly or indirectly correlated with insect stress and immune response (Fig. 9B, C). In Ae. Aegypti, allatostatin interacts with hemocytes to increase humoral immune responses to combat cellular stress and inflammation 65 . FMRFamide signaling in Drosophila increases in response to cellular stress and inflammation to promote stress-induced sleep and cellular recovery to help bolster immune function 66 . Sickie acts directly on Drosophila innate immunity by activating the transcription factor Relish and subsequent transcription of immune defense antimicrobial peptides (AMPs) 61 . Neuropeptide F (NPF) also acts indirectly to modulate stress responses through activation of downstream signaling pathways. In dipterans, NPF independently activates insulin signaling pathways and coordinately the target or rapamycin (TOR) pathway 67 . TOR pathway is an insect nutrient-sensing pathway that regulates cellular metabolism and homeostasis in response to nutrient flux and stress 68,69 . TOR signaling naturally increases with age and increasing cellular dysfunction and stress 70 . With a substantial number of neuropeptides/hormones and their correlating receptors significantly abundant in heads of sporozoite-infected An. gambiae and functioning to maintain healthy mosquito cells and tissues, the decreased abundance of these select four neuropeptides in infected mosquitoes is puzzling. If the "normal' transcriptome profile of naïve mosquitoes reflects changes associated with age-related immune system dysfunction and the loss of cellular homeostasis, or some other age-related pathology, then one could hypothesize that infected mosquitoes might be refractory to those effects.
Immune system. Activation of the immune system in response to infection is energetically costly and likely results in decreased fitness reflected by other life history traits with the potential to impact vectorial capacity. Mosquito anti-Plasmodium immunity is stage specific and involves diverse defense processes, although the extent of immune defense in salivary-gland-stage infected mosquitoes, and especially the head and salivary glands (post sporozoite salivary gland invasion), remains unclear 15 . Since Plasmodium, unlike arboviruses, does not establish infection in the mosquito head, the immunity-related transcriptome changes observed in our study are likely to mostly derive from the combined head (including sensory appendages) and salivary gland transcriptomes. In heads of sporozoite-infected and control sequencing libraries, several components of An. gambiae humoral and cellular, innate immune defenses were identified. An. gambiae Pattern Recognition Receptors (PRRs) include 6 fibrinogen proteins (FBGs), 7 3-glucan binding proteins (3GNBPs), 7 peptidoglycan recognition proteins (PGRPs), 24 leucine-rich repeat immune proteins (LRIMs) including Plasmodium infection responsive leucine-rich proteins, 13 thioester-containing proteins (TEPs), and 2 members of C-type lectin family proteins (CTLs)-AgCTL4 and AgCTLM2 (Fig. 10A, B). A large portion of the PRR transcripts were not present at significantly differential abundances between heads of sporozoite-infected and control mosquitoes; indeed, neither of the transcripts encoding AgCTLM2 or AgCTL4 displayed significant differential abundance.
Anopheline TEPs are antimicrobial and antiparasitic PRRs complement-like factors that drive pathogen clearance through lysis and melanization 71 . Of those, only AgTEP13 (0.34-fold; p-adj = 0.03; AGAP008407) was significantly more abundant in heads of sporozoite-infected mosquitoes, and AgTEP1 (− 1.95; p-adj = < 0.0001; AGAP010815) and AgTEP12 (− 0.94; p-adj = 0.011; AGAP008654) significantly diminished in treatment groups compared with controls (Fig. 10A). While the roles of TEP12 and TEP13 in anopheline immune responses are not currently known, the well-studied TEP1 functions as an enhancer of mosquito phagocytosis 72 . Since TEP1 is incapable of colocalizing and killing Plasmodium sporozoites, its decreased abundance in treatment groups is not unexpected 71,73 . As such, the increased abundance of AgTEP1 in naïve controls likely reflects aging-related immune system dysfunction, as aging decreases phagocytic efficiency and contributes to immunosenescence 46 . LRIMs and APLs are PPRs that mediate innate immune surveillance and response. Transcripts for 20 LRIMs and APLs were downregulated in heads of infected mosquitoes relative to naïve controls, while 4 of the 24 An. gambiae LRIMs were more abundant in treatments. Of the 4 AgLRIMs upregulated, 2 were significantly more abundant, whereas 9 of the 20 LRIMs and APLs were significantly less abundant in treatments than controls (Fig. 10B). These data are largely consistent with previous studies on the transcriptional regulation of LRIMs  www.nature.com/scientificreports/ and APLs in response to Plasmodium, gram-negative/gram-positive bacteria, and fungal challenge 74,75 . While AgLRIM12 (0.72-fold) and AgLRIM5 (0.82-fold; ACOM043147) were both significantly more abundant in heads of infected mosquitoes they are not yet conclusively associated with Plasmodium infection 75,76 ; it is therefore possible the observed differential abundance of AgLRIMs reflect non-immunity responses or potential aging effects that are only partially impacted by P. falciparum infection. Similarly identified AgFBGs, which function in phagocytosis and complement activation, were predominately diminished in treatments relative to naïve controls, with AgFBG (− 1.47-fold; p-adj = 0.015; AGAP011231), AgFBG (− 1.98-fold; p-adj = 0.004; AGAP010869), AgFBG (− 1.29-fold; p-adj = 0.02; AGAP004917) and AgFBG (− 1.29-fold; p-adj = 0.007; AGAP006914) significantly reduced (Fig. 10B) 72 . In Drosophila, aging is associated with increased microbial loads and dysbiosis that trigger increased expression of immune genes and stronger immune activation 77 . Excessive immune activation ultimately results in immune dysfunction and decreased lifespan 78 . In that context, the large number of AgLRIMs, AgAPLs and AgFBGs that are more abundant in naïve control transcriptomes may reflect aging-related dysfunction of the An. gambiae immune system. Accelerated immune system dysfunction in naïve controls also rationalizes the observed AgPGRP transcript abundances identified in replicate sequencing libraries. Our data indicate that transcripts for two AgPGRPs that have been characterized as negative regulators of immune IMD signaling defense, AgPGRP-LB (1.46-fold; p-adj = 0.02; AGAP001212) and AgPGRP-LD (0.69-fold; p-adj = 0.05; AGAP005552: Fig. 10A), are significantly more abundant in heads of sporozoite-infected mosquitoes, suggesting a suppression of the IMD immune pathway 79 . The IMD signaling initiation receptor AgPGRP-LC (− 0.25-fold; p-adj = 0.02; AGAP005203), identified as significantly reduced in treatments compared with controls, further supports the presence of diminished IMD signaling in sporozoite-infected An. gambiae. PGRP immune regulator deficiency results in unregulated IMD signaling causing dysplasia, inflammation, apoptosis, and decreased lifespan 80 . Preservation of immune homeostasis would be advantageous for both Anopheles hosts and Plasmodium pathogens to maintain vectorial capacity. It is possible that An. gambiae challenged with biologically relevant, low-level P. falciparum results in a mild and prolonged engagement of the mosquito immune system that ultimately primes the system for protection from secondary infection and importantly, the negative impact of aging-related immune system degeneration 81 . Indeed, Plasmodium-induced anopheline immunological priming may represent an advantageous selective adaptive pressure that maximizes mosquito host fitness and, consequently, vectorial capacity.
Serine proteases constitute the majority of mosquito immune system moderators that are activated in response to pathogen recognition 72 . While there are hundreds of serine proteases functioning in a variety of mosquito molecular processes, eight CLIP-domain serine proteases (CLIPs)-SPCLIP1, CLIPA2, CLIPA7, CLIPA8, CLIPA14, CLIPA28, CLIPB14, and CLIPB15-mediate P. berghei melanization 72,82,83 . No AgCLIP transcripts were significantly more abundant in heads of sporozoite-infected mosquitoes, although eight were less abundant in treatment groups for which mRNAs of AgCLIPA7 (− 1.34-fold: p-adj = < 0.0001; AGAP011792), AgCLIPA14 (− 1.2-fold; p-adj = < 0.0001; AGAP011788) and AgCLIPA28 (− 0.92-fold: p-adj = 0.02; AGAP010968) were significantly under-represented relative to uninfected controls (Fig. 10C). These data were expected because anopheline CLIPs predominately moderate melanization of rodent Plasmodium ookinetes and are not associated with immune defenses targeting sporozoites or human Plasmodium 83 . Indeed, the significant abundance of AgCLIPs in uninfected middle-aged An. gambiae may result from aging-related immune dysfunction. The absence of elevated levels of AgCLIPs highlights another putative benefit of Plasmodium-mediated immune priming. Another way to moderate immune effector mechanisms such as melanization is performed by Serpins (serine protease inhibitors). Levels of two Serpins, AgSerpin6 (− 0.64-fold; p-adj = < 0.0001; AGAP009212) and AgSerpin10 (− 0.42; p-adj = 0.03; AGAP005246), that are upregulated in response to Plasmodium ookinetes 84,85 , were significantly under-represented compared with uninfected controls (Fig. 10C). As melanization is not expected to be active in An. gambiae carrying salivary-gland-stage P. falciparum infections, there would be no requirement for Serpin moderation, and such genes as Serpin6 and Serpin10 would no longer be required at this timepoint 86 . Likewise, putative immune dysfunction in uninfected controls may impair proper regulation of melanization and, as such, Serpin transcription.
A wide range of immune effectors are downstream of PRR-controlled pathways and other recognition pathways that include anti-microbial peptides (AMPs), reactive oxygen species (ROS), reactive nitrogen species (RNS), prophenoloxidase (PPOs), and apoptosis-and phagocytosis-related genes 72 . In Anopheles, five classes of AMPs including attacins, defensins, diptericins, cecropins, and gambicins have been described 87 . In our data set, three defensins (Def), three cecropins (Cec), one gambicin (Gam) and one attacin (Att) are present (Fig. 10E). Of these, only two defensins and Cec1 were significantly differentially abundant between heads of sporozoiteinfected and uninfected control mosquitoes, with AgDef2 (0.96-fold; p-adj = < 0.0001; AGAP004632), AgDef4 (1.96-fold; p-adj = 0.005; AGAP005416) and AgCec1 (1.78-fold; p-adj = 0.001; AGAP000693) all significantly abundant in treatment groups compared with controls. Interestingly, these shifts plus the significantly increased abundance of AgGNBPB2 (1.46-fold; p-adj = 0.02; AGAP001212; Fig. 10A), a PRR that co-localizes P. falciparum, are the only other indication of an An. gambiae immune response to P. falciparum infection in this study. Of these, only AgCec1 is associated with anti-sporozoite activity. AgDef1 and AgDef4 demonstrate antibacterial activity but have not shown any anti-Plasmodium activity 88 . Interestingly, studies in Ae. aegypti have determined that differential Def abundances can be observed without immune challenge, and often after pathogen clearance, suggesting alternative functions for the AMPs 89 . As such, the role of AgDefs in response to sporozoite infections remains unclear. ROS and RNS are vital effectors of anopheline immune defenses active against malaria parasites 72,90,91 . ROS formation in dipterans is postulated to involve NADPH oxidase (NOX)-based production of superoxide anions that are transformed into H 2 O 2 by superoxide dismutase (SOD) and subsequent production of ROS 90 . NOX is also likely required for heme peroxidase (HPX) and nitric oxide synthase (NOS) activation, both of which potentiate nitration and generations of RNS 91 10E). In that light, it is logical to expect the significant under-representation of transcripts encoding catalase, glutathione peroxidase and thioredoxin peroxidase detoxifying enzymes observed in heads of sporozoite-infected An. gambiae when compared with uninfected controls (Fig. 10E).
In Anopheles and other mosquitoes, CLIPs activate PPOs, the primary effectors mediating melanization 92 . In that light, the lack of differentially abundant AgCLIP transcripts aligns with the observation that AgPPOs were predominately under-represented in heads of sporozoite-infected mosquitoes compared with uninfected controls. AgPPO2 (− 1.70-fold; p-adj = 0.001; AGAP006258), AgPPO5 (− 1.30-fold; p-adj = < 0.0001; AGAM012616), AgPPO6 (− 1.38-fold; p-adj = 0.002; AGAP004977), and AgPPO9 (− 1.50-fold; p-adj = 0.006; AGAP004978) all significantly diminished in An. gambiae carrying salivary-gland-stage sporozoites (Fig. 10). The higher PPO levels in heads of uninfected mosquitoes is surprising and may reflect a higher background of metabolic stress. However, this pattern is consistent with the fact that melanized sporozoite-stage Plasmodium has never been observed. In Ae. aegypti, allatostatin neuropeptide interacts with hemocytes and may have subsequently increased AgPPO expression to address dysbiosis and inflammation resulting from IMD destabilization. This finding is particularly interesting because it correlates mosquito immune system fitness with neuropeptide production and has significantly broader implications for aging-related immune system dysfunction. It also provides a plethora of anopheline molecular mechanisms that putatively synergize with Plasmodium immune priming, thereby providing a significant advantage that ultimately facilitates malaria transmission.
PRRs activate the Toll pathway through cleavage of Spaetzle proteins (SPZs) to activate Dorsal/NK-kB transcription factor (Rel1) and regulate antimicrobial peptides 72 . The transcriptome profile of these components also supports an aging-related immune system deregulation associated with malarial treatments that is importantly not seen in uninfected controls. Transcripts encoding AgRel1 (− 0.44-fold; p-adj = 0.002; AGAM009515) along with one Spaeztle, AgSPZ6 (− 0.80-fold; p-adj = 0.002; AGAP005126), are both significantly under-represented in heads of sporozoite-infected mosquitoes (Fig. 11A). Inasmuch as the Toll pathway has not been linked to defenses against the human malaria parasite P. falciparum 18,95 , it is not surprising that the majority of Toll pathway genes were also significantly under-represented in treatments compared with uninfected controls (Fig. 11A). In aging Drosophila, Toll signaling becomes increasingly dysfunctional, causing inflammation, systemic organismal defects, and metabolic deficiencies that result in shortened lifespans 77 . In that light, decreased activation of Toll in the heads of sporozoite-infected mosquitoes is part of what we propose to be a broadly synergistic Anopheles-Plasmodium paradigm.
As is the case for Toll-mediated signaling, the anopheline IMD immune signaling pathway involves binding of the PRR, PGRP-LC, which triggers activation of NK-kB transcription factor (Rel2) and production of AMPs. Transcripts for all components of IMD signaling, including PGRP-LC, were identified in sequencing libraries, although AgPGRP-LC (− 0.26-fold; p-adj = 0.02; AGAP005203) was the sole transcript significantly diminished in heads of sporozoite-infected mosquitoes (Figs. 10A, 11B). Even the modest differential abundance of AgPGRP-LC, which is a positive regulator of IMD signaling and PGRP expression, provides a rationale for the lack of significantly differential IMD signaling pathway transcripts in sporozoite-infected mosquitoes (Fig. 11B) 72 . In any case, the destabilization of anopheline IMD signaling pathway components would putatively render the immune signaling system a poor target for Plasmodium immunological priming, as the long-lasting immune response would be unpredictable and difficult to adapt 81 .
The anopheline JAK-STAT signaling pathway has been implicated to function in a variety of dipteran developmental processes, as well as in antimicrobial, antiparasitic and antiviral immune defenses 96 . It is negatively regulated by PIAS and suppressor of cytokine signaling (SOCS) 97 and induced by the binding of Unpaired (Upd) to Domeless (Dome) to form a receptor dimer that, in turn, activates Janus Kinase (JAK) or its homolog Hopscotch (Hop) and ultimately STAT which triggers the production of genes regulated by the JAK-STAT pathway and  (Fig. 11B).
As is the case for IMD pathways, anopheline JAK-STAT signaling primarily targets pre-sporozoite stages of Plasmodium, making the lack of JAK-STAT activation in the salivary-gland-stage infected An. gambiae examined here expected 80 . This is particularly evident by the significant abundance of JAK-STAT's primary negative regulator AgSOC in heads of sporozoite-infected mosquitoes. Dipteran SOCs deactivate the Hop-Dome dimer, block STAT activation, and target components of the Hop-Dome dimer complex for degradation 98 . As such, the increased abundance of AgUPD3 in treatment groups may reflect replenishment activity to restore UPD3 levels following P. falciparum midgut invasion or simply maintenance of normal UPD3 equilibrium. In Drosophila, immune system deregulation and mitochondrial degeneration significantly increase JAK-STAT expression in response to growing systemic inflammation and oxidative stress associated with aging. Excessive JAK-STAT activation ultimately leads to a loss of tissue homeostasis and cell death 77 . The active inhibition of JAK-STAT signaling in infected mosquitoes would explain the significantly diminished abundance of AgSTAT2 and AgHop and may also reflect aging-related dysfunction in uninfected An. gambiae controls. As such, the anopheline JAK-STAT signaling pathway may represent a target for Plasmodium priming, reducing the detrimental effects of agingrelated immune deregulation thereby synergistically benefiting both parasite and host.
General metabolism and aging. Metazoan aging can be characterized by the loss of cellular and tissue homeostasis, deterioration of vital biological processes with negative impacts on signaling pathways that regulate nutrient sensing, metabolic allocation, and the maintenance of healthy cells and tissues that together promote mortality 46 . The insect nutrient-sensing pathway TOR regulates cellular metabolism and growth in response to nutrient flux 69 . In Drosophila, increased TOR signaling is associated with aging, and the inhibition of TOR significantly extends fly lifespan and overall viability 99 . Both anopheline TOR transcripts-AgTOR (− 0.44-fold;  (Fig. 12A). Additionally, a transcript for the TOR kinase inhibitor tuberous sclerosis complex 2 (AgTSC2; − 0.46-fold; p-adj = < 0.001; AGAP003445), which in Drosophila inhibits TOR signaling and acts on fat bodies as a mechanism of lifespan expansion, was found to be significantly diminished in treatment groups compared with controls (Fig. 12A) 99 . The AMP-activated kinase (AMPK) signaling pathway maintains insect energy homeostasis by regulating ATP consumption and production; it is activated in response to low intracellular AMP levels, prompting ATP re-accumulation by reducing costly anaerobic processes as well as activating phosphofructokinase (Pfk), thereby increasing glycogenolysis 100 . Interestingly, the transcripts for both AgAMPK (− 0.75-fold; p-adj = < 0.0001; AGAP002686) and AgPfk (− 0.33-fold; p-adj = 0.05; AGAP007642) were significantly diminished in heads of sporozoite-infected An. gambiae compared with non-infected mosquitoes (Fig. 12B). In contrast, AMPK overabundance in An. stephensi is associated with increased innate immune responses and significant alterations in mitochondrial biogenesis 101 . In aging D. melanogaster, increased levels of AMPK improve tissue homeostasis via an anti-inflammatory mechanism to promote longevity 102 . It is possible that, compared with treatment transcriptomes, the age-related increased immune response observed in heads of naïve An. gambiae results in reduced nutrient and energy homeostasis. This hypothesis is also supported by the significantly diminished abundance of hepatocyte nuclear factor 3 (AgHNF3; − 2.17-fold; p-adj = < 0.0001; AGAP001671) and hepatocyte nuclear factor 4 (AgHNF4; − 0.98-fold; p-adj = 0.005; AGAP002155) in heads of sporozoite-infected An. gambiae mosquitoes www.nature.com/scientificreports/ compared with naïve controls. In Ae. Aegypti, RNAi depletion of HNF4 was shown to result in significant downregulation of transcripts encoding triacylglycerol catabolism, which is essential for fat body metabolism 103 . The head transcriptomes of naïve An. gambiae controls provide substantial evidence of a robust endemic sterile inflammatory response as well as decreased mitochondrial function, increased mitochondrial damage and mitophagy. In D. melanogaster, mitochondrial dysfunction is closely associated with decreased mitochondrial transcript levels of genes such as cytochrome oxidase C (COX) and cytochrome B1 (CB) 104 . Indeed, in heads of sporozoite-infected An. gambiae, transcripts for COX and CB were present at very high levels, significantly more abundant than in uninfected controls (Fig. 12C). In general, naïve An. gambiae displayed a significant overabundance of transcripts associated with ubiquitination and autophagy, indicating a higher degree of cell dysfunction than in the treatment group (Fig. 12F). For example, PTEN-induced kinase 1 (PINK1), a serine/threonine kinase that specifically targets damaged mitochondria for autophagy, is significantly diminished in Plasmodium-infected An. gambiae compared with non-infectious mosquitoes (− 0.66-fold; p-adj = 0.0004; AGAP004315) ( Fig. 12B), which may be an indication of a higher degree of mitochondrial dysfunction in naïve aged mosquitoes 105 .
In insects, the forkhead box O (FOXO) transcription factor signaling pathway is an important moderator of protein quality control, stress resistance, and lifespan extension pathways 46 . FOXO functions with its transcriptional target translation initiation factor, 4E-binding protein (4E-BP), to moderate proteostasis in Drosophila by removing damaged/abnormal proteins 106 . While Ag4E-BP was not significantly different, AgFOXO (− 2.36fold; p-adj = < 0.0001; AGAP000662) was significantly diminished in heads of sporozoite-infected An. gambiae compared with naïve controls (Fig. 12D). Similarly, the heat-shock transcription factor (AgHsTF; − 1.00-fold; p-adj = < 0.0001; AGAP029908) and numerous subsequently induced heat shock protein transcripts (Hsps) were all significantly diminished in treatment groups compared with naïve control An. gambiae (Fig. 12E). In dipterans, protein abnormalities resulting from oxidative stress gradually increase with age and are associated with an increase in HsTF and Hsp production to combat injury from abnormal and malformed proteins 107 . Lastly, consistent with these transcriptome shifts, our analyses have also uncovered several transcripts involved in ubiquitin-mediated proteosome protein degradation pathways, the abundance of which are significantly diminished in sporozoite infected mosquitoes (Fig. 12F). When viewed collectively, it is likely these differentially abundant transcripts are indicative of an accumulation of abnormal proteins with detrimental effects in naïve, uninfected mosquitoes that importantly does not occur in age-matched counterparts that harbor salivary-glandstage P. falciparum sporozoites after low intensity, biologically relevant blood-meal infections.

Conclusions
This study is one of the first examinations of the effects of biologically relevant Plasmodium sporozoite-stage infections on anopheline mosquitoes where the impact of Plasmodium infection on vector behavior and physiology significantly drives vectorial capacity. To gain insight into this process in as biologically relevant a manner as possible, we have carried out a set of comprehensive comparative analyses of the transcriptome profile of the head (including sensory appendages) and salivary gland tissues of female anophelines that have undergone low-intensity P. falciparum infections representative of those occurring naturally in disease-endemic regions. To our eyes, the significant reduction in prevalence and other experimental challenges associated with these studies are outweighed by the necessity to examine laboratory-based Plasmodium infections under the most biologically relevant conditions possible. As might be expected, these infections have broad effects on transcript abundance (as a partial proxy for differential gene expression) across a large number of gene families. Here, we have focused on those that impact neuronal function and chemosensation as well as age-related immunity and metabolic homeostasis.
The transcriptome profile shifts we have uncovered reveal a variety of systematic impacts that broadly suggest there are synergistic advantages for mosquitoes harboring salivary gland sporozoites (and therefore likely to be malarious) that align with significant increases in their vectorial capacity (Fig. 13). These notably include an expected increase in chemosensory sensitivity that could reasonably lead to more successful host seeking of infectious mosquitoes, along with a surprisingly distinctive transcriptome alteration that collectively aligns with an anti-aging paradigm that would likely provide a distinct advantage to these mosquitoes. When compared with similarly aged, uninfected mosquitoes, the synergistic advantage of naturally relevant malaria infections provides at least a partial rationale for the enduring persistence of Plasmodium pathogens and, indeed, global malaria.

Materials and methods
Ethics statement. All animal work was conducted in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health (NIH), USA. The protocols and procedures used in this study were approved by the Animal Care and Use Committee of the Johns Hopkins University (Permit Number: MO18H82) and the Johns Hopkins School of Public Health Ethics Committee. Commercial anonymous human blood was used for parasite cultures and mosquito feeding; thus, informed consent was not required, and the use of human blood for these types of experiments is exempt from human subject research.

Mosquito rearing and plasmodium infections.
An. gambiae s.s. (An. gambiae Keele strain) 108 mosquitoes were reared and maintained at 27 °C and 80% humidity with a 12-h light:12-h dark cycle at Johns Hopkins School of Public Health insectary (Baltimore, MD) according to standard procedures 18 . Newly emerged, 0-dayold female An. gambiae were split into nine independent groups (rearing cages; six infected or treatment, and three non-infected or control). All control and treatment groups were given their first and only bloodmeal at 5 days old 109  ) to give a final gametocytemia for each mosquito blood feeding. After feeding, treatment and control groups were independently sorted at 4 °C and unfed and partially fed mosquitoes were removed and disposed. The remaining fully engorged mosquitoes in treatment and control groups were incubated for 18 days in an insectary maintained at 27 °C and 80% relative humidity, with a photoperiod of 12 h light:12 h dark, with dusk and dawn periods of 1 h each at the beginning of each scotophase. All groups had continuous access to sterile distilled water and a 10% sucrose/distilled water solution throughout the duration of the study. The infection of P. falciparum was confirmed at the oocyst stage by sampling at least 8-10 mosquitoes from each of the 6 treatment groups (rearing cages) at 8 days post infection (Table 1). In total, 58 randomly selected mosquitoes were dissected in PBS, stained with 0.1% PBS-buffered mercurochrome (Mil-liporeSigma, Burlington, MA) and examined under a light-contrast microscope (Olympus Life Science, Waltham, MA). At the oocyst stage, P. falciparum infection prevalence of the 58 randomly selected infected mosquitoes was 62%, with a median count of two oocysts per mosquito midgut (SEM = 0.33).

Mosquito dissections and nucleic acid extractions.
After collection, each mosquito in this study was individually stored prior to manual trans-section of heads (including sensory appendages) and thorax including salivary glands, referred to as 'head' or 'head tissues' , from body/carcass which were stored for future use (Fig. 1B). Total RNAs extracted from the heads of established treatment-and control-group An. gambiae mosquitoes were used for Illumina NovaSeq sequencing. Total DNAs were also extracted from the heads of treatment mosquitoes for PCR to validate the P. falciparum infection at sporozoite stage in individual mosquitoes. All head tissues were collected from control-and treatment-group mosquitoes immediately following the 18-day incubation period (Fig. 1). Dissections of control-and treatment-group mosquito heads were performed over 3 days; each group was handled on a different day to prevent cross contamination. All mosquitoes were handled wearing gloves and all dissection equipment washed with glassware detergent, autoclaved, and decontaminated with RNase Away (MilliporeSigma, Burlington, MA). To eliminate potential hemolymph contamination during mosquito dissections, all equipment was sterilized using 90% ethanol and decontaminated with RNase Away (MilliporeSigma, Burlington, MA) in between individual mosquito handling. To maintain tissue integrity, dissections were performed on cold-anesthetized mosquitoes, and work limited to one group-treatment or control. For treatment groups, the head with intact salivary glands from one mosquito was dissected into a distinct numbered and cataloged 1.5-mL test tube containing 200µL Tri Reagent (MilliporeSigma, Burlington, MA). The smaller starting volume of Tri Reagent (MilliporeSigma, Burlington, MA) aided subsequent tissue homogenization. PCR infection validation studies on treatment-group mosquitoes required that all dissected heads remain separated. Prior to dissections, each treatment-group An. gambiae mosquito was randomly assigned an identification number used for labeling of all associated dissected tissues. In total, the heads of 176 treatment mosquitoes were collected into 176 distinct labeled tubes. For control lines, the heads from 30 mosquitoes within a single group were dissected and pooled into one labeled 1.5-mL test tube containing 400 µL Tri Reagent (MilliporeSigma, Burlington, MA). This process was repeated for the two remaining control groups. All sample Figure 13. Schematic depicting transcriptional shifts in An. gambiae following infection with P. falciparum sporozoites compared to naïve controls. Text is weighted to emphasize degree of transcript abundance change. For RNA extractions, sample tubes were processed in small batches of ten, with treatment and control tubes handled separately to prevent cross contamination. A batch of ten sample tubes containing An. gambiae head tissues in Tri Reagent (MilliporeSigma, Burlington, MA) was thawed on ice and the head tissues homogenized first manually using individual sterile, DNA/RNA-free disposable pellet pestles (Fischer Scientific, Waltham, MA) and second, mechanically by sonication (Q700; QSonica, Newtown, CT) to ensure thorough head tissue and cell disruption. Sonication was performed in short 10-to 20-s bursts with samples on ice to minimize nucleic acid shearing. To prevent cross contamination, one pellet pestle was used per sample tube, and the sonicator probe was thoroughly sterilized using 90% ethanol and decontaminated with RNase Away (Millipore-Sigma, Burlington, MA) in between individual tube processing. Following complete tissue homogenization, RNA extractions were performed. An additional 200µL Tri Reagent (MilliporeSigma, Burlington, MA) was added to each treatment sample tube bringing the total volume to 400 µL, equivalent to control sample volumes. Sample tubes were then hand shaken and incubated at room temperature for 5 min. A 100-µL aliquot of UltraPure phenol:chloroform:isoamyl alcohol (25:24:1, v/v; Thermo Fischer, Waltham, MA) was added to each of the tubes, which were then vigorously hand shaken for 1 min and incubated at room temperature for 5 min. For complete phase separation, sample tubes were centrifuged at 13,000 RPM for 15 min at 4 °C. The resulting aqueous phase (containing RNAs) from each sample tube was carefully transferred into a coordinately labeled new RNase-free 1.5-mL test tube containing 100µL 70% ethanol and thoroughly mixed by pipetting. Sample tubes containing the remaining interphase and organic phase were frozen at − 70 °C in preparation for DNA extractions. Total RNAs were extracted from the aqueous phase/ethanol solutions according to the manufacturer's protocol using the PureLink RNA Mini Kit (Thermo Fischer, Waltham, MA). This process was repeated for all 179 treatment and 3 control sample tubes. Total RNA concentrations and purities were measured using NanoDrop (Thermo Fischer, Waltham, MA) and any detected phenol contaminants removed by 1-butanol extraction 111 . Control-sample total RNAs were kept separate to allow for three distinct biological replicates pertinent to RNAseq analyses. For the treatment replicates, after each individual sample was confirmed to contain P. falciparum sporozoites via PCR validation, three distinct pools were prepared for RNAseq analyses. Ultimately, 65 of the 176 treatment heads were positively confirmed for P. falciparum infection and correlating total RNAs pooled into three replicates (21 randomly selected samples per replicate).
DNA extractions were performed on all 179 treatment sample tubes containing remaining interphase and organic phase solutions to perform PCR P. falciparum validation studies. Like RNA extractions, treatment samples were processed in small batches of ten to preserve DNA integrity. A batch of treatment sample tubes was thawed on ice, centrifuged at 13,000 RPM for 5 min at 4 °C, and any remaining aqueous phase was removed and disposed of to prevent RNA contamination of DNA extracts. Then, 200µL back extraction buffer was added to each sample tube and hand mixed by inversion for 3 min. For complete phase separation, sample tubes were centrifuged at 13,000 RPM for 30 min at 4 °C. The resulting aqueous phase (containing DNAs) from each sample was carefully transferred into a coordinately labeled, new sterile 1.5-mL test tube containing 200µL 70% ethanol and thoroughly mixed by pipetting. Total DNAs were extracted from the aqueous phase/ethanol solutions according to the manufacturer's protocol using the DNEasy Blood and Tissue Kit (Qiagen, Germantown, MD). Total DNA concentrations and purities were measured using NanoDrop (Thermo Fischer, Waltham, MA) and any detected phenol contaminants removed by 1-butanol extraction 111 . The total gDNA was finally resuspended at a final concentration of 10 ng/µL and kept at − 20 °C for further PCR experiments.

PCR-based plasmodium infection validation.
PCR experiments were performed to validate the presence of P. falciparum sporozoites within each treatment sample. PCR experiments used total DNAs extracted from dissected heads and salivary glands of individual treatment An. gambiae female mosquitoes (Ref. Mosquito Dissections and Nucleic Acid Extractions). PCR was performed using an MJ Thermal Cycler (Bio-Rad Inc, Hercules, CA) in combination with the DreamTaq Mastermix technology (Thermo Scientific, Waltham, MA) according to the manufacturer's protocol. PCR products were visualized on a 1% agarose gel using the Quickload Purple 1-Kb DNA ladder (New England Biolabs, Ipswich, MA). Primer3 Plus 112 was used to design primer pairs for the P. falciparum circumsporozoite gene (CS) and the An. gambiae ribosomal s7 gene (AgRsp7) used as a control to eliminate the risk of false negatives occurring from putative experimental error 19 . The primer sets used were as follows: All PCR primers were validated and optimized in preliminary studies. PCR experiments were conducted for all treatment samples, with a correlating negative control (blank) and positive P. falciparum infected positive control (Fig. 2). Treatment samples producing CS PCR products were reconfirmed to be positively infected through a second follow-up PCR experiment. www.nature.com/scientificreports/ collected from each of the three treatment replicates. Final mRNAs were fragmented, and cDNA was synthesized, amplified, digested, and purified using Illumina TruSeq chemistry protocols (Illumina, San Diego, CA). In preparation for sequencing, each of the six libraries-three treatment and three controls-was labelled with a 16-bp barcode to identify sequencing data generated from each library [treatment replicate 1 (T1) = TTG CAG  GTT GCA GCGG, treatment T2 = CAA CTG TAT GAG TGCC, treatment T3 = AGT CAA TGT AGA GGCC, control  C1 = ACA AGC TTA CAA CTAT, control C2 = CGT CCT AAC CAT CGGC, and control C3 = TTG AGA ACT CTC  TAAC]. Barcoded treatment and control replicate libraries were hybridized across four lanes of two S4 flow cells for cBot (Illumina, San Diego, CA) cluster generation and sequencing (paired-end, 2 × 150 bp; total 300 cycles). A total of 36 raw data sets of sequencing outputs, 6 per library, were generated using the Illumina software assembler (Illumina, San Diego, CA).