Demographic patterns of human antibody levels to Simulium damnosum s.l. saliva in onchocerciasis-endemic areas: An indicator of exposure to vector bites

Abstract

Background

In onchocerciasis endemic areas in Africa, heterogenous biting rates by blackfly vectors on humans are assumed to partially explain age- and sex-dependent infection patterns with Onchocerca volvulus. To underpin these assumptions and further improve predictions made by onchocerciasis transmission models, demographic patterns in antibody responses to salivary antigens of Simulium damnosum s.l. are evaluated as a measure of blackfly exposure.

Methodology/Principal findings

Recently developed IgG and IgM anti-saliva immunoassays for S. damnosum s.l. were applied to blood samples collected from residents in four onchocerciasis endemic villages in Ghana. Demographic patterns in antibody levels according to village, sex and age were explored by fitting generalized linear models. Antibody levels varied between villages but showed consistent patterns with age and sex. Both IgG and IgM responses declined with increasing age. IgG responses were generally lower in males than in females and exhibited a steeper decline in adult males than in adult females. No sex-specific difference was observed in IgM responses.

Conclusions/Significance

The decline in age-specific antibody patterns suggested development of immunotolerance or desensitization to blackfly saliva antigen in response to persistent exposure. The variation between sexes, and between adults and youngsters may reflect differences in behaviour influencing cumulative exposure. These measures of antibody acquisition and decay could be incorporated into onchocerciasis transmission models towards informing onchocerciasis control, elimination, and surveillance.

Introduction

Female blackflies of the Simulium damnosum sensu lato (s.l.) complex are the predominant vectors of Onchocerca volvulus in Africa. This filarial nematode causes human onchocerciasis, commonly known as river blindness and proposed for elimination of transmission (EOT) in 12 countries by 2030 [1]. Repeat exposure to bites of infective blackflies is a key driver of parasite acquisition, and high biting rates are important determinants of transmission intensity and resurgence following control interventions, particularly those based on mass drug administration (MDA) of ivermectin [2,3]. Therefore, ‘stop-MDA surveys’ and post-treatment surveillance protocols would be improved by monitoring exposure to vector bites in addition to exposure to the parasite and assessment of residual infection in informative age and sex groups depending on the epidemiological setting [4,5]. Annual biting rates (ABR, the number of bites/person/year) are predominantly estimated by performing human landing catches (HLCs) [6], but this is both labour-intensive and often considered unethical. Furthermore, HLC methods likely lead to biased estimates of ABRs as they are typically performed on a few adult males at the riverbank close to blackfly breeding sites, where biting rates are likely to be higher and vector collectors are maximally exposed [7]. Thus, HLCs may not capture the true exposure representative of the community, age or sex groups or their activity patterns. Furthermore, when blackfly densities are high, HLCs may be unable to capture all biting events thereby underestimating the biting rate [8,9]. Applying novel methods to measure individual and population exposure to blackfly vector bites would greatly inform our understanding of exposure patterns to onchocerciasis and inform mathematical transmission modelling of control and elimination [2].

Proteins in the saliva of blood-seeking arthropod vectors of human and veterinary diseases provoke an immunomodulatory response in vertebrates following their exposure to vector bites (e.g. by mosquitoes, sand flies, triatomine bugs, and tsetse flies) [10]. Measurement of host anti-saliva antibodies has proven a useful surrogate marker to monitor individual host biting rates [11–15], seasonal variations in exposure and infection [16,17], and to evaluate vector control interventions [18–22]. Immunoassays to measure human anti-saliva antibody responses to blackfly bites have not been available until recently developed by the authors, specifically to measure human IgG and IgM antibody responses to saliva of S. damnosum s.l. which is the main vector in Africa including the study sites in Ghana [23].

In the present study we report on the demographic patterns revealed by applying these assays to community residents in the Bono East region of Ghana, a savannah setting. This area was under vector control during the Onchocerciasis Control Programme in West Africa (OCP, 1974–2002) [24,25], and is where the first ivermectin community trials were conducted in the late 1980’s [26]. However, transmission persists despite many years of annual MDA, even after adopting biannual MDA in 2010 [27,28]. With the risk of residual transmission leading to possible resurgence of human infections if MDA were to stop [3], the area has been the subject of comprehensive entomological [27,29–31], parasitological [32] and parasite genomic [33] research, and a source of valuable epidemiological and entomological data for parameterizing onchocerciasis transmission models [34–36]. Due to the lack of reliable tools to independently measure human exposure to vector bites, mathematical models have assumed age- and sex-dependent exposure rates based on patterns of human infection with O. volvulus. However, modelling outputs used to inform EOT and surveillance strategies are sensitive to such patterns [4,37]. The aim of this study is to apply our novel human IgG and IgM immunoassays as indicators of individual biting exposure, both between and within (by age and sex) onchocerciasis-endemic communities. The generated demographic patterns will help to scrutinize mathematical modelling assumptions.

Material and methods

Ethical statement

Ethical clearance was obtained from the Council for Scientific and Industrial Research (CSIR) Institutional Review Board (RPN008/CSIR-IRB/2019) in Accra, Ghana. Residents were informed of the objectives of the study and participants provided fully informed written consent; parents or guardians provided consent for <18-year-olds. At the time of the study, ivermectin treatment against onchocerciasis was administered as part of the on-going national biannual MDA programme [38]. Each participant was rewarded with a bar of washing soap, a commercial sachet of malt drink powder (Milo, Nestlé) and a can of condensed milk. Blackfly collections were carried out by local OCP vector collectors following standard HLC techniques [39]. Individual identifiers on data records were anonymized prior to analyses and data storage.

Study area

The study was conducted during the wet season in August 2019 in rural villages in an onchocerciasis-endemic area of the Pru river basin, Pru District, in the Bono East region of Ghana. Fig 1 provides a map of the study area and the location of the study communities, created using the ‘ggplot2’ and the ‘rworldmap’ package in R software [40,41]. The shapefiles for the base layers of the map were retrieved from the Database of Global Administrative Areas (GADM, https://gadm.org/download_country.html) and from DIVA-GIS (http://www.diva-gis.org) under a CC BY license, with permission from GADM and DIVA-GIS.

[Figure omitted. See PDF.]

Fig 1. Map of Ghana indicating study communities in the Pru District, Bono East region.

Shapefiles for the base layers of the map retrieved from GADM (https://gadm.org/download_country.html) and DIVA-GIS (http://www.diva-gis.org) under a CC BY license, with permission from GADM and DIVA-GIS.

https://doi.org/10.1371/journal.pntd.0010108.g001

Recruitment, blood sampling and demographic information

Villages were selected on the criteria that (a) there was evidence of exposure to S. damnosum s.l. vector bites [30,42], (b) they were under current ivermectin treatment as part of the national onchocerciasis control programme [43], (c) the villages comprised >200 residents, and (d) the villages represented variable ABRs based on previous records [30,42]. Five villages were selected; due to the proximity of two of the villages, these were grouped into a single cluster (Asubende and Senyase (ASU/SEN) resulting in four village clusters. These included ASU/SEN, with high biting rates (400–850 flies/person/month), Beposo [BEP] with moderate values (100–350), and Fawoman [FAW] and Ohiampe [OHI] with lower biting rates (25–175 flies/person/month) [30,42].

To obtain blood samples from the village residents, two recruitment approaches were adopted: in villages with <300 residents (ASU/SEN, FAW, and OHI) all residents were invited to participate. In the larger village (BEP) an age/sex stratified sample was identified based on assigning random numbers generated using R software [44], to the pages of the paper-based census records compiled by the MDA programme. The households contained on randomly selected pages were noted and invited to participate until reaching an estimated sample of houses containing 250 individuals as required from statistical calculations (i.e. ~1000 individuals in total across the four village clusters). Children younger than 4 years old were not recruited.

A total of 958 participants were finally recruited, all tested for IgG antibodies, and a subset of 500 individuals tested for IgM antibodies (Table 1). For age stratification purposes, the identified populations in all villages were divided into eight age categories (in years): 5–10, 11–20, 21–30, 31–40, 41–50, 51–60, 61–70, > 71, aiming to recruit 32 residents per age class in each village, split equally between the sexes. In the case of IgM, samples were selected by random number generation assigned to the ordered age-stratified full list of samples and aimed to test around 40–60% of the samples tested by IgG per age class, with some variability in percentage tested in the oldest age group. The numbers tested per age and sex strata per village cluster are shown in Figs A and B and Table C in S1 File. The geographical coordinates, elevation and number of individuals sampled and tested for each immunoglobulin per village cluster are shown in Table 1.

[Figure omitted. See PDF.]

Table 1. Study communities, their geographical coordinates, elevation and numbers of people sampled and tested for IgG and IgM per community in the Pru District, Ghana.

https://doi.org/10.1371/journal.pntd.0010108.t001

Participants in each village/cluster were invited to the local school and assigned a subject/sample identification number written on a personalized card. Two to four ml of blood were collected into Ethylene Diamine Tetra Acetic acid (EDTA) tubes by venipuncture and kept cool for 2–3 hours in an insulated cool box until centrifugation at 2500 rpm for 15 min to separate the plasma. Samples were stored at 7° C until ELISA testing. Metadata on name, age, sex, number of years of residence, house number, name of household head, and history of clinical onchocerciasis were collected.

Immunoassays to measure human exposure to blackfly bites

Enzyme-linked immunosorbent assays (ELISA) previously developed were performed to measure anti-S. damnosum s.l. saliva IgG and IgM human responses [23]. Briefly, host-seeking S. damnosum s.l. females were collected following standard OCP vector collector techniques in one study location (ASU), near the Pru river [39]. All collected flies were stored in a cool box until dissected on the same day. The collected flies were anesthetized in a –20° C freezer for 10 min after which their salivary glands were removed, and aliquots stored in Tris-buffered saline (TBS) (one gland per μl TBS, pH 7.5) at –20° C until further use. Flat-bottom 96-well microtiter plates (ThermoFisher Scientific) were coated with blackfly salivary gland homogenate (SGH) with 0.2 μg SGH/well (for IgG ELISA) or 0.025 μg SGH/well (for IgM ELISA) and incubated overnight at 4° C. The plates were blocked with 6% non-fat dried milk (Bio-Rad) in phosphate-buffered saline (PBS) with 0.05% Tween 20 (PBS-Tw) and incubated with plasma diluted 1/100 (IgG ELISA) or 1/50 (IgM ELISA) in 2% non-fat dried milk. The plates were washed and incubated with peroxidase-conjugated anti-human IgG (1/1,000) or IgM antibody (1/70,000) (Sigma-Aldrich; Bethyl Laboratories, Inc). The ELISA was developed using an orthophenylendiamine (OPD) solution in a phosphate-citrate buffer (pH 5.5) with 0.1% hydrogen peroxide. The reaction was stopped after 5 min with 10% sulfuric acid and the absorbance, optical density (OD) value was measured at 492 nm using a Tecan Infinite M200 microplate reader (Schoeller). Further details of blackfly collection, dissection, and the immunoassays were previously reported [23].

Data standardization

All plasma samples were tested in duplicate. Samples with a coefficient of variation (CoV) of more than 20% were retested. Each plate included a blank control, of which the OD value was subtracted from the sample OD values. A set of two positive (PC) and two negative control (NC) samples were included in each plate to correct for inter-plate variability according to the following formula: Standardized Optical Density (SOD) = ODsample/(average ODPC−average ODNC). Furthermore, a PC sample was titrated in duplicate on three separate plates at seven serial dilutions from 1/50 to 1/3,200. The average of the three log-logit transformed standard curves was used to convert and standardize sample SOD values which are reported below as anti-saliva antibody arbitrary units/ml.

Sample size calculation

Sample sizes were calculated to achieve 90% statistical power (with type I error, α = 0.05) to detect a difference in mean IgG SOD antibody levels between males and females, for an unpaired two-sample effect size DCohen of 0.23 and a variance in SOD of 0.16 [45]. By adding 15% to the calculated sample size as a correction factor for subsequent non-parametric statistical testing, 460 people per sex group were required [46]. This estimate is also based on the reasonable assumption that the differences in the mean antibody responses of the sexes between clusters would be minimal as confirmed by subsequent analyses (test of village cluster × sex interaction term: P>0.243). The equivalent statistical power to detect a difference in IgM responses between sexes with an effect size DCohen of 0.3 and a variance in SOD of 0.16 was 90%. Calculations were made in the R package ‘power’ using a two-sided t-test.

Statistical analyses

Differences between clusters or sex categories were statistically tested by Wilcoxon rank sum or Kruskal-Wallis tests with post-hoc Holm adjustment for multiple comparisons. Changes in immune responses with age and/or sex were tested by fitting generalized linear models (GLM) to the IgG and IgM anti-blackfly saliva antibody values, where a gamma distribution and log-link function gave the best fit by log-likelihood goodness-of-fit statistics. Age, sex and cluster interaction terms were tested, treating age as a continuous variable. Residence duration (years) and age showed high multicollinearity; hence only age was retained in the models. Correlation coefficients (rs) between median IgG and IgM responses by age and sex were estimated by the Spearman’s rank method. All statistical analyses were performed in R software [44], and graphical representations created using the ‘ggplot2’ package in R [40].

Results

Univariate analyses of population antibody response distributions

A total of 958 participants (186–263 per village cluster) between the ages of ≥ 5 to < 96 years old were recruited and sampled for anti-saliva IgG antibodies; demographic characteristics are summarized in Table 2 and in Figs A and B in S1 File. No differences were observed in the age compositions between village clusters or between sexes, nor in the frequencies of participants per age group (Figs A and B in S1 File). Thus, it is unlikely that the variations in age-sex compositions of village clusters were a source of statistical bias.

Univariate analyses of demographic variables detected significant differences between the median IgG responses among village clusters (P<0.001) (Fig 2A), and the median IgG responses of males and females (P<0.05) (Fig 2B). Equivalent differences were not observed in IgM responses (Fig 3); however, both IgG and IgM median antibody responses declined with age (P<0.001) (Figs 2D and 3D). A breakdown of the IgG and IgM antibody distribution with age per individual cluster is visualized in Figs D and E in S1 File.

[Figure omitted. See PDF.]

Fig 2. Human anti-blackfly saliva IgG antibody responses according to cluster, sex, and age.

Boxplots showing the distributions of IgG responses (A) by cluster, (B) by sex, (C) by age, and scatterplot (D) by age and sex (points) showing the best-fit lines (solid lines). IgG levels are shown as units/ml. F: female; M: male. Village clusters: ASU/SEN: Asubende/ Senyase; BEP: Beposo; FAW: Fawoman; OHI: Ohiampe.

https://doi.org/10.1371/journal.pntd.0010108.g002

[Figure omitted. See PDF.]

Fig 3. Human anti-blackfly saliva IgM antibody responses according to cluster, sex and age.

Boxplots showing the distributions of IgM responses (A) by cluster, (B) by sex, (C) by age, and scatterplot (D) by age and sex (points) showing the best-fit lines (solid lines). IgM levels are shown as units/ml. F: female; M: male. Village clusters: ASU/SEN: Asubende/ Senyase; BEP: Beposo; FAW: Fawoman; OHI: Ohiampe.

https://doi.org/10.1371/journal.pntd.0010108.g003

[Figure omitted. See PDF.]

Table 2. Summary statistics of participants and median antibody titres by immunoglobulin and village/cluster.

https://doi.org/10.1371/journal.pntd.0010108.t002

Multivariate analyses of IgG responses according to village cluster, sex and age

Accounting for cluster ID, sex, age, and the age x sex interaction term, multivariate models revealed that IgG responses were significantly lower in BEP than in ASU/SEN (P<0.01), similar between ASU/SEN and FAW (P = 0.406), and higher in OHI than in ASU/SEN (P<0.001) (Figs 4A and 5). Generally, males exhibited lower IgG responses than females (P<0.001) (Figs 4B and 5). Fig 5 illustrates the (exponentiated) regression coefficient estimates of the final model (summarized in Table F in S1 File).

[Figure omitted. See PDF.]

Fig 4. Regression effect plots for all explanatory variables contained in the final model for the IgG responses.

Each plot visualizes the effect of a specific explanatory variable while fixing the others to their reference level. The panels show the variations in responses between village clusters (A), sex (B), age (C), and an interaction between age and sex (D). The vertical lines in panel A and B and the shaded regions around the lines in panel C and D represent the 95% confidence intervals. M = Male; F = Female. Village clusters: ASU/SEN: Asubende/Senyase; BEP: Beposo; FAW: Fawoman; OHI: Ohiampe.

https://doi.org/10.1371/journal.pntd.0010108.g004

[Figure omitted. See PDF.]

Fig 5. Exponentiated regression coefficient estimates for the final model per assay.

The regression coefficient estimates for the final model are shown for IgG responses (A), and IgM responses (B). The circles represent the exponentiated mean coefficient estimate, the value shown above the corresponding circle, and the horizontal lines represent the 95% confidence intervals of the estimate. Significance levels *** P<0.001; ** P<0.01; * P<0.05. The referents were Asubende/Senyase for village clusters, and females for both sex and the age × sex interaction term. Village clusters: BEP: Beposo; FAW: Fawoman, OHI: Ohiampe.

https://doi.org/10.1371/journal.pntd.0010108.g005

IgG anti-saliva antibody responses decreased with increasing age in both sexes (P<0.001) (Figs 4C and 5), but the response in males declined at a faster rate than in females (test of age × sex interaction term; P<0.001) (Figs 4D and 5). Both adult male and female participants ≥18 years old exhibited lower IgG antibody levels than those aged <18 years (Pfemales<0.001; Pmales<0.001) (Fig 5). Males and females in pre-teenage years showed a similar antibody decline with increasing age (age × sex interaction: P = 0.328), whereas in adults the decline was more rapid in males than in females (age × sex interaction: P<0.01) (Fig 4D).

Multivariate analyses of IgM responses according to village cluster, sex or age

Equivalent multivariate analyses of anti-saliva IgM antibody responses for 500 recruits (Table 1) indicated no significant variation between study clusters (PBEP = 0.403; PFAW = 0.301; POHI = 0.632) or sex (P = 0.405) (Figs 5, 6A and 6B). However, there was a general decline in response magnitude with increasing age (P<0.01) (Figs 5 and 6C), but in contrast to IgG, the decline with age was not significantly different between sexes (test of age × sex interaction term: P = 0.726) (Fig 6D). Fig 5 illustrates the (exponentiated) regression coefficient estimates of the final model (summarized in Table F in S1 File).

[Figure omitted. See PDF.]

Fig 6. Regression effect plots for all explanatory variables contained in the final model for the IgM responses.

Each plot visualizes the effect of a specific explanatory variable while fixing the others. The panels show the variations in responses between village clusters (A), sex (B), age (C), and an interaction between age and sex (D). The vertical lines in panel A and B and the shaded regions around the lines in panel C and D represent the 95% confidence intervals. M = Male; F = Female. Village clusters: ASU/SEN: Asubende/Senyase; BEP: Beposo; FAW: Fawoman; OHI: Ohiampe.

https://doi.org/10.1371/journal.pntd.0010108.g006

Correlation between IgG and IgM responses

The correlation between IgG and IgM antibody responses was generally low (rs = 0.10, P<0.05), and not dissimilar considering females or males alone (females: rs = 0.12, P<0.1; males: rs = 0.08, P = 0.21) (Fig F in S1 File).

Discussion

This study measured human antibody responses against blackfly saliva as an indicator of individual bite exposure; an approach that has been validated for several vector–host systems with the exception of blackflies [47,48]. Simuliids are important disease vectors of human and bovine onchocerciasis [49], mansonelliasis caused by Mansonella ozzardi [50], and vesicular stomatitis virus [51], yet they have not received as much attention as sand fly, mosquito or tick vectors. Only recently, we successfully developed novel immunoassays against blackfly saliva and this is the first study to apply these tools to an epidemiological setting. The appropriate use of these assays will also improve our understanding of O. volvulus transmission dynamics and will be used to monitor changes in biting rates and the success of vector control interventions [23]. Heterogenous exposure to blackfly bites is an important determinant of the observed age- and sex-dependent profiles of O. volvulus infection [5]. Evaluating age- and sex-associated exposure patterns using empirically collected data can help inform transmission dynamics models, as these currently use assumed age- and sex-dependent exposure patterns that have, at best, been derived from fitting the models to age- and sex-specific infection profiles [5]. Such assumptions influence the choice of informative age groups for serological monitoring in stop-MDA surveys [4], and have implications for understanding the contribution of currently untreated groups to transmission and morbidity [37], as well as for the design and evaluation of potential prophylactic strategies [52]. Testing such assumptions with independently collected data is crucial for better parameterization of transmission models and improved design of epidemiological surveys and surveillance strategies. Therefore, we explored age-and sex-related patterns of IgG and IgM antibody titres against S. damnosum s.l. saliva across well characterized endemic communities.

Both IgG and IgM antibody levels were high in children and gradually declined with increasing age. Similar trends have been observed in human antibody responses to the saliva of several mosquito species [53–58], posing the question of whether these patterns indicate decreasing exposure to vector bites with age, or increasing immune tolerance and desensitization with persistent or cumulative saliva exposure [59–62]. Investigation of human IgG responses to sand fly saliva supports the latter proposition, with higher anti-saliva antibody responses observed in new compared to long-term residents of a sand fly-endemic region [63]. Similarly, desensitization to salivary antigens was detected in an area colonized for more than 25 years by Aedes mosquitoes, compared to an area where individuals had been exposed for no longer than 5 years [54]. Following that pattern, median IgG antibody responses also tended to be lower in the current study villages where higher S. damnosum s.l. biting rates were previously reported [30,42]. Antibody responses to Anopheles mosquito bites measured after the summer season of high vector abundance were considerably higher than those before the summer season; notwithstanding they consistently declined with increasing host age [54,55]. Interestingly, in that study the decline appeared to be antigen-dependent, as the trend was not detected in antibody responses to a specific recombinantly-expressed protein as opposed to the whole salivary gland homogenate [64–66]. Several immunogenic proteins were recently detected in S. damnosum s.l. saliva of which most were well-known salivary antigens [23]. Future expression of these in recombinant forms may be instructive and increase assay sensitivity. This is especially interesting for IgM as these antibody responses were shown to be less specific than the IgG responses [23].

The median IgG response was lower in males than females, which can be attributed to the greater rate of decline of the IgG response with increasing age in males. This may reflect sex differences in behaviour such as daily habits, occupation, education, or clothing, that influence physical exposure to blood-seeking blackflies. Male occupants of most ages are responsible for agriculture, farming and fishing in the Pru region and may be less well covered by protective clothing, whereas women are more covered, spend more time at home performing domestic duties and/or engage in long-distance trading activities [67,68]. The current study was limited in not recording the daily activities of the participants, though biting blackflies appeared to be ubiquitous throughout the day within villages. That females show greater levels of non-specific innate and adaptive immune responsiveness than males, particularly post-puberty, suggests that hormonal involvement (reviewed in [69]) may contribute to the sex differences observed in this study.

A non-mutually exclusive alternative driving factor behind the observed decline in antibody levels may also be a lower exposure to blackfly bites with increasing age. Such patterns of exposure were predicted by fitting age- and sex-structured onchocerciasis transmission mathematical models to age- and sex-specific profiles of O. volvulus skin microfilariae in another African savannah setting [5]. However, the model fits also indicated that vector exposure of females should increase (rather than decrease) with age. If this is correct, and if women indeed are more intensely exposed as they age, it further supports the case for desensitization with increased long-term biting exposure. Therefore, one next step to better understand exposure patterns with age and sex is to fit dynamic models of antibody acquisition and decay to the (IgG) data obtained here [70,71]. Future studies of the molecular and cellular mechanisms that underly immune tolerance and progressive desensitization to blackfly saliva would also be most informative. Variations in IgG subclass responsiveness are also possible as shown against Anopheles or Aedes saliva [55,72]. Interestingly, human IgG4 amongst bee-keepers was found to be associated with immunotolerance to bee venom [73,74].

Less clear in this study were the age and sex-related trends in IgM responses. Although they also declined with age supporting the immunological desensitization hypothesis, IgM responses are generally shorter lived than IgG responses, hence, likely to be more indicative of recent exposure. A shorter half-life together with a lack of cumulative increase after repeat exposure may partially explain the large number of low IgM responses observed in the village residents, and the lack of correlations between individual host IgG and IgM responses, particularly as the current study was limited to cross-sectional sampling during the high biting season.

Future studies would benefit from quantifying short- and long-term kinetics in individual anti-saliva Ig responses in the context of seasonal fluctuations in vector abundance and distance to vector breeding habitats, to refine our understanding of the link between biting rates and the Ig responses [16,53,54]. In fact, it would be very informative if antibody data generated using our immunoassays could be used in spatial analyses to better understand patterns of vector–human contact with increasing distance from breeding sites. At present, mapping exercises such as those used for Onchocerciasis Elimination Mapping (OEM) collect information on the distribution of breeding sites and vector presence. Our novel tool could complement OEM to identify high-risk locations where exposure to vector bites would provide additional information to seroprevalence surveys to guide start-MDA decisions and identify informative age/sex groups for sampling [75]. Most of the information about the relationship between vector density and distance from breeding sites pertains to African savannah settings (such as those explored here), with less data available to characterize such a relationship in forest and forest-savannah mosaic settings. Therefore, if our anti-vector saliva assays could be combined with (seroprevalence) parasite exposure assays for a range of epidemiological settings, it would be possible to obtain valuable information to help elimination efforts. However, this necessitates the testing and validation of our assays for other species/cytoforms of the S. damnosum complex.

The predominant vector species in the Bono East region are the savannah members of the complex, S. damnosum sensu stricto/S. sirbanum [27,30]. It remains to be established if the anti-S. damnosum s.l. IgG and IgM antibody responses represent a damnosum complex-specific marker, or if member-specific markers of exposure would be more sensitive. We acknowledge that more data need to be collected to ascertain the validity of our immunoassays to reliably measure exposure to vector bites. Regarding empirical approaches, there are no experimental or observational data yet to precisely quantify the relationship between anti-salivary antibody levels and the number of vector bites. As colonizing blackflies in the laboratory is notoriously difficult, collecting field data to understand the kinetics of the antibody responses over a specified time frame would be helpful, particularly in settings with strong seasonality, in which vector biting ceases or greatly decreases for several months during the year. Regarding theoretical approaches, and as mentioned earlier, the exposure profiles predicted by age- and sex-structured onchocerciasis transmission models are broadly consistent with the proposed desensitization hypothesis. If this proves to be the case, we expect vector saliva-naïve children still to respond well to salivary antigens as an indication of continued exposure to vector bites pre- or post-MDA campaigns in areas with no vector control or without major ecological changes affecting vector density. Our assay could then be used to understand potential secular trends in vector biting rates due to anthropogenic change. Ideally, multiplex assays could be developed to test simultaneously for both exposure to vector bites and to parasite antigens.

Conclusion

Serological anti-saliva assays are useful tools to complement information collected by HLCs by measuring human-vector contact and revealing heterogeneities in exposure at the individual and community level that cannot be unravelled by HLCs alone. By novel application to four onchocerciasis-endemic communities, this study successfully evaluated age- and sex-related demographic patterns in blackfly bite exposure. The analyses uncovered the possibility of age- and sex-specific immunotolerance or desensitization to blackfly saliva, likely resulting from cumulative blackfly exposure with age. Concomitant studies of infection levels in humans and flies, vector abundance, and immune responses to blackfly saliva and parasite antigens would greatly help to better understand transmission risk and intensity, and improve parameterization of transmission models with which to inform optimal interventions and surveillance strategies to achieve and protect onchocerciasis elimination.

Supporting information

S1 File. Supplementary Figures and Tables.

Fig A–Proportion of tested individuals. Fig B–Age distribution of the population shown per cluster and sex. Table C–Number of people sampled and tested for IgG and IgM per age group and sex. Fig D–IgG antibody distribution according to age shown per cluster. Fig E–IgM antibody distribution according to age shown per cluster. Table F–Summary of exponentiated regression coefficient estimates. Fig F–Correlation between IgG and IgM responses.

https://doi.org/10.1371/journal.pntd.0010108.s001

(DOCX)

Acknowledgments

Special thanks go to Dr. Francis Veriegh, Frank Aboagye-Twum and Bright Idun for assisting in the blood sampling and the collection of blackfly salivary glands. We gratefully acknowledge the researchers at the Noguchi Memorial Institute for Medical Research in Accra for allowing us access to their laboratory facilities. We also thank all the participants of the study, who were crucial to the execution of the study.

Citation: Willen L, Milton P, Hamley JID, Walker M, Osei-Atweneboana MY, Volf P, et al. (2022) Demographic patterns of human antibody levels to Simulium damnosum s.l. saliva in onchocerciasis-endemic areas: An indicator of exposure to vector bites. PLoS Negl Trop Dis 16(1): e0010108. https://doi.org/10.1371/journal.pntd.0010108

References

1. World Health Organization. Department of Control of Neglected Tropical Diseases. Ending the neglect to attain the Sustainable Development Goals–A road map for neglected tropical diseases 2021–2030. Geneva World Health Organization. 2021 (Editor: Malecele MN). Available: https://www.who.int/publications/i/item/9789240010352 (accessed 11 December 2021).

2. Hamley JID, Milton P, Walker M, Basáñez MG. Modelling exposure heterogeneity and density dependence in onchocerciasis using a novel individual-based transmission model, EPIONCHO-IBM: Implications for elimination and data needs. PLoS Negl Trop Dis. 2019;13: 1–21. pmid:31805049

3. Walker M, Stolk WA, Dixon MA, Bottomley C, Diawara L, Traoré MO, et al. Modelling the elimination of river blindness using long-term epidemiological and programmatic data from Mali and Senegal. Epidemics. 2017;18: 4–15. pmid:28279455

4. Hamley JID, Walker M, Coffeng LE, Milton P, de Vlas SJ, Stolk WA, et al. Structural uncertainty in onchocerciasis transmission models influences the estimation of elimination thresholds and selection of age groups for seromonitoring. J Infect Dis. 2020;221: S510–S518. pmid:32173745

5. Filipe JAN, Boussinesq M, Renz A, Collins RC, Vivas-Martinez S, Grillet ME, et al. Human infection patterns and heterogeneous exposure in river blindness. Proc Natl Acad Sci USA. 2005;102: 15265–15270. pmid:16217028

6. Otabil KB, Gyasi SF, Awuah E, Obeng-Ofori D, Rodríguez-Pérez MA, Katholi CR, et al. The search for an efficient black fly trap for xenomonitoring of onchocerciasis. J Parasitol Res. 2018;2018: 1–10. pmid:30057801

7. Renz A, Wenk P. Studies on the dynamics of transmission of onchocerciasis in a Sudan-savanna area of North Cameroon IV. The different exposure to Simulium bites and transmission of boys and girls and men and women, and the resulting manifestations of onchocerciasis. Ann Trop Med Parasitol. 1987;81: 215–228. pmid:3662664

8. Talom BAD, Enyong P, Cheke RA, Djouaka R, Hawkes FM. Capture of high numbers of Simulium vectors can be achieved with Host Decoy Traps to support data acquisition in the onchocerciasis elimination endgame. Acta Trop. 2021;221: 106020. pmid:34157291

9. Jacobi CA, Enyong P, Renz A. Individual exposure to Simulium bites and intensity of Onchocerca volvulus infection. Parasit Vectors. 2010;3: 1–6. pmid:20051120

10. Arcà B, Ribeiro JM. Saliva of hematophagous insects: a multifaceted toolkit. Curr Opin Insect Sci. 2018;29: 102–109. pmid:30551815

11. Fontaine A, Pascual A, Orlandi-Pradines E, Diouf I, Remoué F, Pagès F, et al. Relationship between exposure to vector bites and antibody responses to mosquito salivary gland extracts. PLoS One. 2011;6: e29107. pmid:22195000

12. Doucoure S, Mouchet F, Cornelie S, DeHecq JS, Rutee AH, Roca Y, et al. Evaluation of the human IgG antibody response to Aedes albopictus saliva as a new specific biomarker of exposure to vector bites. PLoS Negl Trop Dis. 2012;6: e1487. pmid:22363823

13. Vlkova M, Rohousova I, Hostomska J, Pohankova L, Zidkova L, Drahota J, et al. Kinetics of antibody response in BALB/c and C57BL/6 mice bitten by Phlebotomus papatasi. PLoS Negl Trop Dis. 2012;6: e1719. pmid:22802977

14. Caljon G, Van Den Abbeele J, Sternberg JM, Coosemans M, De Baetselier P, Magez S. Tsetse fly saliva biases the immune response to Th2 and induces anti-vector antibodies that are a useful tool for exposure assessment. Int J Parasitol. 2006;36: 1025–1035. pmid:16777113

15. Schwarz A, Sternberg JM, Johnston V, Medrano-Mercado N, Anderson JM, Hume JCC, et al. Antibody responses of domestic animals to salivary antigens of Triatoma infestans as biomarkers for low-level infestation of triatomines. Int J Parasitol. 2009;39: 1021–1029. pmid:19248784

16. Kostalova T, Lestinova T, Sumova P, Vlkova M, Rohousova I, Berriatua E, et al. Canine antibodies against salivary recombinant proteins of Phlebotomus perniciosus: A longitudinal study in an endemic focus of canine leishmaniasis. PLoS Negl Trop Dis. 2015;9: e0003855. pmid:26111018

17. Quinnell RJ, Soremekun S, Bates PA, Rogers ME, Garcez LM, Courtenay O. Antibody response to sand fly saliva is a marker of transmission intensity but not disease progression in dogs naturally infected with Leishmania infantum. Parasit Vectors. 2018;11: 1–12. pmid:29291748

18. Brosseau L, Drame PM, Besnard P, Toto JC, Foumane V, Le Mire J, et al. Human antibody response to Anopheles Saliva for comparing the efficacy of three malaria vector control methods in Balombo, Angola. PLoS One. 2012;7: e44189. pmid:23028499

19. Drame PM, Poinsignon A, Besnard P, Le Mire J, Dos-Santos MA, Sow CS, et al. Human antibody response to Anopheles gambiae saliva: An immuno-epidemiological biomarker to evaluate the efficacy of insecticide-treated nets in malaria vector control. Am J Trop Med Hyg. 2010;83: 115–121. pmid:20595489

20. Gidwani K, Picado A, Rijal S, Singh SP, Roy L, Volfova V, et al. Serological markers of sand fly exposure to evaluate insecticidal nets against visceral leishmaniasis in India and Nepal: A cluster-randomized trial. PLoS Negl Trop Dis. 2011;5: e1296. pmid:21931871

21. Dinesh DS, Das P, Picado A, Davies C, Speybroeck N, Ostyn B, et al. Long-lasting insecticidal nets fail at household level to reduce abundance of sandfly vector Phlebotomus argentipes in treated houses in Bihar (India). Trop Med Int Health. 2008;13: 953–958. pmid:18482197

22. Schwarz A, Juarez JA, Richards J, Rath B, Machaca VQ, Castro YE, et al. Anti-triatomine saliva immunoassays for the evaluation of impregnated netting trials against Chagas disease transmission. Int J Parasitol. 2011;41: 591–594. pmid:21426907

23. Willen L, Basáñez MG, Dvorak V, Veriegh FBD, Aboagye FT, Idun B, et al. Human immune response against salivary antigens of Simulium damnosum s.l.: A new epidemiological marker for exposure to blackfly bites in onchocerciasis endemic areas. PLoS Negl Trop Dis. 2021;15: e0009512. pmid:34157020

24. Hougard J-M, Alley ES, Yaméogo L, Dadzie KY, Boatin BA. Eliminating onchocerciasis after 14 years of vector control: A proved strategy. J Infect Dis. 2001;184: 497–503. pmid:11471108

25. O’Hanlon SJ, Slater HC, Cheke RA, Boatin BA, Coffeng LE, Pion SDS, et al. Model-based geostatistical mapping of the prevalence of Onchocerca volvulus in West Africa. PLoS Negl Trop Dis. 2016;10: e0004328. pmid:26771545

26. Alley ES, Plaisier AP, Boatin BA, Dadzie KY, Remme J, Zerbo G, et al. The impact of five years of annual ivermectin treatment on skin microfilarial loads in the onchocerciasis focus of Asubende, Ghana. Trans R Soc Trop Med Hyg. 1994;88: 581–584. pmid:7992347

27. Lamberton PHL, Cheke RA, Winskill P, Tirados I, Walker M, Osei-Atweneboana MY, et al. Onchocerciasis transmission in Ghana: Persistence under different control strategies and the role of the simuliid vectors. PLoS Negl Trop Dis. 2015;9: e0003688. pmid:25897492

28. Osei-Atweneboana MY, Eng JKL, Boakye DA, Gyapong JO, Prichard RK. Prevalence and intensity of Onchocerca volvulus infection and efficacy of ivermectin in endemic communities in Ghana: a two-phase epidemiological study. Lancet. 2007;369: 2021–2029. pmid:17574093

29. Post RJ, Cheke RA, Boakye DA, Wilson MD, Osei-Atweneboana MY, Tetteh-Kumah A, et al. Stability and change in the distribution of cytospecies of the Simulium damnosum complex (Diptera: Simuliidae) in southern Ghana from 1971 to 2011. Parasit Vectors. 2013;6: 1–11. pmid:23281838

30. Lamberton PHL, Cheke RA, Walker M, Winskill P, Osei-Atweneboana MY, Tirados I, et al. Onchocerciasis transmission in Ghana: Biting and parous rates of host-seeking sibling species of the Simulium damnosum complex. Parasit Vectors. 2014;7: 511. pmid:25413569

31. Cheke RA, Basáñez M, Perry M, White M, Garms R, Obuobie E, et al. Potential effects of warmer worms and vectors on onchocerciasis transmission in West Africa. Philos Trans R Soc B Biol Sci. 2015;370: 20130559. pmid:25688018

32. Frempong KK, Walker M, Cheke RA, Tetevi EJ, Gyan ET, Owusu EO, et al. Does increasing treatment frequency address suboptimal responses to ivermectin for the control and elimination of river blindness? Clin Infect Dis. 2016;62: 1338–1347. pmid:27001801

33. Doyle SR, Bourguinat C, Nana-Djeunga HC, Kengne-Ouafo JA, Pion SDS, Bopda J, et al. Genome-wide analysis of ivermectin response by Onchocerca volvulus reveals that genetic drift and soft selective sweeps contribute to loss of drug sensitivity. PLoS Negl Trop Dis. 2017;11: e0005816. pmid:28746337

34. Plaisier AP, van Oortmarssen GJ, Habbema JDF, Remme J, Alley ES. ONCHOSIM: a model and computer simulation program for the transmission and control of onchocerciasis. Comput Methods Programs Biomed. 1990;31: 43–56. pmid:2311368

35. Basáñez MG, Churcher TS, Grillet ME. Onchocerca-Simulium interactions and the population and evolutionary biology of Onchocerca volvulus. Adv Parasitol. 2009;68: 263–313. pmid:19289198

36. Basáñez MG, Walker M, Turner HC, Coffeng LE, de Vlas SJ, Stolk WA. River blindness: Mathematical models for control and elimination. Adv Parasitol. 2016;94: 247–341. pmid:27756456

37. Hamley JID, Blok DJ, Walker M, Milton P, Hopkins AD, Hamill LC, et al. What does the COVID-19 pandemic mean for the next decade of onchocerciasis control and elimination? Trans R Soc Trop Med Hyg. 2021;115: 269–280. pmid:33515042

38. Turner HC, Osei-Atweneboana MY, Walker M, Tettevi EJ, Churcher TS, Asiedu O, et al. The cost of annual versus biannual community-directed treatment of onchocerciasis with ivermectin: Ghana as a case study. PLoS Negl Trop Dis. 2013;7: e2452. pmid:24069497

39. Walsh JF, Davies JB, Le Berre R, Garms R. Standardization of criteria for assessing the effect of Simulium control in onchocerciasis control programmes. Trans R Soc Trop Med Hyg. 1978;72: 675–676. pmid:734734

40. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York, 2009.

41. South A. rworldmap: A New R package for Mapping Global Data. R J. 2011;3: 35–43.

42. Veriegh FBD. Impacts of semi-annual ivermectin treatment and the spatio-temporal distribution and abundance of Simulium damnosum sensu lato vectors on the transmission of Onchocerca volvulus. PhD Thesis. Kwame Nkrumah University of Science and Technology (KNUST), Ghana. 2020.

43. Biritwum NK, de Souza DK, Asiedu O, Marfo B, Amazigo UV, Gyapong JO. Onchocerciasis control in Ghana (1974–2016). Parasit Vectors. 2021;14: 1–9. pmid:33388087

44. R Core Team (2015). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available: https://www.r-project.org/ (accessed 11 December 2021).

45. Cohen J. Statistical Power Analysis for the Behavioral Sciences. Second edition. Lawrence Erlbaum Associates; 1988.

46. Ryan TP. Sample Size Determination and Power. John Wiley & Sons, Inc.; 2013. https://doi.org/10.1002/9781118439241

47. Andrade BB, Teixeira CR. Biomarkers for exposure to sand flies bites as tools to aid control of leishmaniasis. Front Immunol. 2012;3: 121. pmid:22661974

48. Lestinova T, Rohousova I, Sima M, de Oliveira CI, Volf P. Insights into the sand fly saliva: Blood-feeding and immune interactions between sand flies, hosts, and Leishmania. PLoS Negl Trop Dis. 2017;11: e0005600. pmid:28704370

49. Doyle SR, Armoo S, Renz A, Taylor MJ, Osei-Atweneboana MY, Grant WN. Discrimination between Onchocerca volvulus and O. ochengi filarial larvae in Simulium damnosum (s.l.) and their distribution throughout central Ghana using a versatile high-resolution speciation assay. Parasit Vectors. 2016;9: 536. pmid:27724959

50. Yarzábal L, Basáñez MG, Ramírez-Pérez J, Ramírez A, Botto C, et al. Experimental and natural infection of Simulium sanchezi by Mansonella ozzardi in the Middle Orinoco region of Venezuela. Trans R Soc Trop Med Hyg. 1985;79: 29–33. pmid:3887680

51. Cupp EW, Maré CJ, Cupp MS, Ramberg FB. Biological transmission of vesicular stomatitis virus (New Jersey) by Simulium vittatum (Diptera: Simuliidae). J Med Entomol. 1992;29: 137–140. pmid:1322991

52. Turner HC, Walker M, Lustigman S, Taylor DW, Basáñez MG. Human onchocerciasis: Modelling the potential long-term consequences of a vaccination programme. PLoS Negl Trop Dis. 2015;9: e0003938. pmid:26186715

53. Rizzo C, Ronca R, Fiorentino G, Verra F, Mangano V, Poinsignon A, et al. Humoral response to the Anopheles gambiae salivary protein gSG6: A serological indicator of exposure to Afrotropical malaria vectors. PLoS One. 2011;6: e17980. pmid:21437289

54. Buezo Montero S, Gabrieli P, Montarsi F, Borean A, Capelli S, De Silvestro G, et al. IgG antibody responses to the Aedes albopictus 34k2 salivary protein as novel candidate marker of human exposure to the tiger mosquito. Front Cell Infect Microbiol. 2020;10: 377. pmid:32850479

55. Rizzo C, Ronca R, Lombardo F, Mangano V, Sirima SB, Nèbiè I, et al. IgG1 and IgG4 antibody responses to the Anopheles gambiae salivary protein gSG6 in the sympatric ethnic groups Mossi and Fulani in a malaria hyperendemic area of Burkina Faso. PLoS One. 2014;9: e96130. pmid:24760038

56. Londono-Renteria B, Drame PM, Montiel J, Vasquez AM, Tobón-Castaño A, Taylor M, et al. Identification and pilot evaluation of salivary peptides from Anopheles albimanus as biomarkers for bite exposure and malaria infection in Colombia. Int J Mol Sci. 2020;21: 691. pmid:31973044

57. Montiel J, Carbal LF, Tobón-Castaño A, Vásquez GM, Fisher ML, Londono-Rentería B. IgG antibody response against Anopheles salivary gland proteins in asymptomatic Plasmodium infections in Nariño, Colombia. Malar J. 2020;19: 42. pmid:31973737

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

59. Peng Z, Simons FER. Mosquito allergy: Immune mechanisms and recombinant salivary allergens. Int Arch Allergy Immunol. 2004;133: 198–209. pmid:14963366

60. Feingold BF, Benjamini E, Michaeli D. The allergic responses to insect bites. Annu Rev Entomol. 1968;13: 137–158.

61. Peng Z, Simons FER. A prospective study of naturally acquired sensitization and subsequent desensitization to mosquito bites and concurrent antibody responses. J Allergy Clin Immunol. 1998;101: 284–286. pmid:9500765

62. Mellanby K. Man’s reaction to mosquito bites. Nature. 1946;158: 554. pmid:21001941

63. Mondragon-Shem K, Al-Salem WS, Kelly-Hope L, Abdeladhim M, Al-Zahrani MH, Valenzuela JG, et al. Severity of Old World cutaneous leishmaniasis is influenced by previous exposure to sandfly bites in Saudi Arabia. PLoS Negl Trop Dis. 2015;9: e0003449. pmid:25646796

64. Rizzo C, Lombardo F, Ronca R, Mangano V, Sirima S, Nèbiè I, et al. Differential antibody response to the Anopheles gambiae gSG6 and cE5 salivary proteins in individuals naturally exposed to bites of malaria vectors. Parasit Vectors. 2014;7: 549. pmid:25428638

65. Ribeiro JMC, Arcà B, Lombardo F, Calvo E, Phan VM, Chandra PK, et al. An annotated catalogue of salivary gland transcripts in the adult female mosquito, Aedes aegypti. BMC Genom. 2007;8: 6. pmid:17204158

66. Londono-Renteria BL, Shakeri H, Rozo-Lopez P, Conway MJ, Duggan N, Jaberi-Douraki M, et al. Serosurvey of human antibodies recognizing Aedes aegypti D7 salivary proteins in Colombia. Front Public Health. 2018;6: 11. pmid:29479525

67. World Health Organization & Onchocerciasis Control Programme in West Africa. Socio-demographic study in the Pru Basin. Ouagadougou 1l–15 March 2002. Available: https://apps.who.int/iris/handle/10665/311605 (accessed 11 December 2021).

68. Kim R, van Moorsel T, Kapstein EB. The National and Regional Socio-Economic Impact of Newmont Ghana’s Ahafo Mine. Steward Redqueen Rep. 2013. Available: https://s24.q4cdn.com/382246808/files/doc_downloads/operations_projects/africa/quick_links/Socio-economic-impact-of-Ahafo-Operations.pdf (accessed 11 December 2021).

69. Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016;16: 626–638. pmid:27546235

70. White MT, Griffin JT, Akpogheneta O, Conway DJ, Koram KA, Riley EM, et al. Dynamics of the antibody response to Plasmodium falciparum infection in African children. J Infect Dis. 2014;210: 1115–1122. pmid:24719471

71. Yman V, White MT, Rono J, Arcà B, Osier FH, Troye-Blomberg M, et al. Antibody acquisition models: A new tool for serological surveillance of malaria transmission intensity. Sci Rep. 2016;6: 1–10. pmid:28442746

72. Remoue F, Alix E, Cornelie S, Sokhna C, Cisse B, Doucoure S, et al. IgE and IgG4 antibody responses to Aedes saliva in African children. Acta Trop. 2007;104: 108–115. pmid:17825239

73. Varga EM, Kausar F, Aberer W, Zach M, Eber E, Durham SR, et al. Tolerant beekeepers display venom-specific functional IgG4 antibodies in the absence of specific IgE. J Allergy Clin Immunol. 2013;131: 1419–1421. pmid:23063581

74. García-Robaina JC, De La Torre-Morín F, Vazquez-Moncholi C, Fierro J, Bonnet-Moreno C. The natural history of Apis-specific IgG and IgG4 in beekeepers. Clin Exp Allergy. 1997;27: 418–423. pmid:9146935

75. World Health Organization. Department of Control of Neglected Tropical Diseases. Report of the Third Meeting of the WHO Onchocerciasis Technical Advisory Subgroup Geneva, Switzerland, 26–28 February 2019 (Editor: Cantey P). Available: https://www.who.int/publications/i/item/9789240006638 (accessed 11 December 2021).