1. Introduction
Malaria and dengue fever are two of the most important vector-borne diseases of public health concern in tropical areas around the globe[
1]. Dengue virus (DENV) has periodic outbreaks in Central and South America[
2], with occasional outbreaks in the US, particularly in Texas and Florida[
3,
4]. Importantly, several recent reports have indicated the occurrence of local malaria transmission in the continental US, while this disease remains endemic in Latin America[
5,
6].
Mosquito-borne pathogens are usually deposited in the skin of the vertebrate host along with arthropod saliva during blood-feeding[
7]. Compelling evidence suggests that mosquito saliva induces profound changes in immune responses locally (at the bite site) and systemically, eventually providing a vehicle for pathogen transmission[
8,
9]. Furthermore, salivary secretions are well known to induce significant antibody responses in the vertebrate host. These antibody responses may vary depending on factors such as age, seasonality, and vector abundance[
10,
11]. Several studies suggest that IgG responses to mosquito salivary proteins may serve as surrogate biomarkers for exposure to mosquito bites and an indirect marker for disease risk in travelers and individuals living in endemic areas[12-14].
We and others have previously shown that not all proteins in mosquito saliva are associated with disease risk[
13], and that only a few proteins are suitable as markers for vector-bite exposure[
15,
16]. However, a small number of proteins may be useful in determining susceptibility to infection. For instance, we recently identified a D7 Long (D7L) (AAEL006424) salivary protein from
Ae. aegypti that was capable of physically binding virions and inhibiting DENV infection in vitro and
in vivo [
17]. Interestingly, measurement of IgG antibodies against this D7L protein revealed higher anti-D7L antibody levels in people with active DENV infection compared to febrile DENV uninfected individuals and healthy participants[
17]. In contrast, antibodies against AgBR1 and NeST1, two immunomodulatory mosquito salivary proteins, were higher in healthy individuals compared to DENV-infected individuals[
18]. These data suggest that targeting the correct salivary proteins in vaccine formulations to induce specific antibodies may prevent or modify the course of mosquito-borne infections.
Immune responses to mosquito saliva can vary according to epidemiological factors such as age, sex, seasonality, use of vector control, etc. Interestingly, previous studies found higher IgG antibody levels after mosquito exposure in males compared to females[
19], highlighting a potential sex-dependent incidence of mosquito-transmitted diseases. Such findings have been partially explained by physiological factors, host genetics, and gender-related social determinants, resulting in differences in exposure[
20,
21]. With the increase in global mosquito-borne infections, a better understanding of sex-dependent host responses may be critical to mitigating the negative consequences of vector-pathogen transfer and hormone-related antigen responses.
Currently, most arthropod-borne disease control relies heavily on decreasing human-vector contact through physical devices (i.e., bed nets) or insecticide treatment[
22,
23]. Effective vaccines or drugs against arthropod-borne diseases are scarce[
24,
25]. Importantly, the increase in insecticide and drug resistance calls for the design and implementation of new tools for disease control, and new protocols to measure exposure to arthropod vectors and estimate the potential risk of acquisition of vector-borne diseases, which will guide public health policy.
Recently, an Ae. aegypti peptide, Nterm-34kDa, was evaluated as a quantitative measure of exposure to Aedes bites, and a positive correlation was observed between intensity of exposure, mosquito abundance, and anti-Nterm-34kDa IgG antibody levels[
26] [
27]
. Also, the gSG6-P1 peptide, identified from the An. gambiae SG6 protein[
28]
has been extensively validated as a biomarker of exposure against Anopheles mosquitoes from the subgenus Cellia and Anopheles[
29]
. Both peptides have been used successfully to determine the level of exposure to mosquito bites associated with mosquito control interventions or risk of disease[
15,
30,
31]
.
In tropical areas, people can be exposed to hundreds of mosquito bites per day, and an increase in mosquito abundance is often associated with an increase in vector-borne pathogen transmission[
32]. Importantly, mosquito species are very particular in their feeding behavior. For instance, Anopheles mosquitoes are preferential nocturnal feeders, while Aedes mosquitoes are diurnal biters. Therefore, it is possible that an individual may sustain Aedes mosquito bites during the day, followed by Anopheles mosquito bites at night. To our knowledge, there are currently no studies that describe the effect of exposure to bites from different mosquito species on skin immune responses and pathogen replication. However, prior work suggests that saliva induces skin responses associated with the potential for pathogen establishment [
33]. Also, recent studies suggest that saliva from different mosquito species may have variable impacts on the same pathogen[
34]. Thus, there is a critical need to better understand how contact with salivary proteins from different mosquito species due to sequential exposure to diurnal and nocturnal biters may impact human arbovirus acquisition and anti-viral immune responses.
In this study, our primary goal was to evaluate the levels of antibodies against salivary antigens of different arthropod vectors of human disease in people with either malaria or dengue to assess whether responses against vector saliva are associated with blood parameters leading to severe clinical presentation. Since it is rare that only a single species of mosquito is found in a specific area, and several arthropod-borne diseases are usual in the tropics, we hypothesized that exposure to salivary proteins with different immunomodulatory properties would impact infection and progression to disease. Malaria and dengue fever transmission often co-occur in tropical areas, mainly because of the overlap in the ecological niches preferred by the main vectors of these infections[35-37]. We therefore leveraged our ongoing malaria and dengue surveillance study in Norte de Santander to evaluate exposure to Anopheles and Aedes mosquito saliva and compare these data with blood parameters at diagnosis. To our knowledge, this is the first study exploring a potential correlation between exposure to saliva of non-vectors in clinical presentation of disease and severity.
2. Results
2.1. P. vivax Malaria and Exposure to the Non-Malaria vector Aedes Aegypti
A total of 49 participants with current P. vivax malaria infection, confirmed by microscopy and rapid diagnostic test (RDT), from the areas of Tarra and Tibu in Norte de Santander, Colombia, were included in the study from 2018 to 2019 (Table 1). Mean parasite count was 6,665 parasites/µL (from 420 to 26,480 parasites/µL). We found that antibody levels were not correlated with parasite count or expression of the Pvs25 gene (gametes, ookinetes); however, IgG antibodies against An. albimanus SGE (Spearman correlation r=-0.6099, p=0.0269), Trans-1 (ρ=-0.7510, p=0.0031), gsG6-P1 (ρ=-0.6648, p=0.0132) and Nterm24kDa (Spearman correlation r=-0.7253, p=0.0050) showed a significant negative correlated with the level of expression of Pvs230 (gametocytes) (Table 2).
Table 1.
Description of the study population age in the malaria endemic areas (Tibu and Tarra, Norte de Santander – Colombia), dengue endemic area (Los Patios and Cucuta, Norte de Santander – Colombia); and healthy individuals from Los Patios (Colombia) and Manhattan – Kansas (USA).
Table 1.
Description of the study population age in the malaria endemic areas (Tibu and Tarra, Norte de Santander – Colombia), dengue endemic area (Los Patios and Cucuta, Norte de Santander – Colombia); and healthy individuals from Los Patios (Colombia) and Manhattan – Kansas (USA).
Infection Status |
All age years (range) |
Female (range) |
Male (Range) |
Malaria |
33.4 (1 – 67), n=49 |
36.0 (16 – 53), n=14 |
32.4 (1 – 67), n=35 |
Dengue |
15.3 (1 – 76), n=124 |
17.6 (1 – 76), n=70 |
12.4 (1 – 69), n=54 |
Healthy |
28.9 (2 – 79), n=103 |
27.0 (2 – 79), n=65 |
32.3 (2 – 72), n=38 |
Table 2.
Correlation analysis between IgG antibody responses against each peptide and parasite count/gametocytemia by gender in Plasmodium positive volunteers. Data is presented in Spearman correlation r.
Table 2.
Correlation analysis between IgG antibody responses against each peptide and parasite count/gametocytemia by gender in Plasmodium positive volunteers. Data is presented in Spearman correlation r.
Peptide |
Pvs25 |
Pvs230 |
Parasite count |
All |
|
|
|
Peroxi-P1 |
0.0121 (p=0.9415) |
-0.2842 (p=0.0795) |
-0.0989 (p=0.5492) |
Trans-1 |
0.0715 (p=0.6655) |
-0.1355 (p=0.4107) |
-0.1987 (p=0.2254) |
Trans-2 |
-0.1937 (p=0.2373) |
-0.0644 (p=0.6970) |
-0.0984 (p=0.5511) |
An. albimanus SGE |
-0.0497 (p=0.7638) |
-0.3072 (p=0.0571) |
-0.0538 (p=0.7451) |
gSG6-P1 |
-0.0901 (p=0.5855) |
-0.2913 (p=0.0720) |
-0.1123 (p=0.4960) |
Nterm-34kDa |
-0.0109 (p=0.9473) |
-0.1553 (p=0.3453) |
-0.0314 (p=0.8495) |
Females |
|
|
|
Peroxi-P1 |
-0.0824 (p=0.7890) |
-0.5440 (p=0.0546) |
-0.2418 (p=0.4262) |
Trans-1 |
-0.3989 (p=0.1770) |
-0.7510 (p=0.0031)
|
-0.1761 (p=0.5650) |
Trans-2 |
0.1978 (p=0.5171) |
-0.3187 (p=0.2886) |
-0.0165 (p=0.9574) |
An. albimanus SGE |
-0.2802 (p=0.3538) |
-0.6099 (p=0.0269)
|
-0.1593 (p=0.6031) |
gSG6-P1 |
-0.3132 (p=0.2974) |
-0.6648 (p=0.0132)
|
-0.1593 (p=0.6031) |
Nterm-34kDa |
-0.5000 (p=0.0819) |
-0.7253 (p=0.0050)
|
-0.3077 (p=0.3064) |
Males |
|
|
|
Peroxi-P1 |
0.0510 (p=0.8047) |
-0.1772 (p=0.3866) |
0.0395 (p=0.8480) |
Trans-1 |
0.2113 (p=0.3001) |
0.1217 (p=0.5356) |
-0.1498 (p=0.4651) |
Trans-2 |
0.2205 (p=0.2790) |
0.1022 (p=0.6193) |
-0.0728 (p=0.7236) |
An. albimanus SGE |
0.0168 (p=0.9353) |
-0.1733 (p=0.3971) |
-0.0147 (p=0.9432) |
gSG6-P1 |
-0.0544 (p=0.7919) |
-0.1433 (p=0.4850) |
-0.0332 (p=0.8722) |
Nterm-34kDa |
0.1829 (p=0.3711) |
0.0715 (p=0.7287) |
0.2151 (p=0.2913) |
We next tested the correlation between IgG antibody levels against the salivary proteins with blood parameters such as red blood cell (RBC) count, white blood cell (WBC) count, platelet count, hemoglobin, hematocrit, and parasite count (Table 3). Interestingly, we found that levels of antibodies against the Aedes peptide Nterm34kDa (non-malaria vector) were significantly negatively correlated with RBC count in males, while females presented a significant positive correlation between anti-Nterm34kDa antibodies and WBC count. We also compared antibody levels with age and observed a positive correlation with the level of antibodies against whole An. albimanus SGE, gSG6-P1, and Trans1 in females but not males (Table 4). When measuring odds ratios, we found that malaria patients with normal leucocyte counts were 3.3 times more likely to present high IgG antibodies against the gSG6-P1 peptide (Figure 1A), and people with low parasitemia were 3.3 times more likely to have low antibodies against Trasferrin 2 (Figure 1B).
Figure 1.
Odd ratios comparing levels of IgG antibodies against gSG6-P1 and Transferrin 2 peptide in people with and without leukemia (A) and low and high parasitemia levels (B).
Figure 1.
Odd ratios comparing levels of IgG antibodies against gSG6-P1 and Transferrin 2 peptide in people with and without leukemia (A) and low and high parasitemia levels (B).
Table 3.
Correlation analysis between IgG antibody responses against each peptide and blood parameters by gender in Plasmodium positive volunteers. Data is presented in Spearman correlation r.
Table 3.
Correlation analysis between IgG antibody responses against each peptide and blood parameters by gender in Plasmodium positive volunteers. Data is presented in Spearman correlation r.
Peptide |
Red blood cell count |
White blood cell count |
Platelet count |
Haemoglobin |
Haematocrit |
All |
|
|
|
|
|
Peroxi-P1 |
-0.2660 (p=0.0847) |
-0.0269 (p=0.8642) |
-0.3167 (p=0.0385)
|
-0.1691 (p=0.2782) |
-0.2145 (p=0.1671) |
Trans-1 |
-0.2462 (p=0.1115) |
0.2950 (p=0.0548) |
-0.3288 (p=0.0314)
|
-0.1691 (p=0.2782) |
-0.2445 (p=0.1141) |
Trans-2 |
-0.2651 (p=0.0858) |
0.1180 (p=0.4510) |
-0.3651 (p=0.0161)
|
-0.2651 (p=0.0974) |
-0.3740 (p=0.0135)
|
An. albimanus SGE |
-0.1076 (p=0.4923) |
0.1589 (p=0.3088) |
-0.311 (p=0.0420)
|
-0.1158 (p=0.4597) |
-0.1609 (p=0.3028) |
gSG6-P1 |
-2769 (p=0.0723) |
0.3205 (p=0.0361)
|
-0.2636 (p=0.0877) |
-0.1189 (p=0.4477) |
-0.1300 (p=0.4061) |
Nterm-34kDa |
-0.3288 (p=0.0313)
|
0.2363 (p=0.1271) |
-0.2078 (p=0.1812) |
-0.1402 (p=0.3699) |
-0.1374 (p=0.3797) |
Females |
|
|
|
|
|
Peroxi-P1 |
-0.3714 (p=0.1910) |
0.6497 (p=0.0119)
|
-0.4330 (p=1220) |
-0.1454 (p=0.62000 |
-0.2571 (p=0. 3748) |
Trans-1 |
-0.1958 (p=0.5023) |
0.7522 (p=0.0019)
|
-0.5391 (p=0.0467)
|
0.1235 (p=0.6741) |
-0.0726 (p=0.8052) |
Trans-2 |
-0.3099 (p=0.2809) |
0.6144 (p=0.0194)
|
-0.2308 (p=0.4273) |
-0.1410 (p=0.6307) |
-0.3099 (p=0.2809) |
An. albimanus SGE |
-0.0330 (p=0.9109) |
0.6763 (p=0.0079)
|
-0.4857 (p=0.0783) |
0.1828 (p=0.5316) |
0.0857 (p=0.7708) |
gSG6-P1 |
-0.0857 (p=0.7708) |
0.6188 (p=0.0153)
|
-0.3626 (p=0.2026) |
0.1564 (p=0.5834) |
0.1121 (p=0.7028) |
Nterm-34kDa |
-01560 (p=0.5942) |
0.5923 (p=0.0256)
|
-0.5560 (p=0.0389)
|
-0.0683 (p=0.8166) |
-0.0989 (p=0.7366) |
Males |
|
|
|
|
|
Peroxi-P1 |
-0.3224 (p=0.0880) |
-0.2376 (p=0.2145) |
-0.2266 (p=0.2371) |
-0.2902 (p=0.1268) |
-0.2572 (p=0.1780) |
Trans-1 |
-0.2645 (p=0.1656) |
0.1199 (p=0.5356) |
-0.02303 (p=0.2294) |
-0.2740 (p=0.1504) |
-0.3069 (p=0.1053) |
Trans-2 |
-0.2334 (p=0.2230) |
-0.0712 (p=0.7135) |
-0.3548 (p=0.0589) |
-0.2930 9p=0.1229) |
-0.3972 (p=0.0329)
|
An. albimanus SGE |
-0.1429 (p=0.4595) |
0.0094 (p=0.9615) |
-0.1947 (p=0.3116) |
-0.2169 (p=0.2584) |
-0.2341 (p=0.2216) |
gSG6-P1 |
-0.3300 (p=0.0804) |
0.2192 (p=0.2534) |
-0.1062 (p=0.5834) |
-0.1468 (p=0.4443) |
-0.1417 (p=0.4634) |
Nterm-34kDa |
-0.4171 (p=0.0244)
|
0.0561 (p=0.7727) |
-0.0017 (p=0.9929) |
-0.1594 (p=0.4090) |
-0.1264 (p=0.5135) |
Table 4.
Correlation analysis between IgG antibody responses against each peptide and age by gender in in P. vivax infected patients. Data is presented in Spearman correlation r.
Table 4.
Correlation analysis between IgG antibody responses against each peptide and age by gender in in P. vivax infected patients. Data is presented in Spearman correlation r.
Peptide |
Age |
All |
|
Peroxi-P1 |
-0.0207 (p=0.8880) |
Trans-1 |
0.3270 (p=0.0218)
|
Trans-2 |
0.1596 (p=2733) |
An. albimanus SGE |
0.0471 (p=0.07478) |
gSG6-P1 |
0.2064 (p=0.1547) |
Nterm-34kDa |
0.0828 (p=0.5719) |
Females |
|
Peroxi-P1 |
0.3645 (p=0.2001) |
Trans-1 |
0.6196 (p=0.0181)
|
Trans-2 |
0.3934 (p=0.1641) |
An. albimanus SGE |
0.5934 (p=0.0253)
|
gSG6-P1 |
0.5334 (p=0.0495)
|
Nterm-34kDa |
0.2356 (p=0.4175) |
Males |
|
Peroxi-P1 |
-0.1172 (p=0.5027) |
Trans-1 |
0.2603 (p=0.1310) |
Trans-2 |
0.0793 (p=0.6507) |
An. albimanus SGE |
-0.0771 (p=0.6599) |
gSG6-P1 |
0.1113 (p=0.5243) |
Nterm-34kDa |
0.0417 (p=0.8121) |
2.2. DENV and Exposure to the Non-DENV vector Anopheles Albimanus
A total of 124 DENV-positive volunteers living in Los Patios (n=75), Cucuta (n=24), and Ocana (n=25) between October 2018 and September 2020 were included in this study. Following DENV classification according to WHO guidelines, we included 63 DENV patients with warning signs and 61 DENV patients without warning signs. No significant differences were observed in the level of IgG antibody levels against Nterm34kDa, AnDarApy1, and gSG6-P1 peptides when comparing DENV groups Mann-Whitney test, p>0.05) (Figure 2). However, we observed a significant positive correlation between RBC count and IgG antibody levels against Nterm34kDa (Spearman correlation sp=-0.2107, p=0.0193) and gSG6-P1 (Sρ=-0.1807, p=0.0455). Significant associations between the levels of antibodies against salivary peptides and blood parameters were observed in males but not females. All comparisons are found in Table 5.
Figure 2.
Scatterplot of IgG antibodies against the mosquito salivary peptides Nterm34kDa, AnDarApy1 and gSG6-P1 in people with dengue fever with warning signs (WWS) and dengue without warning signs. Mann-Whitney test, p>0.05.
Figure 2.
Scatterplot of IgG antibodies against the mosquito salivary peptides Nterm34kDa, AnDarApy1 and gSG6-P1 in people with dengue fever with warning signs (WWS) and dengue without warning signs. Mann-Whitney test, p>0.05.
2.3. Healthy Individuals from Endemic and Nonendemic Areas
In 2018, we recruited participants in Manhattan, Kansas for a study focused on evaluating exposure to blood-sucking arthropods. A total of 27 participants donated blood samples in the summer (June – August) and fall (September – October). Our analyses revealed a significant reduction in antibody levels against the Nterm-34kDa, but not for the gSG6-P1 in the Fall (Figure 3), suggesting an association between the intensity of exposure to mosquito bites and IgG antibody levels against mosquito salivary peptides. Next, we compared antibody levels in healthy US participants with levels in healthy individuals living in areas with endemic DENV and malaria. In 2018, samples were collected from 54 healthy volunteers living in houses where a DENV case was reported in Los Patios – Norte de Santander (Colombia). We observed a significant negative correlation between age and IgG antibodies against Nterm34kDa (ρ=0.4182, p=0.000) and gSG6-P1 (ρ=-0.3553, p=0.0003). This significant negative correlation between age and anti-Nterm34kDa remained even when stratifying the data by gender (Table 5). When stratifying data by location, we observed that people from the US demonstrated significant negative correlations between antibody levels against both peptides and age, while healthy individuals from Colombia showed a positive correlation between age and gSG6-P1, but not Nterm-34kDa (Table 6).
Table 5.
Correlation analysis between IgG antibody responses against each peptide and blood parameters by gender. In dengue fever patients. Data is presented in Spearman correlation r.
Table 5.
Correlation analysis between IgG antibody responses against each peptide and blood parameters by gender. In dengue fever patients. Data is presented in Spearman correlation r.
Peptide |
Red blood cell count |
White blood cell count |
Platelet count |
Haemoglobin |
Haematocrit |
All |
|
|
|
|
|
AnDarApy-1 |
0.0955 (p=0.2934) |
0.0277 (p=0.7607) |
-0.0831 (p=0.3609) |
0.0999 (p=0.2715) |
0.1092 (p=0.2294) |
gSG6-P1 |
0.1807 (p=0.0455)
|
-0.0829 (p=0.3617) |
-0.0484 (p=0.5949) |
0.0921 (p=0.3112) |
-0.1097 (p=0.2269) |
Nterm-34kDa |
0.2107 (p=0.0193)
|
-0.0429 (p=0.6375) |
-0.0649 (p=0.4754) |
0.1543 (p=0.0885) |
0.1587 (p=0.0797) |
Females |
|
|
|
|
|
AnDarApy-1 |
-0.0571 (p=0.6390) |
-0.0641 (p=0.5980) |
0.0739 (p=0.5434) |
0.0847 (p=0.4858) |
0.0520 (p=0. 6689) |
gSG6-P1 |
0.0216 (p=0.8591) |
-0.1233 (p=0.3092) |
0.0338 (p=0.7815) |
0.0015 (p=0.9899) |
-0.0156 (p=0.8981) |
Nterm-34kDa |
0.0741 (p=0.5422) |
0.0480 (p=0.6932) |
-0.0153 (p=0.9001) |
0.1431 (p=0.2373) |
0.01078 (p=0.3474) |
Males |
|
|
|
|
|
AnDarApy-1 |
0.2797 (p=0.0425)
|
0.2063 (p=0.1382) |
-0.3052 (p=0.0263) |
0.1138 (p=0.4173) |
0.1681 (p=0.2290) |
gSG6-P1 |
0.4099 (p=0.0023)
|
0.0559 (p=0.6908) |
-0.2038 (p=0.1433) |
0.2074 (p=0.1362) |
0.2408 (p=0.0824) |
Nterm-34kDa |
0.3904 (p=0.0039)
|
0.0304 (p=0.8289) |
-0.1611 (p=0.2493) |
0.1922 (p=0.1680) |
0.2288 (p=0.0994) |
Table 6.
Correlation analysis between IgG antibody responses against each peptide and age. Data is presented in Spearman correlation r.
Table 6.
Correlation analysis between IgG antibody responses against each peptide and age. Data is presented in Spearman correlation r.
Peptide |
Age |
All |
gSG6-P1 |
-0.3533 (p=0.0003)
|
Nterm-34kDa |
-0.4182 (p=0.0000)
|
Correlations by gender |
Females |
|
gSG6-P1 |
-0.3575 (p=0.0035)
|
Nterm-34kDa |
-0.4087 (p=0.0007)
|
Males |
|
gSG6-P1 |
-0.3702 (p=0.0370)
|
Nterm-34kDa |
-0.4105 (p=0.0196)
|
Correlations by location |
Colombia |
gSG6-P1 |
0.2702 (p=0.0481)
|
Nterm-34kDa |
0.0796 (p=0.5675) |
US |
gSG6-P1 |
-0.4772 (p=0.0005)
|
Nterm-34kDa |
-0.4049 (p=0.0039)
|
Figure 3.
IgG antibody levels against the mosquito salivary peptides in healthy participants residents of Kansas (U.S.) in summer 2018 compared to levels in the IgG antibody levels found in volunteers followed during the Fall of 2018 (n=25). Significance was measured by the Mann-Whitney test p<0.05 (*=0.0, **=0.00, ***=0.000 and ns= not significant).
Figure 3.
IgG antibody levels against the mosquito salivary peptides in healthy participants residents of Kansas (U.S.) in summer 2018 compared to levels in the IgG antibody levels found in volunteers followed during the Fall of 2018 (n=25). Significance was measured by the Mann-Whitney test p<0.05 (*=0.0, **=0.00, ***=0.000 and ns= not significant).
3. Discussion
Mosquito saliva is composed of a plethora of molecules that are injected into the skin to counteract host responses and facilitate blood uptake. This saliva also induces production of host antibodies that correlate with the intensity of mosquito exposure. Notably, levels of IgG antibodies against whole or specific salivary proteins have been categorized as a reliable tool to measure exposure to mosquito bites and disease transmission intensity[
16,
38]. In malaria,
Plasmodium gametocytes in vertebrate blood are the infectious stages for the mosquito vector. Interestingly, previous studies showed that the carriage of specific parasite stages influenced host attractiveness to mosquitoes, with people that harbored gametocytes in their blood being more attractive to mosquitoes than those presenting parasitemia with only asexual stages[39-42]. Here, we evaluated gametocyte carriage by measuring the expression of
Plasmodium vivax Pvs25, while Pvs230 was used to determine gametocyte carriage[
43,
44]. Upon analysis of the level of antibodies against mosquito saliva, we observed that expression of Pvs230 was significantly negatively correlated with IgG antibody levels of all mosquito antigens tested in this study, including the non-malaria vector Aedes peptide Nterm34kDa. However, this association was not observed with Pvs25, suggesting that expression of Pvs230 may be associated with exposure to mosquito bites from both
Anopheles and
Aedes mosquitoes. Previous entomological data from this region indicate that Anopheles and Aedes mosquitoes are prevalent in the area[
45,
46]. However, a limitation of this study is that the gSG6-P1 and the Nterm34kDa peptides we evaluated as markers of mosquito exposure have not been extensively validated for their correlation with mosquito abundance[47-49], and we did not collect specific mosquito data from the sites in our study which could be used to confirm exposure to these mosquito species.
One of the objectives of this study was to determine if antibody levels against mosquito saliva were correlated with blood parameters. In the case of the malaria cohort, we observed that the non-malaria vector Aedes peptide Nterm34kDa was negatively associated with RBC count. This finding is notable as one of the signs of severe malaria is anemia, yet we did not observe correlations between any of the Anopheles salivary antigens and RBC count. The observed negative correlation between IgG anti-Nter34kDa and RBC count suggests that individuals with the lowest RBC counts may have the highest exposure to non-malaria vector Aedes mosquitoes or may produce more antibodies against their saliva. Interestingly, a previous study reported that anemia accelerates blood intake by
Ae. aegypti, although it may negatively impact egg production[
50]. Also, prior work has demonstrated that DENV acquisition by
Ae. aegypti was inversely correlated with iron concentration in human serum[
51]. Thus, it is possible that DENV benefits from the
Aedes preference to feed on anemic individuals. Indeed, among our DENV-infected participants, we observed a significant positive correlation between RBC and both Nterm34kDa and gSG6-P1 peptides, suggesting increased exposure to mosquito bites while infected. Moreover, a recent study suggested that infection with DENV increased attractiveness to Aedes mosquitoes[
52], similar to observations made with
Anopheles mosquitoes and
Plasmodium gametocyte carriers[
39,
42].
Mosquito abundance is associated with temperature, humidity, and other environmental factors[53-55]. In temperate and subtropical regions, mosquito abundance drastically changes in the summer vs. fall seasons, with higher mosquito populations in the former[
56,
57]. Previously, we observed a significant decrease in anti-tick antibodies during the fall months[
58]. Interestingly, we found a significant decrease in IgG antibody levels against the Nter34kDa peptide in the fall compared to summer in people living in Kansas. Our results are in agreement with previous studies, which suggest a significant decrease in antibodies against
Ae. aegypti mosquitoes after cessation of exposure[
13]. However, we did not observe a decrease in antibodies against gSG6P1, although prior reports have indicated that antibodies against the gSG6 protein are short-lived[
16], suggesting a potential continued exposure. In Kansas,
An. quadrimaculatus and
An. pseudopunctipennis are abundantly found in the summer, with
An. quadrimaculatus populations receding in late fall[
59]. Studies conducted in Northern Argentina indicate that
An. pseudopunctipennis presents two peaks of abundance, one in spring and another in fall/autumn[
60]. Together, these studies provide a potential explanation for the persistence of antibody levels in samples collected in early fall. Endemicity also appears to affect antibody levels against mosquito antigens in healthy individuals. Specifically, we observed that healthy volunteers living in a non-endemic area for either malaria or DENV (Kansas) presented negative correlations between age and anti-saliva antibodies, while healthy people living in Norte de Santander showed a positive correlation. We found a negative correlation between age and antibodies against Ae. aegypti whole salivary gland extract in healthy individuals living in Norte de Santander in 2015[
10] that could be associated with the development of tolerance as observed in other studies[
19,
38]. However, in our recent study in a population with DENV, we did not observe any associations with age[
61], as shown in other studies[
62], suggesting that several factors may contribute to these discrepancies. Further and larger studies are needed to establish the relationship between age and response to arthropod salivary antigens.
In addition to previous studies showing sex-associated differences in DENV incidence, geographical area also appears to strongly influence the rates of association and occurrence. For instance, a study including DENV cases from at least six Asian countries showed a higher incidence of dengue fever in males[
63], while other studies in Central and South America demonstrated a higher incidence in females[
64,
65]. Herein, we observed that males and females respond differently to mosquito salivary antigens. These results may have implications for pathogen transmission and clinical disease presentation. A major implication of our study is that characteristics, including age, seasonality, and vector control must be considered when assessing serum levels of IgG against
Ae. aegypti salivary proteins as a surrogate of risk of human exposure to mosquito bites and pathogen transmission. Such information is needed in epidemiological studies aimed at control and prevention of mosquito-borne diseases. Furthermore, data on biological sex at the population level could potentially be used to inform calculations of total disease burden in regions where vector-borne diseases have the highest impact.
In conclusion, our data indicate that individuals living in malaria- and dengue-endemic areas are exposed to several mosquito species in a given area. Moreover, exposure to malaria and dengue mosquito vectors may be associated with clinical and immune responses to infection with the reciprocal pathogens transmitted by these specific species, as well as other mosquito species that may not be implicated as vectors of human disease. Our findings suggest that concurrent exposure to multiple disease-carrying mosquito vectors and their salivary proteins with differing immunomodulatory properties could influence the transmission and pathogenesis of malaria, dengue fever, and other vector-borne illnesses. A better understanding of the molecular mechanisms underlying how exposure to multiple and sequential bites from numerous mosquito species influences immune responses and pathogen transmission will advance the development of immune-targeting interventions to reduce disease spread. Moreover, including data on antibody responses against the main mosquito species found in specific areas will provide new knowledge on the interplay of such species, pathogen transmission, and disease severity in humans, which will help to inform the development of more effective vector control and disease prevention efforts.
Author Contributions
conceptualization, BLR, JAM; methodology, BLR; sample collection, JCC, CAP, MUG, LYG-S and LAJ-V; sample testing, MMH and OO; formal analysis, MMH, OO and BLR; investigation, BLR, JAM; resources, BLR and JAM; writing—original draft preparation, BLR, JAM, JGS and BW; draft revisions BLR, JAM, JGS and BW, project administration, BLR and JAM; funding acquisition, BLR and JAM.