The Phlebotomus papatasi systemic transcriptional response to trypanosomatid-contaminated blood does not differ from the non-infected blood meal
Parasites & Vectors volume 14, Article number: 15 (2021)
Leishmaniasis, caused by parasites of the genus Leishmania, is a disease that affects up to 8 million people worldwide. Parasites are transmitted to human and animal hosts through the bite of an infected sand fly. Novel strategies for disease control require a better understanding of the key step for transmission, namely the establishment of infection inside the fly.
The aim of this work was to identify sand fly systemic transcriptomic signatures associated with Leishmania infection. We used next generation sequencing to describe the transcriptome of whole Phlebotomus papatasi sand flies when fed with blood alone (control) or with blood containing one of three trypanosomatids: Leishmania major, L. donovani and Herpetomonas muscarum, the latter being a parasite not transmitted to humans.
Of the trypanosomatids studied, only L. major was able to successfully establish an infection in the host P. papatasi. However, the transcriptional signatures observed after each parasite-contaminated blood meal were not specific to success or failure of a specific infection and they did not differ from each other. The transcriptional signatures were also indistinguishable after a non-contaminated blood meal.
The results imply that sand flies perceive Leishmania as just one feature of their microbiome landscape and that any strategy to tackle transmission should focus on the response towards the blood meal rather than parasite establishment. Alternatively, Leishmania could suppress host responses. These results will generate new thinking around the concept of stopping transmission by controlling the parasite inside the insect.
Leishmaniasis, a disease caused by parasites of the genus Leishmania, is endemic in 85 territories across the globe, with more than 600,000 cases worldwide and a prevalence of 4,000,000 . Leishmania spp. infect vertebrates through the bite of an infected sand fly vector (Diptera: Phlebotominae). The acute form of disease, visceral leishmaniasis (VL), or kala-azar, is fatal in 95% of untreated cases, with up to 50,000 people dying annually though non-fatal infections causing dermatological symptoms (cutaneous leishmaniasis), the most common consequence of infection . Although the ongoing VL elimination program in the Indian subcontinent is proving successful against the most severe clinical forms of VL , the elimination of leishmaniasis will likely require a combination of transmission-blocking strategies and novel treatments. This is especially the case in light of reports of resistance to drugs used to treat human infections [3, 4] as well as resistance to the pesticides used to control vector populations [5,6,7]. However, to develop approaches to blocking transmission, a better understanding is needed of the basic biology that underlines the interactions between parasite and insect vector.
The sand fly responses to blood-feeding have been investigated, with several gene families shown to be transcribed and/or expressed in response to a blood meal , including those encoding digestive enzymes, such as trypsins and chymotrypsins, pathogen recognition molecules and components of the peritrophic matrix, a protective chitinous mesh which lines the midgut after ingestion . However, few sand fly genes or transcripts specifically associated with Leishmania infection have been identified. There is some evidence to suggest that Leishmania are able to suppress host responses to promote survival and infection establishment. Analysis of cDNAs isolated from dissected sand flies Phlebotomus papatasi  and P. perniciosus  midguts revealed that several transcripts which are enriched after receipt of a blood meal are depleted when flies are fed blood containing Leishmania. These included digestive proteases, such as trypsins, as well as peritrophins, which are chitin-binding components of the peritrophic matrix that serves as a temporary barrier to leishmania . Moreover, the influence of Leishmania infection on physiological responses (oviposition and digestion) or traits (longevity) are far from dramatic .
Recently, we described both the host  and parasite  transcriptomes in a trypanosomatid–dipteran insect infection model, namely, Drosophila melanogaster and its natural trypanosomatid parasite Herpetomonas muscarum. We showed that parasite feeding resulted in differential transcription of two nuclear factor-κB (NF-κB) pathways, the Toll and the Immunodeficiency (Imd) signaling pathways, as well as the dual-oxidase (DUOX) pathway and STAT-dependent epithelial stem cell proliferation pathway. We found  that the H. muscarum transcriptome during infection closely resembled that reported for Leishmania major during Phlebotomus duboscqi infection . Transcriptional responses in Drosophila were detected in whole flies, and so even if H. muscarum infection was gastrointestinal, there seemed to be a systemic response involving the gut, the fat body and several secreted neuropeptides, all of which were presumably important in inter-tissue communication .
Given these findings, we believed a comparison of the Drosophila systemic transcriptional responses to those of sand flies during infection would provide valuable knowledge. Common transcriptomic signatures between the two dipteran insects would indicate an evolutionarily conserved response to trypanosomatid immune challenge. Such responses would be of great interest in terms of developing broad-spectrum transmission blocking strategies for trypanosomatid diseases. Conversely, different responses would imply clade-specific host–parasite interactions, with the possibility of potential suppression of the responses described above. However, no comprehensive data were available on the systemic response of the sand fly to Leishmania. Therefore, we sought to add to the body of work already available for sand fly transcriptomic responses to trypanosomatid infection using next generation sequencing (RNA-seq) in whole flies.
Here, we describe the transcriptome of sand fly P. papatasi at three time points corresponding to important stages of trypanosomatid infection: 1 day post blood meal (PBM); following blood-meal digestion and when parasites can be found attached to the midgut epithelium (4 days PBM); when parasites have migrated anteriorly in the gut and are found in the thoracic midgut and the stomodeal valve of the fly (9 days PBM; Fig. 1) . Infections were conducted in the context of both permissive (Leishmania major) and refractory (Leishmania donovani) infections, as well as with monoxenous (infects only insects) trypanosomatid Herpetomonas muscarum, which is not a natural parasite of sand flies. Using this strategy, we hoped to identify host transcriptional signatures associated with permissive and refractory infection outcomes, in addition to identifying evolutionarily conserved host responses as described above.
Our results indicate that there is very little difference between the transcriptomes of flies fed an infected blood meal and those fed a non-infected blood meal. Comparison of early and late time points of blood-only fed flies showed the transcription of genes from several immune pathways—including the Imd, Toll and JAK-STAT signaling pathways. Activation of these responses despite the absence of parasites in the meal may be a pro-active strategy by the sand flies to prevent infection.
Materials and methods
Phlebotomus papatasi maintenance
A laboratory colony of P. papatasi (originating from Turkey) was maintained in the insectary of the Charles University in Prague under standard conditions (26 °C, 60–70% humidity, 14/10 h ligh/dark photoperiod; fed on 50% sucrose) as described previously .
Leishmania donovani (MHOM/ET/2010/GR374), L. major LV561 (LRC-L137; MHOM/IL/1967/Jericho-II) and Herpetomonas muscarum  were cultured in M199 medium (Sigma-Aldrich, St. Louis, MO, USA) containing 10% heat-inactivated foetal bovine serum (FBS; Gibco, Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 1% Basal Medium Eagle vitamins (Sigma-Aldrich), 2% sterile urine, 250 μg/ml amikacin (Amikin; Bristol-Myers Squibb, Princeton Pike, NJ, USA) at 23 °C (L. donovani, L. major) or 28 °C (H. muscarum).
Phlebotomus papatasi infections
Leishmania and H. muscarum promastigotes from log-phase cultures (3–4 days post inoculation) were resuspended in defibrinated and heat-inactivated rabbit blood (LabMediaServis, Jaroměř, Czech Republic) at concentration of 1 × 106 promastigotes/ml, which corresponds to 500–1000 promastigotes per P. papatasi female . Sand fly females (5–9 days old) were infected by feeding on the suspension through a chick-skin membrane (BIOPHARM, Žďár, Czech Republic). Engorged sand flies were maintained under the same conditions as the colony. Each batch was left to develop an infection for 1, 4 or 9 days PBM. In terms of age, on the ninth day PBM sand flies would have been 14–18 days old. Across infections, each batch that we compared had the same age unless we were comparing batches with the same infection. There was no uninfected control for the age of the flies because our question was the difference between blood meal (our uninfected control) and infected blood meal. Age-related signatures would have been contained within the non-infected blood meal. Since systemic sand fly immunity does not become constitutively active during early healthy aging in the absence of infection, we attributed the observed immune activity to receiving a blood meal (infected or not).
RNA extraction and sequencing
Poly-A mRNA was purified from total RNA using oligodT magnetic beads, and strand-specific indexed libraries were prepared using the KAPA Stranded RNA-Seq kit followed by ten cycles of amplification using KAPA HiFi DNA polymerase (KAPA Biosystems, Wilmington, MA, USA). Libraries were quantified and pooled based on post-PCR analysis (Agilent bioanalyzer; Agilent Technologies, Santa Clara, CA, USA), and 75-bp paired-end reads were generated on the Illumina HiSeq v4 sequencing system (Illumina Incs., San Diego, CA, USA) following the manufacturer’s standard sequencing protocols. All raw sequencing reads are available on the European Nucleotide Archive under study accession number PRJEB35592.
Read mapping and differential expression analysis
Reads were mapped to the P. papatasi genome (Ppapl1 v1; Vectorbase)  using the HISAT2 alignment program2 . Those reads that mapped uniquely and in their proper pair were extracted and used to assemble transcripts de novo with the Cufflinks tool (Tuxedo suite) . These newly assembled transcripts were combined with the VectorBase transcript assembly to create a new set of transcripts using the CuffMerge Tool. Both the sequences of the assembled transcripts and the new annotation file (.gtf) are given in the Electronic Supplementary Material (ESM) data files. Reads were then counted against the transcripts generated by the Cufflinks tool using featureCounts . The counts data for the (two) technical replicates for each sample (= each sample was sequenced twice) were collapsed prior importing into R for differential expression analysis (pairwise Wald tests) in DESeq2 . A DESeq dataset was produced from the counts’ matrix using the inbuilt function using the experiment design: design = ~ Batch + Day + Condition. Pairwise Wald tests were used to identified differential transcript levels between samples. Reported P values are adjusted using the Benjamin–Hochberg method (the default in DESeq) to correct for multiple testing.
Analysis performed to conclude that several transcripts lacked conserved domains or that novel transcripts presented conserved domains used the NCBI conserved domains search tool (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) and the CDD v3.18 database. This is a superset including NCBI-curated domains and data imported from the Pfam, SMART, COG, PRK and TIGRFAM databases.
Results and discussion
Read mapping, de novo transcript assembly and differential expression analysis
Sand flies aged 5 to 9 days were fed blood or parasite-containing blood, resulting in four treatment groups of Phlebotomus papatasi females: (i) blood meal only, (ii) blood meal + L. major, (iii) blood meal + L. donovani and (iv) blood meal + Herpetomonas muscarum. Flies were then left to develop an infection or digest their non-infected blood meal for 1, 4 or 9 days (actual age = 18–19 days on ninth day PBM). Across infections, flies were compared at the same time point (same age) unless flies with the same infection were being compared.
RNA was purified from whole sand flies at 1, 4 and 9 days PBM, with the aim to identify any systemic responses over and above tissue-specific signals. Each time point analyzed included 20 flies; thus, 3 (time points) × 20 sand flies were processed per treatment group. This experiment was independently repeated three times for each treatment group so that ultimately 3 (replicates) × 60 (sand flies) =180 sand flies were processed for each treatment group in all three experiments. The goal of our experiment was to determine whether there was a difference between blood meal versus infected blood meal. Therefore, the non-infected blood meal served as the control.
As expected, only L. major was able to establish infection following digestion of the blood meal and defecation at day 4 PBM (Additional file 1: Figure S1). Of these infections, 95% developed into severe L. major infections at 8 days PBM, while the infections by the other two parasites were cleared (Additional file 1: Figure S1). At this time point (4 days PBM), L. major parasites were found in both the anterior and thoracic midgut as well as at the stomodeal valve (SV) (Additional file 2: Figure S2). In mature Leishmania infections, the SV is forced open and becomes blocked with parasites embedded in the promastigote secretory gel, a viscous mixture of phosphoglycans secreted by the parasites. This opening of the SV is essential for colonization of the foregut and transmission by regurgitation .
The resulting reads were sequenced and mapped against the P. papatasi genome (Ppapl1, Vectorbase) . The number of reads generated per sample ranged from 1.08 to 12.05 million reads with 69.7–79.3% of these mapping to the P. papatasi genome in each sample (Additional file 3: Table S1). Upon visual inspection of read mapping using Integrative Genomes Viewer , it appeared that > 20% of reads were mapping to regions that lacked annotated features. To include these potentially novel genes in our analysis we assembled de novo concordantly mapped read pairs (from all samples) into 16,025 transcripts. The assembled transcripts were then merged with the existing annotation of 11,834 transcripts to give a final set of 18,592 unique transcripts (see ESM data files). This represents approximately 97.2 Mb of P. papatasi transcriptome with an average transcript length of 4190 bp. All reads were then counted against the final set of transcripts for differential expression analysis.
Principal component analysis (PCA) showed a high degree of difference between the fly transcriptomes at day 1 PBM and those at day 4 or 9 PBM (Fig. 2), with transcriptomes from days 4 and 9 PBM appearing to be similar. It was also noted that samples did not clearly group in accordance with trypanosomatid feeding status.
Differential expression associated with trypanosomatid presence in the blood meal
A few differentially expressed transcripts specifically associated with trypanosomatids were present in the blood meal (Additional file 4: Table S2). We found no significant difference in transcript abundance between blood-fed and L. major-fed flies at any time point. Furthermore, we found in excess of 12,000 genes for which we rejected the hypothesis that expression had changed by twofold or more in either direction in pairwise comparisons (Wald test) between blood-only fed and blood+trypanosomatid-fed flies (Fig. 3; Table 1).
Differential abundance for some transcripts, however, was observed after H. muscarum and L. donovani feeding compared to blood-only feeding on control flies (Fig. 3d, g, respectively). There were significantly fewer transcripts for the gene PPAI009043, an orthologue to the D. melanogaster signalling protein Rho GTPase activating protein at 54D (RhoGAP54D), in flies fed H. muscarum than in the blood-fed controls at day 1 PBM (log2 fold-change 1.13, P-adj = 0.005). The Aedes aegypti and Anopheles gambiae RhoGAP54D orthologues are upregulated in blood-fed mosquitos compared to sugar-fed controls [23, 24]. Given this and that we did not see this response after Leishmania feeding, we suggest that this transcriptomic response may be H. muscarum specific. The biological significance of reduced RhoGAP54D transcription in this context remains unclear; however, the protein is linked to epithelial morphogenesis during Drosophila development  and so may also play a role in the mature insect gut.
In L. donovani-fed flies there were significantly fewer transcripts for the putative transporter TrpA1 (PPAI004036; log2 fold-change 2.8, P-adj = 0.025) versus blood-only fed flies at 9 days PBM. TrpA1 is more generally associated with chemo- and thermo-sensing [26, 27] in Drosophila; however a study by Du et al.  links TrpA1 to the expulsion of food-borne pathogens by increased defecation and the DUOX pathway (discussed in detail later in text). It can be speculated that a reduction in TrpA1 transcripts after L. donovani feeding may indicate modification of host defensive pathways to promote survival. We also found significantly more CUFF.12679 transcripts (log2 fold-change 16.8, P-adj = 0.0001) in L. donovani-fed flies than in blood-only fed flies. This novel transcript lacks conserved domains or sequence similarity to known dipteran gene transcripts.
Similarily, direct comparisons between trypanosomatid infections yielded few differentially expressed transcripts (Additional file 5: Table S3). At day 1 PBM the only differentially expressed transcript between the three infections was that of trypsin 1 (PPAI010956; P-adj = 0.035) which was twofold enriched in H. muscarum-fed flies compared to those fed L. donovani.
After defecation at around 4 days PBM, it is thought only parasites able to establish in the ectoperitrophic space persisted to develop mature infection . Despite the differences in the infection outcome reported in laboratory infections across the three trypanosomatids , there were few differences in the host transcriptome at this critical time point. Two transcripts were found to be significantly differentially abundant: one corresponding to the PPAI000999 gene and the other a novel transcript CUFF.14170. Both transcripts were found at significantly higher levels (P-adj = 0.04 and 2.27E-09, repectively; log2 fold-changes 4 and 18, respectively) in H. muscarum-fed flies compared to those fed L. donovani. PPAI000999 encodes a protein predicted to bind to chitin (GO:0006030, GO:0008061 and smart00494). The novel transcript CUFF.14170 has no known conserved domains, and BLAST searches against dipteran sequences did not yield any significant hits.
The most variation between the three infections was found at 9 days PBM, where six transcripts were found to be differentially expressed between Leishmania-fed and H. muscarum-fed flies. Compared to sand flies fed H. muscarum, those fed L. donovani had significantly more transcripts for TrpA1 (PPAI004036) and significantly fewer for the putative zinc metalloprotease PPAI010164 and novel transcript CUFF.12679. Sand flies fed L. major had significantly more transcripts for the hypothetical protein PPAI002947. Additionally, feeding with H. muscarum resulted in significantly more CUFF.14170 transcripts, a novel transcript identified in this study which lacks conserved domains, than both Leishmania infections (P-adj = 1.79E−09).
Overall, the above observations suggest that blood-feeding status is the major source of transcriptional variation in these sand flies—and not trypanosomatid infection. As such, we further investigated transcriptomic changes after blood-feeding alone in P. papatasi.
The P. papatasi transcriptome after blood-feeding
Shortly after blood-feeding there were significant changes in transcription that may be a universal response to a blood meal. The transcriptomes at day 1 PBM appeared to be very different to those at 4 (and 9) days PBM, with 12,289 significantly differentially regulated transcripts (Additional file 6: Table S4). However, after defecation of the blood-meal remnants, the transcriptome was comparatively stable with 264 differentially regulated transcripts (4 vs 9 days PBM; Additional file 7: Table S5). Due to the large number of differentially expressed transcripts highlighted by these comparisons we first investigated transcripts whose log2 fold-change was > 4 in either direction between time points. From this subset we were able to focus our analysis on a number of key genes and pathways that are discussed in subsequent sections (Additional file 8: Table S6; Additional file 9: Table S7).
Early transcriptomic responses to ingestion of blood meal are related to digestion, metabolism and immunity
Of the 217 transcripts differentially regulated by > 4-fold between 1 and 4 days PBM, 197 transcripts were found to be comparatively enriched at day 1 PBM and 20 were comparatively enriched at day 4 PBM; 98 of these transcripts did not contain known conserved domains.
Transcripts for putative and known trypsins were one of the most highly represented groups differentially regulated between day 1 and day 4 PBM. We observed upregulation of nine transcripts for putative trypsins and chymotrypsins—including the previously characterized chymotrypsins 1 (PPAI010833), chymotrypsin 3 (PPAI005023) and trypsin 4 (PPAI010456) [8, 29, 30]. We also observed upregulation of transcripts that may represent novel trypsins based on conserved domains and their similarity to other dipteran trypsin/chymotrypsin sequences, as they are not included in the current genome annotation (Ppap v1)  (CUFF.11666, CUFF.9493, CUFF.6542 and the chymotrypsins CUFF.15058, CUFF.16005, CUFF.15086, CUFF.14587, CUFF.12454). In contrast, the transcript putatively encoding for trypsin 1 (PPAI010956) was shown to be enriched at day 4 PBM compared to the earlier time point. The roles of trypsin and chymotrypsin-like serine proteases during blood digestion in hematophagous insects are well characterized, with expression levels varying according to the type of blood meal and the time since the last blood meal. Our findings agree with those reported previously showing upregulation of trypsins 3/4 and chymotrypsin 1 in response to the blood meal, as well as the decrease of trypsin 1 .
In addition to the trypsins themselves, five transcripts whose products are predicted to contain trypsin inhibitor-like domains (PPAI003932, PPAI000270, PPAI000272, PPAI000274, PPAI003557) were also comparatively enriched at day 1 PBM (vs day 4 PBM). It is possible that the corresponding proteins play roles in the regulation of the trypsin 1 as well as other trypsins (e.g. trypsin 2), reported to be downregulated after blood-feeding .
Several transcripts encoding for proteins with predicted serine protease/proteolytic activity, the sequences of which do not resemble trypsins/chymotrypsins, were also comparatively enriched at day 1 PBM. These included two known genes (PPAI009419, PPAI009871) and three novel transcripts (CUFF.6132, CUFF.6133, CUFF.16132). Serine proteases are implicated in several other cellular processes, including innate immune signaling—notably in Toll pathway activation —and the melanization response . The predicted protein for PPAI009419 shares approximately 51% identity with the Culex quinquefasciatus CLIPA15 (also known as masquerade) across its sequence. CLIPA proteases interact with and regulate other CLIPs and with the prophenoloxidases (PPO) involved in melanization [33, 34]. This response produces reactive quinones which then polymerize to form the dark insoluble pigment melanin. These molecules can encapsulate and isolate invading pathogens or toxic compounds. They also locally generate high local levels of cytotoxic reactive oxygen species (ROS) and prevent gas diffusion, starving the invading pathogen of oxygen. In addition to the putative CLIPA transcript, four pro-phenoloxidase transcripts were upregulated in early blood meal (PPO1: PPAI008831, PPAI010450; PPO2: PPAI012836, PPAI012835). These zymogens are the rate-limiting enzymes in the production of melanin. PPO1/2 and CLIPA15 were also found to be upregulated immediately after blood-feeding in Anopheles gambiae , suggesting that this is a conserved response to blood-feeding in dipterans.
We also observed differential transcription of another group of proteins reported to play vital roles in protection against invading pathogens—peritrophins. These core components of the peritrophic matrix have been shown to be a major barrier against infection establishment. Knockdown of Peritrophin 1 (Per1) in P. papatasi results in an approximately 40% increase in Leishmania major load at 48 h after parasite ingestion . In our study, Per1 transcripts were highly enriched at day 1 PBM (vs day 4 PBM) with a log2 fold-change of 9.96. Of the 32 annotated peritrophins in the P. papatasi genome, 14 were found to be significantly differentially regulated between days 1 and 4 PBM (Table 2). The majority of transcripts were comparatively enriched at day 1 PBM; however Per2 and Per28 transcripts were more abundant at later time points. Ramalho-Ortigão et al.  showed that P. papatasi Per1 transcripts were enriched in flies fed a blood meal compared to a sugar meal, while peritrophin 2 (Per2) transcripts were comparatively depleted in blood-fed flies. Additionally, the group showed that transcripts for both Per1 and Per2 were depleted in L. major-infected flies compared to those fed only blood . Our data largely agree with these findings. However, in our study transcript levels were not statistically significantly different between trypanosomatid-fed sand flies and those fed blood only—although we did observe fewer transcripts for Per2 (PPAI009723) in trypanosomatid-fed flies at day 4 PBM (Additional file 10: Figure S3). Other than Per2, the patterns in peritrophin transcript abundance for trypanosomatid-fed flies resembled those of the blood-fed controls.
Additionally, transcripts for another chitin-binding protein, PPAI000188, were significantly more abundant at 4 days PBM than at day 1 PBM. The sequence of PPAI000188 resembles the Lutzomyia longipalpis protein ChiBi (EU124616.1 , 84% protein sequence identity). ChiBi has been shown to be enriched in L. longipalpis fed with blood containing L. infantum chagasi . Its upregulation here in P. papatasi in the absence of Leishmania may indicate this upregulation is a more general response to a blood meal, rather than an infection-specific response.
In addition to trypsins, transcripts of several other groups of genes associated with digestion and nutrient uptake were differentially regulated PBM. Several transcripts for lipid metabolism-associated genes were found to be upregulated at day 1 PBM. In addition, eight transcripts corresponded to known extracellular carboxylic ester hydrolases (PPAI002323, PPAI003061, PPAI003086, PPAI005115, PPAI005116, PPAI005680, PPAI009133, PPAI008993). Similarly, transcripts for a putative sterol transfer protein (PPAI008838), and two paralogous membrane fatty acid desaturase genes (PPAI008098 and PPAI002108) were shown to be comparatively enriched at day 1 PBM. One transcript, CUFF.7417, does not correspond to a known gene; however the transcript showed strong sequence similarity to the extracellular carboxylic ester hydrolases paralogues PPAI005115 and PPAI005116 mentioned above (90% identity). Additionally, CUFF.7417 is immediately downstream of PPAI005115/6 in the genome and as such we propose this represents a previously unknown paralogue.
Four transcripts coded for proteins with solute carrier domains (cl00456). These transcripts encode for the two paralogous sodium-coupled monocarboxylate transporters (SCMTs; PPAI005125 and PPAI007402) and two putative SCMTs (CUFF.14648 and CUFF.14649). The SCMTs are transmembrane proteins that move molecules with a single carboxylate group, including pyruvate and lactate, across the plasma membrane in a proton-dependent manner; they are associated with the insect midgut brush border 
We found two transcripts, CUFF.17209 and CUFF.15972, whose products are predicted to contain the conserved insect allergen-related repeat domain (pfam06757). These transcript sequences also showed similarity to reported cDNAs for P. papatasi microvillar proteins MVP1 and 2, respectively (> 89% identity to mRNA sequences). These proteins were also found previously to be upregulated in sand flies upon ingestion of a blood meal compared to sucrose-fed flies . These transcripts could not be assigned to an annotated gene in the current vector base genome (Ppal1) . The function of these proteins is not well understood although they appear to have a conserved signal peptide at the N-terminus and lack transmembrane domains.
Finally, three olfactory (Or57 [PPAI013155], Or99 [PPAI013290] and the putative protein PPAI002404) and a gustatory receptor orthologous to sweet taste receptors of Drosophila (Gr9 [PPAI010978]) were upregulated at day 1 PBM compared to later time points. It is likely these sensory receptors were involved in sensing and acquisition of the blood meal and that subsequent decreases in their transcript abundances may indicate these sensors were not required after digestion.
The transcriptome after defecation of the blood meal is comparatively stable
The two later time points in this study had similar transcriptomic signatures, with only six transcripts comparatively enriched by > 2-fold at 9 days PBM (vs. 4 days). These transcripts corresponded to two glutamate receptors (PPAI003634, PPAI008275), apoptosis inhibitor survivin (PPAI002284), two histone methyltransferases (PPAI005539, PPAI005538) and a mucin (PPAI009152). Mucins have been implicated in the interaction with Leishmania parasites. Given that several immunity-related transcripts (including peritrophins, mucins and melanization pathway genes) were upregulated, we postulated that a general immune response was triggered upon ingestion of a blood meal and therefore investigated the transcription of the members of the two major innate pathways after a blood meal, namely the Toll and Immunodeficiency (Imd) signaling pathways. Both pathways have been shown to play a role in the response to trypanosomatids [11, 38,39,40,41,42]. We also investigated members of the DUOX and JAK-STAT pathways, both of which have been implicated in the D. melanogaster–H. muscarum interaction . Differential regulation statistics for these transcripts can be found in Additional file 11: Table S8.
Blood ingestion alone is associated with increased innate immune gene transcription
In blood-fed flies, transcripts putatively encoding two early Toll pathway genes (Spätzle and GNBP3) were found to be significantly enriched at day 1 PBM compared to days 4 and 9 PBM (Fig. 4a). One exception was the Spätzle-processing enzyme (SPE), the putative transcript for which was enriched at the latter two time points along with several intracellular Toll pathway components, indicating the possibility of a positive feedback loop back to Spätzle. These trends were broadly consistent in sand flies fed blood only as well as those fed with each of the trypanosomatids. However, only flies fed with blood containing L. major or L. donovani promastigotes had significantly higher levels of transcripts encoding Toll pathway inhibitor Cactus at day 1 PBM compared to day 4 PBM (> 2-fold, P-adj = 0.01 and 5.08E-05, respectively). Cactus transcript abundance was not significantly different between days 1 and 4 PBM in sand flies fed blood or H. muscarum.
A similar pattern emerged for the IMD pathway (Fig. 4b). Transcripts for putative peptidoglycan recognition proteins (PGRPs) were more prevalent at day 1 PBM compared to later time points (Fig. 4b). However, only putative PGRP 2 (CUFF.5670) was found to be statistically significantly enriched (2.23-fold) at day 1 PBM (vs days 4 and 9 PBM). The transcripts putatively encoding IMD, as well as several other proteins downstream of IMD in the pathway were found to be significantly enriched at 4 and 9 days PBM (vs day 1 PBM), including DREDD, TAK1 and IKKβ. We also observed significant enrichment of transcripts putatively encoding negative regulators of the IMD pathway Caspar, dUSP36, Trabid at days 4 and 9 PBM. Interestingly, the IMD transcription factor Relish was not significantly differentially regulated in blood-only-fed flies, while flies fed blood containing L. major or L. donovani promastigotes showed enrichment of putative Relish transcripts at day 1 PBM compared to 4 and 9 days PBM. As such, while there is overall upregulation of IMD pathway transcription with or without trypanosomatids in the blood meal, there may be important differences in the expression levels of the innate effectors the meal regulates when Leishmania are present.
Both Toll and IMD result in the expression of a suite of anti-microbial peptides. Transcripts for these immune effectors were not significantly differentially regulated after blood feeding alone. However, flies fed with blood containing L. major or L. donovani promastigotes were found to have significantly more transcripts for the AMP defensin at day 1 PBM compared to day 4 PBM (2.3 and 1.75 log2 fold-change, respectively). This was not observed in H. muscarum-fed flies. While transcript levels for other anti-microbial peptides did change after trypanosomatid feeding, for example we observed elevated transcript levels for cecropin and attacin in some trypanosomatid infections (Fig. 4c), overall these were not found to be statistically significant changes.
In addition to AMP expression, the IMD pathway can also result in the transcription of the NADPH oxidase, DUOX, through interaction of IMD with MEKK1 . This transmembrane protein is responsible for the production of ROS species in the gut epithelium in response to microbes. We found that the level of DUOX transcripts was significantly higher at days 4 and 9 PBM compared to day 1 PBM in all feeding conditions (log2 fold-change 2.98–3.33; Fig. 5a), with no significant difference in DUOX transcript abundance between 4 and 9 days PBM in any infection condition. Similarly, we saw significant increases in transcripts for genes upstream of DUOX across infection conditions including: the transcription factor ATF2, p38 kinase and MEKK1. As such, induction of DUOX pathway transcription appears to be a generalized response to blood-feeding rather than an infection-specific response.
The JAK–STAT pathway is also associated with the dipteran response to trypanosomatids
Finally, given the association between the JAK-STAT pathway (Fig. 5b), dipteran gut morphology and immunity , particularly in a trypanosomatid infection context , we also investigated the transcription of key components of this pathway after blood-feeding. We observed a higher abundance of putative Upd1 transcripts at day 1 PBM compared to later time points; however this change was only shown to be statistically significant for sand flies fed with blood and L. major where there was a 2.2-fold enrichment of putative Upd1 transcripts. Furthermore, putative transcripts for the JAK-STAT transcription factor STAT92E were twofold enriched in flies in all infection conditions at the two later time points (vs 1 day PBM). We also observed a modest enrichment of transcripts for cytokine Upd2 and the transmembrane receptor Domeless at days 4 and 9 PBM compared to earlier time points (fold-changes 1.19 and 1.74, respectively). The transcription pattern for signaling protein hopscotch resembled that of Domeless; however these transcripts were only found to be statistically significantly enriched in trypanosomatid-fed flies. Together these observations suggested an increase in JAK-STAT signaling a few days after a blood meal in P. papatasi. Further work to investigate if this signalling translates to changes in gut homeostasis, such as the increased stem cell proliferation observed in the Drosophila–Herpetomonas model, will be important. Currently, however, as transcript abundance for STAT92E is enriched in blood only fed controls this response does not appear to be trypanosomatid-specific.
Given the magnitude of the transcriptomic changes associated with blood-feeding alone, and the little variation between blood meals spiked with trypanosomatids that produce very different infections, we speculate that the aforementioned defensive responses are not infection specific. Such a strong response to the blood meal alone is not surprising given the additional stresses associated with the hematophagous habit . The high-risk nutrient attainment method drives the insects to take large volumes of blood at each meal; for example, mosquitoes and tsetse flies expand by up to threefold their pre-meal size after a blood meal [45, 46], which puts enormous mechanical strain on the tissues. In addition to volume of the meal, the content of their meal presents additional problems, such as excess water/ions , toxic compounds  and bacterial expansion in response to the rich meal [49, 50]. Due to the warm-blooded nature of their victim’s blood temperature, the body temperature of blood-feeding arthropods can rapidly (< 60 s) increase by over 10 °C during their meal [51, 52]. All of these factors must also be overcome even in blood meals which do not contain parasites. Additionally, activation of immune pathways following a blood meal (but in the absence of infection) may be a strategy to ‘pre-emptively’ protect the host against pathogens/toxic compounds which may be present in the newly ingested blood meal. Such anticipatory responses have been reported in other hematophagous insects, including the malaria vector Anopheles gambiae .
Moreover, it is known that in sand flies, the blood meal is followed by a decrease in overall gut bacterial diversity  coupled to an increased abundance of aerobic bacteria . It is possible that these changes may mask any effects from the presence of trypanosomatids. However, there was also no significant difference between day 9 PBM L. major and the other day 9 infections. This is important since by day 9 PBM, the blood meal has long been digested and it is only the L. major that is left while the other parasites are cleared. This underlines the non-specificity of the P. papatasi response and implies that for the sand fly, L. major is just another feature in the microbiome landscape. Alternatively, L. major establishment in the midgut could mediate a suppression of host responses. Consistent with this, a recent study comparing RNA-seq data generated from the midgut found just a 1% overall difference between sand flies receiving a blood meal and those receiving a Leishmania-contaminated blood meal . The limitation of our hypothesis is that we have performed our experiments in whole flies and so it could be that there are infection-specific responses in other tissues. More work is needed to verify that what we see at the organismal level is also the case at the level of immune-competent tissues (gut, hemocytes, fat body).
All data generated or analyzed during this study are included in this published article (and its additional information files). Raw data can be found at the European Nucleotide Archive, Project PRJEB35592 (ERP118668).
Courtenay O, Peters NC, Rogers ME, Bern C. Combining epidemiology with basic biology of sand flies, parasites, and hosts to inform leishmaniasis transmission dynamics and control. PLoS Pathog. 2017;13:e1006571.
Sundar S, Singh OP, Chakravarty J. Visceral leishmaniasis elimination targets in India, strategies for preventing resurgence. Expert Rev Anti Infect Ther. 2018;16:805–12.
Ponte-Sucre A, Gamarro F, Dujardin J-C, Barrett MP, Lopez-Velez R, Garcia-Hernandez R, et al. Drug resistance and treatment failure in leishmaniasis: A 21st century challenge. PLoS Negl Trop Dis. 2017;11:e0006052.
Basselin M, Denise H, Coombs GH, Barrett MP. Resistance to pentamidine in Leishmania mexicana involves exclusion of the drug from the mitochondrion. Antimicrob Agents Chemother. 2002;46:3731–8.
Hassan MM, Widaa SO, Osman OM, Numiary MSM, Ibrahim MA, Abushama HM. Insecticide resistance in the sand fly, Phlebotomus papatasi from Khartoum State Sudan. Parasites Vectors. 2012;5:46.
Dhiman RC, Yadav RS. Insecticide resistance in phlebotomine sandflies in Southeast Asia with emphasis on the Indian subcontinent. Infect Dis Poverty. 2016;5:106.
Fawaz EY, Zayed AB, Fahmy NT, Villinski JT, Hoel DF, Diclaro JW. Pyrethroid insecticide resistance mechanisms in the adult Phlebotomus papatasi (Diptera: Psychodidae). J Med Entomol. 2016;53:620–8.
Abrudan J, Ramalho-Ortigão M, O’Neil S, Stayback G, Wadsworth M, Bernard M, et al. The characterization of the Phlebotomus papatasi transcriptome. Insect Mol Biol. 2013;22:211–32.
Dillon RJ, Lane RP. Influence of Leishmania infection on blood-meal digestion in the sandflies Phlebotomus papatasi and P. langeroni. Parasitol Res. 1993;79:492–6.
Dostálová A, Votýpka J, Favreau AJ, Barbian KD, Volf P, Valenzuela JG, et al. The midgut transcriptome of Phlebotomus (Larroussius) perniciosus, a vector of Leishmania infantum: comparison of sugar fed and blood fed sand flies. BMC Genomics. 2011;12:223.
Wang L, Sloan MA, Ligoxygakis P. Intestinal NF-κB and STAT signalling is important for uptake and clearance in a Drosophila–Herpetomonas interaction model. PLoS Genet. 2019;15:e1007931.
Sloan MA, Brooks K, Otto TD, Sanders MJ, Cotton JA, Ligoxygakis P. Transcriptional and genomic parallels between the monoxenous parasite Herpetomonas muscarum and Leishmania. PLoS Genet. 2019;15:e1008452.
Inbar E, Hughitt VK, Dillon LAL, Ghosh K, El-Sayed NM, Sacks DL. The transcriptome of Leishmania major developmental stages in their natural sand fly vector. mBio. 2017;8:e00029-17.
Dostálová A, Volf P. Leishmania development in sand flies: parasite–vector interactions overview. Parasites Vectors. 2012;5:276.
Volf P, Volfova V. Establishment and maintenance of sand fly colonies. J Vector Ecol. 2011;36(Suppl 1):S1-9.
Pruzinova K, Sadlova J, Seblova V, Homola M, Votypka J, Volf P. Comparison of bloodmeal digestion and the peritrophic matrix in four sand fly species differing in susceptibility to Leishmania donovani. PLoS ONE. 2015;10:e0128203.
Giraldo-Calderón GI, Emrich SJ, MacCallum RM, Maslen G, Dialynas E, Topalis P, et al. VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic Acids Res. 2015;43:D707–13.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative Genomics Viewer. Nat Biotechnol. 2011;29:24–6.
Bonizzoni M, Dunn WA, Campbell CL, Olson KE, Dimon MT, Marinotti O, et al. RNA-seq analyses of blood-induced changes in gene expression in the mosquito vector species Aedes aegypti. BMC Genomics. 2011;12:82.
Marinotti O, Calvo E, Nguyen QK, Dissanayake S, Ribeiro JMC, James AA. Genome-wide analysis of gene expression in adult Anopheles gambiae. Insect Mol Biol. 2006;15:1–12.
Greenberg L, Hatini V. Systematic expression and loss-of-function analysis defines spatially restricted requirements for Drosophila RhoGEFs and RhoGAPs in leg morphogenesis. Mech Dev. 2011;128:5–17.
Kim SH, Lee Y, Akitake B, Woodward OM, Guggino WB, Montell C. Drosophila TRPA1 channel mediates chemical avoidance in gustatory receptor neurons. Proc Natl Acad Sci USA. 2010;107:8440 LP – 8445.
Hamada FN, Rosenzweig M, Kang K, Pulver SR, Ghezzi A, Jegla TJ, et al. An internal thermal sensor controlling temperature preference in Drosophila. Nature. 2008;454:217–20.
Du EJ, Ahn TJ, Kwon I, Lee JH, Park J-H, Park SH, et al. TrpA1 Regulates defecation of food-borne pathogens under the control of the Duox pathway. PLoS Genet. 2016;12:e1005773.
Ramalho-Ortigão M, Jochim RC, Anderson JM, Lawyer PG, Pham V-M, Kamhawi S, et al. Exploring the midgut transcriptome of Phlebotomus papatasi: comparative analysis of expression profiles of sugar-fed, blood-fed and Leishmania-major-infected sandflies. BMC Genomics. 2007;8:300.
Ramalho-Ortigao M, Kamhawi S, Rowton E, Ribeiro JMC, Valenzuela J. Cloning and characterization of trypsin and chymotrypsin-like proteases from the midgut of the sandfly vector Phlebotomus papatasi. Insect Biochem Mol Biol. 2003;33:163–71.
Buchon N, Poidevin M, Kwon H-M, Guillou A, Sottas V, Lee B-L, et al. A single modular serine protease integrates signals from pattern-recognition receptors upstream of the Drosophila Toll pathway. Proc Natl Acad Sci USA. 2009;106:12442–7.
An C, Zhang M, Chu Y, Zhao Z. Serine protease MP2 activates prophenoloxidase in the melanization immune response of Drosophila melanogaster. PLoS ONE. 2013;8:e79533–e79533.
Yassine H, Kamareddine L, Chamat S, Christophides GK, Osta MA. A serine protease homolog negatively regulates TEP1 consumption in systemic infections of the malaria vector Anopheles gambiae. J Innate Immun. 2014;6:806–18.
Kan H, Kim C-H, Kwon H-M, Park J-W, Roh K-B, Lee H, et al. Molecular control of phenoloxidase-induced melanin synthesis in an insect. J Biol Chem. 2008;283:25316–23.
Coutinho-Abreu IV, Sharma NK, Robles-Murguia M, Ramalho-Ortigao M. Characterization of Phlebotomus papatasi Peritrophins, and the role of PpPer1 in Leishmania major survival in its natural vector. PLoS Negl Trop Dis. 2013;7:e2132.
Jochim RC, Teixeira CR, Laughinghouse A, Mu J, Oliveira F, Gomes RB, et al. The midgut transcriptome of Lutzomyia longipalpis: comparative analysis of cDNA libraries from sugar-fed, blood-fed, post-digested and Leishmania infantum chagasi-infected sand flies. BMC Genomics. 2008;9:15.
Javed MA, Coutu C, Theilmann DA, Erlandson MA, Hegedus DD. Proteomics analysis of Trichoplusia ni midgut epithelial cell brush border membrane vesicles. Insect Sci. 2019;26:424–40.
Telleria EL, SantAnna MRV, Ortigão-Farias JR, Pitaluga AN, Dillon VM, Bates PA, et al. Caspar-like gene depletion reduces leishmania infection in sand fly host Lutzomyia longipalpis. J Biol Chem. 2012;287:12985–93.
Tinoco-Nunes B, Telleria EL, da Silva-Neves M, Marques C, Azevedo-Brito DA, Pitaluga AN, et al. The sandfly Lutzomyia longipalpis LL5 embryonic cell line has active Toll and Imd pathways and shows immune responses to bacteria, yeast and Leishmania. Parasit Vectors. 2016;9:222.
Boulanger N, Bulet P, Lowenberger C. Antimicrobial peptides in the interactions between insects and flagellate parasites. Trends Parasitol. 2006;22:262–8.
Boulanger N, Lowenberger C, Volf P, Ursic R, Sigutova L, Sabatier L, et al. Characterization of a defensin from the sand fly Phlebotomus duboscqi induced by challenge with bacteria or the protozoan parasite Leishmania major. Infect Immun. 2004;72:7140–6.
Hu C, Aksoy S. Innate immune responses regulate trypanosome parasite infection of the tsetse fly Glossina morsitans morsitans. Mol Microbiol. 2006;60:1194–204.
Ha EM, Lee KA, Park SH, Kim SH, Nam HJ, Lee HY, et al. Regulation of DUOX by the Gαq-phospholipase Cβ-Ca2+ pathway in Drosophila gut immunity. Dev Cell. 2009;16:386–97.
Myllymaki H, Ramet M. JAK/STAT pathway in Drosophila immunity. Scand J Immunol. 2014;79:377–85.
Benoit JB, Denlinger DL. Bugs battle stress from hot blood. Elife. 2017;6:e33035.
Lehane MJ, Lehane MJ, Lehane MJ. The biology of blood-sucking in insects. Cambridge: Cambridge University Press; 2005.
Beyenbach KW, Piermarini PM. Transcellular and paracellular pathways of transepithelial fluid secretion in Malpighian (renal) tubules of the yellow fever mosquito Aedes aegypti. Acta Physiol. 2011;202:387–407.
Sterkel M, Oliveira JHM, Bottino-Rojas V, Paiva-Silva GO, Oliveira PL. The dose makes the poison: nutritional overload determines the life traits of blood-feeding arthropods. Trends Parasitol. 2017;33:633–44.
Wang Y, Gilbreath TM III, Kukutla P, Yan G, Xu J. Dynamic gut microbiome across life history of the malaria mosquito Anopheles gambiae in Kenya. PLoS ONE. 2011;6:e24767.
Volf P, Kiewegova A, Nemec A. Bacterial colonisation in the gut of Phlebotomus duboseqi (Diptera: Psychodidae): transtadial passage and the role of female diet. Folia Parasitol (Praha). 2002;49:73–7.
Lahondere C, Insausti TC, Paim RMM, Luan X, Belev G, Pereira MH, et al. Countercurrent heat exchange and thermoregulation during blood-feeding in kissing bugs. Elife. 2017;6:e26107.
Lahondère C, Lazzari CR. Thermal effect of blood feeding in the telmophagous fly Glossina morsitans morsitans. J Therm Biol. 2015;48:45–50.
Upton LM, Povelones M, Christophides GK. Anopheles gambiae blood feeding initiates an anticipatory defense response to Plasmodium berghei. J Innate Immun. 2015;7:74–86.
Kelly PH, Bahr SM, Serafim TD, Ajami NJ, Petrosino JF, Meneses C, et al. The gut microbiome of the vector Lutzomyia longipalpis is essential for survival of Leishmania infantum. MBio. 2017;8:e01121-e1216.
Coutinho-Abreu IV, Serafim TD, Meneses C, Kamhawi S, Oliveira F, Valenzuela JG. Leishmania infection induces a limited differential gene expression in the sand fly midgut. BMC Genomics. 2020;21:608.
We thank the staff of the DNA Pipelines at Wellcome Sanger Institute for sequencing and generating sequencing libraries.
This work was supported by the European Commission, Horizon 2020 Infrastructure Infravec2 project (https://infravec2.eu). JS and PV were supported by ERD Funds, project CePaViP (CZ.02.1.01/16_019/0000759). MJS and JAC were supported by Wellcome via their core support for the Wellcome Sanger Institute (WSI) through grant 206194. Work in Oxford was supported by a Consolidator grant from the European Research Council (310912 Droso-Parasite, to PL), project grant BB/K003569 from the BBSRC (to PL) and a Wellcome Trust doctoral scholarship (to MAS).
Ethics approval and consent to participate
Not applicable for this study.
Consent for publication
Not applicable for this study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Time-series of infection intensity (% of infected females) for the three trypanosomatids tested. As expected, L. major developed late stage infections while the other two parasites were lost during defecation of blood meal remains.
Time-series of infection localization for the three trypanosomatids used.
Read mapping summaries. This table shows the read mapping information for each sample; for example, number of reads, percentage of read mapped, among others.
Transcripts associated with trypanosomatid presence in the blood meal. This table shows the fold-changes and differential regulation statistics (including P values) for transcripts whose abundance differed between trypanosomatid fed flies and blood-fed control flies.
Transcripts associated with specific trypanosomatids in the blood meal. This table shows the fold changes and differential regulation statistics (including P values) for transcripts whose abundance differed between trypanosomatid infections.
Transcripts significantly differentially regulated between 1 day and 4 days post blood meal (blood-only) in P. papatasi.
Transcripts significantly differentially regulated between 4 day and 9 days post blood meal (blood-only) in P. papatasi.
Transcripts of interest which are differentially regulated between 1 day and 4 days post blood meal (blood-only) in P. papatasi. This is a streamlined version of Additional file 6: Table S4 showing transcripts of interest discussed in the text.
Transcripts of interest which are differentially regulated between 4 days and 9 days post blood meal (blood-only) in P. papatasi. This is a streamlined version of Additional file 6: Table S4 showing transcripts of interest discussed in the text.
Differential regulation statistics for transcripts of dipteran immune pathways of interest (Toll, Imd, DUOX and JAK-STAT) across samples. ns Not significantly differentially regulated.
Log-normalized transcript counts for Peritrophin 2 (Per2) in P. papatasi throughout infection. Error bars show the standard error of the mean.
About this article
Cite this article
Sloan, M.A., Sadlova, J., Lestinova, T. et al. The Phlebotomus papatasi systemic transcriptional response to trypanosomatid-contaminated blood does not differ from the non-infected blood meal. Parasites Vectors 14, 15 (2021). https://doi.org/10.1186/s13071-020-04498-0