Anopheles salivary antigens as serological biomarkers of vector exposure and malaria transmission: A systematic review with multilevel modelling
Abstract
Background:
Entomological surveillance for malaria is inherently resource-intensive and produces crude population-level measures of vector exposure which are insensitive in low-transmission settings. Antibodies against Anopheles salivary proteins measured at the individual level may serve as proxy biomarkers for vector exposure and malaria transmission, but their relationship is yet to be quantified.
Methods:
A systematic review of studies measuring antibodies against Anopheles salivary antigens (PROSPERO: CRD42020185449). Multilevel modelling (to account for multiple study-specific observations [level 1], nested within study [level 2], and study nested within country [level 3]) estimated associations between seroprevalence with Anopheles human biting rate (HBR) and malaria transmission measures.
Results:
From 3981 studies identified in literature searches, 42 studies across 16 countries were included contributing 393 study-specific observations of anti-Anopheles salivary antibodies determined in 42,764 samples. A positive association between HBR (log transformed) and seroprevalence was found; overall a twofold (100% relative) increase in HBR was associated with a 23% increase in odds of seropositivity (OR: 1.23, 95% CI: 1.10–1.37; p<0.001). The association between HBR and Anopheles salivary antibodies was strongest with concordant, rather than discordant, Anopheles species. Seroprevalence was also significantly positively associated with established epidemiological measures of malaria transmission: entomological inoculation rate, Plasmodium spp. prevalence, and malarial endemicity class.
Conclusions:
Anopheles salivary antibody biomarkers can serve as a proxy measure for HBR and malaria transmission, and could monitor malaria receptivity of a population to sustain malaria transmission. Validation of Anopheles species-specific biomarkers is important given the global heterogeneity in the distribution of Anopheles species. Salivary biomarkers have the potential to transform surveillance by replacing impractical, inaccurate entomological investigations, especially in areas progressing towards malaria elimination.
Funding:
Australian National Health and Medical Research Council, Wellcome Trust.
Editor's evaluation
We believe this systematic review on the use of serological data to monitor anopheline mosquito exposure will add to the existing literature and help provide important insight into how these markers may be used to understand malaria transmission.
https://doi.org/10.7554/eLife.73080.sa0Introduction
Sensitive and accurate tools to measure and monitor changes in malaria transmission are essential to track progress towards malaria control and elimination goals. Currently, the gold standard measurement of malaria transmission intensity is the entomological inoculation rate (EIR), a population measure defined as the number of infective Anopheles mosquito bites a person receives per unit of time. EIR is calculated as the human biting rate (HBR; measured at the population level by entomological vector-sampling methodologies [gold standard: human landing catch]) multiplied by the sporozoite index (proportion of captured Anopheles with sporozoites present in their salivary glands). However, estimation of EIR and HBR via entomological investigations is inherently labour and resource intensive, requiring trained collectors, specialised laboratories, and skilled entomologists. Furthermore, these approaches provide a crude population-level estimate of total vector exposure at a particular time and location, precluding investigation of heterogeneity and natural transmission dynamics of individual-level vector–human interactions (Monroe et al., 2020). For example, indoor human landing catches provide poor estimates of outdoor biting and thus total vector exposure (Mathenge et al., 2005). The sensitivity of EIR is further compromised in low transmission settings where the number of Plasmodium-infected specimens detected is low and often zero.
Evaluation of the human antibody response to Anopheles spp. salivary proteins has the potential to be a logistically practical approach to estimate levels of exposure to vector bites at an individual level. Several Anopheles salivary proteins have been shown to be immunogenic in individuals naturally exposed to the bites of Anopheles vectors and have been investigated as serological biomarkers to measure Anopheles exposure (Badu et al., 2012b; Drame et al., 2013a; Drame et al., 2010a; Drame et al., 2010b; Drame et al., 2015; Rizzo et al., 2011a; Rizzo et al., 2011b; Stone et al., 2012; Drame et al., 2012), malaria transmission (Londono-Renteria et al., 2015a; Ya-Umphan et al., 2017; Noukpo et al., 2016), and as an outcome for vector control intervention studies (Drame et al., 2013a; Drame et al., 2010a; Drame et al., 2010b; Noukpo et al., 2016; Idris et al., 2017). However, a major shortcoming of the literature is that studies are largely descriptive and do not quantify the association between entomological and malariometric measures and anti-Anopheles salivary antibody responses. We undertook a systematic review with multilevel modelling to quantify the association between HBR, EIR, and other markers of malaria transmission, with anti-Anopheles salivary antibody responses, and to understand how these associations vary according to transmission setting and dominant Anopheles vectors which can exhibit different biting behaviours. In particular, we were interested in comparing the African context (where Anopheles gambiae and Plasmodium falciparum predominates) to non-African settings (where An. gambiae is absent and where both P. falciparum and Plasmodium vivax are prevalent). This knowledge is pertinent to advance the use of salivary antibody biomarkers as a vector and malaria transmission serosurveillance tool.
Methods
Search strategy and selection criteria
We performed a systematic review with multilevel modelling according to the Meta-analysis of Observational Studies in Epidemiology (MOOSE) and Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines (Moher et al., 2009; Stroup et al., 2000) (Reporting Standards Document). Five databases were searched for published studies investigating antibodies to Anopheles salivary antigens as a biomarker for mosquito exposure or malaria transmission published before 30 June 2020. The protocol (Appendix 1) was registered with PROSPERO (CRD42020185449).
The primary criterion for inclusion in this systematic review was the reporting of estimates of seroprevalence or total levels of immunoglobulin (Ig) in human sera against Anopheles salivary antigens. We considered for inclusion cross-sectional, cohort, intervention, and case–control studies of individuals or populations living in all geographies with natural exposure to Anopheles mosquitoes. Studies that were solely performed in participants not representative of the wider naturally exposed population (i.e. mosquito-allergic patients, soldiers, returned travellers) were excluded.
Measures
Outcomes
The primary outcome of our systematic review was antibodies (seroprevalence or levels, including all Ig isotypes and subclasses) against any Anopheles salivary antigens (full-length recombinant proteins, peptides, and crude salivary extract). Study-reported salivary antibody data was extracted at the most granular level (i.e. for each site; time point), with each observation of seroprevalence or levels included as a study-specific salivary antibody observation. As measurement of antibody levels does not produce a common metric between studies, only values of seroprevalence could be included in multilevel modelling analyses. Therefore, to maximise data, authors of studies that reported only antibody levels were contacted and asked to classify their participants as ‘responders’ or ‘non-responders’ according to seropositivity (antibody level relative to unexposed sera). Studies that provided antibody levels or categorised seropositivity based upon arbitrary cut-offs are included in narrative terms only.
Exposures
The primary exposures of interest were the entomological metrics HBR (average number of bites received per person per night) and EIR (infectious bites received per person per year). Secondary exposures included study-reported prevalence of Plasmodium spp. infection (confirmed by either microscopy, rapid diagnostic test (RDT), or polymerase chain reaction [PCR]) and seroprevalence of antimalarial antibodies against pre-erythrocytic and blood stage Plasmodium spp. antigens. Where exposure estimates were not provided, we attempted to source data from other publications by the authors or used the site geolocation (longitude and latitude) and year to obtain estimates of EIR from the Pangaea dataset (Yamba et al., 2018), P. falciparum rates in 2–10 year olds (PfPR2-10), and dominant vector species (DVS) from the Malaria Atlas Project (MAP; The Malaria Atlas, 2017). Malarial endemicity classes were derived by applying established endemicity cut-offs to MAP PfPR2-10 estimates (Bhatt et al., 2015). For the purposes of the modelling analyses, we defined DVS as where An. gambiae sensu lato (s.l.) was the only DVS, where An. gambiae s.l. was present with additional DVS, or where An. gambiae s.l. was absent. Studies of salivary antigens where exposure variables could not be sourced and data could not be extracted were excluded.
Statistical analysis
Where observations of the seroprevalence of antibodies against the same salivary antigen and exposure of interest were reported in more than one study, generalised linear multilevel modelling (mixed effects, logistic) was used to quantify associations between the exposures of interest and salivary antibody seroprevalence measurements (Song et al., 2019). Random intercepts for study and country were estimated to account for nested dependencies induced from multiple study-specific salivary antibody observations (level 1) from the same study (level 2) and studies from the same country (level 3). Additionally, study-level random slopes for the entomological and malariometric exposure parameters were estimated to model study-specific heterogeneity in the effect of the exposure of interest (HBR/EIR/malaria prevalence/antimalarial antibody seroprevalence). The associations between the various exposures and the different salivary antigens were analysed separately; however, observations of IgG seroprevalence against the recombinant full-length protein (gSG6) and synthetic peptide (gSG6-P1, the one peptide determined in all studies utilising peptides) form of the gSG6 antigen were analysed together.
Potential effect modification of the associations between exposures and anti-Anopheles salivary antibody responses was explored. In analyses quantifying the associations between HBR, as well as EIR, and seropositivity, we included an interaction term with DVS and for vector collection method (human landing catch or other indirect measures, e.g. light traps, spray catches, etc.). For the association between Plasmodium spp. prevalence and seropositivity, interaction terms with malaria detection methodology (light microscopy or PCR) and malarial species (P. falciparum only, or P. falciparum and P. vivax) were estimated.
For the exposure measures (HBR, EIR, malaria prevalence, and antimalarial antibody seroprevalence), the data were log transformed since there were non-linear associations between the exposure measures on the original scale and seroprevalence – supported empirically by superior model fit as indicated by Akaike’s information criterion (AIC) and Bayesian information criterion (BIC) fit indices (Appendix 1—table 1). To aid interpretation, we present our results as a relative increase in the odds of the gSG6 IgG seropositivity for a twofold or, in other words, a 100% relative increase in the exposures. Intraclass correlation coefficients (ICCs) were estimated for country- and study-specific heterogeneity using estimated model variance components. In order to explore the presence of study-level influence in (HBR and EIR) effect estimate modelling, the Generalised Linear Latent and Mixed Models (gllamm) package (Rabe-Hesketh et al., 2000) was used to produce Cooks distance statistics (Cook, 1977) at the study level from the generalised linear multilevel models. A conservative cut-off threshold for Cooks distance (4 /n) was used to guide sensitivity analyses, where studies were excluded, in turn, to assess outlier influence. All statistical analyses were performed using STATA v15.1.
Risk of bias in individual studies
Risk of bias was assessed by one reviewer using the Risk of Bias in Prevalence Studies tool (Hoy et al., 2012). The risk of bias pertains to the reported observations of anti-Anopheles salivary antibody seroprevalence included in the multilevel modelling.
Results
Literature searches identified 158 potentially relevant studies, of which 42 studies were included in the systematic review (Figure 1) and are described in Table 1. From these studies, we extracted n = 393 study-specific observations of anti-Anopheles salivary antibodies determined from antibody measurements in a total of 42,764 sera samples. These studies were performed in 16 countries mostly in hypo- or mesoendemic areas of Africa (32 studies), with a minority performed in South America (four studies), Asia (four studies), and the Pacific (two studies). Studies were classified according to their DVS which reflected the region where the study was conducted. An. gambiae s.l. was a DVS in all African study sites (n = 151 study-specific observations from 23 studies where An. gambiae s.l. was the only DVS and n = 68 from 16 studies where An. gambiae s.l. was present with additional DVS [i.e. Anopheles funestus, Anopheles pharoensis]), with the exception of one study, which together with the 10 non-African studies contributed n = 174 study-specific estimates where An. gambiae s.l. was absent. Most observations came from cross-sectional (n = 191 from 16 studies) or repeated cross-sectional studies (n = 137 from 18 studies), with n = 60 from cohort studies (six studies) and n = 5 from case–control studies (two studies).
The salivary antigen most commonly assessed was An. gambiae salivary gland 6 (gSG6), as a full-length protein (n = 67 from 8 studies) and synthetic peptide (An. gambiae salivary gland 6 peptide 1 [gSG6-P1]; n = 270 from 24 studies). Additional salivary antigens assessed included An. gambiae gSG6-P2 (n = 119 from three studies), recombinant cE5 (n = 15 from two studies), g-5’nuc (n = 3 from one study), and recombinant An. funestus fSG6 (n = 6 from two studies) and f-5’nuc (n = 3 from one study). Seven studies measured antibodies to whole salivary gland extracts (SGE) from An. gambiae (n = 24 from four studies), Anopheles darlingi (n = 5 from two studies), Anopheles albimanus (n = 2 from one study), and Anopheles dirus (n = 3 from one study), while one study assessed antibodies against synthetic peptides of An. albimanus (n = 2) (Table 1). All studies investigated total IgG and only five determined an additional isotype or subclass (Drame et al., 2015; Lawaly et al., 2012; Rizzo et al., 2014a; Rizzo et al., 2014b; Waitayakul et al., 2006). The paucity of studies investigating these latter-mentioned antibody types and Anopheles salivary biomarkers precluded extensive multilevel analyses; instead, we present their associations in narrative terms in Appendix 10. Analyses reported below focus on quantifying the relationships between HBR, EIR, and markers of malaria transmission with total IgG to An. gambiae gSG6. The distributions of exposure observations were: HBR (n = 197 from 24 studies, median: 3.0 bites per person per night, IQR: 0.9–12.1; range: 0–121.4), EIR (n = 60 from 8 studies, median: 7.3 infectious bites received per person per year, IQR: 0–36.4; range: 0–585.6), and Plasmodium spp. prevalence (n = 266 from 22 studies, median: 9.1%; IQR: 4–22%; range: 0–94.6%).
Generalised linear multilevel modelling (mixed effects, logistic) of n = 132 study-specific observations from 12 studies estimated a positive association between Anopheles spp.-HBR (log transformed) and seroprevalence of IgG to An. gambiae gSG6 salivary antigen (Drame et al., 2010a; Drame et al., 2015; Rizzo et al., 2011a; Stone et al., 2012; Drame et al., 2012; Ya-Umphan et al., 2017; Soma et al., 2018; Traoré et al., 2019; Sagna et al., 2013b; Sarr et al., 2012; Ali et al., 2012; Pollard et al., 2019; Figure 2, Appendix 4—table 1). As we have log transformed HBR to account for the non-linear relationship between HBR and log odds of gSG6 IgG seropositivity, we have presented estimated odds ratios for different incremental percent increases in HBR (Figure 2—figure supplement 1). For example, the magnitude of the association was such that a twofold (100% relative) increase in HBR was associated with a 23% increase (OR: 1.23; 95% CI: 1.10–1.37; p<0.001) in the odds of anti-gSG6 IgG seropositivity (Figure 2). Heterogeneity in the effect of HBR on gSG6 across studies was observed (likelihood ratio χ2(1) = 109.25, p<0.001); the 95% reference range of study-specific effects for a twofold increase in HBR ranged from a 12% reduction to a 70% increase in odds (OR: 0.88–1.70). There was no evidence that the association between HBR and gSG6 IgG varied according to vector collection method (human landing catch or other indirect methods; p=0.443) or study design (longitudinal cohort or cross-sectional/repeated cross-sectional; p=0.138). Given the global heterogeneity in the distribution of Anopheles species, we sought to quantify the extent to which the association between An. gambiae gSG6 IgG seropositivity and HBR is moderated by DVS. We observed that the magnitude of the association between An. gambiae gSG6 IgG seropositivity and HBR was greatest in African studies where An. gambiae s.l. was the only dominant vector (p<0.001, Appendix 5); a twofold increase in HBR was associated with a 37% increase (OR: 1.37; 95% CI: 1.19–1.58; p<0.001) in the odds of gSG6 IgG seropositivity compared to an attenuated association for African studies where An. gambiae s.l. was not the only DVS (OR: 1.14 per twofold increase in HBR; 95% CI: 0.98–1.33; p=0.079) and non-African studies where An. gambiae s.l. was absent (OR: 1.05 per twofold increase in HBR; 95% CI: 1.03–1.08; p<0.001). In order to quantify the relationship between gSG6 IgG seroprevalence and HBR, for given HBR values we estimated gSG6 IgG seroprevalence by producing model-based predicted probabilities overall and by DVS (Figure 3). In African studies where An. gambiae s.l. is the only DVS, predicted seroprevalence of An. gambiae gSG6 ranged from 21% (95% CI: 0–45%) to 86% (95% CI: 67–100%) for an HBR of 0.1–100 bites per person per night, respectively (Figure 3, Figure 3—figure supplement 1).
A positive association was also found between seroprevalence of anti-gSG6 IgG antibodies and EIR in analysis of n = 38 study-specific observations from eight studies (Figure 4, Appendix 6) [Rizzo et al., 2011b; Ya-Umphan et al., 2017; Soma et al., 2018; Ali et al., 2012; Ambrosino et al., 2010; Perraut et al., 2017; Pollard et al., 2019; Badu et al., 2012b]. For a twofold increase in EIR, the odds of anti-gSG6 IgG seropositivity increased by 11% (OR: 1.11; 95% CI: 1.05–1.17; p<0.001), with heterogeneity in the study-specific effects (95% reference range: 1.00–1.24; likelihood ratio χ2(1) = 15.02, p<0.001). There was no evidence of effect modification by either vector collection method (p=0.095) or DVS (p=0.080) on the association between seroprevalence of anti-gSG6 IgG and EIR.
Similar positive associations were also found between anti-gSG6 IgG levels, HBR, and EIR in 11 studies [Drame et al., 2015; Drame et al., 2012; Stone et al., 2012; Soma et al., 2018; Ali et al., 2012; Rizzo et al., 2011a; Rizzo et al., 2011b; Poinsignon et al., 2010b; Poinsignon et al., 2008a; Charlwood et al., 2017; Sagna et al., 2013b] and 3 studies [Rizzo et al., 2011b; Ya-Umphan et al., 2017; Ali et al., 2012], respectively, but 7 studies showed no association between HBR and levels of IgG to gSG6 [Drame et al., 2010a; Ya-Umphan et al., 2017; Traoré et al., 2018; Traoré et al., 2019; Sarr et al., 2012; Poinsignon et al., 2009; Pollard et al., 2019].
The association between anti-gSG6 IgG seroprevalence and population-level prevalence of Plasmodium spp. infection was investigated. Generalised linear multilevel modelling (mixed effects, logistic) of n = 212 from 14 studies that measured Plasmodium spp. prevalence contemporaneously in their study [Badu et al., 2012b; Rizzo et al., 2011b; Ya-Umphan et al., 2017; Soma et al., 2018; Koffi et al., 2015; Traoré et al., 2019; Badu et al., 2015; Sagna et al., 2013b; Sarr et al., 2012; Perraut et al., 2017; Drame et al., 2010a; Idris et al., 2017; Proietti et al., 2013; Kerkhof et al., 2016] showed that for a twofold increase in the prevalence of Plasmodium spp. infection the odds of gSG6 IgG seropositivity increased by 38%, although the confidence intervals were wide (OR: 1.38; 95% CI: 0.89–2.12; p=0.148) and heterogeneity in the study-specific effects was observed (95% reference range: 0.30–6.37; likelihood ratio χ2(1) = 235.5, p<0.001) (Figure 5 and Appendix 7). In the association between gSG6 IgG seropositivity and Plasmodium spp. infection, there was no evidence for a moderating effect of Plasmodium spp. detection method (light microscopy or PCR, p=0.968), or species (African studies with P. falciparum versus non-African studies where P. falciparum and P. vivax are co-prevalent, p=0.538).
Additionally, 14 studies reported observations of anti-gSG6 IgG levels and the prevalence of Plasmodium spp. infections measured contemporaneously in their study. The median anti-gSG6 IgG antibody levels increased with increasing Plasmodium spp. prevalence in six of these studies (Drame et al., 2010a; Ya-Umphan et al., 2017; Idris et al., 2017; Poinsignon et al., 2010b; Sarr et al., 2012; Kerkhof et al., 2016), or in Plasmodium spp.-infected compared to non-infected individuals (Londono-Renteria et al., 2015a; Montiel et al., 2020), but showed no association in eight studies (Rizzo et al., 2011b; Soma et al., 2018; Koffi et al., 2015; Traoré et al., 2018; Traoré et al., 2019; Badu et al., 2015; Sagna et al., 2013b; Poinsignon et al., 2009). Furthermore, we also investigated associations with serological measures of malaria exposure and found that for a twofold increase in pre-erythrocytic and blood stage antigen seroprevalence there was a 2.19-fold (OR: 2.19; 95% CI: 1.18–4.04; p=0.013) and 41% to 5.69-fold (OR range: 1.41–5.69; p range: <0.001 to 0.523) increase in the odds of anti-gSG6 IgG seropositivity, respectively (Appendix 8).
To give epidemiological context, we estimated anti-gSG6 seroprevalence by producing model-based predicted probabilities by malarial endemicity class (a categorical variable derived by applying established cut-off values for the PfPR2-10 extracted from MAP). Generalised linear multilevel modelling (mixed effects, logistic) on 297 study-specific salivary antibody observations from 22 studies shows that the estimated anti-gSG6 IgG seroprevalence is higher for the higher endemicity classes (eliminating malaria: 20% [95% CI: 8–31%]; hypoendemic: 34% [95% CI: 19–49%]; mesoendemic: 52% [95 CI: 35–68%]; hyperendemic settings: 47% [95% CI: 27–64%]; holoendemic: 78% [95% CI: 67–90%]; p<0.001; Table 2). Interactions with DVS or region (Africa/non-Africa) could not be explored due to collinearity with malaria endemicity class. Therefore, in addition using Bayes best linear unbiased predictions (BLUPs) we estimated country-specific gSG6 IgG seroprevalence from an intercept-only multilevel model fitted to 301 study-specific salivary antibody observations from 22 studies. It showed that IgG seroprevalence to An. gambiae gSG6 was lowest in countries in the Pacific region where An. gambiae is absent (Vanuatu [31%] and Solomon Islands [32%]) and highest in countries where An. gambiae is a DVS (Benin [72%] and Burkina Faso [65%]; Appendix 9).
Assessments of internal and external study validity revealed there was a moderate risk of selection bias (Appendix 2) due to the study-specific inclusion criteria of populations at higher risk of malaria which contributed gSG6 seroprevalence observations. Sensitivity analyses exploring potential study-level outlier influence on the estimated associations between anti-gSG6 IgG seroprevalence, HBR, and EIR showed no evidence of bias (effect estimates for each sensitivity analysis were consistent with model estimates overall) for studies identified as exhibiting potential influence (HBR: n = 6; EIR: n = 6).
Discussion
This systematic review and multilevel modelling analysis provides the first quantification of a positive non-linear association between seroprevalence of An. gambiae gSG6 IgG antibodies and HBR and demonstrated that its magnitude varied with respect to the DVS present in the area. Importantly, this review identified a paucity of studies conducted outside of Africa, as well as investigating salivary antigens representing different Anopheles spp. and antigenic targets. gSG6 antibodies were positively associated with the prevalence of Plasmodium spp. infection as well as established epidemiological measures of malaria transmission: malaria endemicity class and EIR. Overall, our results demonstrate that antibody seroprevalence specific for Anopheles spp. salivary antigens has the potential to be an effective measure of vector exposure and malaria transmission at the population and, potentially, individual level.
An. gambiae gSG6 IgG seropositivity increased with increasing HBR, although these increases had diminishing impact on An. gambiae gSG6 IgG seropositivity at higher levels of HBR (approximately greater than two bites per person per night). In our study, 17 studies performed across Africa (Angola, Benin, Burkina Faso, Cote d’Ivoire, and Senegal) and the Asia Pacific (Cambodia, Myanmar, and the Solomon Islands) reported an HBR < 2, demonstrating the applicability of gSG6 as a biomarker of HBR across a broad range of malaria-endemic regions. We also observed that the association was strongest in areas where An. gambiae s.l. was the only DVS (i.e. concordant An. gambiae species-specific HBR with An. gambiae gSG6 antibodies). Associations, albeit weaker, were also observed between discordant species-specific HBR and gSG6, most likely because the An. gambiae SG6 gene shares moderate sequence identity with vector species that are dominant in other regions (Africa: 80% An. funestus; Asia: 79% Anopheles stephensi and Anopheles maculatus; 54% An. dirus; Pacific: 52.5% Anopheles farauti), and is absent from the DVS of the Americas (An. albimanus and An. darlingi) (Arcà et al., 2017). The generalisability of An. gambiae gSG6 IgG as a biomarker of exposure to other Anopheles spp. may therefore be limited. However, our review also identified a paucity of studies investigating additional salivary antigenic targets and Anopheles species not present in Africa. The identification of novel salivary antigens that are species-specific will be valuable in quantifying exposure to the other Anopheles vectors that share limited identity with An. gambiae SG6 (such as An. farauti and An. dirus), as well as Anopheles spp. which lack SG6 (as done for An. albimanus and An. darlingi; Londono-Renteria et al., 2020a; Londono-Renteria et al., 2020b). An Anopheles species-specific serological platform could advance vector surveillance by more accurately capturing exposure to DVS in the South American and Asia Pacific regions which exhibit diverse biting behaviours and vector competence (DVS typically bite outdoors during the night and day, respectively; The Malaria Atlas, 2017; Sinka et al., 2012; Sinka et al., 2010; Trung et al., 2005; Herrera et al., 2015; Chaumeau et al., 2018), as well as the increasing threat of urban malaria from An. stephensi in Africa (Takken and Lindsay, 2019; Sinka et al., 2020).
This review demonstrated that the prevalence of Anopheles salivary antibodies increased with increasing prevalence of Plasmodium spp. infection (although confidence intervals were wide and we observed heterogeneity in the effect between studies) as well as established epidemiological measures of malaria transmission: malaria endemicity class and EIR. Anti-salivary antibodies, such as SG6 IgG, may therefore have the potential to serve as a proxy measure for receptivity of a population to sustain malaria transmission. Their application could be particularly relevant in pre-elimination areas, or non-endemic areas under threat of imported malaria, where Anopheles salivary antibodies are more readily detectable than parasites; salivary antibodies were predicted to be prevalent (20%) in areas defined as eliminating malaria (<1% PfPR2-10). Furthermore, if SG6 IgG seroprevalence can be effectively combined with a measurement of the sporozoite index, salivary antibodies as a marker of HBR could help overcome sensitivity limitations of EIR in low transmission areas. Additional measures could include estimates of malaria prevalence or serological biomarkers that are species- or life stage-specific (e.g. Plasmodium spp. pre-erythrocytic antigens as biomarkers for recent parasite inoculation). Indeed, positive associations between antibodies specific for Plasmodium spp. pre-erythrocytic and blood stage antigens with gSG6 were demonstrated in analyses of data from diverse malaria-endemic areas. Serological tools combining salivary antigens with antigens specific for the different Plasmodium spp. could be easy to employ and complement malaria surveillance programmes. These tools may be particularly useful in the Asia Pacific, a region of relatively low malaria transmission with goals of elimination, but the highest burden of P. vivax malaria where blood stage infection can be caused by relapses from dormant liver stages. In these areas, parasite prevalence may therefore overestimate ongoing malaria transmission, making vector surveillance tools essential to informing elimination strategies in the Asia Pacific and other regions where P. vivax is endemic.
The gold standard entomological measures HBR and EIR provide crude population-level estimates of vector and malaria exposure that are specific in space and time and preclude investigation of individual-level heterogeneity and natural transmission dynamics. Our study demonstrated that salivary biomarkers measured at the individual level, such as gSG6 IgG, can be used to quantify total vector exposure at the population level, without requiring laborious entomological experiments. However, validating an individual-level serological measure, which demonstrates considerable individual-level variation, against the imperfect population-level gold standards of HBR and EIR is challenging and reflected in the variation in study-specific estimates in the association between gSG6 IgG and HBR in modelling analyses. However, the accuracy of salivary antibodies to measure individual-level exposure to Anopheles bites is yet to be validated; literature searches identified no studies investigating this association at the individual level. Without detailed measurements of individual-level vector exposure, or a detailed knowledge of the half-life of Anopheles salivary antibodies post biting event, the true accuracy of salivary antibodies, such as SG6 IgG, to measure individual-level HBR remains unknown. This knowledge is particularly pertinent where Anopheles salivary biomarkers might be applied to assess the effectiveness of a vector control intervention or used to measure temporal changes in malaria transmission; particularly in areas or populations where there is considerable heterogeneity in individual-level risk of Anopheles exposure (e.g. unmeasured outdoor biting due to occupational exposure for forest workers; Sandfort et al., 2020).
The broad nature of our inclusion and quality criteria was a key strength of our systematic review, which aimed to provide a comprehensive analysis of all Anopheles salivary biomarkers and determine their associations with entomological and malariometric measures of transmission. However, this review has two main limitations. First, despite the inclusive nature, assessment of the external validity of the review revealed a moderate risk of bias; some studies exhibited a high risk of selection bias as they were performed in specific high-risk populations not representative of the overall population (i.e. children only). This is accounted for to some degree by specification of a random effect (i.e. intercept) for study, which accounts for unmeasured study-specific factors that may introduce study-specific measurement error to measurement of the outcome. Second, with respect to internal validity, there may be potential selection bias introduced by the exclusion of studies reporting zero HBR (7 observations from three studies; Rizzo et al., 2011b; Pollard et al., 2019; Sagna et al., 2013b), EIR (22 observations from three studies; Ya-Umphan et al., 2017; Rizzo et al., 2011b; Soma et al., 2018), and malaria prevalence (15 observations from three studies; Idris et al., 2017; Sagna et al., 2013b; Kerkhof et al., 2016) estimates, given we modelled the log of these factors. However, adding a small constant (e.g. 0.001) to a zero value to permit modelling of a log estimate can also introduce considerable bias (i.e. seemingly small differences between values become very large on the log scale). In light of this, we also chose to provide estimates of association and gSG6 IgG seroprevalence according to a selected range of epidemiologically relevant hypothetical HBRs (no widely accepted HBR classification exists in the literature) and according to widely accepted, discrete, endemicity classes according to MAP estimates (which permitted inclusion of all studies) to provide epidemiological context. However, there is the potential for misclassification of malarial endemicity class derived from geospatially extracted MAP predictions of PfPR2-10 which increase in uncertainty in areas with scarce data. Similarly, we used MAP vector occurrence data to inform DVS categories for 7 (out of 42) studies. Cross-referencing these 7 studies with a 2017 updated database for African vectors (using data for the nearest neighbouring village) identified 10 discrepant datapoints from 3 studies (from a total of 28 datapoints from 7 studies) (Snow, 2017). Any misclassification events may cause us to underestimate the standard error in the effect of malaria endemicity class and DVS on gSG6 IgG.
Conclusions
In order to advance progress towards malaria elimination, the World Health Organization has called for innovative tools and improved approaches to enhance vector surveillance and monitoring and evaluation of interventions (World Health Organization, 2017). Our systematic review has provided evidence that Anopheles salivary antibodies are serological biomarkers of vector and malaria exposure, by quantifying their positive association with Anopheles-HBR and established epidemiological measures of malaria transmission. These salivary biomarkers have the potential to replace crude population-level estimates of entomological indices with a precise and scalable tool that measures Anopheles vector exposure at the individual level. This approach could be expanded into a serosurveillance tool to assess the effectiveness of vector control interventions, define heterogeneity in malaria transmission, and inform efficient resource allocation that would ultimately accelerate progress towards elimination.
Appendix 1
Supplementary methodology
Search strategy
We performed a systematic review with multilevel modelling of the published literature according to the MOOSE guidelines (Stroup et al., 2000) and the PRISMA specifications (Moher et al., 2009). The protocol was registered with PROSPERO (CRD42020185449).
The electronic databases PubMed, Scopus, Web of Science, African Index Medicus, and the Latin American and Caribbean Health Sciences Literature (LILACS) were searched for studies published before 30 June 2020 investigating Anopheles salivary antigens as a biomarker for mosquito exposure or malaria transmission. Search terms were as follows: Anophel* AND saliva* AND (antibod* OR sero* OR antigen OR marker* OR biomarker* OR gSG6* OR gSG* OR SG* OR cE5). The reference lists of included studies were screened for additional studies, and Google Scholar was used to identify additional works by key authors. No formal attempt was made to identify unpublished population studies as it would have required significant description of the design, methods, and analysis used in these studies, and a review of ethical issues.
Selection criteria
The primary criterion for inclusion in this systematic review was the reporting of observations of seroprevalence or total levels of Ig antibodies (including all isotypes and subclasses) in human sera against recombinant or synthetic peptide Anopheles salivary antigens. We considered for inclusion cross-sectional studies, cohort studies, intervention studies, and case–control studies of individuals or populations (including sub-populations) living in all geographies with natural exposure to Anopheles mosquitoes. Studies that were solely performed in participants not representative of the wider population (i.e. mosquito-allergic patients, soldiers, returned travellers) were excluded. The minimum quality criteria for inclusion in this review were antibody detection performed using enzyme-linked immunosorbent assay (ELISA), multiplex or Luminex assays.
The exposure variables of interest included entomological and malariometric parameters, including (i) HBR, defined as the number of bites received per person per unit of time; (ii) EIR, defined as the number of infectious bites per person per unit of time, calculated as the HBR multiplied by the sporozoite index; (iii) estimates of malaria prevalence; and (iv) population-level seroprevalence estimates against Plasmodium spp. malarial antigens. To ensure HBR estimates were given for the same unit of time (bites per person per night), biting rates given per week were divided by 7, and biting rates given per month we multiplied by 12 and divided by 365. Similar approaches were employed to ensure consistent units for EIR (infectious bites per person per year). Plasmodium spp. infections had to be confirmed by either microscopy, RDT, or molecular methods (PCR). Plasmodium spp. diagnosis was included for all Plasmodium spp. combined and the species level if provided. Where exposure estimates were not provided, we attempted to source data from other publications by the authors or used the site geolocation and year to obtain estimates of EIR from the Pangaea dataset (Yamba et al., 2018). P. falciparum rates in 2–10 year olds (globally, 2000–2017) and DVS from the MAP (The Malaria Atlas, 2017). Studies of salivary antigens where exposure variables could not be sourced and data that could not be extracted were excluded.
Selection of studies
One author performed database searches and screened reference lists to identify possible studies. One author screened studies against inclusion criteria, with discussion and input from a second reviewer.
Approaches to include all available studies
The authors of any studies that did not contain relevant information on the study design, populations, eligibility criteria, or key study data were contacted and relevant data requested. Authors were contacted via an initial email detailing the precise nature of the systematic review and the data required. If the authors did not reply to three email requests or were unable to provide relevant data, the studies were deemed to insufficiently meet inclusion/quality criteria and were excluded. As measurement of antibody levels does not produce a common metric between studies, authors were asked to classify their participants as ‘responders’ or ‘no-responders’ according to seropositivity (antibody level relative to unexposed sera) within each study to allow comparisons of seroprevalence between studies (Cutts et al., 2020; Cutts et al., 2014; Fowkes et al., 2010). Studies that were only able to provide antibody levels or categorised seropositivity based upon arbitrary cut-offs were excluded from multilevel modelling analyses and included in narrative terms. Where the salivary antibody response and exposure variable were measured in the same population and reported in multiple publications, the study with the largest sample size was included, otherwise the earliest study was included.
Data extraction
Data were extracted using a data collection form by one reviewer. Any data that was provided at the sub-population level was extracted at the lowest level, that is, if a study was performed across multiple sites, and an estimate for both salivary antibody seroprevalence/levels and the exposure of interest is given for each site, it was included the site level, rather than an aggregated level.
Measures
Outcomes
The primary outcome of interest of our systematic review was the reported antibody response (both seroprevalence and levels of all Ig subclasses and isotypes) to Anopheles salivary antigens. Multilevel modelling analyses were performed where the seroprevalence of antibodies against the same antigen and the exposure of interest were reported in more than one study.
Exposures
The primary exposures of interest included in the multilevel modelling analyses were the HBR and EIR, a measure of the average number of bites received per person per night and infectious bites received per person per year, respectively. Secondary exposures assessed include the prevalence of any Plasmodium spp. infection (including P. falciparum only, P. vivax only, or untyped infections). Additional secondary exposures include the P. falciparum infection rate in 2–10 year olds extracted from MAP, as well as the seroprevalence of antimalarial antibodies against pre-erythrocytic and blood stage antigens.
Clinical and methodological heterogeneity were explored using prespecified variables to minimise spurious findings. Variables considered for inclusion were study design (cohort, cross-sectional, repeated cross-sectional), DVS, study participants (adults only, children only, adults and children), preparation of salivary antigen (recombinant full-length protein, synthetic peptide), malaria detection methodology (light microscopy, RDT, PCR), and entomological vector collection methodology (human landing catch, light traps, and spray catches).
Statistical analysis
Where there were sufficient data to pool observations of the same exposure and outcome measures, generalised linear multilevel modelling was used to undertake analyses quantifying associations between the exposures of interest and salivary antibody seroprevalence measurements. Models were generalised through use of the logit link function and binomial distribution (statistical notation for HBR model shown below as Equation 1). Seroprevalence was modelled in binomial form as the number of individuals seropositive to the total sample size. A three-level random effects model with a nested framework was used to account for dependency in the data, with random intercepts for country (level 3) and study (level 2) estimated. Hence, level 1 units represented multiple salivary antibody observations within a study induced by the study design (i.e. multiple time points, sites, age categories). Additionally, study-level random slopes for entomological and malariometric exposures were estimated to permit the effects to vary across studies. Model structure was determined empirically through likelihood ratio tests (p<0.05), with the exception of country at the third,level which was included in HBR and EIR analyses to estimate country-specific seroprevalence estimates of anti-salivary antibodies. The associations between the various exposures and the different salivary antigens were analysed separately; however, observations of IgG seroprevalence against the recombinant full-length protein (gSG6) and synthetic peptide (gSG6-P1, the one peptide determined in all studies utilising peptides) form of the gSG6 antigen were analysed together, with a fixed term for antigen construct considered for inclusion in the model. Of note, gSG6 peptide 2 (gSG6-P2) was excluded from being analysed with gSG6 and gSG6-P1 as the two studies that reported anti-gSG6-P2 IgG seroprevalence also reported the seroprevalence of anti-gSG6-P1 IgG and only one could be included. Potential effect modification of the associations between the exposures of interest and the anti-Anopheles salivary antibody responses was explored and undertaken by estimating interaction terms for DVS (An. gambiae s.l. only, An. gambiae s.l. and other DVS, or An. gambiae s.l. absent) and for vector collection method (human landing catch or other indirect measures, e.g. light traps, spray catches, etc.). For the association between Plasmodium spp. prevalence and gSG6 IgG seropositivity, interaction terms for malaria detection methodology (light microscopy or PCR), and malarial species type (P. falciparum only, or P. falciparum and P. vivax) were estimated. Other variables considered for inclusion in adjusted models were study design, participant, and salivary antigen construct; however, these variables showed no association with anti-gSG6 IgG and were thus excluded.
AIC and BIC fit indices were used to determine the best-fitting functional forms for the association between log odds of gSG6 IgG seropositivity and HBR, EIR, and Plasmodium spp. prevalence – linear, log, quadratic, and cubic functions were fitted, with a log transformation exhibiting superior model fit (Appendix 1—table 1). To aid interpretation, we present our results as a relative increase in the odds of the gSG6 IgG seropositivity for a twofold (100% relative) increase in the exposures. Additional relative percent changes in HBR and EIR are also presented.
Empirical Bayes BLUPs were used to estimate the probability of gSG6 IgG seropositivity in the average study and country, which is equivalent to an estimated gSG6 IgG seroprevalence. In order to maximise the number of included studies in our modelling, we predicted anti-gSG6 seroprevalence according to endemicity class, derived by applying established endemicity cut-offs to PfPR2-10 estimates (Bhatt et al., 2015) extracted from MAP using site year and geolocation (if MAP data unavailable endemicity as stated in study). ICCs and 95% reference ranges were estimated for country-, study-, and slope-specific heterogeneity (where appropriate) using estimated model variance components.
Statistical notation for the generalised linear multilevel model (mixed effects, logistic) used to estimate the association between An. gambiae gSG6 IgG seropositivity and HBR
The model can be formally written as
where
~ N(0, ), ~ N(0, ) and ~ N(0, ), (2)
where is the vector of model covariates, is the model constant and represents the log odds (probability) of gSG6 IgG seropositivity for a log HBR of zero, is the fixed effect for log HBR for country j and study i, is the random effect (i.e. intercept) for between-country heterogeneity in probability of gSG6 IgG seropositivity, is the random effect (i.e. intercept) for between-study heterogeneity in probability of gSG6 IgG seropositivity, and is the random effect (i.e. coefficient) for between-study heterogeneity in the effect of log HBR.
Risk of bias in individual studies
For cross-sectional, cohort or intervention studies, selection bias was assessed by reviewing the studies’ inclusion and exclusion criteria. Any case–control studies or studies that presented salivary antibody data stratified by malaria infection status were included in narrative terms only. Risk of bias was assessed by one reviewer using the Risk of Bias in Prevalence Studies tool (Hoy et al., 2012). The risk of bias pertains to the reported observations of anti-Anopheles salivary antibody seroprevalence included in the multilevel modelling.
Appendix 2
Risk of bias assessment
Risk of bias was assessed for each study by one independent reviewer using the Risk of Bias in Prevalence Studies tool (Hoy et al., 2012). This tool comprises 10 items and a summary assessment to assess the external validity (selection and non-response bias) and internal validity (measurement bias) of the study’s seroprevalence observations. The risk of bias pertains to the reported observations of anti-Anopheles salivary antibody seroprevalence included in the multilevel modelling.
With regard to external validity, seven of the studies included in the review were performed in specific populations (i.e. children only) that were not representative of the national population and were deemed to be at high risk of selection bias. Only seven studies included some form of random sampling, and frequently insufficient detail was provided on the sampling frame; as such most studies were included as high risk of selection bias. Furthermore, no studies reported participant response rate, and as such were indicated as high risk of non-response bias.
In terms of internal validity, all studies had an acceptable case definition, with the same mode of data collection, a valid instrument, and an acceptable prevalence period, so were all deemed to be of low risk. However, 12 studies did not include a denominator, instead only reporting the study sample size and prevalence estimate, and were included as high risk.
Overall, due to the specific nature of some of the sample populations for which these prevalence observations are given (i.e. children only) and as participant non-response rate is not given, we conclude that there is a moderate risk of study bias. According to the Risk of Bias in Prevalence Studies tool (Hoy et al., 2012), this implies that future research is likely to have an impact on our confidence in the prevalence observations.
Appendix 3
Reasons for study exclusion
Appendix 4
Association between gSG6 IgG seropositivity and HBR
Appendix 5
Association between gSG6 IgG seropositivity and HBR, moderated by dominant vector species
Appendix 6
Association between gSG6 IgG seropositivity and EIR
Appendix 7
Association between gSG6 IgG seropositivity and malaria prevalence
Appendix 8
Association between gSG6 IgG seropositivity and antimalarial antibody seroprevalence
Antibodies against P. falciparum pre-erythrocytic stage antigens
The pooled analysis of 159 study-specific observations from eight studies showed that a twofold increase in PfCSP IgG seropositivity was associated with a 2.19-fold (OR: 2.19; 95% CI: 1.18–4.04; p=0.013) increase in odds of anti-gSG6 IgG seropositivity (Stone et al., 2012; Ya-Umphan et al., 2017; Koffi et al., 2015; Koffi et al., 2017; Ambrosino et al., 2010; Perraut et al., 2017; Proietti et al., 2013; Kerkhof et al., 2016). Furthermore, we observed that gSG6 IgG levels increased with increasing PfCSP IgG seroprevalence in four studies (Ya-Umphan et al., 2017; Koffi et al., 2015; Koffi et al., 2017; Kerkhof et al., 2016), with another study contributing only a single estimate (Stone et al., 2012).
Antibodies against P. falciparum blood stage antigens
Furthermore, we observed a twofold increase PfAMA1 IgG seroprevalence was associated with a 2.47-fold (OR: 2.47; 95% CI: 2.25–2.71; p<0.001) increase in odds of gSG6 IgG seropositivity based upon 62 study-specific observations from eight studies (Stone et al., 2012; Ya-Umphan et al., 2017; Koffi et al., 2015; Koffi et al., 2017; Perraut et al., 2017; Yman et al., 2016; Proietti et al., 2013; Idris et al., 2017). A similar association was observed for PfMSP119 IgG, with twofold increase in seroprevalence associated with 2.49-fold (OR: 2.49; 95% CI: 1.21–5.12; p=0.014) increase in odds of gSG6 IgG seropositivity. This association was derived from 163 study-specific observations from 10 studies (Stone et al., 2012; Ya-Umphan et al., 2017; Koffi et al., 2015; Koffi et al., 2017; Perraut et al., 2017; Proietti et al., 2013; Yman et al., 2016; Badu et al., 2015; Kerkhof et al., 2016; Idris et al., 2017). Analysis of 47 study-specific observations from three studies indicated that a twofold increase in PfMSP2 IgG seroprevalence was associated with a 41% (OR: 1.41; 95% CI: 1.21–1.65; p<0.001) increase in odds of gSG6 IgG seropositivity (Ya-Umphan et al., 2017; Perraut et al., 2017; Yman et al., 2016). While 17 study-specific observations from two studies showed a twofold increase in PfMSP3 IgG seroprevalence was associated with a 2.66-fold (OR: 2.66; 95% CI: 2.36–3.00; p<0.001) increase in odds of gSG6 IgG seropositivity (Stone et al., 2012; Yman et al., 2016).
The pooled analysis of 128 study-specific observations from five studies showed that a twofold increase in PfGLURP IgG seroprevalence was associated with a 3.05-fold (OR: 3.05; 95% CI: 2.58–3.61; p<0.001) increase in odds of gSG6 IgG seropositivity (Koffi et al., 2015; Koffi et al., 2017; Ambrosino et al., 2010; Perraut et al., 2017; Kerkhof et al., 2016). And 18 study-specific observations from five studies indicated that a twofold increase in P. falciparum schizont extract IgG seropositivity was associated with a 5.69-fold (OR: 5.69; 95% CI: 0.03–1188.69; p=0.523) increase in odds of gSG6 IgG seropositivity (Idris et al., 2017; Koffi et al., 2015; Sarr et al., 2012; Perraut et al., 2017; Koffi et al., 2017).
We observed that increasing seroprevalence of IgG antibodies against PfAMA1 saw increased levels of anti-gSG6 IgG in three studies (Idris et al., 2017; Koffi et al., 2015; Koffi et al., 2017), but no association in another (Ya-Umphan et al., 2017). The levels of gSG6 IgG increased with increasing PfMSP119 IgG seroprevalence in three studies (Idris et al., 2017; Koffi et al., 2015; Badu et al., 2015), but showed no association in three other studies (Ya-Umphan et al., 2017; Koffi et al., 2017; Kerkhof et al., 2016). No association between gSG6 IgG levels and MSP2 IgG seroprevalence was observed in one study (Ya-Umphan et al., 2017). PfGLURP IgG seroprevalence and gSG6 IgG antibody levels were reported in three studies, with one study reporting increased levels (Koffi et al., 2015), one study reporting no association (Kerkhof et al., 2016), and one study reporting decreased levels of anti-gSG6 IgG with increasing anti-PfGLURP seroprevalence (Koffi et al., 2017). One study showed increasing gSG6 IgG levels with increasing P. falciparum schizont extract IgG, while three other studies showed no association (Koffi et al., 2015; Koffi et al., 2017; Sarr et al., 2012). Of note, one study provided a single seroprevalence estimate of antibodies against PfAMA1, PfMSP119, and PfMSP3 so no relationships can be drawn (Stone et al., 2012).
Antibodies against P. vivax antigens
In pooled analyses of 115 study-specific observations from two studies (Idris et al., 2017; Kerkhof et al., 2016), we observed that a twofold increase in the seroprevalence of PvAMA1 was associated with a 3.87-fold (OR: 3.87; 95% CI: 3.46–4.32; p<0.001) increase in the odds of anti-gSG6 IgG seropositivity. Furthermore, in 103 study-specific observations from two studies (Idris et al., 2017; Kerkhof et al., 2016), a twofold increase in PvMSP119 IgG seroprevalence was associated with a 2.37-fold (OR: 2.37; 95% CI: 2.26–2.50; p<0.001) increase in the odds of anti-gSG6 IgG seropositivity. However, neither study showed an association between the levels of gSG6 IgG and the seroprevalence of PvAMA1 and PvMSP119 IgG (Idris et al., 2017; Kerkhof et al., 2016).
Appendix 9
Country and study-specific predicted probability of gSG6 IgG seropositivity
In order to obtain estimates of gSG6 IgG seroprevalence for each country and study, an intercept-only three-level random-effects logistic regression was fitted to 301 study-specific observations from 22 studies. The predicted probability of gSG6 IgG seropositivity was calculated at the country level (Appendix 9—figure 1), indicating that the seroprevalence was lowest in the Pacific region (Vanuatu [31%] and Solomon Islands [32%]) and highest in Benin (72%) and Burkina Faso (65%). Furthermore, the predicted probability of gSG6 IgG seropositivity was calculated at the study level (Appendix 9—figure 2), indicating that the seroprevalence was lowest in Ambrosino et al., 2010 (13%) and highest in Drame et al., 2015 (91%).
Appendix 10
Association between alternative salivary biomarkers and exposures of interest
Our systematic review identified a paucity of studies that assessed the relationship between our exposures of interest and most alternate Anopheles salivary biomarkers (i.e. non-An. gambiae gSG6 IgG), thus preventing the estimation of a pooled association. The exceptions being that we observed that a twofold increase in HBR was associated with a 12% increase (OR: 1.12; 95% CI: 1.02–1.24; p=0.017) in odds of anti-An. funestus fSG6 IgG seropositivity (six study-specific observations from two studies; Rizzo et al., 2011a; Ali et al., 2012; Appendix 10—table 1), as well as a 12.97-fold (OR: 12.97; 95% CI: 10.95–15.36; p<0.001) and 4.04-fold (OR: 4.04; 95% CI: 3.60–4.54; p<0.001) increase in odds of anti-gSG6-P2 IgG seropositivity associated with a twofold increase in seroprevalence of PfCSP and PfGLURP IgG, respectively (115 and 116 study-specific observations from two studies, respectively; Ambrosino et al., 2010; Kerkhof et al., 2016, Appendix 10—tables 2 and 3). The associations between exposures of interest and the additional salivary biomarkers are further discussed in narrative terms in below.
Human biting rate
In addition to the increased odds of An. funestus fSG6 seropositivity with increasing HBR, the majority of studies reported a positive association between HBR and the seroprevalence and levels of anti-gSG6-P1 IgM (Drame et al., 2015), the levels of gSG6-P2 IgG (Poinsignon et al., 2008a), the seroprevalence and levels of anti-cE5 IgG (Rizzo et al., 2014a), the levels of anti-fSG6 IgG (Rizzo et al., 2011a; Ali et al., 2012), the seroprevalence and levels of anti-f5’nuc IgG (Ali et al., 2012), and the median levels of anti-An. gambiae SGE IgG and IgG4 (Drame et al., 2010b; Lawaly et al., 2012; Remoue et al., 2006). One study reported similar median levels of anti-gSG6 IgG1 across populations and time points, whilst reporting that anti-gSG6 IgG4 titre increased with increasing HBR in one of the populations, but not in the other (Rizzo et al., 2014b). Similarly, there was no consistent association between HBR and the levels of anti-cE5 IgG (Marie et al., 2015), levels of anti-An. gambiae SGE IgE (Lawaly et al., 2012) and the seroprevalence and levels of anti-g5’nuc IgG (Ali et al., 2012).
Entomological inoculation rate
Ali et al., 2012 reported higher seroprevalence and levels anti-fSG6 IgG and anti-f5’nuc IgG with increasing EIR, while anti-g5’nuc IgG seroprevalence and levels were not associated with EIR. An additional study reported gSG6-P2 IgG seroprevalence estimates of 0% for three sites, irrespective of EIR (Ambrosino et al., 2010).
Malaria prevalence
Two studies showed that increased Plasmodium spp. prevalence was associated with higher median levels of anti-An. gambiae SGE IgG (Drame et al., 2010b; Brosseau et al., 2012), while another study showed different anti-An. gambiae SGE IgG levels for very similar prevalence of malaria and slightly lower levels of anti-An. gambiae SGE IgE and IgG4 for the time point with greater malaria prevalence (Lawaly et al., 2012). Kerkhof et al., 2016 showed increasing levels of anti-gSG6-P2 IgG for higher prevalence of any Plasmodium spp. infection, while Londono-Renteria et al., 2020a showed lower levels of IgG antibodies against TRANS-P1, TRANS-P2, PEROX-P1, PEROX-P2, and PEROX-P3 in the site with higher PCR-confirmed malaria prevalence. Additionally, several case-controlled studies, and two cross-sectional study, reported median antibody levels stratified by malaria infection status. These studies show higher levels of anti-An. darlingi SGE IgG (Andrade et al., 2009), anti-An. gambiae SGE IgG (Remoue et al., 2006), anti-An. dirus SGE IgG and IgM (Waitayakul et al., 2006), and IgG antibodies against SGEs of two Colombian strains of An. albimanus in Plasmodium spp.-infected individuals compared to non-infected (Montiel et al., 2020). While Montiel et al., 2020 observed no association between anti-An. darlingi SGE IgG levels and infection status.
Antimalarial antibody seroprevalence
Our multilevel modelling indicated that there were 12.97-fold (OR: 12.97; 95% CI: 10.95–15.36; p<0.001) and 4.04-fold (OR: 4.04; 95% CI: 3.60–4.54; p<0.001) increase in odds of anti-gSG6-P2 IgG seropositivity associated with a twofold increase in the seroprevalence of PfCSP and PfGLURP IgG, respectively (Ambrosino et al., 2010; Kerkhof et al., 2016; Appendix 10—table 2 and Appendix 10—table 3). However, we observed weak positive associations between the levels of IgG antibodies against gSG6-P2 peptide and the seroprevalence of IgG antibodies against PfMSP119, PfGLURP and PvMSP119, but no association with PfCSP or PvAMA1 (Kerkhof et al., 2016).
Data availability
The current manuscript is a systematic review with multilevel modelling of study level data. The constructed dataset and associated code used for analyses are available at https://github.com/ellenakearney/Anopheles_salivary_biomarker_systematic_review (copy archived at swh:1:rev:eadb7ab1cbf93c730fb463f43e1b13b9ae16ddd3). Data for 'Dominant malaria vector species globally, 2010' (Sinka et al., 2012) and 'Plasmodium falciparum parasite rate in 2-10 year olds globally, 2000-2017' (Weiss et al., 2019) can be accessed from the Malaria Atlas Project Explorer (https://malariaatlas.org/explorer/#/).
-
GitHubID Anopheles_salivary_biomarker_systematic_review. Anopheles salivary biomarkers review.dta.
References
-
Malaria transmission in two rural communities in the forest zone of GhanaParasitology Research 108:1465–1471.https://doi.org/10.1007/s00436-010-2195-1
-
Induced immunity against the mosquito Anopheles stephensi: reactivity characteristics of immune seraMedical and Veterinary Entomology 13:53–64.https://doi.org/10.1046/j.1365-2915.1999.00143.x
-
Haematophagous arthropod saliva and host defense system: a tale of tear and bloodAnais Da Academia Brasileira de Ciencias 77:665–693.https://doi.org/10.1590/s0001-37652005000400008
-
Biomarkers for susceptibility to infection and disease severity in human malariaMemorias Do Instituto Oswaldo Cruz 106 Suppl 1:70–78.https://doi.org/10.1590/s0074-02762011000900009
-
An updated catalogue of salivary gland transcripts in the adult female mosquito, Anopheles gambiaeThe Journal of Experimental Biology 208:3971–3986.https://doi.org/10.1242/jeb.01849
-
Detection of immunogenic proteins from Anopheles sundaicus salivary glands in the human serumRevista Da Sociedade Brasileira de Medicina Tropical 48:410–416.https://doi.org/10.1590/0037-8682-0185-2015
-
Identification of antigenic proteins from salivary glands of female Anopheles maculatus by proteomic analysisAsian Pacific Journal of Tropical Biomedicine 6:924–930.https://doi.org/10.1016/j.apjtb.2016.08.012
-
Immune interactions between mosquitoes and their hostsParasite Immunology 28:143–153.https://doi.org/10.1111/j.1365-3024.2006.00805.x
-
Immunoglobulin G antibody profiles against Anopheles salivary proteins in domestic animals in SenegalJournal of Medical Entomology 48:691–693.https://doi.org/10.1603/me10183
-
Characterization of Aedes communis, Aedes aegypti and Anopheles stephensi mosquito saliva antigens by immunoblottingInternational Archives of Allergy and Immunology 112:169–174.https://doi.org/10.1159/000237450
-
Function and evolution of a mosquito salivary protein familyThe Journal of Biological Chemistry 281:1935–1942.https://doi.org/10.1074/jbc.M510359200
-
An insight into the sialome of Anopheles funestus reveals an emerging pattern in anopheline salivary protein familiesInsect Biochemistry and Molecular Biology 37:164–175.https://doi.org/10.1016/j.ibmb.2006.11.005
-
Immunological aspects of the immune response induced by mosquito allergensInternational Archives of Allergy and Immunology 165:271–282.https://doi.org/10.1159/000371349
-
Spatial repellents and malaria transmission in an endemic area of Cambodia with high mosquito net usageMalariaWorld Journal 8:11.
-
Detection of influential observation in linear regressionTechnometrics: A Journal of Statistics for the Physical, Chemical, and Engineering Sciences 19:15.https://doi.org/10.2307/1268249
-
Characterization of immunogenic proteins in Anopheles gambiae salivary glands and their potential use as a marker of exposure to malariaThe American Journal of Tropical Medicine and Hygiene 79:320.
-
Spatio-temporal dynamics of asymptomatic Malaria: bridging the gap between annual Malaria resurgences in a sahelian environmentThe American Journal of Tropical Medicine and Hygiene 97:1761–1769.https://doi.org/10.4269/ajtmh.17-0074
-
Mosquito-borne diseases and omics: salivary gland proteome of the female aedes aegypti mosquitoOmics : A Journal of Integrative Biology 21:45–54.https://doi.org/10.1089/omi.2016.0160
-
Salivary gland transcriptome analysis during Plasmodium infection in malaria vector Anopheles stephensiInternational Journal of Infectious Diseases 13:636–646.https://doi.org/10.1016/j.ijid.2008.07.027
-
The Anopheles gambiae transcriptome - a turning point for malaria controlInsect Molecular Biology 26:140–151.https://doi.org/10.1111/imb.12289
-
Biomarkers of vector bites: Arthropod immunogenic salivary proteins in vector-borne diseases controlGeneral Methods in Biomarker Research and Their Applications 1:1177–1205.https://doi.org/10.1007/978-94-007-7696-8
-
Human antibody response to Anopheles gambiae saliva: A new immuno-epidemiological marker to evaluate the effectiveness of insecticides treated nets (ITNs)The American Journal of Tropical Medicine and Hygiene 79:358.
-
Human antibody response to Anopheles gambiae saliva: an immuno-epidemiological biomarker to evaluate the efficacy of insecticide-treated nets in malaria vector controlThe American Journal of Tropical Medicine and Hygiene 83:115–121.https://doi.org/10.4269/ajtmh.2010.09-0684
-
Human antibody response to Anopheles salivary gSG6-P1 peptide: New immuno-epidemiological tool for evaluating the efficacy of insecticides treated nets (ITNs) in malaria vector controlThe American Journal of Tropical Medicine and Hygiene 83:260.https://doi.org/10.4269/ajtmh.2010.83.211
-
New salivary biomarkers of human exposure to malaria vector bitesAnopheles Mosquitoes: New Insights into Malaria Vectors 1:755–795.https://doi.org/10.5772/3392
-
Immunological targeting of critical insect AntigensAmerican Entomologist 48:150–162.https://doi.org/10.1093/ae/48.3.150
-
Toward a catalog for the transcripts and proteins (sialome) from the salivary gland of the malaria vector Anopheles gambiaeThe Journal of Experimental Biology 205:2429–2451.https://doi.org/10.1242/jeb.205.16.2429
-
The immunomodulatory factors of bloodfeeding arthropod salivaParasite Immunology 22:319–331.https://doi.org/10.1046/j.1365-3024.2000.00309.x
-
Prospects for malaria elimination in Mesoamerica and HispaniolaPLOS Neglected Tropical Diseases 9:e0003700.https://doi.org/10.1371/journal.pntd.0003700
-
The innate and adaptive response to mosquito saliva and Plasmodium sporozoites in the skinAnnals of the New York Academy of Sciences 1342:37–43.https://doi.org/10.1111/nyas.12661
-
Assessing risk of bias in prevalence studies: modification of an existing tool and evidence of interrater agreementJournal of Clinical Epidemiology 65:934–939.https://doi.org/10.1016/j.jclinepi.2011.11.014
-
Proteomics of anopheles vectors of malariaTrends in Parasitology 34:961–981.https://doi.org/10.1016/j.pt.2018.08.009
-
Characterization of human IgE and mouse IgG1 responses to allergens in three mosquito species by immunoblotting and ELISAInternational Archives of Allergy and Immunology 126:206–212.https://doi.org/10.1159/000049515
-
Epidemiological consequences of immune sensitisation by pre-exposure to vector salivaPLOS Neglected Tropical Diseases 11:10.https://doi.org/10.1371/journal.pntd.0005956
-
Members of the salivary gland surface protein (SGS) family are major immunogenic components of mosquito salivaThe Journal of Biological Chemistry 286:40824–40834.https://doi.org/10.1074/jbc.M111.280552
-
Immunological consequences of arthropod vector-derived salivary factorsEuropean Journal of Immunology 41:3396–3400.https://doi.org/10.1002/eji.201190075
-
At the interface between parasite and host: the salivary glands of the African malaria vector Anopheles gambiaeParassitologia 48:573–580.
-
The Anopheles gambiae salivary protein gSG6: an anopheline-specific protein with a blood-feeding roleInsect Biochemistry and Molecular Biology 39:457–466.https://doi.org/10.1016/j.ibmb.2009.04.006
-
Antibody response against Anopheles albimanus (Diptera: Culicidae) salivary protein as a measure of mosquito bite exposure in HaitiJournal of Medical Entomology 47:1156–1163.https://doi.org/10.1603/me09240
-
Long-lasting permethrin-impregnated clothing protects against mosquito bites in outdoor workersThe American Journal of Tropical Medicine and Hygiene 93:869–874.https://doi.org/10.4269/ajtmh.15-0130
-
Identification and Pilot Evaluation of Salivary Peptides from Anopheles albimanus as Biomarkers for Bite Exposure and Malaria Infection in ColombiaInternational Journal of Molecular Sciences 21:E691.https://doi.org/10.3390/ijms21030691
-
Antibody Responses Against Anopheles darlingi Immunogenic Peptides in Plasmodium Infected HumansFrontiers in Cellular and Infection Microbiology 10:455.https://doi.org/10.3389/fcimb.2020.00455
-
A brief review on biomarkers and proteomic approach for malaria researchAsian Pacific Journal of Tropical Medicine 8:253–262.https://doi.org/10.1016/S1995-7645(14)60327-8
-
Inhibition of the complement system by saliva of Anopheles (Nyssorhynchus) aquasalisInsect Biochemistry and Molecular Biology 92:12–20.https://doi.org/10.1016/j.ibmb.2017.11.004
-
Allergy to hematophagous arthropods bitesCurrent Dermatology Reports 3:6–12.https://doi.org/10.1007/s13671-013-0065-7
-
Operational assessment of long-lasting insecticidal nets by using an anopheles salivary biomarker of human-vector contactThe American Journal of Tropical Medicine and Hygiene 95:1376–1382.https://doi.org/10.4269/ajtmh.15-0541
-
Novel salivary gland allergens from tropical mosquito species and IgE reactivity in allergic patientsThe World Allergy Organization Journal 13:100099.https://doi.org/10.1016/j.waojou.2020.100099
-
Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae)Biological Journal of the Linnean Society 97:854–866.https://doi.org/10.1111/j.1095-8312.2009.01226.x
-
Immunoblot analysis of salivary allergens in 10 mosquito species with worldwide distribution and the human IgE responses as these allergensThe Journal of Allergy and Clinical Immunology 101:498–505.https://doi.org/10.1016/s0091-6749(98)70357-4
-
Immune responses to mosquito saliva in 14 individuals with acute systemic allergic reactions to mosquito bitesThe Journal of Allergy and Clinical Immunology 114:1189–1194.https://doi.org/10.1016/j.jaci.2004.08.014
-
Mosquito allergy: immune mechanisms and recombinant salivary allergensInternational Archives of Allergy and Immunology 133:198–209.https://doi.org/10.1159/000076787
-
Mosquito allergy and mosquito salivary allergensProtein Peptide Lett 14:975–981.https://doi.org/10.2174/092986607782541088
-
Mosquito salivary gland antigens identified by circulating human antibodiesArchives of Dermatology 125:219–222.https://doi.org/10.1001/archderm.1989.01670140071012
-
Mosquito biting modulates skin response to virus infectionTrends in Parasitology 33:645–657.https://doi.org/10.1016/j.pt.2017.04.003
-
Functional analyses yield detailed insight into the mechanism of thrombin inhibition by the antihemostatic salivary protein cE5 from Anopheles gambiaeThe Journal of Biological Chemistry 292:12632–12642.https://doi.org/10.1074/jbc.M117.788042
-
Human IgG response to Anopheles gambiae salivary proteins as an immuno-epidemiological marker of exposure to malaria vector bitesThe American Journal of Tropical Medicine and Hygiene 79:218.
-
Development of a new biomarker of exposure to Anopheles bites based on human antibody responses to salivary proteins: From the concept to the applicationsThe American Journal of Tropical Medicine and Hygiene 83:184.
-
First attempt to validate the gSG6-P1 salivary peptide as an immuno-epidemiological tool for evaluating human exposure to Anopheles funestus bitesTropical Medicine & International Health 15:1198–1203.https://doi.org/10.1111/j.1365-3156.2010.02611.x
-
Development of a biomarker of exposure to Anopheles bites based on human antibody response to salivary proteins: From concept to application in the fieldPathogens and Global Health 107:455.
-
Evaluation of the antibody response to Anopheles salivary antigens as a potential marker of risk of malariaTransactions of the Royal Society of Tropical Medicine and Hygiene 100:363–370.https://doi.org/10.1016/j.trstmh.2005.06.032
-
Role of arthropod saliva in blood feeding: sialome and post-sialome perspectivesAnnual Review of Entomology 48:73–88.https://doi.org/10.1146/annurev.ento.48.060402.102812
-
The Anopheles gambiae cE5, a tight- and fast-binding thrombin inhibitor with post-transcriptionally regulated salivary-restricted expressionInsect Biochemistry and Molecular Biology 42:610–620.https://doi.org/10.1016/j.ibmb.2012.04.008
-
Sialogenins and immunomodulators derived from blood feeding parasitesToxins and Hemostasis 1:131–152.https://doi.org/10.1007/978-90-481-9295-3
-
Epidemiological applications of assessing mosquito exposure in a malaria-endemic areaArthropod Vector: Controller of Disease Transmission, Vector Saliva-Host Pathogen Interactions2017 2:209–229.https://doi.org/10.1016/B978-0-12-805360-7.00012-5
-
Spatial risk of urban exposure to Anopheles and Aedes mosquito bites in Africa using salivary antibody-based biomarkersThe American Journal of Tropical Medicine and Hygiene 99:47.
-
Differential acquisition of human antibody responses to Plasmodium falciparum according to intensity of exposure to Anopheles bitesTransactions of the Royal Society of Tropical Medicine and Hygiene 106:460–467.https://doi.org/10.1016/j.trstmh.2012.05.006
-
A global map of dominant malaria vectorsParasites & Vectors 5:69.https://doi.org/10.1186/1756-3305-5-69
-
Salivary gland proteome of the human malaria vector, Anopheles campestris-like (Diptera: Culicidae)Parasitology Research 112:1065–1075.https://doi.org/10.1007/s00436-012-3233-y
-
Increased threat of urban malaria from anopheles stephensi mosquitoes, AfricaEmerging Infectious Diseases 25:1431–1433.https://doi.org/10.3201/eid2507.190301
-
Evaluation of malaria urban risk using an immuno-epidemiological biomarker of human exposure to anopheles bitesThe American Journal of Tropical Medicine and Hygiene 98:1353–1359.https://doi.org/10.4269/ajtmh.17-0231
-
Behavioural heterogeneity of Anopheles species in ecologically different localities in Southeast Asia: a challenge for vector controlTropical Medicine & International Health 10:251–262.https://doi.org/10.1111/j.1365-3156.2004.01378.x
-
Mosquito saliva alone has profound effects on the human immune systemPLOS Neglected Tropical Diseases 12:e0006439.https://doi.org/10.1371/journal.pntd.0006439
-
Anopheles gambiae circumsporozoite protein-binding protein facilitates plasmodium infection of mosquito salivary glandsThe Journal of Infectious Diseases 208:1161–1169.https://doi.org/10.1093/infdis/jit284
-
Use of an anopheles salivary biomarker to assess malaria transmission risk along the Thailand-Myanmar borderThe Journal of Infectious Diseases 215:396–404.https://doi.org/10.1093/infdis/jiw543
-
Anopheles salivary biomarker as a proxy for estimating Plasmodium falciparum malaria exposure on the Thailand–Myanmar borderThe American Journal of Tropical Medicine and Hygiene 99:350–356.https://doi.org/10.4269/ajtmh.18-0081
-
New markers in Anopheles gambiae salivary glands after Plasmodium berghei infectionVector Borne and Zoonotic Diseases 13:119–127.https://doi.org/10.1089/vbz.2012.0964
Article and author information
Author details
Funding
National Health and Medical Research Council (1134989)
- Julie A Simpson
- Freya JI Fowkes
National Health and Medical Research Council (1166753)
- Freya JI Fowkes
National Health and Medical Research Council (1196068)
- Julie A Simpson
Australian Government (Australian Government Research Training Program Scholarship)
- Ellen A Kearney
Wellcome Trust (220211)
- Victor Chaumeau
Victorian State Government (Operational Infrastructure Support Program received by the Burnet Institute.)
- Ellen A Kearney
- Paul A Agius
- Julia C Cutts
- Freya JI Fowkes
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Franck Remoue, Anne Poinsignon, Bruno Arcà, Vincent Corbel, Richard Paul, André Sagna, Kingsley Badu, Berlin Londono-Renteria, Jacques Derek Charlwood, William Stone, Chris Drakeley, Karen Kerkhof, Sylvie Manguin, and Yunita Armiyanti for responding to requests for further information/data for the systematic review.
Copyright
© 2021, Kearney et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,220
- views
-
- 209
- downloads
-
- 10
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
- Genetics and Genomics
Alzheimer’s disease (AD) is a complex degenerative disease of the central nervous system, and elucidating its pathogenesis remains challenging. In this study, we used the inverse-variance weighted (IVW) model as the major analysis method to perform hypothesis-free Mendelian randomization (MR) analysis on the data from MRC IEU OpenGWAS (18,097 exposure traits and 16 AD outcome traits), and conducted sensitivity analysis with six models, to assess the robustness of the IVW results, to identify various classes of risk or protective factors for AD, early-onset AD, and late-onset AD. We generated 400,274 data entries in total, among which the major analysis method of the IVW model consists of 73,129 records with 4840 exposure traits, which fall into 10 categories: Disease, Medical laboratory science, Imaging, Anthropometric, Treatment, Molecular trait, Gut microbiota, Past history, Family history, and Lifestyle trait. More importantly, a freely accessed online platform called MRAD (https://gwasmrad.com/mrad/) has been developed using the Shiny package with MR analysis results. Additionally, novel potential AD therapeutic targets (CD33, TBCA, VPS29, GNAI3, PSME1) are identified, among which CD33 was positively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. TBCA and VPS29 were negatively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. GNAI3 and PSME1 were negatively associated with the main outcome traits of AD, as well as with LOAD, but had no significant causal association with EOAD. The findings of our research advance our understanding of the etiology of AD.
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.