Monitoring mosquito nuisance for the development of a citizen science approach for malaria vector surveillance in Rwanda
Malaria Journal volume 20, Article number: 36 (2021)
Many countries, including Rwanda, have mosquito monitoring programmes in place to support decision making in the fight against malaria. However, these programmes can be costly, and require technical (entomological) expertise. Involving citizens in data collection can greatly support such activities, but this has not yet been thoroughly investigated in a rural African context.
Prior to the implementation of such a citizen-science approach, a household entomological survey was conducted in October–November 2017 and repeated one year later in Busoro and Ruhuha sectors, in southern and eastern province of Rwanda, respectively. The goal was to evaluate the perception of mosquito nuisance reported by citizens as a potential indicator for malaria vector hotspots. Firstly, mosquito abundance and species composition were determined using Centers for Disease Control and Prevention (CDC) light traps inside the houses. Secondly, household members were interviewed about malaria risk factors and their perceived level of mosquito nuisance.
Tiled roofs, walls made of mud and wood, as well as the number of occupants in the house were predictors for the number of mosquitoes (Culicidae) in the houses, while the presence of eaves plus walls made of mud and wood were predictors for malaria vector abundance. Perception of mosquito nuisance reported indoors tended to be significantly correlated with the number of Anopheles gambiae sensu lato (s.l.) and Culicidae collected indoors, but this varied across years and sectors. At the village level, nuisance also significantly correlated with An. gambiae s.l. and total mosquito density, but only in 2018 while not in 2017.
Perception of mosquito nuisance denoted in a questionnaire survey could be used as a global indicator of malaria vector hotspots. Hence, involving citizens in such activities can complement malaria vector surveillance and control.
Malaria remains a public health concern in Rwanda despite the gains made in malaria reduction in the past decades . Since 2012, malaria has increased every year, thereby impeding the progress made up to 2011 . From 2012 to 2016, the country reported an eight-fold increase in malaria cases. Additionally, malaria-related deaths increased from 325 in 2012 to 663 in 2016. The eastern and southern parts of the country have been the most afflicted regions. The increase in malaria cases has been observed in all districts, including districts that were previously defined as being at pre-elimination phase. This increase was observed across all age groups, suggesting the entire population is at risk of acquiring a malaria infection .
In Rwanda, malaria vector surveillance is carried out monthly in 12 sentinel sites across the country. It aims to monitor key parameters of malaria vectors, and provides entomological data, such as entomological inoculation rate (EIR), to guide the planning of vector control interventions [3, 4]. At present, the monitoring of malaria mosquito density is combined with malaria incidence and helps to determine the spatio-temporal spread of infections. However, active mosquito surveillance is conducted in only 3% (12/416) of the sectors of Rwandan territory, while in the remaining 97%, mosquito monitoring is not implemented . The current mosquito surveillance approach is challenging to be implemented in all areas of the country due to the inadequate local capacity in entomology and the high costs related to extending vector surveillance countrywide, hence alternative approaches are desired .
Passive surveillance through citizen-science is a tool to track mosquito presence and spread, such as of Aedes albopictus, a potential vector of dengue, Zika and chikungunya viruses [6, 7]. Studies carried out in The Netherlands and Spain showed that citizens provided scientifically valuable information through questionnaires and through sending mosquito samples to laboratories in charge of mosquito identification. This can ease cost constraints for mosquito surveillance [6, 8]. Mosquito nuisance reported by citizens via a questionnaire in combination with actual mosquito samples collected by citizens revealed the presence or absence of two known biotypes of Culex pipiens and their hybrids, which can be important vectors of West Nile virus .
In the current study, prior to the implementation of a wider citizen-science programme for malaria vector surveillance in Ruhuha, Rwanda, mosquito species composition was determined using a conventional mosquito trapping method in two rural sectors by means of two cross-sectional surveys performed in 2017 and 2018. In addition, factors that could explain the observed spatial distribution of mosquito species collected were analysed. Lastly, it was investigated whether perceived mosquito nuisance reported by the participants could provide an indication of potential (malaria) vector hotspots, especially in areas where entomological surveillance for malaria vectors is not implemented.
Household and entomological surveys were carried out in Ruhuha and Busoro sectors, located respectively in Bugesera district (eastern province), and in Nyanza district (southern province) in Rwanda (Fig. 1). The choice of the study sites was based on the large number of malaria cases reported in Ruhuha since 2012, and the long-term working relationship with the local health centre. Ruhuha sector covers 54 km2 and is sub-divided into 35 villages . An estimated 24,000 people are living in more than 5100 households (HHs) . The area is a high malaria-endemic zone . Busoro is a sector covering 74 km2 and is sub-divided into 41 villages. It has a total population of approximately 34,000 people living in 8,000 HHs . Irrigated rice fields are the main type of land use in both sectors. Both sectors differ from each other by the area of irrigated rice fields which is 40% larger in Busoro (178 hectares or 2.4% of the land surface) in comparison with Ruhuha (93 hectares or 1.7% of the land surface) [13, 14]. Additionally, Ruhuha is located near the shores of Lake Cyohoha South . The wetlands potentially serve as favourable places for mosquito breeding.
Household and entomological surveys were conducted in Ruhuha and Busoro for three weeks in October–November 2017 and were repeated in a modified form in the same period in 2018. Data from the household survey were coupled with entomological data collected for both sectors and years. In addition to the demographic characteristics, factors were defined that could explain the variation in mosquito abundance. Additionally, the relationship between perceived mosquito nuisance experienced by the citizens in their house and mosquitoes collected using CDC light traps in Busoro and Ruhuha was assessed.
The studied HHs were selected among those that were part of a larger household survey . Six villages per sector were selected by simple random sampling. In Busoro sector, the villages of Kireranyana, Gikombe, Muhindo, Karambi, Runazi and Rucyamo were selected while in Ruhuha sector, Kagasera, Kibaza, Kiyovu, Kamweru, Rusenyi, and Mubano villages were selected for the study (Fig. 1). At the village level, lists of HHs were provided by the village leaders, and a systematic random sampling was used to draw a sample of HHs to be visited. As a result, 30 to 31 HHs were selected for each village for the large household survey. From these HHs, 6–8 HHs were selected for the household survey. In 2017, 82 HHs were thus selected randomly (42 from Busoro and 40 from Ruhuha) and considered for the study. In 2018, the same survey was repeated by selecting randomly and directly from the list, excluding HHs during the previous survey in 2017. In other words, HHs interviewed in 2017 were different from those from 2018. In 2018, 84 HHs were selected, 42 from Busoro and 42 from Ruhuha. If household members were absent at the time of the interview, the interviewers progressed to the next selected house. Only residents above 18 years of age and who consented to participate in the study were interviewed. The house occupied by the interviewed resident was selected for the mosquito collection after obtaining informed consent.
For the household survey conducted in 2017 and 2018, perception of mosquito nuisance experienced indoors and per season was assessed by the citizens in Busoro and Ruhuha. Six to 8 household surveys were carried out per day. The average interview time was 45 min up to one hour per household. The questions were originally written in English and translated into the local language (Kinyarwanda), and prior to the survey translated back into English. The questionnaire was pre-tested in a pilot study of 10 HHs selected randomly in the neighbouring Mareba sector, eastern province, for its consistency, and then revised consequently. Field data collectors were trained in addressing the questionnaire before conducting the survey. Only participants who consented were interviewed in Kinyarwanda. All translations were made by professional translators including members of the project team and cross-checked by native speakers. The translators were asked to review and cross-check the items and identify any problems in wording, terminology, ability to understand, and relevance.
Questionnaire data were collected in electronic forms using Open Data Kit (ODK) Collect set-up [12, 16]. The questionnaire encompassed questions with closed- and open-ended questions. The questionnaire contained different sections with questions related to demographic data of the participants, perception of mosquito nuisance and its seasonality, knowledge on malaria, and to vector and malaria control practices. For statistical analysis, data on demographic characteristics, household characteristics and mosquito nuisance perceived indoors by the respondents were used.
Demographic data included gender, age, marital status, education level, occupation, Ubudehe category (a community-based social categorization of household and dependents into different groups based on their income) , size of the household, bed net ownership, and mobile phone ownership of the participants. Other characteristics, such as house features (type of wall (mud/clay or wood), type of floor (cement or mud), type of roof (iron or tile sheets), presence of eaves, and livestock ownership (species owned and location where they were kept) were also included in the questionnaire.
Perception of mosquito nuisance
Respondents were asked to answer the question whether they experienced mosquito nuisance in their environment. Participants who experienced mosquito nuisance were requested to scale the level of nuisance experienced in the house on a 6-point Likert scale (0 = no nuisance, 1 = very little to 5 = very much). They were also requested to indicate when (rainy or dry seasons) they experienced mosquito nuisance.
Mosquito collection by CDC light traps
Mosquitoes were collected in 2017 and 2018 in 166 selected houses in Ruhuha (82) and Busoro (84) among the interviewed participants using a miniature CDC light trap. The CDC light traps were set up in the bedroom and hung at the foot end of the bed, with the shield of the trap at 1.5 m from the floor . The traps were set up at 18:00, and the owner of the room was instructed to put off the light of the trap and tie the bag connected to the collecting cup at 06:00 the next morning to avoid mosquitoes escaping from the traps. After their collection from the traps, mosquitoes were stored in labelled petri dishes before morphological identification in the laboratory.
Mosquitoes were identified using standard morphological identification keys for anophelines and culicines [19, 20]. Mosquitoes were scored as unfed or blood fed. The mosquitoes were then stored under silica gel in labelled Eppendorf tubes with the codes of the respondent interviewed and kept in an envelope with the name of the village under cool conditions at the central Laboratory of Entomology for identification of Plasmodium falciparum infection and further molecular species identification.
Blood meal identification
For samples collected in 2017, each mosquito abdomen was ground in 100 µl of phosphate-buffered saline (PBS), and then filled up to 1 ml PBS. Blood meals were identified by direct enzyme-linked immunosorbent assay (ELISA) using antihost (IgG) conjugate against human and cow proteins in a single-step assay . The non-reacting samples were then tested subsequently using goat IgG. ELISA results were read visually . The anthropophilic rate was determined as the proportion of mosquitoes that fed exclusively on human blood among all fed mosquitoes.
The head and thorax of all female Anopheles gambiae sensu lato (s.l.) collected in 2017 using CDC light traps were used to test for the presence of circumsporozoite protein (CSP) of P. falciparum using ELISA. A sample with an optical density (OD) value above the cut-off (cut-off = 2 × mean OD of 7 negative samples) was considered positive . The sporozoite rate was calculated as the number of mosquitoes infected with P. falciparum sporozoites divided by the total number of mosquitoes processed.
Molecular species identification
For molecular identification of the sibling species of An. gambiae s.l., a random sample of 9% (n = 233/2514) of the total number of An. gambiae s.l. collected during the two years and from both sectors were identified using the rDNA-polymerase chain reaction (PCR) assay . If the initial PCR testing failed to amplify a sample, then the PCR analysis was repeated once or twice until successful amplification was achieved. If a sample could not be identified after three rounds of PCR, it was scored as unknown .
Household survey data were imported from ODK software into Microsoft Excel (2016) and checked for consistency in the values and answers. They were electronically loaded onto a central server for backup, translated and coupled with data from the entomological survey. Statistical analysis was undertaken in SPSS 23.0 (SPSS Inc, Chicago, IL, USA), and included the calculation of frequencies and Chi-squares statistics. Bivariate analysis of correlation between the dependent variable (number of An. gambiae s.l. or Culicidae) and the independent variables (number of members of the household, house structural features (floor, wall, roof), species of animal kept in the house and presence of eaves was performed to determine predictors that could explain the mosquito abundance indoors. Only predictors that had a screening significance lower than 0.1 were then considered for the final models. For this purpose, generalized linear models (GLM, negative binomial with log link) were used. Besides house structural features, household size and livestock ownership, other factors included in the GLM were sector (Ruhuha/Busoro) and the year of study (2017/2018). All data on mosquito collections were entered into Excel to calculate the sporozoite rate and human blood index and the various mosquito species identified were summarized as proportions. Furthermore, Spearman correlation analysis was used to analyse the relationship between mosquito nuisance reports by the respondents and number of mosquitoes and species collected in the respondents’ houses.
One-hundred and sixty-six respondents were enrolled both in the household survey and mosquito collection, out of 167 (n = 166/167, 99%) requested respondents to participate. Table 1 provides an overview of the demographic characteristics of the participants. Overall, both sectors were similar in their demographics. In both sectors more than half of the respondents were female (64%). The average age of the respondents in Busoro was 44 years (n = 84), and 41 years (n = 82) in Ruhuha. Almost 76% of the respondents were schooled (no category) while the remainder (24%) was unschooled. Respondents from Busoro were more highly educated (81%) than those from Ruhuha (71%). Most of the respondents were farmers (93%), followed by self-employed (1%), private officer (1%), student (1%), and unemployed (4%). However, more respondents in Busoro than in Ruhuha were farmers (98 vs 89%). In the latter sector, more individuals were unemployed. There was no significant difference when comparing Ubudehe categories (a community-based social categorization of household and dependents into different groups based on their incomes) between the two sectors. Most of the participants owned a bed net (74%), but there was a significant difference in bed net ownership between Busoro and Ruhuha (81 vs 66%). Almost half of the participants owned a mobile phone (46%). There was no significant difference in mobile phone ownership when comparing both sectors.
Most of the participants had mud floor houses (83%) with mud walls (59%), closed eaves (65%), and iron sheet roofing (69%) (Table 2). There were no significant differences between both sectors for house features (eaves presence, floor, wall features), except for iron sheeting (100 in Ruhuha vs 38% in Busoro) (Table 2).
Overall, there was a significant difference in livestock ownership when comparing both sectors (P = 0.029). Of the respondents, 72% owned at least one species of livestock (cows, pigs, poultry, rabbits, goats or sheep; Table 3). However, participants in Busoro owned more livestock than participants in Ruhuha (80 vs 65%). Goats were the most often owned livestock. Poultry ownership differed significantly between Busoro and Ruhuha (P = 0.003): respondents in Busoro owned more poultry indoor than Ruhuha.
Mosquito species composition by CDC light traps
From the mosquitoes collected in 2017 and 2018, 74% (n = 7370) were collected in 2017 (Additional file 1) and 26% (n = 2595) in 2018 (Additional file 2). Of all mosquitoes, 74.2% were morphologically identified as culicines and 25.2% as anophelines. Among female mosquitoes collected, 25.6% (n = 2210) were fed and 77.8% (n = 7755) were unfed. Of the total anophelines collected (n = 2514), 94.2% was An. gambiae s.l. Other Anopheles species collected included Anopheles brohieri (0.2%), Anopheles funestus (1.2%), Anopheles maculipalpis (1.1%), Anopheles pharoensis (0.6%), Anopheles rufipes (0.5%), and Anopheles ziemanni (2.2%). Busoro recorded 72.8% of the total An. gambiae s.l. collected, the remaining 27.2% was collected in Ruhuha. Anopheles funestus (n = 31), An. pharoensis (n = 14) and An. ziemanni (n = 56) were the most frequently encountered other human-biting Anopheles species identified in Busoro, while An. brohieri (n = 2),) An. rufipes (n = 13), and An. maculipalpis (n = 5) were Anopheles species known as non-human-biting. In Ruhuha, An. pharoensis (n = 1) was another human biting malaria vector, and An. brohieri (n = 2), An. maculipalpis (n = 22) were the non-human-biting anophelines. The malaria vector An. funestus was not collected in Ruhuha.
From 7451 culicines, 99.7% (n = 7425) consisted of Culex spp. and 0.3% (n = 26) of Mansonia spp. The highest proportion of culicines identified was Cx. quinquefasciatus (99.6%), with 89.2% (n = 6647) from Busoro and 10.4% (n = 778) from Ruhuha, respectively.
Of 1,046 Anopheles spp. collected, 100 were selected randomly. Ninety-eight specimens were An. gambiae s.l., while the other two were An. maculipalpis and An. rufipes, respectively. Of all 100 blood-fed specimens, 82 (82%) had fed on a single host (human, goat or bovine), while 5 (5%) had fed on mixed hosts and 13 were unspecified for the antigens assayed, suggesting that these An. gambiae s.l. had fed on other hosts than humans, goats or cattle (Table 4). For the remaining 87 Anopheles, 58% (n = 58) of the An. gambiae s.l. were engorged with human blood and 1% (n = 1) An. gambiae s.l. was engorged with blood of both human and goat origin.
Plasmodium falciparum infection rates
Of the 1,046 mosquitoes tested by ELISA, P. falciparum CSP antigen was detected in 14 out of 1013 (1.3%), 11 out of 971 tested (1.1%) An. gambiae s.l. and 3 out of 42 tested An. ziemanni (7.1%). The overall sporozoite rate of anopheline species was 1.3%. The infection rate was 0.5% (3/573) in Busoro and 2.3% (11/473) in Ruhuha.
Molecular identification of members of the Anopheles gambiae complex
Of the 9.4% (236/2514) of the An. gambiae s.l. selected for sibling species identification from 2017 and 2018, 145 (61.4%) were identified as An. gambiae sensu stricto (s.s.) and 74 (31.4%) as Anopheles arabiensis (Additional file 3). Sixteen (6.8%) were not amplified and one sample was contaminated (0.4%).
Factors explaining mosquito abundance
After bivariate analysis, different species kept in the houses, such as poultry, rabbit, goat, and sheep (Table 3), did not show a statistically significant correlation (P < 0.1) with the dependent variables. Therefore, they were not selected for GLM analysis. The house features that included the materials used for the construction of the roof and the wall, and the numbers of occupants in the house, were predictors for the number of mosquitoes (Culicidae) in the houses. The predictors roof, wall and size of the household had a statistically significant effect on the number of mosquitoes indoors, while the floor composition did not contribute statistically to a difference (Table 5).
There were significantly more Culicidae in Busoro than in Ruhuha, and the incidence rate ratio for 2017 was 2.5 times that of 2018 (P < 0.001; Table 5). Houses with tiled roofs were more exposed to mosquitoes than houses with an iron roof (P = 0.002; Table 5). Likewise, houses with walls made of mud and wood had a larger number of mosquitoes (Culicidae) than houses with walls made with mud (P < 0.001; Table 5).
Similarly, for the total number of only An. gambiae s.l., there were more females collected in Busoro than in Ruhuha. However, the year effect was not significant (Table 5). Houses with closed eaves and mud walls were more likely to not have malaria vectors resting inside the house than houses with open eaves or walls made with wood and mud.
Perceived mosquito nuisance per sector
In total, 96% (n = 159, both 2017 and 2018) of the respondents reported to have experienced at least some mosquito nuisance. There was a significant difference between both sectors in the proportion of respondents that reported mosquito nuisance (P = 0.050) (Table 6). Of those that did perceive nuisance, 6% (n = 10) experienced “very little”, 25% (n = 43) experienced “little”, 17% experienced “some”, 32% experienced “much” (n = 52, 32%) and 16% reported “very much” nuisance. There was no significant difference between the two sectors when comparing the level of mosquito nuisance experienced indoors (P = 0.177; Table 7).
Temporal variation in perceived mosquito nuisance
When asked in what season respondents perceive mosquito nuisance, most respondents reported to experience “very much” nuisance during the long rainy season (March–May; Fig. 2b). Interestingly, more respondents in Busoro (n = 73) perceived mosquito nuisance in this season than in Ruhuha (n = 54). Respondents perceived “little” to “very much” nuisance during the short rainy season (September until November; Fig. 2a). Respondents mostly perceived “some” to “little” nuisance during the small rainy season (December-February; Fig. 2c) and mostly “very little” to “little” in the big dry season (June until August; Fig. 2d).
Spatial variation in perceived mosquito nuisance
Within each sector, there was substantial spatial variation at village level in the degree of mosquito nuisance perceived indoors (Fig. 3). In Busoro, the highest perceived mosquito nuisance was from Runazi (x̅ = 3.9), Rucyamo (x̅ = 3.8), and Kireranyana (x̅ = 3.6), followed by Gikombe (x̅ = 3.1), Muhindo (x̅ = 2.9) and Karambi (x̅ = 2.5) (Fig. 3). In Ruhuha, 82 households from Kamweru (x̅ = 3.5) and Rusenyi (x̅ = 3.6) reported to experience much nuisance in their houses while the remaining participants from Kibaza (x̅ = 2.9), Kiyovu (x̅ = 2.9), and Mubano (x̅ = 3.1) reported having experienced some nuisance. Interestingly, respondents from Kagasera reported having experienced little mosquito nuisance (x̅ = 1.8) during the period of study.
Relationship between mosquito abundance and perceived indoor mosquito nuisance
First, a Spearman correlation coefficient was computed to assess the relationship between the number of mosquitoes (total Culicidae and total An. gambiae s.l.) collected and the perceived nuisance level experienced in the houses by the participants. When data from both years and sites were aggregated, there was a moderate but significant correlation between the total Culicidae and perceived nuisance indoors (rs = 0.316, n = 166, P < 0.001). Similarly, there was a moderate, but significant correlation between the number of An. gambiae s.l. collected indoors and the perceived mosquito nuisance (rs = 0.281, n = 166, P < 0.001).
When the data were analysed per year and per sector for Culicidae, there was a strong spatial difference. No significant correlation was observed between total mosquito numbers and perceived nuisance indoors in Busoro for both years (2017: rs = 0.055, n = 42, P = 0.731 and 2018: rs = 0.200, n = 42, P = 0.204; Fig. 4c, d), whereas there were moderate, but significant correlations between the total mosquito numbers and perceived mosquito nuisance in Ruhuha for both years of study (2017: rs = 0.389, n = 40, P = 0.013 and 2018: rs = 0.305, n = 42, P = 0.049; Fig. 4a, b; Table 8).
For An. gambiae s.l. per sector and year, there was a strong temporal difference. There were significant correlations in both sectors in 2017 (Busoro: rs = 0.37, n = 42, P = 0.016; Ruhuha: rs = 0.45, n = 40, P = 0.004; Fig. 5a, c). However, in 2018, these significant correlations between perceived nuisance and An. gambiae s.l. were absent for both sectors (Busoro: rs = 0.119, n = 42, P = 0.452; Ruhuha: rs = 0.251, n = 40, P = 0.109; Fig. 5b, d; Table 8).
When data were analysed one spatial level lower, i.e., by village, there were significant correlations between average nuisance level reported and An. gambiae s.l. (rs = 0.607, P = 0.002, Fig. 6a), as well as between average nuisance level and total mosquitoes (rs = 0.528, P = 0.008, Fig. 6b). When analysed per year separately, the correlations were strong and significant for 2018 (black dots, Fig. 6), but not significant for the data from 2017 (grey dots, Fig. 6).
The mosquito nuisance derived from the questionnaires revealed a significant nuisance caused by mosquitoes, notably from the Culicidae. The findings show that there was a high mosquito density including Culex species and mainly An. gambiae s.l. for both years in both sectors. A higher density of mosquitoes and malaria vectors were partially explained by the house construction materials for roof and walls, as well as by the presence of eaves and number of occupants in the house. These factors also contribute to a high level of perceived mosquito nuisance, although correlations were space (sector) and time (year) dependent. The presence of P. falciparum-infected mosquitoes in the houses, as shown for the 2017 data, contributes to the risk of contracting malaria. The results are supported by several other studies that show that individuals who live in rural areas in poorly constructed houses are exposed to more mosquito bites, and hence to an increased risk of malaria transmission [26,27,28,29,30,31]. In a study conducted in Kenya, houses made of wood and mud exhibited a significant effect on mosquito abundance in the houses . Factors such as closed eaves reduced rates of house entry by anopheline mosquitoes compared to fully open eaves as was also demonstrated in other studies from Tanzania, The Gambia and Kenya [33,34,35,36,37].
Spatio-temporal variations in mosquito abundance were observed between Busoro and Ruhuha. In Busoro, mosquitoes, including numbers of An. gambiae s.l., were collected in significantly higher numbers in comparison with Ruhuha. This spatial difference may be explained by the fact that although both sectors have wetlands, Busoro is characterized by a larger area (178 ha) than Ruhuha (93 ha) that is dedicated to rice irrigation . The difference was strongly influenced by collections from one of the villages (Rucyamo), which is the village closest to the irrigated fields (Fig. 1) and which contributed to 65% of all mosquitoes collected. This village was also the village with the next highest level of mosquito nuisance reported. Living near the rice field, the chance of more mosquitoes and to experience more nuisance was higher because the wetlands provide good mosquito breeding habitat [38, 39].
In the present study, the predominant sibling species of the An. gambiae complex was An. gambiae s.s. for both years and sectors combined (66%). This dominance of An. gambiae s.s. is similar to a study conducted in one site near Kigali City in 2007 by President’s Malaria Initiative-Rwanda, in which it was reported that An. gambiae s.s. accounted for 93.6% of the total of 157 An. gambiae s.l. examined by PCR, while An. arabiensis accounted for the remaining 6.4% . However, the finding was contrary to a study where the characterization of An. gambiae s.l. from 10 sentinel sites revealed that the predominant sibling species was An. arabiensis (83%) . Although An. funestus was recognized as the dominant Anopheles species in previous studies from Rwanda [42, 43], this species was collected only in Rucyamo, the village closest to the more permanently inundated wetlands and rice fields, which are ideal habitats for this species.
For 2017, the P. falciparum infection rate in Busoro was higher compared to the infection rate in Ruhuha. Other anopheline species, though collected in low numbers, should not be neglected in strategies for malaria control and elimination, because they can transmit other mosquito-borne diseases, such as the Babanki virus (BBKV) that is transmitted by An. brohieri, as found in Senegal . Albeit in low numbers, both An. pharoensis and An. ziemanni can transmit P. falciparum as observed in studies from Cameroon, Ethiopia, Guinea Bissau, Tchad, and Kenya, respectively [38, 39, 45,46,47]. Although Culex spp. have not been incriminated as vector of disease in Rwanda, the species caused a high burden of nuisance in Cameroon . Considering this group of species as potential vector for other vector-borne diseases will be important in the framework of integrated vector management. Reducing their numbers would substantially reduce mosquito nuisance experienced, and thus enhance community involvement in uptake of vector control for malaria prevention .
The respondents experienced mosquitoes as a significant problem in their daily life, especially during the rainy season that lasts from March to May. The increase in perception of mosquito nuisance over the seasons corresponds to the increase of vector density at the start of the long rainy season, which is mainly due to an increase in wetland area providing suitable habitats for larval development of mosquitoes. The percentage of participants that reported much nuisance was higher in Busoro (87%) than in Ruhuha (67%). Living in the vicinity of marshlands increased the chance to experience higher mosquito abundance and hence higher nuisance level as observed from participants from Busoro and Ruhuha [13, 14]. In addition, the presence of blocked ditches produced by exploitation of sand for house construction was noted. Such ditches are known to be artificial mosquito breeding sites [50, 51].
The correlations found between nuisance and number of mosquitoes can be explained as larger numbers of mosquitoes collected indoors will result in more biting activity and, hence a higher level of mosquito nuisance. The results suggest that levels of perceived mosquito nuisance are in some way indicative of mosquito densities indoors. Consequently, it could be argued that perceived mosquito nuisance in the peridomestic area can be used as a global indicator for malaria transmission risk. Such data can be obtained by filling out a questionnaire indicating the level of nuisance expressed on a Likert scale. A study conducted in Algeria demonstrated that perception of citizens can help to identify occurrence of Aedes albopictus in a residential neighbourhood in Bir-Khadem . This helped to put in place vector control measures that could prevent the propagation of Ae. albopictus to other areas and to avoid the massive use of insecticides for vector control, which could ultimately lead to insecticide resistance .
It should be noted that when considering each sector separately, perceived mosquito nuisance was significantly correlated to the numbers of An. gambiae when data from both years were added together, while for Culicidae, perceived mosquito nuisance was correlated to the number of Culicidae for Ruhuha but not for Busoro, even after adding datasets of both study years together. The reasons why there was no correlation for Busoro remains unclear. This may be explained by the fact that in Busoro mosquito densities were more extreme than in Ruhuha and that variation in nuisance by these high densities could no longer be caught in the level of perceived mosquito nuisance.
We have shown that data on mosquito nuisance can indicate malaria vector abundance and hence identify malaria transmission risk. An important next step, however, is to scale-up the approach, validate it in other settings and possibly integrate it in vector surveillance efforts. ICT-based technologies, such as apps and web-based platforms, can greatly support such initiatives, as evidenced by the success of the Mosquito Alert application in Spain . This requires investment costs upfront (ICT, human resources etc.), but these may largely be paid off by boosting vector surveillance data collection, as well as by building durable partnerships between public health authorities and citizens .
Poor housing construction significantly led to increased malaria vector density and thus possibly malaria risk in rural Rwanda. This suggests that good house construction needs to be considered as one of the vector control strategies that can be provided for poor populations. At the largest scale in this study, i.e., if data for years and sectors are combined, the relationships between the level of perceived mosquito nuisance and mosquito density at family and species level were clearly shown. Perception of mosquito nuisance denoted in a questionnaire survey could be used as an indicator of mosquito abundance and, hence, for An. gambiae s.l. occurrence. Involving citizens in reporting the level of mosquito nuisance can contribute to improved malaria vector surveillance and control.
Availability of data and materials
Data supporting the conclusions of this article are included within the article. The datasets used and/or analysed during the present study are available from the corresponding author upon reasonable request.
Centers for Disease Control
Entomological inoculation rate
Open Data Kit
Phosphate Buffered Saline solution
Enzyme-Linked ImmunoSorbent Assay
Polymerase chain reaction
Generalized linear model
President’s malaria initiative
Karema C, Aregawi MW, Rukundo A, Kabayiza A, Mulindahabi M, Fall IS, et al. Trends in malaria cases, hospital admissions and deaths following scale-up of anti-malarial interventions, 2000–2010. Rwanda Malar J. 2012;11:236.
MOH. Revised National Malaria Contingency Plan 2016–2020. Kigali; 2017.
Tusting LS, Bousema T, Smith DL, Drakeley C. Measuring changes in Plasmodium falciparum transmission. Precision, accuracy and costs of metrics. Adv Parasitol. 2014;84:151–208.
Hakizimana E, Karema C, Munyakanage D, Githure J, Mazarati JB, Tongren JE, et al. Spatio-temporal distribution of mosquitoes and risk of malaria infection in Rwanda. Acta Trop. 2018;182:149–57.
MOH. National Strategic Plan for Integrated Vector Management (2013 - 2017). Kigali, Rwanda; 2013.
Palmer JRB, Oltra A, Collantes F, Delgado JA, Lucientes J, Delacour S, et al. Citizen science provides a reliable and scalable tool to track disease-carrying mosquitoes. Nat Commun. 2017;8:916.
Heym EC, Schröder J, Kampen H, Walther D. The nuisance mosquito Anopheles plumbeus (Stephens, 1828) in Germany—a questionnaire survey may help support surveillance and control. Front Public Health. 2017;5:278.
Vogels CBF, van de Peppel LJJ, van Vliet AJH, Westenberg M, Ibañez-Justicia A, Stroo A, et al. Winter activity and aboveground hybridization between the two biotypes of the west nile virus vector Culex pipiens. Vector-Borne Zoonotic Dis. 2015;15:619–26.
Ingabire CM, Alaii J, Hakizimana E, Kateera F, Muhimuzi D, Nieuwold I, et al. Community mobilization for malaria elimination: application of an open space methodology in Ruhuha sector. Rwanda Malar J. 2014;13:167.
Kateera F, Ingabire CM, Hakizimana E, Kalinda P, Mens PF, Grobusch MP, et al. Malaria, anaemia and under-nutrition: three frequently co-existing conditions among preschool children in rural Rwanda. Malar J. 2015;14:440.
Rulisa S, Kateera F, Bizimana JP, Agaba S, Dukuzumuremyi J, Baas L. Malaria prevalence, spatial clustering and risk factors in a low endemic area of eastern Rwanda: a cross sectional study. PLoS ONE. 2013;8:e69443.
Asingizwe D, Poortvliet PM, Koenraadt CJM, Van VAJH, Ingabire CM, Mutesa L, et al. Role of individual perceptions in the consistent use of malaria preventive measures : mixed methods evidence from rural Rwanda. Malar J. 2019;18:270.
Musabe BL. Crop Intensification program as a tool of poverty reduction in Rwanda.: University of Rwanda; 2012.
Hakizimana E. Integrated vector management (IVM) as a tool for community empowerment towards malaria elimination in Rwanda. Wageningen University; 2019.
Rwanyiziri G, Rugema J. Climate Change Effects on Food Security in Rwanda: Case study of wetland rice production in Bugesera district. Rwanda J. 2013;1:35–51.
Raja A, Tridane A, Gaffar A, Lindquist T, Pribadi K. Android and ODK based data collection framework to aid in epidemiological analysis. J Public Health Inform. 2014;5:228.
Fenny AP, Yates R, Thompson R. Social health insurance schemes in Africa leave out the poor. Int Health. 2018;10:1–3.
Mboera LEG, Kihonda J, Braks MAH, Knols BGJ. Influence of centers for disease control light trap position, relative to a human-baited bednet, on catches of Anopheles gambiae and Culex quinqufasciatus in Tanzania. Am J Trop Med Hyg. 1998;59:595–6.
Gillies MT, Coetzee M. A Supplement to the Anophelinae of Africa South of the Sahara. Publ South African Inst Med Res. 1987;55:63.
Service MW. Mosquito ecology: field sampling methods. New York: Elsevier; 1993.
Beier JC, Perkins PV, Wirtz RA, Koros J, Diggs D, Gargan TP, et al. Bloodmeal identification by direct enzyme-linked immunosorbent assay (ELISA), tested on Anopheles (Diptera: Culicidae) in Kenya. J Med Entomol. 1988;25:9–16.
Beier JC, Koros JK. Visual assessment of sporozoite and bloodmeal ELISA samples in malaria field studies. J Med Entomol. 1991;28:805–8.
Dotson EM, Linser P, Koekemoer L, Benedict MQ. Methods in Anopheles. 2nd Edition. Atlanta; 2011.
Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am Soc Trop Med Hyg. 1993;49:520–9.
Minakawa N, Mutero CM, Githure JI, Beier JC, Yan G. Spatial distribution and habitat characterization of Anopheline mosquito larvae in western Kenya. Am J Trop Med Hyg. 1999;61:1010–6.
Snyman K, Mwangwa F, Bigira V, Kapisi J, Clark TD, Osterbauer B, et al. Poor housing construction associated with increased malaria incidence in a cohort of young ugandan children. Am J Trop Med Hyg. 2015;92:1207–13.
Kateera F, Mens PF, Hakizimana E, Ingabire CM, Muragijemariya L, Karinda P, et al. Malaria parasite carriage and risk determinants in a rural population: a malariometric survey in Rwanda. Malar J. 2015;14:16.
Tuyishimire J. Spatial modelling of malaria risk factors in Ruhuha. Rwanda J Ser - Life Nat Sci. 2013;1:17.
Jatta E, Jawara M, Bradley J, Jeffries D, Kandeh B, Knudsen JB, et al. How house design affects malaria mosquito density, temperature, and relative humidity : an experimental study in rural Gambia. Lancet Planet Health. 2015;2:e498-508.
Kaindoa EW, Finda M, Kiplagat J, Mkandawile G, Nyoni A, Coetzee M, et al. Housing gaps, mosquitoes and public viewpoints: a mixed methods assessment of relationships between house characteristics, malaria vector biting risk and community perspectives in rural Tanzania. Malar J. 2018;17:298.
Charlwood JD, Pinto J, Ferrara PR, Sousa CA, Ferreira C, Gil V, et al. Raised houses reduce mosquito bites. Malar J. 2003;2:45.
Zhou G, Munga S, Minakawa N, Githeko AK, Yan G. Spatial relationship between adult malaria vector abundance and environmental factors in western kenya highlands. Am J Trop Med Hyg. 2007;77:29–35.
Ogoma SB, Lweitoijera DW, Ngonyani H, Furer B, Tanya L, Mukabana WR, et al. Screening mosquito house entry points as a potential method for integrated control of endophagic filariasis, arbovirus and malaria vectors. PLoS Negl Trop Dis. 2010;4:e773.
Njie M, Dilger E, Lindsay SW, Kirby MJ. Importance of eaves to house entry by anopheline, but not culicine, mosquitoes. J Med Entomol. 2009;46:505–10.
Menger DJ, Omusula P, Wouters K, Oketch C, Carreira AS, Durka M, et al. Eave screening and push-pull tactics to reduce house entry by vectors of malaria. Am J Trop Med Hyg. 2016;94:868–78.
Mburu MM, Juurlink M, Spitzen J, Moraga P, Hiscox A, Mzilahowa T, et al. Impact of partially and fully closed eaves on house entry rates by mosquitoes. Parasit Vectors. 2018;11:383.
Ondiba IM, Oyieke FA, Ong GO, Olumula MM, Nyamongo IK, Estambale BBA. Malaria vector abundance is associated with house structures in Baringo County. Kenya PLoS One. 2018;13:e0198970.
Tabue RN, Nem T, Atangana J, Bigoga JD, Patchoke S, Tchouine F, et al. Anopheles ziemanni a locally important malaria vector in Ndop health district, north west region of Cameroon. Parasit Vectors. 2014;7:262.
Sanford MR, Cornel AJ, Nieman CC, Dinis J, Marsden CD, Weakley AM, et al. Plasmodium falciparum infection rates for some Anopheles spp. from Guinea-Bissau, west Africa. F1000Res. 2014;3:243.
Lansana K. Rwanda entomology monitoring: Technical report presidential malaria initiative (PMI), Kigali-Rwanda. 2008.
Hakizimana E, Karema C, Munyakanage D, Iranzi G, Githure J, Tongren JE, et al. Susceptibility of Anopheles gambiae to insecticides used for malaria vector control in Rwanda. Malar J. 2016;15:582.
Munyantore S. Historique de la lutte antipaludique au Rwanda. Rev Med Rwandaise. 1989;21:14–28.
Loevinsohn ME. Climatic warming and increased malaria incidence in Rwanda. Lancet. 1994;343:714–8.
Diagne N, Fontenille D, Konate L, Faye O, Lamizana MT, Legros F, et al. Les Anopheles du Senegal. Bull Soc Path Exot. 1994;87:267–77.
Hinzoumbé CK, Péka M, Antonio- C, Donan-gouni I, Awono-ambene P, Samè- A, et al. Malaria vectors and transmission dynamics in Goulmoun, a rural city in south-western Chad. BMC Infect Dis. 2009;9:71.
Kibret S, Lautze J, Boelee E, Mccartney M. How does an Ethiopian dam increase malaria ? Entomological determinants around the Koka reservoir. Trop Med Int Health. 2012;17:1320–8.
Kamau L, Mulaya N, Vulule JM. Evaluation of potential role of Anopheles ziemanni in malaria transmission in western Kenya. J Med Entomol. 2006;43:774–6.
Nchoutpouen E, Talipouo A, Djiappi-tchamen B, Djamouko- L, Kopya E, Ngadjeu CS, et al. Culex species diversity, susceptibility to insecticides and role as potential vector of lymphatic filariasis in the city of Yaounde Cameroon. PLoS Negl Trop Dis. 2019;13:e0007229.
Ingabire CM, Rulisa A, Van Kempen L, Muvunyi C, Koenraadt CJM, Van Vugt M, et al. Factors impeding the acceptability and use of malaria preventive measures: implications for malaria elimination in eastern Rwanda. Malar J. 2015;14:136.
Castro MC, Kanamori S, Kannady K, Mkude S, Killeen GF, Fillinger U. The importance of drains for the larval development of lymphatic filariasis and malaria vectors in dares salaam, United Republic of Tanzania. PLoS Negl Trop Dis. 2010;4:e693.
De Silva PM, Marshall JM. Factors contributing to urban malaria transmission in sub-saharan Africa: a systematic review. J Trop Med. 2012;2012:819563.
Benallal KE, Garni R, Bouiba L, Harrat Z. First Detection of Aedes (Stegomyia) albopictus (Diptera: Culicidae ) in Algiers, the capital city of Algeria. J Arthropod-Borne Dis. 2019;13:420–5.
Ashepet MC, Jacobs L, van Oudheusden M, Huyse T. Wicked solution for wicked problems: citizen science for vector-borne disease control in Africa. Trends Parasitol. 2020;S1471–4922(20):30289.
We are grateful to the Ruhuha and Busoro health centres’ leadership and research assistants for their support and participation in the study. Special thanks to Jeanine Loonen for her technical assistance and Domina Asingizwe who participated in field activities. We also thank the lay community members for their willingness to participate in the study.
This publication is part of the project “Environmental Virtual Observatories for Connective Action (EVOCA)”, project duration 2016–2020, which is funded by Wageningen University, The Netherlands through its Interdisciplinary Research and Education Fund (INREF) (2100710609).
Ethics approval and consent to participate
The ethical approval was guaranteed to the study (Approval Notice: No 405/ CMHS/IRB/2016) by the Institutional Review Board of the College of Medicine and Health Sciences, University of Rwanda.
Consent for publication
Participants agreed on the publication of findings reported in this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Mosquito species collected using CDC light traps in selected villages in Busoro and Ruhuha sector, Rwanda (2017)
Mosquito species collected using CDC light traps in selected villages in Busoro and Ruhuha sector, Rwanda (2018)
Members of the Anopheles gambiae complex found among samples of Anopheles gambiae s.l. tested from Ruhuha and Busoro sector, Ruhuha
About this article
Cite this article
Murindahabi, M.M., Takken, W., Misago, X. et al. Monitoring mosquito nuisance for the development of a citizen science approach for malaria vector surveillance in Rwanda. Malar J 20, 36 (2021). https://doi.org/10.1186/s12936-020-03579-w