Antibody Responses Against Anopheles darlingi Immunogenic Peptides in Plasmodium Infected Humans

Introduction: Malaria is still an important vector-borne disease in the New World tropics. Despite the recent decline in malaria due to Plasmodium falciparum infection in Africa, a rise in Plasmodium infections has been detected in several low malaria transmission areas in Latin America. One of the main obstacles in the battle against malaria is the lack of innovative tools to assess malaria transmission risk, and the behavioral plasticity of one of the main malaria vectors in Latin America, Anopheles darlingi.

Methods: We used human IgG antibodies against mosquito salivary gland proteins as a measure of disease risk. Whole salivary gland antigen (SGA) from Anopheles darlingi mosquitoes was used as antigen in Western blot experiments, in which a ~65 kDa protein was visualized as the main immunogenic band and sent for sequencing by mass spectrometry. Apyrase and peroxidase peptides were designed and used as antigens in an ELISA-based test to measure human IgG antibody responses in people with different clinical presentations of malaria.

Results: Liquid chromatography–mass spectrometry revealed 17 proteins contained in the ~65 kDa band, with an apyrase and a peroxidase as the two most abundant proteins. Detection of IgG antibodies against salivary antigens by ELISA revealed a significant higher antibody levels in people with malaria infection when compared to uninfected volunteers using the AnDar_Apy1 and AnDar_Apy2 peptides. We also detected a significant positive correlation between the anti-peptides IgG levels and antibodies against the Plasmodium vivax and P. falciparum antigens PvMSP1 and PfMSP1. Odd ratios suggest that people with higher IgG antibodies against the apyrase peptides were up to five times more likely to have a malaria infection.

Conclusion: Antibodies against salivary peptides from An. darlingi salivary gland proteins may be used as biomarkers for malaria risk.

Introduction

Malaria is one of the most important parasitic diseases of the tropics and subtropics. Despite global decreases in malaria cases, Latin America has experienced an increase in cases since 2015. Colombia is among the top five Latin American countries with the highest numbers of malaria cases reported annually (Chaparro-Narvaez et al., 2016). In 2019 the number of malaria cases officially reported was 78,513 with Chocó, Nariño, Cordoba, and Antioquia, the departments with the highest number of malaria cases (31.4, 20.4, 11.6, and 9.3%, respectively) (Instituto Nacional de Salud, 2019). Several factors contribute to this phenomenon. First, more than 80% of the Colombian territory is suitable for the transmission of malaria due to the presence of a diverse number of mosquito species capable of serving as vectors of Plasmodium spp. (Rodriguez et al., 2011). Second, at least three species of Plasmodium cause infection in Colombia, with Plasmodium vivax as the main species, followed by Plasmodium falciparum and Plasmodium malariae (Arevalo-Herrera et al., 2012). Forced migration and precarious socioeconomic conditions also favor the spread and prevalence of Plasmodium spp. (Martens and Hall, 2000; Jitthai, 2013; Rodrigues et al., 2018).

Globally, most efforts to prevent malaria transmission focus on P. falciparum, the deadliest malaria parasite species. However, P. vivax is also a significant cause of morbidity and mortality (Chaparro-Narvaez et al., 2016; Geleta and Ketema, 2016; Zain Ul et al., 2016), and improved methods of monitoring its transmission are needed. Malaria transmission in Colombia occurs year-round, with peaks typically occurring between February and June (Ruiz et al., 2006). Malaria incidence is highly associated with changes in climatic conditions (i.e., temperature, elevation, and humidity) (Ruiz et al., 2006; Cassab et al., 2011). More importantly, ~87% of the reported cases come from people living in rural areas with limited access to health care (Rodriguez et al., 2011). Frequently, research about malaria transmission in these areas is restricted due to the lack of infrastructure to support research, including limited access to electricity.

Currently, the gold standard to determine malaria transmission intensity is the Entomological Inoculation Rate (EIR), which represents the number of infectious mosquito bites a person receives per night. This provides an indirect estimation of the risk of acquiring the disease. One of the challenges of this technique is the need to detect parasites in mosquitoes that are collected in the field. Mosquito trapping systems can be biased and overrepresent the most abundant species in a specific area. In addition, current methods are unable to pinpoint the exact individuals who have been bitten or whether there are differences in the number of bites different people receive. Thus, current surveillance efforts lack individualized, deployable, and cost-effective tools to measure human-vector contact and improve measurement of malaria transmission risk. Such tools would be particularly valuable in assessing the efficacy of vector control interventions in protecting people from malaria vector bites.

Our previous studies demonstrated that IgG antibodies against salivary proteins of major disease vectors correlate with the clinical presentation of dengue fever and malaria (Londono-Renteria et al., 2010; Londoño-Rentería et al., 2015a; Montiel et al., 2020a), and can be used to evaluate vector control interventions (Londono-Renteria et al., 2015c). In fact, IgG antibodies against a salivary peptide (gSG6-P1) from a major African malaria vector has been validated as a reliable marker to measure exposure to several Anopheles spp. (Poinsignon et al., 2009; Badu et al., 2012; Sagna et al., 2013). However, recent studies suggest that the gSG6 gene is not present in the subgenus Nyssorhynchus, which is the most prevalent American Anopheles subgenus (Arca et al., 2017). Although, previous studies have shown that the gSG6-P1 may be useful in areas of the New World where species belonging to subgenus other than Nyssorhynchus (Londono-Renteria et al., 2015b; Montiel et al., 2020a), there is a need to detect and validate novel biomarkers for exposure to these Anopheles spp. There are at least 40 Anopheles species in Colombia (Montoya-Lerma et al., 2011). The main vectors in the country are Anopheles nuneztovari, Anopheles albimanus and Anopheles darlingi—all three species within the subgenus Nyssorhynchus (Montoya-Lerma et al., 2011). Recently, we described several peptides that could be used to determine exposure to Anopheles albimanus bites, with antibodies against transferrin and salivary peroxidase proteins associated with the clinical presentation of malaria (Montiel et al., 2020a).

In this study, we evaluated the immunogenicity of salivary proteins in An. darlingi salivary gland antigen (SGA) and designed four peptides to measure IgG responses in people living in malaria endemic areas. The objective was to evaluate whether IgG antibodies against these peptides could be reliable markers to detect the risk of disease. An. darlingi is one of the most prolific mosquito vectors of malaria in Central and South America (Mirabello and Conn, 2006; Villarreal-Trevino et al., 2015; Prado et al., 2019), known for its behavioral plasticity and highly anthropophilic behavior. This mosquito species bite indoors and outdoors (endophagic/exophagic) and display exophilic/endophilic opportunistic behaviors (Moutinho et al., 2011; Reinbold-Wasson et al., 2012; Moreno et al., 2015). Our study suggests that IgG antibodies against An. darlingi salivary proteins are useful to assess malaria risk in areas where this species is endemic. This is the first study describing An. darlingi salivary peptides as potential tools to evaluate immune responses against mosquito saliva with the aim to improve accurate estimation of malaria transmission intensity, which is imperative for directing control efforts and predicting the effects of such interventions.

Materials and Methods

Ethical Considerations and Human Sample Collection

Serum samples from 179 participants were collected by active and passive case detection from November 2016 to October 2017 in Turbo and El Bagre as part of a malaria study in Colombia. These represent two malaria endemic areas with differing malaria incidence as measured by the Annual Parasite Index (API: confirmed cases during 1 year/population under surveillance ×1,000): the API in Turbo was 1.22 in 2016 and 0.77 in 2017, while malaria incidence in El Bagre was higher, with an API of 24.7 in 2016 and 21.3 in 2017 (Londoño-Rentería et al., 2015a; Londono-Renteria et al., 2020; Montiel et al., 2020b). Participant characteristics are described in Table 1. No concurrent entomological data was collected for this study. However, these field sites were selected based on a previous study reporting entomological data, where more than 96% of mosquitoes collected in Turbo were An. albimanus, while more than 70% of the mosquitos in El Bagre were An. darlingi (Gutierrez et al., 2009). Human sample collection and testing protocols were approved by the IRB of Kansas State University (IRB #8956 and IBC #1206) and the Medical School at the Universidad de Antioquia in Medellín, Colombia (Record 011 dated 28 July 2016).

TABLE 1

Table 1. Description of the serum samples by malaria infection status and location.

Detection of Malaria Parasites

On-site malaria diagnosis was done by microscopy, then confirmed by nested PCR (nPCR) at the Universidad de Antioquia (Colombia) following methods described previously (Londoño-Rentería et al., 2015a; Londono-Renteria et al., 2020). Based on these results, samples were classified as malaria positive or malaria negative, based on the presence/absence of parasites or a positive PCR result. In case of a disagreement, the PCR results was taken as the final diagnosis.

An. darlingi Mosquito Rearing

All An. darlingi were reared in the insectary at the Naval Medical Research Unit No. 6 (NAMRU-6) in Iquitos, Peru as described previously (Villarreal-Trevino et al., 2015). Larvae were reared in a dedicated larvae room (26.8 ± 0.7°C and 76.1 ± 6.3% relative humidity, and 12-h light:12-h dark photoperiod) and adults were maintained in a dedicated adult insectary room (25.9 ± 0.8°C and 69.7 ± 5.7% relative humidity, and 12-h light:12-h dark photoperiod) (Villarreal-Trevino et al., 2015). Eight to ten- day old adult mosquitoes were used for dissections of salivary glands. These mosquitoes were blood fed with chicken blood at an age of 3–5 days post-emergence, and provided with 10% sucrose solution ad libitum.

Salivary Gland Antigen (SGA) Preparation

Female An. darlingi mosquitoes were cold-anesthetized, washed in 70% ethanol, and placed in PBS (pH 7.2) for salivary gland dissection. Salivary glands were pooled (100 mosquitoes per pool for a total of two pools) and placed in 1X PBS. To prepare the antigen extract, the pool of salivary glands was allowed to freeze at −80°C and thaw at 4°C four times to induce cell rupture and release of proteins; the resulting SGA was kept in PBS at −80°C until use. Protein concentration was determined using a Thermo Scientific NanoDrop™ (Thermo Fisher Scientific, Wilmington, DE, USA) (Londono-Renteria et al., 2013, 2020; Montiel et al., 2020a).

Immunoblotting

To identify immunogenic proteins in An. darlingi SGA, pooled sera from ten uninfected subjects were used (five subjects from Turbo and five subjects from El Bagre) were used in the Immunoblotting experiments performed according to Londono-Renteria et al. (2010) (Londono-Renteria et al., 2010). In brief, 2 μg/well of SGA was separated in a 4–15% polyacrylamide gel (BioRad) at 100 V and transferred to a nitrocellulose membrane (Turbo PVDF, 0.45-m pore size, 8.5 by 13.5 cm, Biorad) in a Transblot Turbo following the standard 30 min transference protocol. Membranes were blocked for 1 h with blocking buffer (2% milk in PBST) and incubated with pooled human sera (1/250) at 4°C overnight. Each membrane was washed five times with wash solution (1X PBST) and incubated with HRP-conjugate anti human-IgG, in a 1/1,000 dilution during 1 h at 37°C. Color development was obtained with HRP chromogenic substrate TMB (Novex, Invitrogen). A pre-stained molecular weight marker of 16.5Ð*210 kDa (BioRad) was used for the estimation of the size of the proteins.

Liquid Chromatography/Mass Spectroscopy (LC/MS)

A band of 65 kDa was identified by human serum in the immunoblot. Gel bands corresponding to the immunogenic bands were excised and sent in duplicate to the National Institutes of Health (NIH/NIAID) for identification by mass spectrometry (MS). Briefly, the samples were reduced in solution containing 50 mM HEPES (pH 8.0) 10% acetonitrile and 5 mM DTT for 40 min at 37°C. After cooling to room temperature, the samples were mixed at 15 mM with iodoacetamide. After 15 min of alkylation, 200 ng of trypsin was added and the samples were incubated at 37°C for 15 h in a final volume of 40μl. The solution was evaporated to near dryness under a vacuum at 50°C. Twenty-five microliters of 0.1% trifluoroacetic acid was added and the pH was adjusted to 2.5 with the addition of 10% trifluoroacetic acid. Samples with an estimated protein content of <2 μg were desalted and concentrated with C18 μ Zip Tips. Samples containing up to 10μg were desalted with C18 OMIX 10 solid phase extraction tips. The digests were eluted with 0.1% TFA, 50% acetonitrile and dried under vacuum. The peptides were dissolved in 12 μl 0.1% formic acid, 3% acetonitrile which was used as the injection solvent. Digested peptides were subjected to the LC/MS analysis using an Orbitrap Fusion mass spectrometer (ThermoFisher Scientific, West Palm Beach, FL) connected with EASY nLC 1,000 liquid chromatography (LC) system. Nano-LC was carried out with a 5 μL injection onto a PepMap 100 C18 3-μm trap column (2 cm, ID 75 μm) and a 2 μm PepMap RSLC C18 column (25 cm, ID 75 μm), both from ThermoFisher Scientific. The LC was operated at a 300 μL/min flow rate with a 100-min linear gradient from 100% solvent A (0.1% formic acid, and 99.9% water) to 40% solvent B (0.1% formic acid, 20% water, and 79.9% acetonitrile) followed by a column wash. A standard data-dependent acquisition performed with a full MS spectrum was obtained by the Orbitrap for m/z 400–2,000 at the resolution of 120,000 with EASY-IC calibration. The precursor ions, with charges from two to eight, were selected, isolated (1.6 m/z window), fragmented by CID, then scanned by the Ion Trap. Survey scans were performed every 2 s and the dynamic exclusion was enabled for 30 s.

Acquisitions were searched against the NCBI-nr (05/2019) proteome and the cRAP.fasta database (theGPM.org) using PEAKS v10 (Bioinformatics Solutions Inc, Ontario, Canada) using a semi-tryptic search strategy with tolerances of 6 ppm for MS and 0.5 Da for MS/MS and carbamidomethylation of cysteine as a fixed modification and oxidation of methionine as a dynamic modification allowing for two missed cleavages. Peptides were filtered with a 0.5% FDR using a decoy database approach and a two spectral matches/peptide requirement.

Anopheles Immunogenic Peptide Selection and Plasmodium Antigens

Immunogenic peptides were designed from the most abundant identified proteins: apyrase (Protein Accession # ETN63669.1) and peroxidase (Protein Accession # ETN66035.1) (Table 2). Proteins were analyzed for the presence of signal peptide (signal P) sequence using the SignalP 5.0 server and for sequence homologies (at least 50% identity with E-value 1 × 10–5) to Anopheles and other major culicid disease vector species using online BLAST program. Protein sequence and structure analyses were then performed for the top Anopheles-specific proteins using the Protean 3D package of the DNASTAR software (DNASTAR Inc., Madison, WI, USA) (Londono-Renteria et al., 2020). The analysis of linear B-cell epitopes, antigenic regions, flexible regions (turns), hydrophobic regions, stability and/or charge density were conducted by importing the FASTA files from NCBI of sequences of interest into SVMTriP (Center for Plant Innovation, University of Nebraska, Lincoln, NE, USA). Highly antigenic, stable, flexible, charged and <50% hydrophobic regions that were specific to Anopheles mosquitoes (>50% homology) were selected using the software for 18–22 amino acid peptides (Table 3). The AnDar_PeroX1, AnDar_PeroX2, AnDar_Apy1, and the AnDar_Apy2 peptides of interest were sent for synthesis to GenScript (Piscataway, NJ, USA).

TABLE 2

Table 2. List of proteins identified in an ~65 kDa protein band from the western-blot using a pool of serum from study participants.

TABLE 3

Table 3. An. darlingi peptide sequences from the apyrase and the peroxidase/oxidase proteins identified by mass spectrometry.

Two recombinant proteins representing the Plasmodium Merozoite Surface Protein from P. falciparum (Pf-MSP) (Montiel et al., 2020a) and P. vivax (Pv-MSP) (Fitzgerald, USA) were used to evaluate antibody responses against the parasite.

ELISA Testing

ELISA conditions were standardized as published elsewhere (Londono-Renteria et al., 2010) and the test was used to measure IgG antibodies against total SGA and individual peptides in human samples. Briefly, 96-well ELISA plates (UltraCruz^® ELISA Plate, Santa Cruz Biotechnology, Dallas, TX) were coated with 100 μL/well of 0.5 μg/ml of An. darlingi SGA or 2 μg/ml of each peptide in 1X PBS and incubated overnight at 4°C. Plates were blocked for 1 h at room temperature with 2% milk in 1X PBST (blocking buffer) and incubated with 100 μL/well of a 1/100 serum dilution in blocking buffer at 37°C for 2 h. Plates were washed three times with wash solution (1 × PBS and 0.1% Tween) and incubated with 100 μL/well of horseradish peroxidase (HRP)-conjugated goat anti-human IgG (1:1,000), antibodies (Abcam, Cambridge, UK) at 37°C for 1.5 h. Colorimetric development was carried out using tetra-methyl-benzidine (Abcam) as a substrate and incubated for 2 min before the reaction was stopped using 2 N sulfuric acid. Absorbance was measured at 450 nm. Each sample was tested in duplicate. Three controls were included in each plate: (1) control blank: two wells without SGA to control for non-specific induction of color for any of the reagents used in the test; (2) negative control: two wells with SGA but without human serum to control for any non-specific color induction of the coating antigen; and (3) positive control to control for plate to plate variation and normalize OD values. To determine a cut off value of exposure to Plasmodium antigens and An. darlingi, we tested serum from 12 volunteers living in the US with no history of travel to South America. The OD cut off values to establish exposure were set at 0.167 for Pf-MSP1, 0.153 for Pv-MSP1 and 0.396 for An. darlingi SGA.

Data Analysis

After testing for normality, we concluded that the data did not meet the normality requirement. Consequently, we performed non-parametrical statistical tests. Specifically, we used the Mann-Whitney U-test to compare two independent groups (i.e., malaria positive vs. negative) and the Kruskal-Wallis test to compare more than two independent groups (i.e., negative, P. vivax malaria and P. falciparum malaria). Spearman Correlation Coefficients were calculated to measure the strength of association between IgG antibodies against salivary proteins and IgG antibodies against Plasmodium antigens as well as age. Odd ratios and Fisher's exact test were used to assess significance of risk calculation. OD values were categorized as low (less than the median OD) or high (equal or above the median OD) for each peptide: AnDar_PeroX1 (median OD = 0.579), AnDar_PeroX2 (median OD = 0.6612), AnDar_Apy1 (median OD = 0.4410) and the AnDar_Apy2 (median OD = 0.4286). All tests were considered significant with a p-value < 0.05.

Results

Immunogenic Proteins in An. darlingi SGA

Western blot testing revealed four major immunogenic bands (~250, ~62, ~28, and ~8 kDa) (Figure 1). The most prominent band was a 65 kDa band that was excised and sent for sequencing. Seventeen proteins were identified by LC-MS (Table 2). The most abundant proteins found were a ~63 kDa apyrase (Accession# ETN63669.1) and a 65 kDa oxidase/peroxidase (Accession# ETN66035.1). Two peptides were designed from each of these proteins (AnDar_Apy1, AnDar_Apy2, AnDar_PeroX1, and AnDar_PeroX2) to evaluate immune responses to these specific salivary proteins in people with active malaria infections. The specific sequence of each peptide is described in Table 3. Apyrases as well as salivary peroxidases are ubiquitous proteins. The An. darlingi apyrase described in the current study has a 63% identity with an apyrase of Anopheles gambiae (AGAP011971), and 49% with an apyrase from Aedes aegypti (AAEL006347) and an apyrase of Cx. quinquefasciatus (CPIJ011010). No significant similarity was found with An. albimanus apyrases. The An. darlingi salivary oxidase/peroxidase has a 87.8% identity with the An. albimanus peroxidase (AAD22196.1). The An. darlingi oxidase/peroxidase also presents a 53% similarity with An. gambiae (AGAP010735) while similarity with Ae. aegypti (AAEL000507) and Cx. quinquefasciatus (CPIJ017579) was 48%. However, the selected immunogenic peptides align with identities between 56 and 85% with these proteins (Supplementary Table 1).

FIGURE 1

Figure 1. Immunoblot of An. darlingi SGA. Pooled serum from uninfected individuals exposed to An. darlingi bites was used to identify immunogenic proteins. (A) Silver stain of whole SGA. (B) Immunoblot.

Antibody Responses Against Salivary Proteins

Serum from a total of 107 females and 72 males with ages between 2 and 83 years old (Age average = 25.1 years old) were included in this study. No significant differences were detected when comparing IgG antibody levels against any of the antigens from either Turbo or El Bagre samples with the exception of AnDarApy2 which was significantly higher in El Bagre (Mann-Whitney test p = 0.0097) (Figure 2). When comparing IgG antibody levels among groups by their infection status (malaria negative vs. malaria positive), we did not detect significant differences between the groups from either field site using the whole SGA as antigen (Mann-Whitney test p > 0.05) (Figure 3). However, when using the apyrase peptides, we found that antibodies against AnDar-Apy2 were significantly higher in malaria positive patients in both Turbo (Mann-Whitney test, p = 0.0001) and El Bagre (Mann-Whitney test, p < 0.0001), while AnDar-Apy1 was higher in malaria positive samples only in El Bagre (Mann-Whitney test, p < 0.0001). Interestingly, the peroxidase-derived peptides, only showed significant differences between infected and uninfected samples from El Bagre, where AnDar-PeroX1 was significantly higher in malaria infected people (Mann-Whitney test, p = 0.0497). When comparing the antibody levels regardless of the field site, both of the apyrase peptides and AnDarPeroX1 were significantly higher in malaria infected samples.

FIGURE 2

Figure 2. Comparison of IgG antibody levels against An. darlingi SGA (A) and peptides from Apyrase (B) and Peroxidase (C) between the study sites (Turbo and El Bagre). ns denotes not significant, p < 0.05 (* = 0.01 to 0.05, ** = 0.001 to 0.01, and *** = 0.0001 to 0.001).

FIGURE 3

Figure 3. Comparison of IgG antibody levels against An. darlingi salivary antigens in samples with malaria infections and uninfected samples from Turbo (A–C), El Bagre (D–F) and both sites combined (G–I). ns denotes not significant, p < 0.05 (* = 0.01 to 0.05, ** = 0.001 to 0.01, and *** = 0.0001 to 0.001).

Antibody Responses Against Plasmodium and Mosquito Salivary Proteins

Spearman correlation analysis showed significant associations between three of the salivary proteins (AnDar_Apy1, AnDar_Apy2, and AnDar_PeroX2) and either PvMSP1 or PfMSP1 (Spearman rank correlation test, p < 0.05). When divided by malaria infection status, a significant positive association between IgG antibodies against AnDarApy2 and AnDarPeroX2 was maintained with both PvMSP1 and PfMSP1 in both malaria infected and uninfected individuals (Table 4). No association with either Plasmodium antigen was detected when comparing IgG antibodies against AnDarPeroX1. When comparing whether there were significant differences in the IgG antibody levels against any salivary antigens and the species of Plasmodium causing infection, we found that antibodies against both apyrases peptides were significantly higher in people with P. vivax infection (Mann-Whitney test, p < 0.0001) and P. falciparum [Mann-Whitney test, p = 0.0011 (AnDar_Apy1), p = 0.0005 (AnDar_Apy2)] when compared against IgG antibodies in uninfected participants. In the case of IgG antibodies against the peroxidase peptides, we found only significant higher antibodies against AnDar_Perox1 in P. vivax infection (Mann-Whitney test, p = 0.0114 (P. vivax), p = 0.1285 (P. falciparum), while antibodies against AnDar_Perox2 were significantly higher only when comparing P. falciparum infected against uninfected (Mann-Whitney test, p = 0.6464 (P. vivax), p = 0.0131 (P. falciparum). No differences were found with SGE antibodies (Mann-Whitney test, p > 0.05) (Figure 4).

TABLE 4

Table 4. Correlation analysis between IgG antibodies against An. darlingi salivary antigens and IgG antibodies against P. falciparum, P. vivax showing samples grouped by infection status.

FIGURE 4

Figure 4. Comparison of IgG antibody levels against An. darlingi apyrase peptides (A,B), peroxidase peptides (C,D), and SGA (E), in samples from individuals infected with either P. falciparum or P. vivax. ns denotes not significant, p < 0.05 (* = 0.01 to 0.05, ** = 0.001 to 0.01, and *** = 0.0001 to 0.001).

Malaria Risk and IgG Antibodies Against An. darlingi Salivary Proteins

Odd ratios were calculated to determine the risk of malaria infection based on the levels of IgG antibodies against whole SGA and the individual salivary peptides. Our analysis indicated that samples with high levels of antibodies against apyrase peptides had 3.1 (AnDar_Apy1) and 5.1 (AnDar_Apy2) times higher probability of being positive for malaria than samples with lower antibody levels (Fisher's exact test, p < 0.0001). Significant odd ratios were not observed when comparing IgG levels of SGA, AnDar_Perox1 or AnDarPerox2 (Table 5).

TABLE 5

Table 5. Odds ratio calculations to measure risk of malaria infection based on level of IgG antibodies against An. darlingi salivary proteins.

Discussion

An. darlingi is one of the most important malaria vectors in the Americas, and is associated with transmission of both P. falciparum and P. vivax (de Arruda et al., 1986; Flores-Mendoza et al., 2004). Some data suggest that infected An. darlingi tend to bite more outdoors and in the early evening, which reduces the effectiveness of bed nets (Moreno et al., 2015; Newell et al., 2018). The (EIR) is the current gold standard to estimate malaria transmission intensity. However, this technique is labor intensive and requires sensitive methods to accurately detect Plasmodium infection in field-captured mosquitoes. In addition, it is impossible to track individual human exposure to infective bites and estimation of risk that can be later used to determine local transmission dynamics. Tools capable of monitoring human-vector contact are needed to more accurately estimate disease risk. Here, we evaluated the IgG responses against An. darlingi salivary proteins and showed through ELISA that peptides against a salivary apyrase and peroxidase are suitable antigens to estimate the risk of suffering a malaria infection.

In an effort to identify candidate salivary peptides from An. darlingi salivary proteins, an immunogenic ~65 kDa protein band was identified by antibodies in serum collected from study participants. LC-MS revealed 17 salivary proteins, with apyrase and a peroxidase as the most predominant proteins. Most of the identified proteins have previously been associated with blood feeding. Salivary apyrases hydrolyze ATP and ADP inhibiting platelet aggregation. They are associated with host-seeking behavior and the ability of mosquitoes to find blood vessels (Calvo et al., 2006). Infection status with human pathogens influences mosquito behavior and changes the content of salivary gland and salivary secretions (Chisenhall et al., 2014a,b). A previous study suggested that apyrase activity decreases as sporozoite load increases, rendering the infected mosquitoes more likely to probe and re-feed after an infectious first blood meal (Thievent et al., 2019). Each probing event increases the likelihood of injecting sporozoites into the skin, which increases the chance of infection. Salivary proteins are also injected during each probing event (Ribeiro, 2000; Chisenhall et al., 2014a), which may explain the high antibody responses against salivary proteins in people with active malaria infection. Interestingly, immunization with mosquito saliva leads to the development of anti-apyrase antibodies, which can block its enzymatic activity. However, a decrease in feeding or probing time was not observed in immunized mice even in the presence of higher titers of anti-apyrase antibodies (Mathews et al., 1996). Previous studies have shown changes in salivary gland physiology of infected mosquitoes (Cotama et al., 2013; Pinheiro-Silva et al., 2015). It is unclear if infection with either P. vivax or P. falciparum changes the expression profile of apyrase or peroxidase proteins in the salivary gland of An. darlingi mosquitoes. This information would be helpful to determine a potential function of these proteins in malaria transmission by An. darlingi mosquitoes.

Besides the apyrase and the salivary peroxidase described above, other important salivary proteins that had been previously reported were also detected in this study. First, two members of the SG1 family were detected: the putative salivary protein SG-1B along with a SG1-like salivary protein. The SG-1B appears to be uniquely expressed in the female salivary glands, while the SG1-like can also be expressed in males (Arca et al., 2005). Recently, a member of the SG1 family, TRIO protein, was evaluated as a potential vaccine candidate against Plasmodium infection since immunization against these proteins induced partial protection against Plasmodium berghei infection. Protection was further improved when animals were co-immunized with saliva and parasite proteins (Dragovic et al., 2018). Second, a long form D7 protein was also detected. Our previous studies demonstrated the immunogenicity of this family of proteins in Ae. aegypti saliva (Londono-Renteria et al., 2018). The D7 protein family has been identified as a member of an odorant binding protein (OBP) superfamily known to bind biogenic amines such as histamine, playing an important role as mediators of inflammation and vascular permeability, and facilitating blood feeding (Calvo et al., 2006). Finally, a 30 kDa salivary antigen family protein was identified. These proteins are exclusively found in the salivary glands of adult female mosquitoes and previous work suggested that they may inhibit platelet aggregation (Yoshida and Watanabe, 2006; Calvo et al., 2007b). A venom allergen 5 was also identified, and this family of proteins is found ubiquitously in animals and plants. Their specific function in Anopheles mosquitoes is still unknown, but it is thought to function either in the suppression of the host immune system, or as an anticoagulant in a wide range of hematophagous arthropods (Ribeiro and Francischetti, 2003; Calvo et al., 2007a).

It is important to note that whole salivary gland proteins contained in the SGA have been successfully used as antigen to detect IgG antibodies in vertebrate serum (Waitayakul et al., 2006; Badu et al., 2012). However, salivary gland dissection is a tedious process requiring skilled personnel (Coleman et al., 2007) and a continuous source of mosquito salivary glands. Not all mosquitoes are easily colonizable, and although using salivary glands from field mosquitoes can offer a better approximation of the antigens that people are actually exposed to (Andrade et al., 2009), it can be challenging to collect them in sufficient numbers and maintain a proper cold chain and protein stabilization before they are processed for antigen preparation (Fontaine et al., 2011a). Thus, the use of recombinant proteins or peptides is advantageous, particularly when consistency and reproducibility are desired. In that regard, we selected four peptides from the two main immunogenic proteins found after sequencing of the ~65 kDa band, an apyrase and a salivary peroxidase. The ELISA results revealed that when using the AnDar_Apy1, AnDar_Apy2, and AnDar_PeroX1 peptides as antigens, the IgG antibody levels were significantly higher in samples with active malaria infection as compared to uninfected samples. These results are consistent with our recent study where we observed that malaria-infected samples had significantly higher antibody levels against An. albimanus salivary peptides (Londono-Renteria et al., 2020).

Previous studies associated this increase in antibodies with a potentially higher exposure to mosquito bites (Fontaine et al., 2011b). Moreover, this increase in IgG antibodies may also be the result of saliva-dependent activation of immune cells leading to the production of cytokines such as IL4 and IL10 (Vogt et al., 2018). IL4 is involved in the activation of B cells and their differentiation to plasma cells (Granato et al., 2014) and recent studies suggest that IL10 promotes IgG4 antibodies (Jeannin et al., 1998), which is the main antibody subtype that recognizes arthropod salivary antigens (Brummer-Korvenkontio et al., 1994; Cardenas et al., 2019). Further studies are needed to determine the mechanisms involved in the differential antibody profiles between people residing under the same conditions but with different malaria infection status.

In contrast to a previous study in Brazil (Andrade et al., 2009), we did not find significant differences between malaria-infected and uninfected samples when using SGA. One of our previous studies showed the possibility of discrete differences in the responses against Anopheles SGA from two different mosquito colony strains, with those differences potentially associated with differences in geographical origins as well as the time of colonization (Londono-Renteria et al., 2020; Montiel et al., 2020a). Another study detected important differences in salivary gland content from wild and colony maintained arthropods from the same species (Maldonado-Ruiz et al., 2019). Consequently, it is possible that the difference between this current study and the one performed in Brazil may reside in the source of the SGA and its components. Specifically, in this study, we used SGA from a recently colonized strain of An. darlingi in Iquitos, Peru, while the previous study used SGA from field-caught mosquitoes in Brazil (Andrade et al., 2009).

Regardless, analysis of the immunogenic proteins among mosquito strains may reveal shared immunogenic proteins that represent suitable candidates to measure risk of disease irrespective of the geographical area or the strains. Thus, it was important to identify the evolutionarily conserved regions within immunogenic proteins as we seek to develop widely applicable biomarkers of malaria risk. Only AnDar_Apy2 was significantly higher in El Bagre, the area with historically higher abundance of An. darlingi relative to Turbo (Gutierrez et al., 2009, 2010; Rosero et al., 2013); all other peptides did not show differences between the two sites. The lack of difference in antibody levels against most of the peptides between the two study sites may be explained by the 56 and 80% identity with other major salivary proteins representing a source of cross-reactivity. However, since we did not include concurrent mosquito collection data, we cannot confirm the concurrent An. darlingi population density at the time of this study. Thus, the results presented here are more sensitive at correlating antibodies against salivary proteins with presence of malaria infection rather than An. darlingi bite intensity.

Previous studies suggested that there is a significant association between the level of IgG antibodies against the most prominent Plasmodium antigens and the development of human immunity against malaria symptoms (Doolan et al., 2009; Rodriguez-Barraquer et al., 2018). Among these antigens, MSP-1 is considered as one of the more important in the development of immunity against malaria and has been investigated as a potential vaccine candidate. Immunity against malaria is thought to be acquired through chronic and sustained exposure to Plasmodium antigens (Rodriguez-Barraquer et al., 2018). However, recent evidence suggests that immune responses against salivary proteins may also contribute to the development of immunity against disease (Manning et al., 2018). Our current study revealed that the level of IgG antibodies against AnDar_Apy 2 and the AnDar_PeroX2 were significantly associated with IgG levels against both PvMSP1 and PfMSP1 in all study groups, suggesting a positive association between intensity of exposure to malaria antigens and both peroxidase and apyrase salivary proteins. Furthermore, when measuring the risk of malaria infection through odd ratio calculations, a significant increase in risk was only observed in antibodies against the apyrase peptides. Previous study showed an association between the level of antibodies against PfMSP1 and the only currently validated salivary biomarker from the An. gambiae salivary peptide, gSG6-P1 (Badu et al., 2012; Ya-Umphan et al., 2018). Our results suggest that AnDarApy1 and AnDarApy2 may be more suitable biomarkers to estimate malaria risk in areas where An. darlingi mosquitoes are endemic. Also interesting is that the correlation between Plasmodium antigens and the apyrase peptides salivary was also significant in the uninfected population. It is important to remember that infected individuals present significantly higher IgG antibodies. People in endemic areas may encounter a significant number of bites containing infective sporozoites throughout their lives. Furthermore, in endemic areas for P. vivax, as the one where this study was carried out, is very difficult to differentiate a recent infection from a recrudescence, and people may suffer several malaria infections in a year. Since this study did not included follow up of patients, we cannot pinpoint the time of the most recent infection. However, our study showed a significant higher antibody levels against both apyrases peptides and the AnDar_Perox1 in P. vivax infected populations suggesting the relevance of these peptides to measure risk of infection with the most prevalent species circulating in this area. Further studies are urgently needed to determine the persistence of antibodies against the salivary proteins and correlate with previous/recent malaria episodes.

Although further studies including a greater number of serum samples from different geographical areas are needed to validate these peptides as biomarkers of both malaria risk and human-vector contact, the current study represents an important advance in evaluating immune responses against mosquito salivary peptides in Latin America.

Conclusion

Peptides derived from the An. darlingi salivary apyrase and peroxidases are suitable candidate markers for measuring human IgG antibody responses that may be associated with the risk of malaria infections.

Data Availability Statement

All datasets generated for this study in the sequencing data are included in the article and the Supplementary Material. Other raw data will be available upon request sent to blondono@ksu.edu.

Ethics Statement

The studies involving human participants were reviewed and approved by Kansas State University and University of Antioquia. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

BL-R, AL, HV, GV, EC, MF, and MC: study design and data analysis, writing, and reviewing. JM and AT-C: sample collection. KE-V, LR, and MB: mosquito rearing and salivary gland dissection. BL-R and JM: ELISA testing. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the USDA National Institute of Food and Agriculture, hatch multistate project 1021430 and the USAID-CDC - 200-2018-M-00025.

Disclaimer

The views expressed in this article reflect the results of research conducted by the authors and do not necessarily reflect the official policy or position of the Department of the Navy, Department of Defense, the Centers for Disease Control and Prevention, the Department of Health and Human Services, nor the U.S. Government.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was possible thanks to the support provided by the Kansas State University Department of Entomology and the U.S. Naval Medical Research Unit No. 6 (NAMRU-6). Special thanks to the NAMRU-6 Iquitos insectary staff for assisting in An. darlingi rearing. We are grateful to the Ministerio de Agricultura y Riego de Peru, Direccion General Forestal y de Fauna Silvestre for permission to conduct research studies on laboratory reared An. darlingi under the auspices of Resolucion Directoral No. 132-2019-MINAGRI-SERFOR-DGGSPFFS. The authors also want to thank the residents of Turbo and El Bagre for participating in our study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00455/full#supplementary-material

Supplementary Table 1. Immunogenic peptides alignment between An. darlingi and other vectors of human disease.

References

Andrade, B. B., Rocha, B. C., Reis-Filho, A., Camargo, L. M., Tadei, W. P., Moreira, L. A., et al. (2009). Anti-Anopheles darlingi saliva antibodies as marker of Plasmodium vivax infection and clinical immunity in the Brazilian amazon. Malar. J. 8:121. doi: 10.1186/1475-2875-8-121

PubMed Abstract | CrossRef Full Text | Google Scholar

Arca, B., Lombardo, F., Struchiner, C. J., and Ribeiro, J. M. (2017). Anopheline salivary protein genes and gene families: an evolutionary overview after the whole genome sequence of sixteen Anopheles species. BMC Genomics 18:153. doi: 10.1186/s12864-017-3579-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Arca, B., Lombardo, F., Valenzuela, J. G., Francischetti, I. M., Marinotti, O., Coluzzi, M., et al. (2005). An updated catalogue of salivary gland transcripts in the adult female mosquito, Anopheles gambiae. J. Exp. Biol. 208(Pt 20), 3971–3986. doi: 10.1242/jeb.01849

PubMed Abstract | CrossRef Full Text | Google Scholar

Arevalo-Herrera, M., Quinones, M. L., Guerra, C., Cespedes, N., Giron, S., Ahumada, M., et al. (2012). Malaria in selected non-amazonian countries of Latin America. Acta Trop. 121, 303–314. doi: 10.1016/j.actatropica.2011.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Badu, K., Siangla, J., Larbi, J., Lawson, B. W., Afrane, Y., Ong'echa, J., et al. (2012). Variation in exposure to Anopheles gambiae salivary gland peptide (gSG6-P1) across different malaria transmission settings in the western Kenya highlands. Malar. J. 11:318. doi: 10.1186/1475-2875-11-318

PubMed Abstract | CrossRef Full Text | Google Scholar

Brummer-Korvenkontio, H., Lappalainen, P., Reunala, T., and Palosuo, T. (1994). Detection of mosquito saliva-specific IgE and IgG4 antibodies by immunoblotting. J. Allergy Clin. Immunol. 93, 551–555. doi: 10.1016/S0091-6749(94)70066-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Calvo, E., Dao, A., Pham, V. M., and Ribeiro, J. M. (2007a). An insight into the sialome of Anopheles funestus reveals an emerging pattern in anopheline salivary protein families. Insect Biochem. Mol. Biol. 37, 164–175. doi: 10.1016/j.ibmb.2006.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Calvo, E., Mans, B. J., Andersen, J. F., and Ribeiro, J. M. (2006). Function and evolution of a mosquito salivary protein family. J. Biol. Chem. 281, 1935–1942. doi: 10.1074/jbc.M510359200

PubMed Abstract | CrossRef Full Text | Google Scholar

Calvo, E., Tokumasu, F., Marinotti, O., Villeval, J. L., Ribeiro, J. M., and Francischetti, I. M. (2007b). Aegyptin, a novel mosquito salivary gland protein, specifically binds to collagen and prevents its interaction with platelet glycoprotein VI, integrin alpha2beta1, and von willebrand factor. J. Biol. Chem. 282, 26928–26938. doi: 10.1074/jbc.M705669200

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardenas, J. C., Drame, P. M., Luque-Burgos, K. A., Berrio, J. D., Entrena-Mutis, E., Gonzalez, M. U., et al. (2019). IgG1 and IgG4 antibodies against Aedes aegypti salivary proteins and risk for dengue infections. PLoS ONE 14:e0208455. doi: 10.1371/journal.pone.0208455

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassab, A., Morales, V., and Mattar, S. (2011). Climatic factors and cases of dengue in Monteria, Colombia: 2003–2008. Rev. Salud Publica 13, 115–128. doi: 10.1590/S0124-00642011000100010

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaparro-Narvaez, P. E., Lopez-Perez, M., Rengifo, L. M., Padilla, J., Herrera, S., and Arevalo-Herrera, M. (2016). Clinical and epidemiological aspects of complicated malaria in Colombia, 2007–2013. Malar. J. 15:269. doi: 10.1186/s12936-016-1323-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Chisenhall, D. M., Christofferson, R. C., McCracken, M. K., Johnson, A. M., Londono-Renteria, B., and Mores, C. N. (2014a). Infection with dengue-2 virus alters proteins in naturally expectorated saliva of Aedes aegypti mosquitoes. Parasit. Vectors 7:252. doi: 10.1186/1756-3305-7-252

PubMed Abstract | CrossRef Full Text | Google Scholar

Chisenhall, D. M., Londono, B. L., Christofferson, R. C., McCracken, M. K., and Mores, C. N. (2014b). Effect of dengue-2 virus infection on protein expression in the salivary glands of Aedes aegypti mosquitoes. Am. J. Trop. Med. Hyg. 90, 431–437. doi: 10.4269/ajtmh.13-0412

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, J., Juhn, J., and James, A. A. (2007). Dissection of midgut and salivary glands from Ae. aegypti mosquitoes. J. Vis. Exp. 5:228. doi: 10.3791/228

PubMed Abstract | CrossRef Full Text | Google Scholar

Cotama, S., Dekumyoy, P., Samung, Y., and Lek-Uthai, U. (2013). Salivary glands proteins expression of Anopheles dirus a fed on Plasmodium vivax- and Plasmodium falciparum-infected human blood. J. Parasitol. Res. 2013:535267. doi: 10.1155/2013/535267

PubMed Abstract | CrossRef Full Text | Google Scholar

de Arruda, M., Carvalho, M. B., Nussenzweig, R. S., Maracic, M., Ferreira, A. W., and Cochrane, A. H. (1986). Potential vectors of malaria and their different susceptibility to Plasmodium falciparum and Plasmodium vivax in northern Brazil identified by immunoassay. Am. J. Trop. Med. Hyg. 35, 873–881. doi: 10.4269/ajtmh.1986.35.873

PubMed Abstract | CrossRef Full Text | Google Scholar

Doolan, D. L., Dobano, C., and Baird, J. K. (2009). Acquired immunity to malaria. Clin. Microbiol. Rev. 22, 13–36. doi: 10.1128/CMR.00025-08

CrossRef Full Text | Google Scholar

Dragovic, S. M., Agunbiade, T. A., Freudzon, M., Yang, J., Hastings, A. K., Schleicher, T. R., et al. (2018). Immunization with AgTRIO, a protein in Anopheles saliva, contributes to protection against Plasmodium infection in mice. Cell Host Microbe 23, 523–535.e525. doi: 10.1016/j.chom.2018.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Flores-Mendoza, C., Fernandez, R., Escobedo-Vargas, K. S., Vela-Perez, Q., and Schoeler, G. B. (2004). Natural Plasmodium infections in Anopheles darlingi and Anopheles benarrochi (Diptera: culicidae) from eastern Peru. J. Med. Entomol. 41, 489–494. doi: 10.1603/0022-2585-41.3.489

PubMed Abstract | CrossRef Full Text | Google Scholar

Fontaine, A., Pascual, A., Diouf, I., Bakkali, N., Bourdon, S., Fusai, T., et al. (2011a). Mosquito salivary gland protein preservation in the field for immunological and biochemical analysis. Parasit. Vectors 4:33. doi: 10.1186/1756-3305-4-33

PubMed Abstract | CrossRef Full Text | Google Scholar

Fontaine, A., Pascual, A., Orlandi-Pradines, E., Diouf, I., Remoue, F., Pages, F., et al. (2011b). Relationship between exposure to vector bites and antibody responses to mosquito salivary gland extracts. PLoS ONE 6:e29107. doi: 10.1371/journal.pone.0029107

PubMed Abstract | CrossRef Full Text | Google Scholar

Geleta, G., and Ketema, T. (2016). Severe malaria associated with Plasmodium falciparum and P. vivax among children in pawe hospital, Northwest Ethiopia. Malar. Res. Treat. 2016:1240962. doi: 10.1155/2016/1240962

PubMed Abstract | CrossRef Full Text | Google Scholar

Granato, A., Hayashi, E. A., Baptista, B. J., Bellio, M., and Nobrega, A. (2014). IL-4 regulates bim expression and promotes B cell maturation in synergy with BAFF conferring resistance to cell death at negative selection checkpoints. J. Immunol. 192, 5761–5775. doi: 10.4049/jimmunol.1300749

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutierrez, L. A., Gomez, G. F., Gonzalez, J. J., Castro, M. I., Luckhart, S., Conn, J. E., et al. (2010). Microgeographic genetic variation of the malaria vector Anopheles darlingi root (Diptera: culicidae) from Cordoba and Antioquia, Colombia. Am. J. Trop. Med. Hyg. 83, 38–47. doi: 10.4269/ajtmh.2010.09-0381

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutierrez, L. A., Gonzalez, J. J., Gomez, G. F., Castro, M. I., Rosero, D. A., Luckhart, S., et al. (2009). Species composition and natural infectivity of anthropophilic Anopheles (Diptera: culicidae) in the states of Cordoba and Antioquia, Northwestern Colombia. Mem. Inst. Oswaldo Cruz 104, 1117–1124. doi: 10.1590/S.0074-02762009000800008

PubMed Abstract | CrossRef Full Text | Google Scholar

Instituto Nacional de Salud (2019). Periodo Epidemiológico XIII. Bogota: Instituto Nacional de Salud.

Google Scholar

Jeannin, P., Lecoanet, S., Delneste, Y., Gauchat, J. F., and Bonnefoy, J. Y. (1998). IgE versus IgG4 production can be differentially regulated by IL-10. J. Immunol. 160, 3555–3561.

PubMed Abstract | Google Scholar

Jitthai, N. (2013). Migration and malaria. Southeast Asian J. Trop. Med. Public Health 44(Suppl 1), 166–200.

PubMed Abstract | Google Scholar

Londono-Renteria, B., Cardenas, J. C., Cardenas, L. D., Christofferson, R. C., Chisenhall, D. M., Wesson, D. M., et al. (2013). Use of anti-Aedes aegypti salivary extract antibody concentration to correlate risk of vector exposure and dengue transmission risk in Colombia. PLoS ONE 8:e81211. doi: 10.1371/journal.pone.0081211

PubMed Abstract | CrossRef Full Text | Google Scholar

Londoño-Rentería, B., Cárdenas, J. C., Giovanni, J. E., Cárdenas, L., Villamizar, P., Rolón, J., et al. (2015a). Aedes aegypti anti-salivary gland antibody concentration and dengue virus exposure history in healthy individuals living in an endemic area in Colombia. Biomedica 35, 572–581. doi: 10.7705/biomedica.v35i4.2530

PubMed Abstract | CrossRef Full Text | Google Scholar

Londono-Renteria, B., Drame, P. M., Montiel, J., Vasquez, A. M., Tobon-Castano, A., Taylor, M., et al. (2020). Identification and pilot evaluation of salivary peptides from Anopheles albimanus as biomarkers for bite exposure and malaria infection in Colombia. Int. J. Mol. Sci. 21:691. doi: 10.3390/ijms21030691

PubMed Abstract | CrossRef Full Text | Google Scholar

Londono-Renteria, B., Drame, P. M., Weitzel, T., Rosas, R., Gripping, C., Cardenas, J. C., et al. (2015b). An. gambiae gSG6-P1 evaluation as a proxy for human-vector contact in the Americas: a pilot study. Parasit. Vectors 8:533. doi: 10.1186/s13071-015-1160-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Londono-Renteria, B., Patel, J. C., Vaughn, M., Funkhauser, S., Ponnusamy, L., Grippin, C., et al. (2015c). Long-lasting permethrin-impregnated clothing protects against mosquito bites in outdoor workers. Am. J. Trop. Med. Hyg. 93, 869–874. doi: 10.4269/ajtmh.15-0130

PubMed Abstract | CrossRef Full Text | Google Scholar

Londono-Renteria, B. L., Eisele, T. P., Keating, J., James, M. A., and Wesson, D. M. (2010). Antibody response against Anopheles albimanus (Diptera: culicidae) salivary protein as a measure of mosquito bite exposure in Haiti. J. Med. Entomol. 47, 1156–1163. doi: 10.1603/ME09240

PubMed Abstract | CrossRef Full Text | Google Scholar

Londono-Renteria, B. L., Shakeri, H., Rozo-Lopez, P., Conway, M. J., Duggan, N., Jaberi-Douraki, M., et al. (2018). Serosurvey of human antibodies recognizing Aedes aegypti D7 salivary proteins in Colombia. Front. Public Health 6:111. doi: 10.3389/fpubh.2018.00111

PubMed Abstract | CrossRef Full Text | Google Scholar

Maldonado-Ruiz, L. P., Montenegro-Cadena, L., Blattner, B., Menghwar, S., Zurek, L., and Londono-Renteria, B. (2019). Differential tick salivary protein profiles and human immune responses to lone star ticks (Amblyomma americanum) from the wild vs. a laboratory colony. Front. Immunol. 10:1996. doi: 10.3389/fimmu.2019.01996

PubMed Abstract | CrossRef Full Text | Google Scholar

Manning, J. E., Morens, D. M., Kamhawi, S., Valenzuela, J. G., and Memoli, M. (2018). Mosquito saliva: the hope for a universal arbovirus vaccine? J. Infect. Dis. 218, 7–15. doi: 10.1093/infdis/jiy179

PubMed Abstract | CrossRef Full Text | Google Scholar

Martens, P., and Hall, L. (2000). Malaria on the move: human population movement and malaria transmission. Emerg. Infect. Dis. 6, 103–109. doi: 10.3201/eid0602.000202

PubMed Abstract | CrossRef Full Text | Google Scholar

Mathews, G. V., Sidjanski, S., and Vanderberg, J. P. (1996). Inhibition of mosquito salivary gland apyrase activity by antibodies produced in mice immunized by bites of Anopheles stephensi mosquitoes. Am. J. Trop. Med. Hyg. 55, 417–423. doi: 10.4269/ajtmh.1996.55.417

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirabello, L., and Conn, J. E. (2006). Molecular population genetics of the malaria vector Anopheles darlingi in Central and South America. Heredity 96, 311–321. doi: 10.1038/sj.hdy.6800805

PubMed Abstract | CrossRef Full Text | Google Scholar

Montiel, J., Carbal, L. F., Tobon-Castano, A., Vasquez, G. M., Fisher, M. L., and Londono-Renteria, B. (2020a). IgG antibody response against Anopheles salivary gland proteins in asymptomatic Plasmodium infections in Narino, Colombia. Malar. J. 19:42. doi: 10.1186/s12936-020-3128-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Montiel, J., Zuluaga, L. M., Aguirre, D. C., Segura, C., Tobon-Castano, A., and Vasquez, A. M. (2020b). Microscopic and submicroscopic Plasmodium infections in indigenous and non-indigenous communities in Colombia. Malar. J. 19:157. doi: 10.1186/s12936-020-03226-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Montoya-Lerma, J., Solarte, Y. A., Giraldo-Calderon, G. I., Quinones, M. L., Ruiz-Lopez, F., Wilkerson, R. C., et al. (2011). Malaria vector species in Colombia: a review. Mem. Inst. Oswaldo Cruz 106(Suppl. 1), 223–238. doi: 10.1590/S0074-02762011000900028

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreno, M., Saavedra, M. P., Bickersmith, S. A., Lainhart, W., Tong, C., Alava, F., et al. (2015). Implications for changes in Anopheles darlingi biting behaviour in three communities in the peri-iquitos region of Amazonian Peru. Malar. J. 14:290. doi: 10.1186/s12936-015-0804-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Moutinho, P. R., Gil, L. H., Cruz, R. B., and Ribolla, P. E. (2011). Population dynamics, structure and behavior of Anopheles darlingi in a rural settlement in the Amazon rainforest of Acre, Brazil. Malar. J. 10:174. doi: 10.1186/1475-2875-10-174

PubMed Abstract | CrossRef Full Text | Google Scholar

Newell, I., Wiskin, C., Anthoney, J., Meza, G., and de Wildt, G. (2018). Preventing malaria in the Peruvian Amazon: a qualitative study in Iquitos, Peru. Malar. J. 17:31. doi: 10.1186/s12936-018-2177-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinheiro-Silva, R., Borges, L., Coelho, L. P., Cabezas-Cruz, A., Valdes, J. J., do Rosario, V., et al. (2015). Gene expression changes in the salivary glands of Anopheles coluzzii elicited by Plasmodium berghei infection. Parasit. Vectors 8:485. doi: 10.1186/s13071-015-1079-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Poinsignon, A., Cornelie, S., Ba, F., Boulanger, D., Sow, C., Rossignol, M., et al. (2009). Human IgG response to a salivary peptide, gSG6-P1, as a new immuno-epidemiological tool for evaluating low-level exposure to Anopheles bites. Malar. J. 8:198. doi: 10.1186/1475-2875-8-198

PubMed Abstract | CrossRef Full Text | Google Scholar

Prado, C. C., Alvarado-Cabrera, L. A., Camargo-Ayala, P. A., Garzon-Ospina, D., Camargo, M., Soto-De Leon, S. C., et al. (2019). Behavior and abundance of Anopheles darlingi in communities living in the Colombian Amazon riverside. PLoS ONE 14:e0213335. doi: 10.1371/journal.pone.0213335

PubMed Abstract | CrossRef Full Text | Google Scholar

Reinbold-Wasson, D. D., Sardelis, M. R., Jones, J. W., Watts, D. M., Fernandez, R., Carbajal, F., et al. (2012). Determinants of Anopheles seasonal distribution patterns across a forest to periurban gradient near iquitos, Peru. Am. J. Trop. Med. Hyg. 86, 459–463. doi: 10.4269/ajtmh.2012.11-0547

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribeiro, J. M. (2000). Blood-feeding in mosquitoes: probing time and salivary gland anti-haemostatic activities in representatives of three genera (Aedes, Anopheles, Culex). Med. Vet. Entomol. 14, 142–148. doi: 10.1046/j.1365-2915.2000.00227.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribeiro, J. M., and Francischetti, I. M. (2003). Role of arthropod saliva in blood feeding: sialome and post-sialome perspectives. Annu. Rev. Entomol. 48, 73–88. doi: 10.1146/annurev.ento.48.060402.102812

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodrigues, P. T., Valdivia, H. O., de Oliveira, T. C., Alves, J. M. P., Duarte, A., Cerutti-Junior, C., et al. (2018). Human migration and the spread of malaria parasites to the New World. Sci. Rep. 8:1993. doi: 10.1038/s41598-018-19554-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez, J. C., Uribe, G. A., Araujo, R. M., Narvaez, P. C., and Valencia, S. H. (2011). Epidemiology and control of malaria in Colombia. Mem. Inst. Oswaldo Cruz 106(Suppl. 1), 114–122. doi: 10.1590/S0074-02762011000900015

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez-Barraquer, I., Arinaitwe, E., Jagannathan, P., Kamya, M. R., Rosenthal, P. J., Rek, J., et al. (2018). Quantification of anti-parasite and anti-disease immunity to malaria as a function of age and exposure. eLife 7:e35832. doi: 10.7554/eLife.35832.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosero, D. A., Naranjo-Diaz, N., Alvarez, N., Cienfuegos, A. V., Torres, C., Luckhart, S., et al. (2013). Colombian Anopheles triannulatus (Diptera: culicidae) naturally infected with Plasmodium spp. ISRN Parasitol. 2013:927453. doi: 10.5402/2013/927453

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiz, D., Poveda, G., Velez, I. D., Quinones, M. L., Rua, G. L., Velasquez, L. E., et al. (2006). Modelling entomological-climatic interactions of Plasmodium falciparum malaria transmission in two Colombian endemic-regions: contributions to a national malaria early warning system. Malar. J. 5:66. doi: 10.1186/1475-2875-5-66

PubMed Abstract | CrossRef Full Text | Google Scholar

Sagna, A. B., Gaayeb, L., Sarr, J. B., Senghor, S., Poinsignon, A., Boutouaba-Combe, S., et al. (2013). Plasmodium falciparum infection during dry season: IgG responses to Anopheles gambiae salivary gSG6-P1 peptide as sensitive biomarker for malaria risk in Northern Senegal. Malar. J. 12:301. doi: 10.1186/1475-2875-12-301

PubMed Abstract | CrossRef Full Text | Google Scholar

Thievent, K., Zilio, G., Hauser, G., and Koella, J. C. (2019). Malaria load affects the activity of mosquito salivary apyrase. J. Insect. Physiol. 116, 10–16. doi: 10.1016/j.jinsphys.2019.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Villarreal-Trevino, C., Vasquez, G. M., Lopez-Sifuentes, V. M., Escobedo-Vargas, K., Huayanay-Repetto, A., Linton, Y. M., et al. (2015). Establishment of a free-mating, long-standing and highly productive laboratory colony of Anopheles darlingi from the Peruvian Amazon. Malar. J. 14:227. doi: 10.1186/s12936-015-0733-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogt, M. B., Lahon, A., Arya, R. P., Kneubehl, A. R., Spencer Clinton, J. L., Paust, S., et al. (2018). Mosquito saliva alone has profound effects on the human immune system. PLoS Negl. Trop. Dis. 12:e0006439. doi: 10.1371/journal.pntd.0006439

PubMed Abstract | CrossRef Full Text | Google Scholar

Waitayakul, A., Somsri, S., Sattabongkot, J., Looareesuwan, S., Cui, L., and Udomsangpetch, R. (2006). Natural human humoral response to salivary gland proteins of Anopheles mosquitoes in Thailand. Acta Trop. 98, 66–73. doi: 10.1016/j.actatropica.2006.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ya-Umphan, P., Cerqueira, D., Cottrell, G., Parker, D. M., Fowkes, F. J. I., Nosten, F., et al. (2018). Anopheles salivary biomarker as a proxy for estimating Plasmodium falciparum malaria exposure on the thailand-myanmar border. Am. J. Trop. Med. Hyg. 99, 350–356. doi: 10.4269/ajtmh.18-0081

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, S., and Watanabe, H. (2006). Robust salivary gland-specific transgene expression in Anopheles stephensi mosquito. Insect. Mol. Biol. 15, 403–410. doi: 10.1111/j.1365-2583.2006.00645.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zain Ul, A., Qadeer, A., Akhtar, A., and Rasheed, A. (2016). Severe acute respiratory distress syndrome secondary to Plasmodium vivax malaria. J. Pak. Med. Assoc. 66, 351–353.

PubMed Abstract | Google Scholar