- Research
- Open access
- Published:
Mosquito breeding water parameters are important determinants for Microsporidia MB in the aquatic stages of Anopheles species
Parasites & Vectors volume 17, Article number: 509 (2024)
Abstract
Background
Microsporidia MB disrupts Plasmodium development in Anopheles mosquitoes, making it a possible biocontrol tool for malaria. As a tool for vector/disease control, its ecological distribution and the factors that determine their occurrence must be defined. We investigated the frequency of Microsporidia MB in Anopheles mosquitoes across selected sites in northern and southern Ghana, as well as the physicochemical parameters of mosquito breeding water that are associated with the occurrence of the fungus, by fitting regression models.
Methods
A non-column extraction method was used to extract DNA from the abdomens of 4255 adult Anopheles mosquitoes that emerged from larvae and pupae collected between August and October of 2021 and 2022. Detection of Microsporidia MB was achieved using quantitative PCR (qPCR), while mosquito species were molecularly identified using short interspersed nuclear elements (SINE), restriction fragment length polymorphism (RFLP) methods, and the ANOSPP algorithm.
Results
Overall Microsporidia MB distribution was 2.2% (92/4255). Male mosquitoes exhibited a higher frequency of infections and had a predicted probability of infection that was 85% higher than that of females. Sites in Ghana's Savannah zone had the highest Microsporidia MB distribution (68.5%). Biochemical oxygen demand in mosquito breeding water was estimated to be positively associated with and significantly predicts Microsporidia MB in mosquitoes with an accuracy of 94%. Increasing ammonium ion concentrations reduced the chances of finding Microsporidia MB-positive mosquitoes. According to our data, all Anopheles mosquitoes, including minor species such as An. squamosus, An. pretoriensis and An. rufipes, had equal probability of Microsporidia MB infection.
Conclusions
These results provide preliminary information on micro-ecological factors that potentially support the sustainability of Microsporidia MB infection in mosquitoes during their aquatic life stages. It will be important, therefore, to explore the impact of strategies for larval source management on these factors to ensure that the symbiont's persistence during the host's aquatic stages may not be adversely affected should it be used as an integrated approach for mosquito/disease control.
Graphical abstract

Background
Novel malaria control measures using mosquito symbionts that target the transmission of the Plasmodium parasite in malaria vectors are being explored. Candidate anti-parasitic micro-organisms for disease control are required to be generationally transferrable and environmentally restricted to the target hosts and present little or no fitness costs to the vector host, among other characteristics [1,2,3]. Microsporidia MB, a fungus, has been identified to naturally infect major malaria anopheline vectors in sub-Saharan Africa [4,5,6] and is associated with reduced Plasmodium oocysts in infected Anopheles arabiensis [5]. Although Microsporidia MB has not yet been tested for its effects on Plasmodium in Anopheles gambiae, this fungus has been reported to be avirulent in this Anopheles species and has transmission characteristics [7] similar to those observed in An. arabiensis [4]. This makes Microsporidia MB a plausible candidate for malaria transmission-blocking in anopheline mosquitoes. Fundamentally, it appears to satisfy most, if not all of the preferred characteristics of an ideal micro-organism for biocontrol of disease.
Larval source management, monitoring programs and an understanding of the ecological determinants of conducive environments for mosquito sustainment remain integral parts of mosquito and disease control. Ultimately, removing larval breeding sites, particularly those created by human activities such as urban farming, would reduce the incidence of mosquitoes and their associated diseases [8,9,10]. Due to the global rise in temperatures, the increasing incidence of rainfall events resulting from climate change [11, 12] and changes in landscapes due to urbanization [13,14,15,16], larval sites continue to emerge in new environments, and mosquitoes are adapting to new niches [17, 18]. Regardless of these changes, larval source mapping and larviciding significantly contributed to disease control [19,20,21,22]. Introducing paratransgenic micro-organisms to mosquitoes at larval breeding sites has also been proposed as a potential management strategy [23]. The underlying rationale is that mosquito larvae obtain a substantial part of their gut microbiota from the water in which they breed and can transfer these micro-organisms from one developmental stage to another [24]. In addition, the physicochemical properties of the mosquito larval breeding water impact the microbial diversity at breeding sites, mosquito larval abundance and body size [25,26,27,28,29].
The ecological distribution of Microsporidia MB is still not clearly understood. In Kenya, where this fungus was first described in An. arabiensis, frequencies were observed to peak some weeks after a rainfall event [5]. The authors of this study suggested that Microsporidia MB transmission between mosquitoes was influenced by environmental factors [5]. Other microsporidian species have been particularly observed to be more prevalent at sites where irrigated rice fields are present [30], and Microsporidia MB has also been reported near irrigation sites [5]. Whether this is mainly because of the permanent availability of breeding sites, or there are ecological factors that influence its presence, remains undefined. One potential route for using Microsporidia MB is by spore dispersal into the environment where mosquitoes would be infected by ingestion [3]; this strategy necessitates a clear understanding of the environmental requirements to explore such opportunities for control. Thus, in this study, we investigated the frequency and distribution of Microsporidia MB among Anopheles mosquitoes from Ghana, as well as the association of these with selected physicochemical properties in Anopheles larval breeding sites. We used regression models to identify those factors that significantly influence the predictability of Microsporidia MB infection in mosquitoes.
Methods
Study areas
The study was carried out in 10 areas spread across the six agro-ecological zones of Ghana (Fig. 1): Zebilla (Sudan Savannah), Sawla (Guinea Savannah), Nkoranza (Forest-Savannah Transition zone), Tarkwa and Prestea (Rain Forest), Kade and Aveyime (Semi-deciduous), Ada, Afife and Dodowa (Coastal Savannah). Most of the sites were selected out of convenience because they were visited periodically for on-going malaria surveillance activities. A few were also chosen because of previous detection of Microsporidia MB in archived mosquito DNA samples [6]. Two rice irrigation sites were purposely included to evaluate whether Microsporidia MB infection correlated with these fields, as suggested from the results of a previous study [5].
Mosquito sampling
Anopheles larvae and pupae were collected from breeding sites within the study areas. Sampling was conducted in two consecutive years: August to September 2021 and September to October 2022. These months were within periods following the peak rainfall in Ghana (Additional file 1: Figure S1). The amount of time required per collection period to obtain appreciable numbers of Anopheles-positive breeding sites for the study varied across the study areas, and we spent between 2 and 12 days to identify 5 to 13 breeding sites per area per collection year. The GPS co-ordinates of each breeding site was recorded using the KoBoCollect app. Anopheles larvae and pupae were collected by the standard dipping method using a 350-mL standard dipper or ladle and then transferred into separate transparent plastic containers with lids that were labeled according to breeding site. In field-improvised insectaries or in the Vector Labs at the Noguchi Memorial Institute for Medical Research (NMIMR), the contents of the plastic containers were transferred into labeled breeder cups and reared to adults. The emerged adult mosquitoes were initially identified based on basic morphology of the leg and wing band patterns using established keys [31]; once identified, they were euthanized by freezing at − 20 °C and stored in 1.5-ml tubes with 70% ethanol. Female and male mosquitoes were separated into different tubes for each breeding site.
Measuring water physicochemical parameters
At each breeding site, selected physicochemical parameters of the mosquito breeding water were determined prior to collecting immature mosquitoes. Water temperature, pH and dissolved oxygen (DO) level were measured in situ using a thermometer (Fisherbrand traceable Lollipop thermometer, Thermo Fisher Scientific, USA), pH meter (Ohaus pH meter starter pen; OHAUS Corp., Parsippany, NJ, USA) and dissolved oxygen meter (Setsczy), respectively, and recorded. After collecting mosquitoes from a site, we collected 2 L of the mosquito breeding water into plastic bottles and stored the water on ice while in the field. These water samples were refrigerated upon return to the laboratory at NMIMR. The water samples were submitted to the Water Research Institute of the Council for Scientific and Industrial Research, Ghana to measure other water parameters, such as biochemical oxygen demand (BOD), turbidity, salinity, nitrate-nitrogen (NO3-N), phosphate-phosphorus (PO4-P), sulfate (SO4), ammonium-nitrogen (NH4-N) and manganese (Mn), copper (Cu), calcium (Ca) and zinc (Zn) ions. These parameters were selected based on the reported germination requirements of some microsporidians [32,33,34] and were limited due to the availability of funding.
DNA extraction from mosquitoes
Adult mosquitoes were retrieved from storage and surface-sterilized in 5% bleach, 70% ethanol and 1× sterile phosphate-buffered saline (PBS). Each mosquito was sectioned into two parts, namely the head + thorax and the abdomen, respectively, with each part placed separately into 1.5-mL tubes. This strategy was adopted to reduce mosquito DNA contamination that would result from using the whole mosquito body because the mosquito abdomen contains the reproductive tissues and midgut where most Microsporidia MB are localized [4, 5, 7]. Total DNA was extracted following an adapted non-column protocol using the DNeasy Blood and Tissue DNA reagents (Qiagen, Hilden, Germany). Mosquito abdomens were rinsed briefly in sterile 1× PBS to remove any excess ethanol. The PBS was then discarded, and 1-mm diameter glass beads and 180 µL of ATL lysis buffer were added to the tubes. A blank (tube without a mosquito sample) was included with each set of samples being extracted to check for cross-contamination. The samples were homogenized in a Fisherbrand™ Bead Mill 24 Homogenizer (Thermo Fisher Scientific, Waltham, MA, USA) and then centrifuged at 12,000 g for 20 min. After centrifugation, the supernatant of the homogenized samples was carefully transferred into new tubes, 20 µL of proteinase K (10 mg/mL)was added and the tubes were incubated at 56 °C overnight. A 200-µL aliquot of buffer AL was added to each sample, and the samples mixed by vortexing and incubated at 70 °C for 10 min. The samples were then centrifuged at 12,000 g for 15 min and supernatants carefully transferred into new tubes. A 300-µL aliquot of ice-cold isopropanol was added to each tube, and the tubes mixed, incubated at − 40 °C for 1 h and then centrifuged at 12,000 g for 30 min at 4 °C. The supernatant of each tube was subsequently discarded, and the DNA pellet was washed with washing buffers AW1 and AW2 for 5 min each. The DNA pellet was dried at 50 °C for 2–4 min, reconstituted in 50 µL AE buffer and stored at − 80 °C for future use.
Identification of mosquito species
Molecular identification of mosquito species was performed with extracted DNA using PCR–restriction fragment length polymorphism (RFLP) [35] followed by short interspersed nuclear elements (SINE) PCR [36] to further distinguish between An. gambiae and Anopheles coluzzii. PCR products were run in 2% agarose gels stained with SYBR Safe and visualized on an ultraviolet (UV) transilluminator. All samples that could not be identified to species using the methods above were submitted to the ANOSPP [37, 38] platform at the Wellcome Sanger Institute for identification.
Detection of Microsporidia MB
Mosquito DNA was pooled (5–10 mosquitoes) in equal volumes according to breeding site for the detection of Microsporidia MB. When a pool tested positive for the presence of the fungus, all samples from that pool were tested individually to enable the frequency of infection to be estimated accurately. Using established Microsporidia MB-specific primers [5], Microsporidia MB DNA was amplified in a 20-µL reaction volume consisting of 10 µL Luna Universal qPCR Master Mix (NEB, Ipswich, MA, USA; catalog no. M3003), 0.5 µLeach of 10 µM forward and reverse primers, 1 µL of DNA template and 8 µL of nuclease-free water. The real-time PCR (RT-PCR) cycling conditions consisted of an initial denaturation at 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 1 min, annealing at 62 °C for 90 s and extension at 72 °C for 60 s, with a final elongation at 72 °C for 5 min. Melting curves were generated at temperatures ranging from 65 °C to 95 °C. Since RT-PCR was used qualitatively for detection, we used the amplification of the Anopheles ribosomal S7 gene with the S7F (5′-TCCTGGAGCTGGAGATGAAC-3′) and S7R (5′-GACGGGTCTGTACCTTCTGG-3′) primers as a DNA quality control. Positive controls used in the reaction were Microsporidia MB DNA obtained from Jeremy Herren’s laboratory in KEMRI. Samples were considered to be positive for Microsporidia MB if their melting curve matched that of the positive control.
Data analyses
Data were analyzed using custom scripts in R statistical software [39]. Chi-square tests were used to assess the independence of mosquito proportions to Microsporidia MB status, study sites and mosquito sex. Breeding sites from each study area were regarded as replicates. The median proportions of mosquito species identified and Microsporidia MB frequencies were compared with the Kruskal–Wallis test. All categorical factors, i.e. site, species and mosquito sex, were fitted into a generalized linear model (GLM) in a binary outcome formula to identify those that significantly correlated with the probability of detecting Microsporidia MB. A binomial distribution was fitted using all the explanatory variables singly and then as interaction terms.
Principal component analysis (PCA) biplots were used to show the correlation of water physiochemical parameters with Microsporidia MB infection. Multivariate binary logistic regression models were fitted to determine the odds of Microsporidia MB positivity given the parameters tested. Four models were fitted and evaluated as follows: model 1 = all water parameters; model 2 = reduced from model 1, using only parameters that showed significance or marginally significant (< 0.1) contribution in model 1 following the analysis of variance (ANOVA) test; models 3 and 4 were expanded forms of model 1 and 2, respectively, where ‘site’ was included to allow downstream analyses of predicted probability of Microsporidia MB at each study area. All models were compared with the ‘tab_model’ function in the sjPlot R package [40], and the best model was selected based on the Akaike Information Criteria (AIC) instead of the pseudo-R2 values because of the different number of parameters in each model, and to control for overfitting due to complexity. Receiver operating characteristic (ROC) curves were generated to show the accuracy of the selected models. Predicted probabilities of Microsporidia MB infection were estimated for each study site as a function of significant predictors in each selected model.
Results
Anopheles species distribution
A total of 4255 adult mosquitoes were obtained from rearing of the immature stages that were collected from breeding sites (Additional file 2: Table S1). Mosquito counts did not statistically differ between the sampling years (Chi-square test; χ2 = 20, df = 16, P = 0.22), but their proportions per site (relative to the total number of mosquitoes collected) were different between study sites (Kruskal test; χ2 = 66.44, df = 9, P < 0.0001). Nine confirmed and distinct Anopheles species were identified (Fig. 2), which included the minor species An. coustani, An. pharoensis, An. pretoriensis, An. rufipes and An. squamosus, together accounting for 2.9% of the mosquitoes analyzed. The predominant An. gambiae sensu stricto (An. gambiae s.s.) (38.6%), An. coluzzii (37.4%) and An. arabiensis (5.4%) were mostly identified through SINE PCR and RFLP while the other species were distinguished with the ANOSPP algorithm. About 3.7% of the samples were determined to be An. gambiae sensu lato (An. gambiae s.l.), and 1.6% were identified broadly as belonging to the genus Anopheles but could not be further distinguished. About 10% of mosquito samples remained unidentified. Zebilla, which is in the Sudan Savannah eco-zone in the north of the country, had the highest species richness, while Dodowa, in the Coastal Savannah in the south of the country, was the least diverse (Fig. 2). Mosquito species distribution, however, showed no association with ecological zones (Kruskal test; χ2 = 10.81, df = 5, P = 0.05).
Frequencies of Anopheles mosquitoes across study sites. Each bar represents the total percentage of mosquitoes from the collection attributed to a specified site. Stacked bars are colored according to species and show the percentage of Anopheles spp. contributed by each site. The study sites are arranged from the northernmost to the southernmost regions (Fig. 1)
Microsporidia MB frequency and distribution
Overall Microsporidia MB infection was 2.2% (92/4255) and was detected in six groups of mosquitoes, including the unspecified Anopheles species and mosquito samples that were unidentified. Approximately 53% of the infected mosquitoes were An. gambiae s.s. and 32% were An. coluzzii (Fig. 3a). Only one An. arabiensis (1/92) was positive for Microsporidia MB, representing 1.1% of the total number of mosquitoes infected with the fungus. Despite this, Microsporidia MB positivity did not show significant association with mosquito species (Chi-square; χ2 = 12.27, df = 11, P = 0.34).
Microsporidia MB positivity was associated with the study site (Chi-square; χ2 = 43.83, df = 9, P < 0.0005) and ecological zone (Chi-square; χ2 = 39.66, df = 5, P < 0.0005), suggesting that infections were more likely to be found in mosquitoes from some sites than from others. However, the data showed no differences in the observed Microsporidia MB frequencies among sites and ecological zones (Kruskal test; P > 0.05), implying that higher numbers of mosquitoes are probably required to observe significant differences in Microsporidia MB infections between sites. A higher number of males were Microsporidia MB positive (59% males vs 40% females; χ2 = 8.432, P = 0.004) (Fig. 3b).
Prediction of Microsporidia MB positivity in study sites, mosquito sex and species
Given the study site, mosquito sex and species, the model that best explained the occurrence of Microsporidia MB was the one that included study site and mosquito sex as multivariate non-interaction terms (Additional file 1: Table S2). In this model, study site and mosquito sex, but not mosquito species, were significant factors that correlated with Microsporidia MB positivity (Additional file 1: Table S3). Male mosquitoes demonstrated the highest odds ratio (OR) for Microsporidia MB infection (OR 1.85, 95% confidence interval [CI] 1.22, 2.85; P = 0.004), with a probability that was approximately 85% higher than in females (Table 1). With reference to Zebilla, where 12% of the total Microsporidia MB abundance was recorded, Sawla was the only site with significant odds for finding infected mosquitoes [OR 3.78; 95% CI 2.04, 7.68; P < 0.001), with 275% higher probability. The model was 71% accurate (Fig. 4a), and predicted that overall, a male mosquito had a 3.23% chance of Microsporidia MB infection while a female had 1.61%. Each study site, however, showed varied probabilities, with the highest in Sawla and the probability almost non-existent in Nkoranza, Dodowa and Aveyime (Fig. 4b; Additional file 1: Figure S2).
Assessment of the selected model for explaining Microsporidia MB status based on study site and mosquito sex. A Receiver operating characteristic (ROC) curve showing area under the curve (AUC), which depicts the accuracy of the model. B The predicted probabilities of detecting Microsporidia MB infection in mosquitoes from sampled sites, showing that the odds for Microsporidia MB positivity status are higher in males across all sites
Physicochemical parameters of the mosquito breeding water among study sites
A few measured water parameters were below the detectable limits, and to allow plotting and further analyses, such values were replaced with the next smallest positive number (Additional file 3: Table S4). All the parameters showed significant differences among the study sites except, biochemical oxygen demand (BOD), ammonium-nitrogen (NH4-N) and manganese (Mn) (Fig. 5). Nkoranza differed most with Ada and Afife, being only similar in the PO4-P, DO concentrations and temperature (Additional file 4: Table S5). Salinity was the most variable parameter, followed by sulfate (SO4), with the latter being generally higher as one moved from the north towards the south (Fig. 5).
Boxplots of estimated water parameters in the study sites showing the distribution of the physicochemical parameters among breeding sites at each study area. The horizontal line depicts the median value. P-values were calculated with the Kruskal-Wallis test to compare the median among study sites for each parameter. Significant P-values are shown in bold. BOD, Biochemical oxygen demand; DO, dissolved oxygen; NH4-N, ammonium-nitrogen; NO3-N, nitrate-nitrogen; PO4-P, phosphate-phosphorus; SAL, salinity; SO4, sulfate; TURB, turbidity
Correlation of Microsporidia MB positivity with physicochemical parameters of the mosquito breeding water
We visualized the correlation between the water parameters and the proportion of Microsporidia MB-infected mosquitoes using PCA biplots. The first two axes explained approximately 65% of the total variance observed in the data, with Ca ions, SO4 and salinity contributing most to the first principal component (Fig. 6a). BOD, PO4-P, NH4-N, DO, Mn and temperature were positively correlated with Microsporidia MB. The remaining parameters tested, with the exception of pH, SO4, Ca ions and salinity, had a negative relationship with the presence of the fungus (Fig. 6A).
Results from the best selected model that explains probability of detecting Microsporidia MB-infected mosquitoes based on physicochemical parameters of the mosquito breeding water. A Principal coordinates analysis (PCA) biplot of the correlation between Microsporidia MB and the quantitative parameters. Arrows show the direction of each variable in relation to the principal components (PC), and the lengths of the arrows represent the strength of the contribution of the variable to the PC. Arrows close together are positively correlated with each other, while those opposite each other are negatively correlated. Variables at a 90º angle have no correlation. B Receiver Operating Characteristic (ROC) curve shows the area under the curve (AUC). C Predicted probabilities of detecting Microsporidia MB infection in mosquitoes from sampled sites based on arbitrary BOD and NH4-N. BOD, Biochemical oxygen demand; DO, dissolved oxygen; NH4-N, ammonium-nitrogen; NO3-N, nitrate-nitrogen; PO4-P, phosphate-phosphorus; SO4, sulfate
Based on the best model (Additional file 1: Table S6) which incorporates turbidity, PO4-P, SO4, NH4-N, BOD, Mn, Cu, Zn, DO and site as explanatory variables for Microsporidia MB status, the odds for detecting Microsporidia MB increases by a factor of 1.04 with every unit increase in BOD [OR (95% CI) = 1.04 (1.02, 1.10); P = 0.02] and decreases with increasing units of NH4-N (Table 2). With 94% accuracy (Fig. 6B), the model predicts varied probabilities for Microsporidia MB infection at the sampled sites which appears dependent on a site-related interaction between the BOD and NH4-N (Fig. 6c; Additional file 4: Table S5).
Discussion
Microsporidia MB infection among mosquitoes is reported to be generally low [5,6,7] with seasonal variations that shows increased distribution among mosquitoes after rainfall [5]. As rainfall also influences the range and dynamics of mosquito breeding sites, identifying micro-environmental factors in mosquito breeding water that impact the sustainability of Microsporidia MB among mosquitoes could contribute to a clearer understanding of how this micro-organism can be used effectively in mosquito-borne disease control. Our data shows that Microsporidia MB positivity is not a random occurrence, but it is associated with the physicochemical properties of mosquito breeding sites. Among the parameters tested, BOD was found to have a significant positive relationship with Microsporidia MB and, depending on the breeding site, Mn ion and ammonium-nitrogen NH4-N could also be marginally inversely involved. In agreement with previous reports, we also found male mosquitoes to have a higher probability of Microsporidia MB infection than females irrespective of the study area, independent of the Anopheles species.
Microsporidians are largely obligate intracellular parasites, which implies that they need a host cell to survive. Their spores are environmentally resistant and metabolically inactive until germination and are the only microsporidian life stage that has been found to be viable outside of a host cell [41]. Microsporidians, therefore, require hosts to be in their optimal growth environment to ensure their sustainability, without necessarily using the resources available in the host’s surroundings [34]. In the present study, we found that BOD significantly raised the odds of detecting Microsporidia MB-infected larvae. BOD is a measure of the amount of oxygen available to micro-organisms during the decomposition of organic matter; as such, BOD is an important parameter for assessing organic pollution in water. While high BOD can harm aquatic life due to a reduced DO level, BOD has also been found to be linked to Anopheles larvae abundance [27, 42]. High organic content is essential for larval nutrition and protects larvae against UV radiation [43]. Together, these findings suggest that there is an increased chance of detecting Microsporidia MB-infected larvae in environments that are very conducive for Anopheles breeding. Such sites are likely to be those that are densely populated with Anopheles larvae. Given the generally low distribution of infection in Ghana [6], this implication is not surprising. Relatively more mosquitoes need to be collected to detect Microsporidia MB infection, which is more probable where large numbers of mosquitoes are found. BOD levels above 305 mg/l are, however, detrimental to Anopheles larvae [44], and such sites are not preferred by ovipositing Anopheles female mosquitoes in nature [45]. The BOD levels measured in the present study were generally well below this experimental value, with an average BOD level that was more than twofold higher at Microsporidia MB-positive sites (79.4 mg/l) than at Microsporidia MB-negative sites (32.8 mg/l). It will be useful to titrate the limits of BOD for mosquito abundance and Microsporidia MB detection to better inform its use for surveillance and evaluation of control activities.
According to our regression model, NH4-N showed decreased odds for Microsporidia MB-infected mosquitoes. Ammonium (NH4) gets into the mosquito aquatic environment through both anthropological and mosquito physiological processes. Mosquito larvae excrete ammonium into their breeding water; consequently, overcrowding would lead to increased NH4 ion concentrations to levels that become toxic to the larval population [46, 47]. Survival of larvae also decreases significantly with extraneous addition of high concentrations of NH4 to the breeding water [48] such as, for example, through agricultural fertilizers [49]. This could partly explain why some study sites, particularly those where agricultural activities were observed, such as in Aveyime and Afife, would have low mosquito numbers. However, these sites were sampled in only one collection period, which limited the total number of mosquitoes obtained in the study. Both NH4 and SO4 have been shown to correlate with increased abundance of An. arabiensis and Culex spp. in rice fields [49] and nitrogen has been shown to correlate with Culex sp. in agricultural wastewater [50]. Our results could not corroborate this association because An. arabiensis accounted for a low proportion of our catch, and those that were found were mostly collected at non-rice field sites. If the estimated water parameters are associated with Microsporidia MB as a direct function of mosquito abundance, then Microsporidia MB frequency would be expected to increase with increasing NH4-N and SO4 levels. In addition to increased levels of NH4, high SO4 also prevents larval development of Anopheles mosquitoes [45, 51]; in the present study, this was demonstrated by a negative correlation between Microsporidia MB status and SO4. In fact, SO4 estimates were on average approximately 1253% higher in rice fields (694.7 mg/l) than in non-agricultural sites and approximately 96% higher in Microsporidia MB-negative sites (156.7 mg/). Microsporidia MB positivity also correlated positively with NH4-N, but the probability of Microsporidia MB decreased with increasing levels of this ion, indicating that NH4-N could be a good indicator for mosquito abundance and the presence of Microsporidia MB. The levels of NH4-N also need to be clearly defined to identify the range of concentrations that are conducive for sustaining Microsporidia MB. According to our physicochemical parameter model, it appears that with similar BOD and NH4-N across sites, there are other ecological factors that may not have been evaluated in this study that would explain the varied probabilities of finding Microsporidia MB-infected mosquitoes at these sites. Extensive studies at agricultural sites could provide more insight into their association with the abundance of dominant Anopheles species, such as An. gambiae, An. coluzzii and An. funestus, and Microsporidia MB infectivity. Such studies are needed as the results of the present study are contrary to what has previously been shown for An. arabiensis in Kenya [5]; NH4, BOD and SO4 concentrations at these sites could also have negative implications on the use of the fungal symbiont to disease control.
Microsporidia MB has so far been isolated only from Anopheles mosquito species and, in particular, in dominant malaria-transmitting vectors, including An. funestus [4,5,6]. Plasmodium infections have also been detected in other members of the Anopheles gambiae complex, such as An. pharoensis, An. rufipes, An. squamosus and An. coustani, all of which are mostly exophilic [52,53,54,55]. While the contribution of these latter mosquito species to the epidemiology of malaria in a geographical area may be negligible in most malaria endemic regions due to their relatively low numbers, these minor vectors are also affected by control activities that are primarily designed to target the major vectors. The wide range of mosquitoes identified with the improved ANOSPP algorithm has allowed us to show, for the first time, that these minor species are equally likely to be infected with Microsporidia MB. Our results suggest that the species of mosquito is not a significant determining factor for Microsporidia MB positivity and that the probability of Microsporidia MB positivity is simply higher in any Anopheles male mosquito. A fruitful line of research would be to investigate the frequency of Microsporidia MB infections in geographical locations where these vectors are found in high numbers to confirm that these vectors can also be infected with the fungus. The results of such studies will further expand the range at which Microsporidia MB-Anopheles research can be extended for a more inclusive Microsporidia MB-mediated control strategy [3].
Conclusions
The implications for mosquito populations naturally infected with Microsporidia MB are highly promising for future control mechanisms involving this endosymbiont. Our study used regression models to show that a number of biotic and abiotic factors can be incorporated into surveillance models for mosquito larval source management, thus maximizing the detection of Microsporidia MB-infected populations in a geographical location. Such models will also be useful in assessing the feasibility and field dynamics of using Microsporidia MB for control, especially with increasing chemical use in agricultural activities [56].
Availability of data and materials
The dataset(s) supporting the conclusions of this article is(are) included within the article (and its additional file(s)).
References
Gao H, Cui C, Wang L, Jacobs-Lorena M, Wang S. Mosquito microbiota and implications for disease control. Trends Parasitol. 2020;36:98–111.
Wang S, Ghosh AK, Bongio N, Stebbings KA, Lampe DJ, Jacobs-Lorena M. Fighting malaria with engineered symbiotic bacteria from vector mosquitoes. Proc Natl Acad Sci USA. 2012;109:12734–9.
Bukhari T, Pevsner R, Herren JK. Microsporidia: a promising vector control tool for residual malaria transmission. Front Trop Dis. 2022;3:957109.
Nattoh G, Maina T, Makhulu EE, Mbaisi L, Mararo E, Otieno FG, et al. Horizontal transmission of the symbiont Microsporidia MB in Anopheles arabiensis. Front Microbiol. 2021;12:647183.
Herren JK, Mbaisi L, Mararo E, Makhulu EE, Mobegi VA, Butungi H, et al. A microsporidian impairs Plasmodium falciparum transmission in Anopheles arabiensis mosquitoes. Nat Commun. 2020;11:2187.
Akorli J, Akorli EA, Tetteh SNA, Amlalo GK, Opoku M, Pwalia R, et al. Microsporidia MB is found predominantly associated with Anopheles gambiae s.s. and Anopheles coluzzii in Ghana. Sci Rep. 2021;11(1):18658. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-021-98268-2.
Nattoh G, Onyango B, Makhulu EE, Omoke D, Angango LM, Kamau L, et al. Microsporidia MB in the primary malaria vector Anopheles gambiae sensu stricto is avirulent and undergoes maternal and horizontal transmission. Parasit Vectors. 2023;16:335.
Matthys B, Vounatsou P, Raso G, Tschannen AB, Becket EG, Gosoniu L, et al. Urban farming and malaria risk factors in a medium-sized town in Cote d’Ivoire. Am J Trop Med Hyg. 2006;75:1223–31.
Yadouléton A, N’Guessan R, Allagbé H, Asidi A, Boko M, Osse R, et al. The impact of the expansion of urban vegetable farming on malaria transmission in major cities of Benin. Parasit Vectors. 2010;3:118.
Klinkenberg E, McCall P, Wilson MD, Amerasinghe FP, Donnelly MJ. Impact of urban agriculture on malaria vectors in Accra. Ghana Malar J. 2008;7:1–9.
Anoopkumar AN, Aneesh EM. A critical assessment of mosquito control and the influence of climate change on mosquito-borne disease epidemics. Environ Dev Sustain. 2022;24:8900–29.
Sewe MO, Sjödin H, Rocklöv J, Tompkins AM, Casallas A, Colón-González FJ, et al. Projecting the risk of mosquito-borne diseases in a warmer and more populated world: a multi-model, multi-scenario intercomparison modelling study. Lancet Planet Health. 2021;5:404–18.
Ferraguti M, Martínez-de la Puente J, Roiz D, Ruiz S, Soriguer R, Figuerola J. Effects of landscape anthropization on mosquito community composition and abundance. Sci Rep. 2016;6:29002.
Labbo R, Fandeur T, Jeanne I, Czeher C, Williams E, Arzika I, et al. Ecology of urban malaria vectors in Niamey. Repub Niger Malar J. 2016;15:314.
Antonio-Nkondjio C, Defo-Talom B, Tagne-Fotso R, Tene-Fossog B, Ndo C, Lehman LG, et al. High mosquito burden and malaria transmission in a district of the city of Douala. Cameroon BMC Infect Dis. 2012;12:275.
Keating J, Macintyre K, Mbogo CM, Githure JI, Beier JC. Characterization of potential larval habitats for Anopheles mosquitoes in relation to urban land-use in Malindi. Int J Health Geogr. 2004;3:9.
Afrane YA, Githeko AK, Yan G. The ecology of Anopheles mosquitoes under climate change: case studies from the effects of deforestation in East African highlands. Ann N Y Acad Sci. 2012;1249:204–10.
Sinka ME, Pironon S, Massey NC, Longbottom J, Hemingway J, Moyes CL, et al. A new malaria vector in Africa: predicting the expansion range of Anopheles stephensi and identifying the urban populations at risk. Proc Natl Acad Sci USA. 2020;117:24900–8.
Antonio-Nkondjio C, Sandjo NN, Awono-Ambene P, Wondji CS. Implementing a larviciding efficacy or effectiveness control intervention against malaria vectors: Key parameters for success. Parasit Vectors. 2018;11:57.
García GA, Fuseini G, Mba Nlang JA, Nsue Maye VO, Bela NR, Wofford RN, et al. Evaluation of a multi-season, community-based larval source management program on Bioko Island, Equatorial Guinea. Front Trop Dis. 2022;3:846955.
Jahir A, Kahamba NF, Knols TO, Jackson G, Patty NFA, Shivdasani S, et al. Mass trapping and larval source management for mosquito elimination on Small Maldivian Islands. Insects. 2022;13:805.
Dongus S, Nyika D, Kannady K, Mtasiwa D, Mshinda H, Fillinger U, et al. Participatory mapping of target areas to enable operational larval source management to suppress malaria vector mosquitoes in Dar es Salaam Tanzania. Int J Health Geogr. 2007;6:37.
Ratcliffe NA, Furtado Pacheco JP, Dyson P, Castro HC, Gonzalez MS, Azambuja P, et al. Overview of paratransgenesis as a strategy to control pathogen transmission by insect vectors. Parasit Vectors. 2022;15:112.
Lindh JM, Borg-Karlson K, Faye I. Transstadial and horizontal transfer of bacteria within a colony of Anopheles gambiae (Diptera: Culicidae) and oviposition response to bacteria-containing water. Acta Trop. 2008;107:242–50.
Hessou-Djossou D, Djègbè I, Ahadji-Dabla KM, Nonfodji OM, Tchigossou G, Djouaka R, et al. Diversity of larval habitats of Anopheles mosquitoes in urban areas of Benin and influence of their physicochemical and bacteriological characteristics on larval density. Parasit Vectors. 2022;15:207.
Seal M, Chatterjee S. Combined effect of physico-chemical and microbial quality of breeding habitat water on oviposition of malarial vector Anopheles subpictus. PLoS ONE. 2023;18:e0282825.
Kinga H, Kengne-Ouafo JA, King SA, Egyirifa RK, Aboagye-Antwi F, Akorli J. Water physicochemical parameters and microbial composition distinguish Anopheles and Culex mosquito breeding sites: Potential as ecological markers for larval source surveillance. J Med Entomol. 2022;59:1817–26.
Akeju VA, Adeniyi Olusi T, Adepeju S-O. Effect of physicochemical parameters on Anopheles mosquitoes larval composition in Akure North Local Government area of Ondo State, Nigeria. JOBAZ. 2022;83:34.
Ouédraogo WM, Toé KH, Sombié A, Viana M, Bougouma C, Sanon A, et al. Impact of physicochemical parameters of Aedes aegypti breeding habitats on mosquito productivity and the size of emerged adult mosquitoes in Ouagadougou City Burkina Faso. Parasit Vectors. 2022;15:478.
Service MW. Mortalities of the immature stages of species B of the Anopheles gambiae complex in Kenya: comparison between rice fields and temporary pools, identification of predators, and effects of insecticidal spraying. J Med Entomol. 1977;13:535–45.
Coetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:70.
Hayman JR, Southern TR, Nash TE. Role of sulfated glycans in adherence of the microsporidian Encephalitozoon intestinalis to host cells in vitro. Infect Immun. 2005;73:841–8.
El-Taweel HA, Tolba MM, Sadaka HA, El-Zawawy LA, Osman MM. Zinc PVA versus potassium dichromate for preservation of microsporidian spores of human origin. Parasitol Res. 2012;111:689–94.
Keohane EM, Weiss LM. The structure, function, and composition of the microsporidian polar tube. In: Wittner M, Weiss L, editors. The microsporidia and microsporidiosis. Washington: ASM Press; 1994. p. 196–224.
Fanello C, Santolamazza F, Della Torre A. Simultaneous identification of species and molecular forms of the Anopheles gambiae complex by PCR-RFLP. Med Vet Entomol. 2002;16:461–4.
Santolamazza F, Mancini E, Simard F, Qi Y, Tu Z, Della TA. Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar J. 2008;7:163.
Boddé M, Makunin A, Ayala D, Bouafou L, Diabaté A, Ekpo UF, et al. High-resolution species assignment of Anopheles mosquitoes using k-mer distances on targeted sequences. Elife. 2022;11:e78775.
Boddé M, Makunin A, Teltscher F, Akorli J, Andoh NE, Bei AK, et al. Improved species assignments across the entire Anopheles genus using targeted sequencing. Front Genet. 2024;15:1456644.
R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016. http://www.r-project.org
Lüdecke D. sjPlot: Data visualization for statistics in social science. R package; 2024. https://cran.r-project.org/web/packages/sjPlot/sjPlot.pdf
Keeling PJ, Fast NM. Microsporidia: biology and evolution of highly reduced intracellular parasites. Annu Rev Microbiol. 2024;56:93–116.
Ononamadu CJ, Datit JT, Imam AA. Insecticide resistance profile of Anopheles gambiae mosquitoes: a study of a residential and industrial breeding sites in Kano Metropolis, Nigeria. Environ Health Insights. 2020;14:1178630219897272.
Berry NL, Overholt EP, Fisher TJ, Williamson CE. Dissolved organic matter protects mosquito larvae from damaging solar UV radiation. PLoS ONE. 2020;15:e0240261.
Tchuinkam T, Mpoame M, Make-Mveinhya B, Simard F, Lélé-Defo E, Zébazé-Togouet S, et al. Optimization of breeding output for larval stage of Anopheles gambiae (Diptera: Culicidae): prospects for the creation and maintenance of laboratory colony from wild isolates. Bull Entomol Res. 2011;101:259–69.
Omondi C, Nonoh J, Ntabo R. Characterization of Anopheles funestus larval habitats in Fiyoni, Kwale County, Kenya: insights into malaria vector ecology and control. Eur J Ecol. 2023;9:40–51.
Cochran DG. Excretory systems. In: Blum MS, editor. Fundamentals of insect physiology. New York: John Wiley and Sons; 1985. p. 91–138.
Clements AN. The biology of mosquitoes, volume 1: development, nutrition and reproduction. CABI: Wallingford; 1992.
Walker ED. Toxicity of sulfide and ammonium to Aedes triseriatus larvae (Diptera: Culicidae) in water-filled tree holes and tires. J Med Entomol. 2016;53:577–83.
Mutero CM, Nganga PN, Wekoyela P, Githure J, Konradsen F. Ammonium sulphate fertiliser increases larval populations of Anopheles arabiensis and culicine mosquitoes in rice fields. Acta Trop. 2004;89:187–92.
Sanford MR, Chan K, Walton WE. Effects of inorganic nitrogen enrichment on mosquitoes (Diptera: Culicidae) and the associated aquatic community in constructed treatment wetlands. J Med Entomol. 2005;42:766–76.
Oyewole IO, Momoh OO, Anyasor GN, Ogunnowo AA, Ibidapo CA, Oduola OA, et al. Physico-chemical characteristics of Anopheles breeding sites: impact on fecundity and progeny development. Afr J Environ Sci Tech. 2009;3:447–52.
Mwangangi JM, Muturi EJ, Muriu SM, Nzovu J, Midega JT, Mbogo C. The role of Anopheles arabiensis and Anopheles coustani in indoor and outdoor malaria transmission in Taveta District, Kenya. Parasit Vectors. 2013;6:114.
Stevenson JC, Simubali L, Mbambara S, Musonda M, Mweetwa S, Mudenda T, et al. Detection of Plasmodium falciparum infection in Anopheles squamosus (Diptera: Culicidae) in an area targeted for malaria elimination, Southern Zambia. J Med Entomol. 2016;53:1482–7.
Carrara GC, Petrarca V, Niang M, Coluzzi M. Anopheles pharoensis and transmission of Plasmodium falciparum in the Senegal River delta West Africa. Med Vet Entomol. 1990;4:421–4.
Saili K, de Jager C, Sangoro OP, Nkya TE, Masaninga F, Mwenya M, et al. Anopheles rufipes implicated in malaria transmission both indoors and outdoors alongside Anopheles funestus and Anopheles arabiensis in rural south-east Zambia. Malar J. 2023;22:1–12.
Philbert A, Lyantagaye SL, Nkwengulila G. Farmers’ pesticide usage practices in the malaria endemic region of North-Western Tanzania: implications to the control of malaria vectors. BMC Public Health. 2019;19:1456.
Acknowledgements
We thank the staff of Professor Samuel Dadzie’s lab at NMIMR for accommodating us in their field collection activities across the country. It would have been difficult to travel to certain sites without this support. We also thank Dr. Jeremy Herren, KEMRI for sharing Microsporidia MB DNA samples which we used as the positive control for our assay.
Funding
This study was funded by ANTI-VeC African Anopheles Symbiont Survey grant award to JA (AV/AASS/009) through Open Philanthropy Project (University of Glasgow). The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.
Author information
Authors and Affiliations
Contributions
JA conceived the project, sought funding and coordinated the execution. NEA, RKE, CD and SO conducted field sampling. NEA, EAA, RKE, SNAT, RMP, DBS and SKDO conducted mosquito dissections. NEA, EAA, RKE, RMP, DBS and SKDO perfomed molecular detection of Microsporidia MB. EAA, NEA and JA analyzed the results and drafted the manuscript. SKD supported fieldwork planning, and logistics. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1:
Monthly average rainfall information for Ghana in 2021 (blue) and 2022 (green) obtained from ClimateEngine.org (https://www.climateengine.org/). Figure S2: Predicted probabilities of Microsporidia MB infection in study sites according to the selected model status ~ site + sex. Table S2: Comparison of logistic regression models fitted to identify significant categorical factors (study site, mosquito species and sex) correlating with Microsporidia MB positivity. Table S3: Results for test of statistical relationship between study site, mosquito species and sex as interaction terms for Microsporidia MB occurrence. Table S6: Comparison of outputs for models on water abiotic parameters association with Microsporidia MB.
Additional file 2: Table S1:
Summary of mosquito collected and Microsporidia MB relative abundance per site.
Additional file 3: Table S4:
Estimated water parameters for selected breeding sites at the study areas.
Additional file 4: Table S5:
Results of Dunn test with Benjamini-Hochberg correction for pairwise comparisons, following Kruskal Wallis test of variance for water parameters among the study sites.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Akorli, E.A., Andoh, N.E., Egyirifa, R.K. et al. Mosquito breeding water parameters are important determinants for Microsporidia MB in the aquatic stages of Anopheles species. Parasites Vectors 17, 509 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-024-06596-9
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-024-06596-9