A comparative assessment of adult mosquito trapping methods to estimate spatial patterns of abundance and community composition in southern Africa
Parasites & Vectors volume 12, Article number: 462 (2019)
Assessing adult mosquito populations is an important component of disease surveillance programs and ecosystem health assessments. Inference from adult trapping datasets involves comparing populations across space and time, but comparisons based on different trapping methods may be biased if traps have different efficiencies or sample different subsets of the mosquito community.
We compared four widely-used trapping methods for adult mosquito data collection in Kruger National Park (KNP), South Africa: Centers for Disease Control miniature light trap (CDC), Biogents Sentinel trap (BG), Biogents gravid Aedes trap (GAT) and a net trap. We quantified how trap choice and sampling effort influence inferences on the regional distribution of mosquito abundance, richness and community composition.
The CDC and net traps together collected 96% (47% and 49% individually) of the 955 female mosquitoes sampled and 100% (85% and 78% individually) of the 40 species or species complexes identified. The CDC and net trap also identified similar regional patterns of community composition. However, inference on the regional patterns of abundance differed between these traps because mosquito abundance in the net trap was influenced by variation in weather conditions. The BG and GAT traps collected significantly fewer mosquitoes, limiting regional comparisons of abundance and community composition.
This study represents the first systematic assessment of trapping methods in natural savanna ecosystems in southern Africa. We recommend the CDC trap or the net trap for future monitoring and surveillance programs.
Adult mosquito sampling is a key component of mosquito surveillance [1,2,3,4], but trapping success may vary across studies due to differences in trapping methods. Different traps vary in their ability to catch certain species and life stages [5,6,7,8]. For example, the dominant species attracted with light-baited traps may differ from those attracted to traps baited with carbon dioxide or live hosts . Additionally, sampling conditions such as the number of nights over which sampling occurs and weather conditions may also influence trapping success , with some species and traps potentially more affected than others. This variation in trapping outcome may not limit inference on the presence or absence of common species (e.g. the information used in global risk maps [11, 12]). However, it does limit inference based on comparing patterns of diversity or abundance across space or time . Given that these comparisons are required to evaluate the ecological or anthropogenic drivers of mosquito populations and disease risk, choice of trapping methods presents challenges and opportunities for optimizing sampling efforts.
Previous studies have evaluated trapping methods in Europe, North America and South America [5,6,7,8, 13] while studies in southern Africa remain relatively limited. The mosquito fauna (Diptera: Culicidae) of southern Africa consists of over 216 species, many of which are endemic to the region . Additional studies evaluating trapping methods are needed to evaluate if traps designed for other locations and species perform equally well in the species-rich communities in southern Africa. For example, a recent comparison of four traps in Germany found that the Biogents Sentinel trap (BG trap) collected the highest number of individuals in urban and snowmelt forest environments where Culex species predominate. In contrast, the Centers for Disease Control miniature light trap (CDC trap) collected the most mosquitoes in floodplain environments where Aedes vexans predominate, suggesting that the preferred trapping device may vary by habitat and community composition . Qualitative evaluations in southern Africa suggest different traps likely collect different subsets of the mosquito community . However, the two studies that systematically evaluated trapping methods in southern Africa focused only on house-based trapping in residential areas [16, 17].
An evaluation of trapping methods is also needed to influence the design of future entomological and pathogen surveillance efforts . Such efforts could provide baseline information for public health interventions by identifying hotspots with a high abundance of vector species [19, 20]. There were more than 29,000 confirmed cases of malaria in southern Africa in 2017 (Botswana, Namibia, South Africa, Lesoto, Swaziland) ; key malaria vectors include members of the Anopheles gambiae complex (An. arabiensis, An. gambiae, An. merus) [22,23,24] and the An. funestus complex . Mosquito-borne livestock and wildlife infections are also a concern, as outbreaks of West Nile virus, Rift Valley fever, Sindbis and Wesselsbron occur [26, 27]. Key vectors for these viral infections include, Aedes caballus, Ae. circumluteolus, Ae. dentatus, Ae. juppi, Ae. mcintoshi, Ae. ochraceus, Culex pipiens, Cx. quinquefasciatus, Cx. univittatus and Cx. theileri (based on suspected or known prime vectors in Africa, reviewed elsewhere ). Previous work has characterized the distribution , ecological drivers  and consequences for malaria risk  of key anopheline species. Understanding the distribution and drivers of non-Anophelinae mosquito species in southern Africa could facilitate informed management of a broader range of mosquito-borne disease or the development of vector control programmes that target multiple infections [15, 31].
In this study, we assessed four commonly-used adult mosquito trapping methods for their use in natural savanna ecosystems. Because surveillance programmes may have multiple aims (e.g. entomological surveys, pathogen surveillance or invasive/nuisance mosquito species control) we compared traps designed for a range of purposes with a range or attraction methods. We assessed the Centers for Disease Control miniature light trap (CDC) and a net trap due to their historic success for entomological surveys in a range of systems in southern Africa [9, 15, 16]. These traps use incandescent light + CO2 and CO2, respectively, to attract host-seeking vectors. We assessed the Biogents Sentinel-2 trap (BG) because of its worldwide success for both general entomological surveys  and targeted vector surveillance (e.g. Ae. aegypti, Ae. albopictus, Cx. quinquefasciatus, Cx. pipiens [10, 13]). It uses visual cues, CO2 and a lure to attract mosquitoes searching for hosts or a place to rest. We additionally compare the Biogents gravid Aedes trap (GAT) that is designed to collect gravid female mosquitoes, particularly container-breeding Aedes species, Ae. aegypti and Ae. albopictus [32, 33]. Gravid mosquitoes are more likely to be carrying pathogens, making the GAT trap potentially useful for pathogen surveillance or for joint pathogen and entomological surveillance efforts when used in combination with other traps.
We assessed these traps following two objectives. We compared estimates of abundance and community composition among all four traps (CDC, BG, GAT, net). We also evaluated traps designed for entomological surveillance (CDC, BG, net) by quantifying the importance of trapping method for inferring the regional distribution of mosquito abundance, mosquito community composition and the distribution of key disease vectors.
We sampled 16 sites within four regions of Kruger National Park (KNP), South Africa (Fig. 1; between 22°31ʹ and 25°31ʹS, 30°45ʹ and 32°00ʹE). We focus on KNP because it is a hotspot of mosquito diversity , it is in a region of southern Africa with regular outbreaks of mosquito-borne disease (e.g. [34, 35]), and it includes a sentinel site for mosquitoes and pathogen surveillance ). Sampling sites within KNP were chosen to cover the geographical extent of the park and to capture the variability in rainfall, geology and vegetation types within KNP . This region experiences summer rainfall between November and April, and we selected our sampling to occur from March to April 2017, when mosquito populations are generally high. We sampled four out of the 22 management regions within the park: Malelane, Satara, Shingwedzi and Punda Maria. We additionally sampled in Skukuza during a pilot organizational and training week where traps were not directly compared. Within each region, we sampled at four sites, which were selected based on multiple criteria (Fig. 1). The primary selection criterion was to sample from water bodies that represented diverse types of wetlands (temporary ponds, permanent ponds, rivers). Additional criteria stipulated that the water bodies were at least 2 km away from one another to avoid sampling mosquitoes from adjacent water bodies and within 25 km from one another for sampling logistics. These distances are justified based on mean mosquito dispersal rates, which range from 35 m to 1.4 km depending on the species and habitat [38,39,40,41]. Although dispersal farther than 2 km is possible, it is uncommon  and we do not expect it to influence our results because only one pair of sites in Shingwedzi was located within 1 km of each other due limited surface water availability in the area.
Trapping and identification
We trapped at each site within a region for four consecutive nights and moved sequentially between regions each week in a random order (Skukuza: 20–23 March 2017; Malelane: 27–30 March 2017; Satara: 3–6 April 2017; Shingwedzi: 17–20 April 2017; Punda Maria: 24–27 April 2017). The four traps we used were the Centers for Disease Control miniature light trap with an incandescent light (Bioquip Inc, Rancho Dominguez, USA), a Biogents Sentinel 2 trap (Biogents AG, Regensburg, Germany), the Biogents gravid Aedes trap (Biogents AG) and the net trap. Although the net trap is not commercially available, it is easily and inexpensively made from netting and poles (see [15, 31]). For consistency in sampling, we set up the traps at a similar distance away from the water body (15–25 m) and approximately 30–50 m away from each other. Our aim was to estimate variation due to trap position, so we rotated the position of all traps after two nights of trapping. This study design allowed us to quantify how position within a site influences catch rates and compare if traps vary in their sensitivity to position and weather conditions (temperature, relative humidity and wind speed).
Our trap-use protocol was based on preliminary trapping work and expertise from mosquito surveillance in southern Africa [15, 31]. Specifically, we equipped the CDC trap, BG trap and net trap with a closed container containing 200–400 g of dry ice for a CO2 bait. To ensure the dry ice would last overnight, we reduced sublimation by wrapping the dry ice in multiple layers of newspaper and ensuring the container was closed except for 2 small holes. We placed the dry ice containers inside the BG traps, at the center of the net trap, and hanging with the CDC trap. Additionally, we equipped the BG traps with BG-Lures (Biogents AG). The GAT trap was baited with water following prior field trials suggesting water and hay-infusion bated GAT traps were equally successful . We chose to use water to ensure consistently-baited traps across the study because although trap entry is not driven by chemical cues in the infusion [44, 45], oviposition behavior is highly variable across infusion types, infusion concentrations and fermentation period . To set the trap, we placed a single piece of long-lasting insecticidal net inside the center of the trap, which was coated in 4.8% alphacypermethrin (Biogents AG). In addition to the spatial locations defined above, we hung the CDC light trap on a branch so it was approximately 1 m from the ground, and we placed the BG and GAT traps on the ground. We ensured the traps were placed close to the vegetation, but not directly under it, which has been shown to improve trapping success . If available, we placed the net trap in slight shade and pulled the netting down to leave a gap of approximately 10 cm between the bottom of the net and the ground.
All traps were monitored twice a day at before-dusk (16:00–17:00 h) and dawn (6:00–7:00 h). The GAT trap was monitored during each visit but was collecting mosquitoes continuously because it does not require batteries. All other traps were set up before dusk and emptied again at dawn. The timing of this sample collection was used because the net trap requires clearing at dawn (mosquitoes may leave the trap after sunrise) and the CDC and BG traps require batteries, which we recharged during the day. This timing focuses our collection on species active between 16:00 h and the following 6:00 h, potentially missing species active just after dawn, the implications of which are discussed below. Due to rainfall, we excluded one trapping night at all sites in Satara and Punda Maria from all analyses.
Our sampling design (16 sites; 3–4 sampling nights each; 4 traps) resulted in 224 total trap-nights. Comparing traps at the same site resulted in 40 BG-CDC comparisons, 40 BG-GAT comparisons, 40 BG-net comparisons, 50 CDC-GAT comparisons and 50 CDC-net comparisons. The differences in sample size are due to both organization errors (11/224, or 5% of trapping nights) and losses due to wildlife interference (13/224, or 6% of trapping nights). Specifically, three BG traps were destroyed over the study by hyena, likely due to the BG-lure, which is specially designed to mimic human sent. Damage to traps from animals is an inherent feature of sampling in wildlife areas, although we took care to minimize the effects of our trapping. Because the consequences of these challenges have not previously been quantified and should be considered in future sample size calculations, we have provided detailed notes on trapping rates and trap-wildlife interactions in Additional file 1: Table S1. Our analyses are conservative in that they only compare data from the 220 successfully-collected, non-damaged trapping nights.
Directly after emptying the traps, we stored the mosquitoes on dry ice in the field and at − 20 °C until identification. For identification of Anophelinae mosquitoes, we used identification literature from Gilles & Coutzee . All anophelines were identified to species except for members of An. funestus complex and the An. gambiae complex, referred to herein as An. gambiae (s.l.), which require molecular methods for identification [38, 47]. For identification of the Culicinae, we used the key by Jupp . Species were identified independently in duplicate by coauthors. Because species identification in the Aedes vexans complex (Ae. vexans, Ae. hirsutus, Ae. fowleri, Ae. durbanensis, Ae. natronius) and the Aedes dentatus complex (Ae. dentatus, Ae. subdentatus Ae. pachyurus, Ae. bevisi, Ae. cumminsii) were inconsistent at the species level, we aggregated them to the species complex level (Additional file 1: Table S2).
Temperature and rainfall within KNP follow a north-south gradient, with the highest values in the southwest . We monitored temperature, relative humidity and wind speed using a Kestrel 3000 handheld weather meter (NK Inc., Boothwyn, PA, USA). We calculated the median temperature, relative humidity, and wind speed across sites within a region for each morning at the time of collection (Additional file 1: Table S3). All environmental variables were standardized for analysis, by centering and dividing by one standard deviation.
Assessing how the number of traps and trapping nights influences mosquito richness
We assessed how observed species richness (the total number of unique species) saturates with sampling effort. We aggregated the data successively over 1, 2, 3 or 4 nights and calculated the cumulative proportion of species identified with an increasing proportion of nights. We also evaluated each trap’s ability to estimate richness by comparing richness estimated in pairs of non-damaged trap types at each site aggregated across all sampling nights using a Wilcoxon signed-rank test and a Bonferroni correction for multiple comparisons. This analysis accounts for trap losses by excluding comparisons from damaged traps.
Assessing whether trap type influences estimates of mosquito abundance and inference on the regional patterns of abundance
To test for differences in abundance, we compared the number of mosquitoes collected per night between pairs of non-damaged traps. First, we quantified the relationship between trap type and abundance with a Wilcoxon signed-rank test and a Bonferroni correction for comparisons among the 4 traps. Then, we quantified the influence of trap type for inferences on the regional patterns of abundance using linear regressions with Poisson errors and a log-link function. The regression analyses assessed the relationship between a trap’s nightly mosquito counts with region of the park and weather conditions (wind speed, temperature, relative humidity). It also included random effect for site (Additional file 2: Table S4). We were additionally interested in the variance due to trap position within a site, estimated with a random effect for trap position. Although we were unable to estimate this parameter in the regression due to the number of damaged traps, our comparisons remain valid because the variation in counts due to trap position was smaller than the variation due to trap type or region. We fitted the regression model separately to data from the CDC trap, the net trap, the BG trap and data from aggregating abundance across all traps at a site. Because no individuals were collected in the GAT trap on most nights, we did not fit the model to data from the GAT trap alone. For each dataset, we conducted model selection using backward selection based on Akaike information criteria with a correction for small sample sizes (AICc) and select the model with the lowest AICc value. We fitted all regression models in R  using the lme4 package .
Assessing whether trap type influences estimates of mosquito community composition and inference on the regional patterns of mosquito communities
To assess if different traps provide different estimates of community composition, we first evaluated if certain species were particularly attracted to one trap over the other. We assessed species-specific trap bias by calculating the difference in the number of individuals for each species sampled between each pair of traps collected at a site over all nights of trapping. Because traps were paired at each site, we tested for differences between the traps using a Wilcoxon signed-rank test and a Bonferroni correction for multiple comparisons. We only compared the species-specific trap bias of common species, defined as being observed in the dataset more than three times. Because 23 common species were compared (k = 23), significant differences between traps occur when P-values are less than P = 0.05/23. We assessed trap bias for rare species visually.
To test for differences in community composition among traps, we used a non-parametric analysis of similarities analysis (ANOSIM), visualized potential differences with non-metric multidimensional scaling (NMDS) and quantified the influence of trap type for inferences on the regional patterns of community composition with hierarchical clustering. For all analyses, we calculated Bray-Curtis dissimilarity matrices based on the trap-specific (BG, CDC, net) abundances of all taxa within a site aggregated across sampling nights. The ANOSIM analysis is a non-parametric test for differences in mosquito communities among traps that compares the ranks of Bray-Curtis dissimilarity measures from samples collected from the same vs different traps . To visualize this, we created an ordination of traps and sites in mosquito community space for each region of the park. Before all ordinations, we applied a Wisconsin transformation followed by a square root transformation to the species matrices, which standardizes by species maxima and reduces the influence of highly abundant taxa, respectively . The ordinations converged on a stable two-dimensional solution, based on stress values. We conducted all community analyses in R using the vegan package .
Describing regional patterns of disease vectors
We additionally describe how known prime vectors for West Nile virus (Cx. pipiens, Cx. quinquefasciatus, Cx. theileri, Cx. univittatus), Rift Valley fever (Ae. dentatus, Ae. mcintoshi, Ae. ochraceus), Sindbis (Cx. univittatus), and Wesselsbron (none found) are distributed across regions [22, 27]. Additional known prime vectors for these infections were not found in the study (Ae. caballus, Ae. circumluteolus, Ae. juppi). Chikungunya and dengue fever outbreaks are less common in South Africa , but we describe the distribution of their vector, Ae. aegypti , because additional known prime vectors in Africa were not found (Ae. africanus, Ae. albopictus, Ae. cordellieri, Ae. furcifer, Ae. luteocephalus, Ae. neoafricanus, Ae. taylori). We focus on the known prime vectors for a conservative description, but additional mosquito species are considered suspected vectors (reviewed in ). We used the numbers of Cx. pipiens complex to approximate the numbers of Cx. pipiens and Cx. quinquefasciatus, Cx. univittatus complex for Cx. univittatus, and Ae. dentatus complex for Ae. dentatus. Although this assumption may hide epidemiologically important variation, it is a valid approximation for these infections because multiple members of the complex are known or suspected vectors. However, we did not plot the distribution of malaria vectors because the most abundant species of An. gambiae (s.l.) in KNP, An. quadriannulatus, is not a malaria vector .
We collected 955 female mosquitoes, 946 (99.1%) of which were identified to the level of species or species complex. The most common species included members of the Cx. univittatus complex, Ae. vexans complex and Cx. pipiens complex. We also collected over 50 An. gambiae (s.l.), An. pretoriensis and Cx. theileri females. Mosquito data are provided in Additional file 3: Table S5.
Mosquito communities can be characterized with the net and CDC traps after multiple trapping nights
Based on all traps together, mosquito community richness was sensitive to the number of sampling nights (Fig. 2; Additional file 4: Figure S1). Taken across all sites in Malelane, 89% (24/27) of the total number of unique species identified after four trapping nights were collected by night 2 and 96% (26/27) were collected by night 3. In Shingwedzi, 79% (19/24) of the species were collected by night 2 and 96% (23/24) by night 3. In Punda Maria, 53% (8/15) of the species were collected by night 2 and 87% (13/15) by night 3. These percentages overestimate the percent of richness captured because species accumulation curves suggest more than four nights of sampling are required to estimate total species richness (Fig. 2). The CDC and net trap together sampled all of the mosquito species captured (range across regions, 93–100%), while the BG and GAT trap captured far fewer species (Fig. 2). Estimates of richness in the BG trap were lower than in the net or CDC trap (Wilcoxon signed-rank test: BG vs net, n = 40, V = 136, P < 0.001; BG vs CDC, n = 40, V = 132.5, P < 0.001) but did not differ between the net or CDC trap (n = 50, V = 51, P = 0.940).
Daily mosquito abundance estimates were comparable between the net and CDC traps but not the BG and GAT trap
Most individuals were collected in the CDC or the net trap, while the BG and GAT trap captured far fewer individuals (Fig. 3a). Together, the CDC and net trap sampled 96% of the individuals collected (range across regions, 94–99%; Additional file 1: Table S1). The net trap collected a mean of 8.6 females per night, the CDC trap collected a mean of 7.4, the BG trap collected a mean of 0.7, and the GAT trap collected a mean on 0.1 females (median: 5.5, 4, 0, 0). Paired by site, estimates of mosquito abundance based on the BG trap were lower than estimates from the net and CDC trap (BG vs CDC, n = 40, V = 22, P < 0.001; BG vs net, n = 40, V = 12, P < 0.001; Fig. 3b), while estimates based on the GAT trap were the lowest (BG vs GAT n = 40, V = 153, P < 0.001). However, estimates from the net and CDC trap did not significantly differ (CDC vs net, n = 50, V = 602.5, P = 0.12). When data from traps were aggregated across sites and nights within a region, trap choice continues to have a large influence on mosquito abundance (Additional file 1: Table S1).
Regression analyses showed that trap choice also influences comparisons of abundance between regions. All traps identified similar trends, with the highest numbers captured in Malelane and the lowest numbers captured in Punda Maria (Fig. 3c, d). This spatial pattern was significantly different for models fit to the CDC data (Additional file 2: Table S4). All traps were influenced by weather conditions, particularly temperature and wind speed (Additional file 2: Table S4). We collected higher numbers in warm, low wind conditions (Fig. 3d). The net trap also collected higher numbers in low relative humidity conditions. Results based on counts aggregated from multiple traps were similar to results based on counts from either the CDC or net trap alone (Fig. 3c; Additional file 2: Table S4).
Mosquito community composition was consistent between the net and CDC trap but not the BG trap
Community composition was similar between the net and CDC trap, but not for the BG trap (ANOSIM, overall: R = 0.126, P = 0.04; pairwise: net vs CDC, P = 0.894; net vs BG, P = 0.023; CDC vs BG, P = 0.009). NMDS ordinations of traps in species space reflect this relationship although some variation across regions of the park does exist (Additional file 4: Figure S2, Table S6). However, we did not find evidence for species-specific bias between any of the traps, suggesting that differences in community composition in the BG trap are driven by the relatively lower abundance collected in the trap (Fig. 4). For common species collected in the net vs CDC trap, the number of individuals collected was not significantly different between the traps (Fig. 4a; no hypothesis tests for individual species were significant). Rare species also did not show trap biases and include Ae. aerarius, Ae. metallicus, Ae. unidentatus, An. maculipalpis, An. ziemanni, Cx. antennatus, Cx. bitaenorhynchus, Cx. nebulosus, Lutzia tigripes, Mansonia africana and Uranotaenia balfouri (Additional file 4: Figure S3). For comparisons with the BG trap, the net and CDC traps both collected higher numbers of individuals, but there were no species or genus shifts driving this pattern (Fig. 4b, c; Additional file 4: Figure S4).
Although trap choice influences estimates of community composition (e.g. community richness, ordinations), hierarchical cluster analysis suggests that may be less important for comparisons between regions (Fig. 5; Additional file 4: Figure S5). Regional patterns in mosquito communities were consistent across trap types, with samples from Malelane and Satara being more similar than samples from Shingwedzi. Mosquito communities estimated from the CDC and net trap were clustered by region (Additional file 4: Figure S5), indicating that communities within a region were more closely related to each other than communities between regions regardless of the net or CDC trap. In contrast, samples from the BG trap were clustered separately from the samples from the net and CDC trap (Fig. 5). We therefore describe the regional patterns of disease vectors based on data from all traps together (Fig. 6). Rift Valley fever vectors were most common in Satara while West Nile virus and Sindbis vectors were more common in Shingwedzi and Punda Maria. We did not find known vectors for Wesselsbron (Ae. caballus, Ae. circumluteolus) within the park.
Mosquito trapping is used for disease surveillance, biodiversity surveys and nuisance-reduction. In light of these multiple, non-exclusive aims, this study compared traps based on four potential goals: collecting large numbers of mosquitoes; estimating mosquito community composition, including invasive, vector or rare species; characterizing the spatial patterns of abundance; and characterizing the spatial patterns of community composition. We expected trade-offs among these aims, for example with traps specializing on one vector to be less successful in estimating community richness (and vice versa, e.g. ). In contrast to this expectation, our results indicated that the net and CDC trap consistently performed best across multiple outcomes.
The CDC and net trap collected higher abundances and more unique species compared to the BG and GAT traps. These differences are based on comparisons in 40 or more trapping nights from 16 sites where both traps were deployed (Additional file 1: Table S1). This result is different from trap comparison studies in Europe , the USA (BG , GAT ) and South America (BG , GAT ), where both the BG and GAT traps have been shown to perform well. One reason for their relatively low success within KNP could be the diversity and types of species present within the park. The GAT trap has been designed to capture container-breeding species, such as Ae. aegypti , and the BG trap performs well in sampling Ae. aegypti, Ae. albopictus and Cx. pipiens [7, 58]. Although Ae. aegypti was present in our study and previous studies , they were relatively rare (Additional file 1: Table S2). The BG trap’s sampling efficiency for the two most-common species complexes in our dataset, members of the Ae. vexans complex and Cx. univitatus complex, is either low  for the Ae. vexans complex or previously unknown for the Cx. univitatus complex. Future work might find improved catch rates for the BG trap by modifying the type of lure used or extending the sampling period to 9:30 h to catch Aedes species that bite at dusk and dawn (peak Ae. aegypti activities between 15:30–19:30 h and 5:30–9:30 h; ). An alternative reason for the low success of the BG trap could be that it is negatively influenced by the presence of alternative refugia and oviposition sites . Habitat heterogeneity is likely to be higher within KNP compared to other, primarily urban or suburban environments where mosquito traps have been evaluated ([13, 55, 56] but see ). Additional comparisons in urban environments in southern Africa are needed to distinguish these two hypotheses.
By providing a detailed comparison between the net and CDC traps, our results suggest that these traps provide similar estimates of community richness and community diversity. A key aim of this work was to evaluate how choice of trap influences spatial estimates of mosquito abundance and mosquito community composition. As a result, our sites were selected to sample a diversity of breeding habitats, some of which were known to have high numbers of mosquitoes and others which we anticipated to have lower numbers. This approach is a better approximation to how the traps would be deployed for surveillance compared to targeting only sites with high numbers of mosquitoes. Despite this, the CDC trap collected 7.4 females per night while the net trap collected 8.6 females per night, similar to previous collection efforts in KNP that report aggregated data, with a median of 31 (range 17–116) females across three net traps and a median of 19 (range 19–33) in three CDC traps . Based on this data, the consistent patterns of community composition and largely consistent patterns of abundance between traps suggest that comparisons between studies using these two methods are possible.
This work sets the stage for follow-up studies comparing mosquito abundance and community composition across space or time. We recommend that such a follow-up study should use the CDC or net trap across multiple nights or sites within a region. More traps and sampling will be needed to characterize spatial patterns of abundance because we observed significant variation among sites and relatively low abundances per trap. The choice between the CDC and net trap should also consider other features of the traps, such as ease of use and specimen quality. For morphological identification, the net trap has the advantage that mosquito specimens can be collected with minimal damage, which makes them easier to identify . It also has no motorized or battery-powered parts, making it difficult to break and straightforward to mend. Therefore, studies requiring precise species identification such as biodiversity assessments or studies conducted in remote locations may prefer the net trap. However, for sampling large numbers of sites or sites in remote locations, the CDC trap has an advantage because the timing of when traps need to be cleared is more flexible compared to the net trap, which has to be cleared at sunrise. These practical considerations may mean that the CDC trap is better suited for large, comparative studies.
After assessing four different mosquito trapping methods in a natural savanna ecosystem, we recommend the net trap, the CDC trap or their combined use for outdoor mosquito surveillance in southern Africa. These traps performed well based on four evaluation criteria: collecting large numbers of mosquitoes; estimating mosquito community composition, including vector or rare species; characterizing the spatial patterns of abundance; and characterizing the spatial patterns of community composition. This suggests they are appropriate for both biodiversity surveys and vector surveillance. As such, this study provides a valuable proof-of-principle for characterizing the spatial patterns of non-vectors as well as vectors for multiple diseases.
Akaike information criteria with a correction for small sample sizes
analysis of variance
Biogents sentinel trap
Centers for Disease Control miniature light trap
Biogents gravid Aedes trap
Kruger National Park
Meyer Steiger DB, Ritchie SA, Laurance SGW. Mosquito communities and disease risk influenced by land use change and seasonality in the Australian tropics. Parasites Vectors. 2016;9:387.
Janko MM, Irish SR, Reich BJ, Peterson M, Doctor SM, Mwandagalirwa MK, et al. The links between agriculture, Anopheles mosquitoes, and malaria risk in children younger than 5 years in the Democratic Republic of the Congo: a population-based, cross-sectional, spatial study. Lancet Planet Health. 2018;2:e74–82.
Diuk-Wasser MA, Brown HE, Andreadis TG, Fish D. Modeling the spatial distribution of mosquito vectors for West Nile virus in Connecticut, USA. Vector Borne Zoonotic Dis. 2006;6:283–95.
Galardo AKR, Zimmerman RH, Lounibos LP, Young LJ, Galardo CD, Arruda M, et al. Seasonal abundance of anopheline mosquitoes and their association with rainfall and malaria along the Matapí River, Amapá, Brazil. Med Vet Entomol. 2009;23:335–49.
Dennett JA, Vessey NY, Parsons RE. A comparison of seven traps used for collection of Aedes albopictus and Aedes aegypti originating from a large tire repository in Harris County (Houston), Texas. J Am Mosq Control Assoc. 2004;20:342–9.
Hoel DF, Kline DL, Allan SA. Evaluation of six mosquito traps for collection of Aedes albopictus and associated mosquito species in a suburban in north central Florida. J Am Mosq Control Assoc. 2009;25:47–57.
Lühken R, Pfitzner WP, Börstler J, Garms R, Huber K, Schork N, et al. Field evaluation of four widely used mosquito traps in Central Europe. Parasites Vectors. 2014;7:268.
L’Ambert G, Ferre J-B, Schaffner F, Fontenille D. Comparison of different trapping methods for surveillance of mosquito vectors of West Nile virus in Rhone Delta, France. J Vector Ecol. 2012;37:269–75.
Jupp PG, McIntosh BM, Nevill EM. A survey of the mosquito and Culicoides faunas at two localities in the Karoo region of South Africa with some observations on bionomics. Onderstepoort J Vet Res. 1980;47:1–6.
Crepeau TN, Healy SP, Bartlett-Healy K, Unlu I, Farajollahi A, Fonseca DM. Effects of Biogents Sentinel Trap field placement on capture rates of adult Asian tiger mosquitoes, Aedes albopictus. PLoS ONE. 2013;8:e60524.
Kraemer MUG, Sinka ME, Duda KA, Mylne AQN, Shearer FM, Barker CM, et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae albopictus. Life. 2015;4:e08347.
Hay SI, Sinka ME, Okara RM, Kabaria CW, Mbithi PM, Tago CC, et al. Developing global maps of the dominant Anopheles vectors of human malaria. PLoS Med. 2010;7:e1000209.
Farajollahi A, Kesavaraju B, Price DC, Williams GM, Healy SP, Gaugler R, et al. Field efficacy of BG-Sentinel and industry-standard traps for Aedes albopictus (Diptera: Culicidae) and West Nile virus surveillance. J Med Entomol. 2009;46:919–25.
Jupp PG. Mosquitoes of southern Africa. Hartebeespoort: Ekogilde Publishers; 1996.
Cornel AJ, Lee Y, Almeida APG, Johnson T, Mouatcho J, Venter M, et al. Mosquito community composition in South Africa and some neighboring countries. Parasites Vectors. 2018;11:331.
Sikaala CH, Killeen GF, Chanda J, Chinula D, Miller JM, Russell TL, et al. Evaluation of alternative mosquito sampling methods for malaria vectors in Lowland South-East Zambia. Parasites Vectors. 2013;6:91.
Govella NJ, Chaki PP, Mpangile JM, Killeen GF. Monitoring mosquitoes in urban Dar es Salaam: Evaluation of resting boxes, window exit traps, CDC light traps, Ifakara tent traps and human landing catches. Parasites Vectors. 2011;4:40.
Fernandes JN, Moise IK, Maranto GL, Beier JC. Revamping mosquito-borne disease control to tackle future threats. Trends Parasitol. 2018;34:359–68.
Maharaj R, Moonasar D, Baltazar C, Kunene S, Morris N. Sustaining control: lessons from the Lubombo spatial development initiative in southern Africa. Malar J. 2016;15:409.
Sharp BL, Kleinschmidt I, Streat E, Maharaj R, Barnes KI, Durrheim DN, et al. Seven years of regional malaria control collaboration Mozambique, South Africa, and Swaziland. Am J Trop Med Hyg. 2007;76:42–7.
World Health Organization. Global Health Observatory (GHO) data. Malaria country profiles. 2017. https://www.who.int/malaria/publications/country-profiles/2017/profile_zaf_en.pdf?ua=1. Accessed 20 Nov 2018.
Burke A, Dandalo L, Munhenga G, Dahan-Moss Y, Mbokazi F, Ngxongo S, et al. A new malaria vector mosquito in South Africa. Sci Rep. 2017;7:43779.
Coetzee M, Craig M, le Sueur D. Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000;16:74–7.
Coetzee M, Hunt RH, Wilkerson R, Torre AD, Coulibaly MB, Besansky NJ. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3619:246–74.
Coetzee M, Fontenille D. Advances in the study of Anopheles funestus, a major vector of malaria in Africa. Insect Biochem Mol Biol. 2004;34:599–605.
Jupp P. Mosquitoes as vectors of human disease in South Africa. S Afr Fam Pract. 2005;47:68–72.
Braack L, Gouveia de Almeida AP, Cornel AJ, Swanepoel R, de Jager C. Mosquito-borne arboviruses of African origin: review of key viruses and vectors. Parasites Vectors. 2018;11:29.
Coetzee M, Hunt RH, Braack L, Davidson G. Distribution of mosquitoes belonging to the Anopheles gambiae complex, including malaria vectors, south of latitude 15°S. S Afr J Sci. 1993;89:227–31.
Charlwood JD. Some like it hot: a differential response to changing temperatures by the malaria vectors Anopheles funestus and An. gambiae s.l.. PeerJ. 2017;5:e3099.
Kabaghe AN, Chipeta MG, Gowelo S, Mburu M, Truwah Z, McCann RS, et al. Fine-scale spatial and temporal variation of clinical malaria incidence and associated factors in children in rural Malawi: a longitudinal study. Parasites Vectors. 2018;11:129.
Hammami P, Tran A, Kemp A, Tshikae P, Kgori P, Chevalier V, et al. Rift Valley fever vector diversity and impact of meteorological and environmental factors on Culex pipiens dynamics in the Okavango Delta, Botswana. Parasites Vectors. 2016;9:434.
Ritchie SA, Buhagiar TS, Townsend M, Hoffmann A, van den Hurk AF, McMahon JL, et al. Field validation of the gravid Aedes trap (GAT) for collection of Aedes aegypti (Diptera: Culicidae). J Med Entomol. 2014;51:210–9.
Johnson BJ, Ritchie SA, Fonseca DM. The state of the art of lethal oviposition trap-based mass interventions for arboviral control. Insects. 2017;8:5.
Beechler BR, Bengis R, Swanepoel R, Paweska JT, Kemp A, van Vuren PJ, et al. Rift Valley fever in Kruger National Park: do buffalo play a role in the inter-epidemic circulation of virus? Transbound Emerg Dis. 2015;62:24–32.
Silal SP, Barnes KI, Kok G, Mabuza A, Little F. Exploring the seasonality of reported treated malaria cases in Mpumalanga, South Africa. PLoS ONE. 2013;8:e76640.
Dahan-Moss Y, Jamesboy E, Gilbert A, Burke A, Kaiser M, Mwamba M, et al. Malaria vector surveillance report, South Africa, January–December 2018. 2018. http://www.nicd.ac.za/wp-content/uploads/2019/05/MALARIA-VECTOR-SURVEILLANCE-REPORT-SOUTH-AFRICA-JANUARY-DECEMBER-2018.pdf. Accessed 16 July 2019.
Pickett STA, Cadenasso ML, Benning TL. Biotic and abiotic variability as key determinants of savanna heterogeneity at multiple spatiotemporal scales. In: du Toit JT, Rogers KH, Biggs HC, editors. The Kruger experience: ecology and management of savanna heterogeneity. Washington: Island Press; 2003. p. 22–40.
Ciota AT, Matacchiero AC, Kilpatrick AM, Kramer LD. The effect of temperature on life history traits of Culex mosquitoes. J Med Entomol. 2014;51:55–62.
Hamer GL, Anderson TK, Donovan DJ, Brawn JD, Krebs BL, Gardner AM, et al. Dispersal of adult Culex mosquitoes in an urban West Nile virus hotspot: a mark-capture study incorporating stable isotope enrichment of natural larval habitats. PLoS Negl Trop Dis. 2014;8:e2768.
Muir LE, Kay BH. Aedes aegypti survival and dispersal estimated by mark-release-recapture in northern Australia. Am J Trop Med Hyg. 1998;58:277–82.
Medeiros MCI, Boothe EC, Roark EB, Hamer GL. Dispersal of male and female Culex quinquefasciatus and Aedes albopictus mosquitoes using stable isotope enrichment. PLoS Negl Trop Dis. 2017;11:e0005347.
Verdonschot PFM, Besse-Lototskaya AA. Flight distance of mosquitoes (Culicidae): a metadata analysis to support the management of barrier zones around rewetted and newly constructed wetlands. Limnologica. 2014;45:69–79.
Johnson BJ, Hurst TP, Quoc HL, Unlu I, Freebairn C, Faraji A, et al. Field comparisons of the gravid Aedes trap (GAT) and BG-Sentinel trap for monitoring Aedes albopictus (Diptera: Culicidae) populations and notes on indoor GAT collections in Vietnam. J Med Entomol. 2017;54:340–8.
Trexler JD, Apperson CS, Schal C. Laboratory and field evaluations of oviposition responses of Aedes albopictus and Aedes triseriatus (Diptera: Culicidae) to oak leaf infusions. J Med Entomol. 1998;35:967–76.
Ponnusamy L, Xu N, Böröczky K, Wesson DM, Ayyash LA, Schal C, et al. Oviposition responses of the mosquitoes Aedes aegypti and Aedes albopictus to experimental plant infusions in laboratory bioassays. J Chem Ecol. 2010;36:709–19.
Gillies T, Coetzee M. A supplement of the Anophelinae of Africa south of the Sahara (Afrotropical Region). In: Publications of the South African Institute for Medical Research. Johannesburg: South African Institute for Medical Research; 1987. p. 1–143.
Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66:804–11.
Venter FJ, Scholes RJ, Eckhardt HC. The abiotic template and its associated vegetation pattern. In: du Toit JT, Rogers KH, Biggs HC, editors. The Kruger experience: ecology and management of savanna heterogeneity. Washington: Island Press; 2003. p. 83–129.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014.
Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Anderson MJ, Walsh DCI. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: what null hypothesis are you testing? Ecol Monogr. 2013;83:557–74.
Legendre P, Gallagher ED. Ecologically meaningful transformations for ordination of species data. Oecologia. 2001;129:271–80.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community Ecology Package. R package version 2.5-3. 2018. https://cran.r-project.org/package=vegan. Accessed Mar 2019.
Brown R, Hing CT, Fornace K, Ferguson HM. Evaluation of resting traps to examine the behaviour and ecology of mosquito vectors in an area of rapidly changing land use in Sabah, Malaysian Borneo. Parasites Vectors. 2018;11:346.
Burkett DA, Kelly R, Porter CH, Wirtz RA. Commercial mosquito trap and gravid trap oviposition media evaluation, Atlanta, Georgia. J Am Mosq Control Assoc. 2004;20:233–8.
Degener CM, Eiras ÁE, Ázara TMF, Roque RA, Rösner S, Codeço CT, et al. Evaluation of the effectiveness of mass trapping with BG-Sentinel traps for dengue vector control: a cluster randomized controlled trial in Manaus, Brazil. J Med Entomol. 2014;51:408–20.
Eiras AE, Buhagiar TS, Ritchie SA. Development of the gravid Aedes trap for the capture of adult female container-exploiting mosquitoes (Diptera: Culicidae). J Med Entomol. 2014;51:200–9.
Meeraus WH, Armistead JS, Arias JR. Field comparison of novel and gold standard traps for collecting Aedes albopictus in northern Virginia. J Am Mosq Control Assoc. 2008;24:244–8.
Ndenga BA, Mutuku FM, Ngugi HN, Mbakaya JO, Aswani P, Musunzaji PS, et al. Characteristics of Aedes aegypti adult mosquitoes in rural and urban areas of western and coastal Kenya. PLoS ONE. 2017;12:e0189971.
Ball TS, Ritchie SR. Sampling biases of the BG-Sentinel trap with respect to physiology, age, and body size of adult Aedes aegypti (Diptera: Culicidae). J Med Entomol. 2010;47:649–56.
Loecher M, Ropkins K. RgoogleMaps and loa: unleashing R graphics power on map tiles. J Stat Softw. 2015;63:1–18.
We thank South African National Parks (SANParks) for their support in conducting this study in Kruger and providing us with mosquito specimens to help with identification. We thank Purvance H. Shikwambana and a great field team: Vicky Beckers, Nina Haver, Louie Krol, Skhumbuza Mdletshe, Karabo Moloi, Nondumiso Myataza and Gijs van Nes. We thank Dr Alan Kemp, Professor A. Paulo Gouveia de Almeida and Dr Anthony Cornel for their assistance with identification.
This study was supported by the Gratama Fund from the University of Leiden (Grant Number 2016.08) to MS, the Uyttenboogaart-Eliasen foundation for comparative entomology to EEG (SUB.2016.12.08) and the RCN-IDEAS travel grant to EEG.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The number of females collected by trap and site. The mosquito column (mosq) indicates the number of females collected; the traps column of trapping nights represented. Table S2. The number of females collected of each species by region of the park. Table S3. Summary of weather conditions sampled within each region of the park. Numbers are displayed as median and range in parentheses.
Additional methods and results for the regression analysis. Table S4. Model parameters, estimates, standard error (SE) and hypothesis tests for the Poisson regression analyses in Fig. 3.
Data on the number of species collected.
. Descriptive results comparing species-specific shifts in mosquito communities collected in the net and CDC trap. Species more commonly collected in a trap are listed if 5 more were collected in that trap after all sampling days at the site. Figure S1. The apparent richness (number of unique species) and diversity for sites within each region. Figure S2. Non-metric multidimensional scaling ordinations of trap differences in mosquito communities in Malelane, Satara, Shingwedzi and Punda Maria. Figure S3. Species-specific trap preferences for the net vs CDC trap difference based on rare species not displayed in Fig. 3. Dots represent the difference in the number of mosquitoes collected in the net vs the CDC trap based on the total number of mosquitoes sampled across nights at each site. Figure S4. The net trap and the CDC trap caught higher numbers of mosquitoes (Fig. 3) and this pattern was not driven by any species or genus-specific trap bias (left figures) but by variation in the total number of the species collected (right figures). Figure S5. Dendrogram of species composition based on Bray-Curtis dissimilarity and the hierarchical clustering algorithm.
About this article
Cite this article
Gorsich, E.E., Beechler, B.R., van Bodegom, P.M. et al. A comparative assessment of adult mosquito trapping methods to estimate spatial patterns of abundance and community composition in southern Africa. Parasites Vectors 12, 462 (2019). https://doi.org/10.1186/s13071-019-3733-z
- Community composition
- Kruger National Park
- South Africa
- Trap bias