Skip to main content

Rapid classification of epidemiologically relevant age categories of the malaria vector, Anopheles funestus

Abstract

Background

Accurately determining the age and survival probabilities of adult mosquitoes is crucial for understanding parasite transmission, evaluating the effectiveness of control interventions and assessing disease risk in communities. This study was aimed at demonstrating the rapid identification of epidemiologically relevant age categories of Anopheles funestus, a major Afro-tropical malaria vector, through the innovative combination of infrared spectroscopy and machine learning, instead of the cumbersome practice of dissecting mosquito ovaries to estimate age based on parity status.

Methods

Anopheles funestus larvae were collected in rural south-eastern Tanzania and reared in an insectary. Emerging adult females were sorted by age (1–16 days old) and preserved using silica gel. Polymerase chain reaction (PCR) confirmation was conducted using DNA extracted from mosquito legs to verify the presence of An. funestus and to eliminate undesired mosquitoes. Mid-infrared spectra were obtained by scanning the heads and thoraces of the mosquitoes using an attenuated total reflection–Fourier transform infrared (ATR–FT-IR) spectrometer. The spectra (N = 2084) were divided into two epidemiologically relevant age groups: 1–9 days (young, non-infectious) and 10–16 days (old, potentially infectious). The dimensionality of the spectra was reduced using principal component analysis, and then a set of machine learning and multi-layer perceptron (MLP) models were trained using the spectra to predict the mosquito age categories.

Results

The best-performing model, XGBoost, achieved overall accuracy of 87%, with classification accuracy of 89% for young and 84% for old An. funestus. When the most important spectral features influencing the model performance were selected to train a new model, the overall accuracy increased slightly to 89%. The MLP model, utilizing the significant spectral features, achieved higher classification accuracy of 95% and 94% for the young and old An. funestus, respectively. After dimensionality reduction, the MLP achieved 93% accuracy for both age categories.

Conclusions

This study shows how machine learning can quickly classify epidemiologically relevant age groups of An. funestus based on their mid-infrared spectra. Having been previously applied to An. gambiae, An. arabiensis and An. coluzzii, this demonstration on An. funestus underscores the potential of this low-cost, reagent-free technique for widespread use on all the major Afro-tropical malaria vectors. Future research should demonstrate how such machine-derived age classifications in field-collected mosquitoes correlate with malaria in human populations.

Graphical Abstract

Background

Despite significant investments in malaria control and research, there were still an estimated 249 million malaria cases and 619,000 deaths in 2021 globally, a significant majority of which occurred in sub-Saharan Africa [1]. Other than the poor economic conditions and weak health systems, the continued high burden of malaria in Africa is attributable to key biological threats, notably malaria parasite resistance to drugs [2,3,4], vector resistance to insecticides [5, 6], increasing occurrence of malaria parasites evading detection by rapid diagnostic tests [7,8,9,10,11], and disruptions from major disease outbreaks such as Ebola and COVID-19 [12,13,14]. Effective vector control, primarily with insecticide-treated nets (ITNs) and indoor residual spraying (IRS), has been the most important component of malaria control in Africa [15]. However, its continued effectiveness requires active innovation to address the current threats, and improved understanding of the major vector species in different settings.

Anopheles funestus is one of the four main malaria vector species in sub-Saharan Africa, the others being An. gambiae, An. arabiensis and An. coluzzii, and it is also one of the most widespread [16,17,18,19]. Anopheles funestus is particularly important in East and Southern Africa, where it is becoming the dominant malaria vector. For example, in parts of Tanzania, An. funestus is reported to be responsible for 86–97% of all new malaria infections [17, 20,21,22]. Its dominance is due to multiple factors, including (i) being highly anthropophilic, and thus preferring to bite humans over other vertebrates [17, 23], (ii) being highly endophilic, i.e. preferring to bite inside human dwellings rather than outside [24], (iii) having significantly higher survival rates than other species [25], (iv) being resistant to commonly used insecticides [17, 18, 26] and (v) preferentially breeding in perennial habitats with year-round productivity [27]. Given its importance and dominance in malaria transmission systems, vector surveillance programs in the respective countries should be designed with special attention to this vector species.

Besides evaluating biting densities and Plasmodium infection rates, accurately determining the age and survival of An. funestus is crucial for monitoring transmission dynamics and assessing the effectiveness of vector control interventions such as ITNs and IRS. Dissection of mosquito ovaries is still the main entomological technique for estimating the age of vector populations [28]. The dissections are usually performed under light microscopes to assess the reproductive history, specifically the parity status, of the mosquitoes. This involves observing whether the ovaries contain coiled tracheolar skeins (indicating non-parous mosquitoes) or stretched-out tracheoles (indicating parous mosquitoes). Non-parous mosquitoes are considered young in this case, whereas parous mosquitoes are considered old and may carry the malaria parasites, having had multiple blood-feedings [28]. Unfortunately, these dissections tend to be laborious and time-consuming, especially when dissecting large numbers of mosquitoes, and are impractical on a large scale.

Furthermore, the reliability of mosquito dissections is limited by their reproductive history. For instance, a female mosquito can have more than one blood meal but still not oviposit, a scenario known as gonotrophic discordance or pre-gravid blood meal [29]. Moreover, since the gonotrophic cycles of Anopheles mosquitoes can be as short as 2–3 days under optimal climatic conditions [30, 31], it is possible for parous mosquitoes to be relatively young, and in rare cases for nulliparous mosquitoes to be several days old due to the scarcity of blood meals (e.g. when ITN coverage and usage is high). Therefore, using parity alone to distinguish between epidemiologically distinct age categories of adult mosquitoes, especially in the context of malaria transmission, which requires 10–14 days of incubation [32], is not always realistic.

All these concerns suggest the need for alternative age-grading techniques that are easy to perform cheaply at scale and can provide accurate representations of epidemiologically important mosquito age categories and populations. The alternative mosquito age-grading methods currently include the analysis of cuticular hydrocarbon patterns using gas chromatography [33] and gene transcription [34,35,36]. Near-infrared spectroscopy (NIRS) (12,500 cm−1 to 4000 cm−1 frequencies) [37], which involves passing infrared light through a mosquito sample to measure the absorbance or reflectance of the organic compound functional groups, has also been used to estimate ages for various mosquito species of both laboratory-reared and wild-collected mosquitoes [38,39,40,41,42,43,44].

More recently, mid-infrared spectroscopy (MIRS) has been used to predict and estimate mosquito age, recording the biochemical composition of mosquito samples at longer wavelength frequencies [45,46,47]. In addition, machine learning (ML) techniques, including convolutional neural networks, have been utilized to differentiate MIRS spectra associated with distinct mosquito ages and species in both laboratory and wild mosquitoes [46, 47]. The infrared-based systems have so far been successful for various applications on three of the four main African malaria vectors (i.e. An. gambiae sensu stricto, An. arabiensis and An. coluzzii [46]), but have yet to be demonstrated for An. funestus. The goal of this study was therefore to test whether a similar ML-MIRS approach could classify adult female An. funestus mosquitoes derived from wild-caught larvae into two epidemiologically relevant age categories: young (0–9 days old, too young to have mature Plasmodium sporozoites in their salivary glands) and old (10 days or older, potentially carrying mature Plasmodium sporozoites given the right climatic conditions), factoring in a parasite incubation period of 10–14 days.

Methods

Mosquito collection

Third and fourth instar mosquito larvae were collected from known aquatic habitats of An. funestus in five different villages in south-eastern Tanzania, namely Tulizamoyo (8.3669°S, 36.7336°E), Kilisa (8.3721°S, 36.5584°E), Lupiro (8.38333°S, 36.66667°E), Ikwambi (7.9833°S, 36.8184°E) and Ruaha (8.9068°S, 36.7185°E). The larvae were transported to the vector biology laboratory (VectorSphere) at the Ifakara Health Institute for further rearing. The larvae were kept in water from their natural breeding habitats and were fed TetraMin® fish food.

Once they pupated, the pupae were separated from the larvae and placed in emergence cages. The emergent adult mosquitoes were maintained at 26–28 °C, 70–85% relative humidity and a 12:12 h light/dark photoperiod, on a 10% sugar solution diet.

Mosquito preservation and scanning

Female adults were collected and individually preserved according to their age, from 1 to 16 days old. A total of 2084 mosquitoes were collected. The female mosquitoes were killed using chloroform and subsequently stored in separate 1.5-ml microcentrifuge tubes containing silica gel for desiccation. The heads and thoraces of the individual female mosquitoes were scanned using an attenuated total reflection–Fourier-transform infrared spectrometer (ATR FT-IR) to obtain mid-infrared spectra with a resolution of 2 cm−1 at 4000–400 cm−1 frequencies as described previously, complete with background spectral calibration [45, 48, 49]. For each sample, 16 sample scans were averaged to obtain the primary output spectrum [46].

Mosquito identification

Although the field collections had been performed in known An. funestus habitats, it was necessary to confirm the identity of the mosquitoes and eliminate any unwanted species. This was accomplished primarily by morphology-based taxonomy using keys of Afro-tropical Anopheles [50] but was complemented by polymerase chain reaction (PCR) identification to sort between sibling species in the An. funestus group. Wild An. funestus complex DNA was extracted from the two legs of adult female mosquitoes. The two legs of an individual An. funestus mosquito were placed separately in 1.5-ml microcentrifuge tubes, followed by 20 µl of TE (Tris-EDTA) buffer, and incubated at 95 °C for 15 min. PCR was then used to differentiate An. funestus from other sibling species, using species-specific primers targeting the non-coding internal transcribed spacer (ITS2) region using the protocol described by Koekmoer et al. [51]. The PCR reaction was performed in a 25 µl volume, consisting of a PCR mixture of 2.5 µl 10× reaction buffer, 25 mM MgCl2, 10 pmol/µl of each primer, 8 mM of each dNTP, 5 units of thermostable Taq DNA polymerase and 3 µl of DNA template. The PCR products were analysed by electrophoresis in 2.5% agarose gel stained with classic view DNA dye for visualization of DNA bands. Only An. funestus mosquitoes were considered for further analysis, and any other species were discarded.

Machine learning

Mosquito spectra with low intensity, abnormal background or atmospheric interferences (with water vapor and carbon dioxide) were discarded [45]. The data from the remaining spectra (N = 2084) were processed and analysed in Python 3.9 using scikit-learn [52] and TensorFlow 2.0 [53, 54]. The data were rescaled using the StandardScaler algorithm, with a mean of 0 and a standard deviation of 1.

Using the StratifiedShuffleSplit algorithm, the dataset was split into training (n = 1875) and test/unseen (n = 209) sets. To train the supervised ML models, An. funestus ages were used as training labels. Anopheles funestus, ranging from 1 to 16 days old, were divided into two epidemiologically relevant age categories, taking into consideration the incubation period of malaria parasites of 10–14 days [32]. The first group included An. funestus that were between 1 and 9 days old and were considered young and incapable of transmitting malaria (i.e. non-infectious age group). The second group included An. funestus that were between 10 and 16 days old and were considered old enough to be capable of transmitting malaria given the right environmental conditions (i.e. potentially infectious).

Multiple standard ML classifiers, including k-nearest neighbours (KNN), logistic regression (LR), support vector machine (SVM), random forest (RF) and extreme gradient boosting (XGBoost), were compared to determine which model predicted the data with the highest classification accuracy. The best-performing model was further optimized by fine-tuning its hyperparameters. The top 100 spectral features (wavenumbers) with the most influence on the model predictions were identified and utilized to reduce the dimensionality of the spectra data, followed by retraining of the best ML classifier.

Moreover, two multi-layer perceptron (MLP) models were trained by reducing the dimensionality of the spectra data using different inputs: (1) the top 100 features extracted from the best-performing ML classifier, and (2) principal components using the scikit-learn library. Both MLP models had six fully connected layers, each containing 500 neurons, to enable the model to learn from the network's weights, as demonstrated previously [47]. To prevent overfitting, a dropout layer with a rate of 0.5 was used, and early stopping was implemented when the validation loss could no longer improve after 400 iterations [55, 56]. The model performance was evaluated using k-fold cross-validation (k = 5) to ensure an unbiased assessment of the standard ML and MLP models, as described previously [37].

To assess the ability of the optimized models to identify all positive instances and avoid false negatives, the recall score (i.e. sensitivity or true positive rate) was estimated as the ratio of correctly age-classified An. funestus to the total number of An. funestus in the respective age category in the dataset. Moreover, to measure the ability of the models to avoid false positives, the precision score (i.e. the positive predictive value) was estimated as the ratio of correctly age-classified An. funestus to the total number of predicted positive instances of the respective age categories. Lastly, we calculated the F1 score, which balances both precision and recall scores by giving equal weight to both measures. This score provides a single value that represents the overall performance of the model in terms of its ability to correctly classify positive and negative cases. A higher F1 score signifies better model performance, where a maximum value of 1 represents flawless precision and recall.

Results

Predicting An. funestus age classes using standard ML models

In the initial comparison of standard ML models, XGBoost emerged as the best classifier with the highest prediction accuracy and lowest standard deviation, achieving 84% accuracy (Fig. 1A). After optimizing the parameters, the XGBoost model was able to classify spectra that were previously unseen with an overall accuracy of 87%. It achieved accuracy of 89% and 84% for young (1–9 days old) and old (10–16 days old) An. funestus females, respectively (Fig. 1B). The recall scores (i.e. sensitivity or true positive rates) of this model were 0.89 and 0.84 for the young and old mosquitoes, respectively, while its precision scores (i.e. the positive predictive value) were 0.87 for both age categories (Table 1).

Fig. 1
figure 1

Machine learning prediction of An. funestus age classes. A Comparison of standard ML classifiers in predicting An. funestus age classes; KNN k-nearest neighbours, LR logistic regression, SVM support vector machine, RF random forest, XGBoost gradient boosting, MLP multilayer perceptron. B Confusion matrix for predicting the age class of An. funestus using XGBoost on an unseen dataset, results for the ML trained with all spectral features

Table 1 Precision, recall and F1 score of XGBoost and multi-layer perceptron (MLP) models for predicting age categories of An. funestus

From the initial XGBoost model, we identified the spectral features that were most important for the prediction. This analysis aimed to reduce the number of training features and enhance the accuracy of the model during retraining (Fig. 2A). When the XGBoost classifier was retrained with the top 100 features, the classification accuracy increased to 89%, correctly predicting young and old An. funestus females with 92% and 85% accuracy, respectively (Fig. 2B).

Fig. 2
figure 2

A Relative importance of XGBoost features that have the most influence in predicting the age classes of An. funestus. B Confusion matrix for predicting the age class of An. funestus using XGBoost on an unseen dataset; the results for the ML retrained with important features/wavenumbers (n = 100) identified by the initial XGBoost model

Prediction of An. funestus age classes using MLP models

We explored the possibility of improving the accuracy by training the MLP classifier using the important wavenumbers (n = 100) identified in the XGBoost predictions. As a result, the MLP achieved an improved accuracy of 94.5% in the unseen test data (Fig. 3A), correctly distinguishing between young and old An. funestus females with accuracy of 95% and 94%, respectively (Fig. 3B).

Fig. 3
figure 3

A MLP training and validation accuracy for An. funestus age classes as training time increases (epoch; number of iterations over the entire dataset during the training process, i.e. seconds/iterations). Confusion matrix for predicting the age class of An. funestus. B The results for the MLP trained with important features/wavenumbers (n = 100) identified by XGBoost. C The results for the MLP method trained with eight principal components

Lastly, in a previous study, we presented evidence that employing principal component analysis (PCA) with eight components effectively reduces the dimensionality of the spectral data [47]. This reduction in dimensionality not only preserved a substantial portion of the data variability but also mitigated overfitting while enhancing the signal-to-noise ratio. By utilizing a reduced set of features, we trained the MLP model to improve its predictive performance [47]. In the present study, when PCA was utilized to reduce the dimensionality of the spectra data, the MLP classifier achieved an overall accuracy of 93% for both young and old An. funestus mosquitoes (Fig. 3C).

Discussion

Anopheles funestus mosquitoes are currently the major vector of malaria transmission in Tanzania, accounting for over 80% of malaria transmission [17, 20,21,22]. Anopheles funestus tends to have better survival rates [25], and is generally a slow-growing mosquito, which adds to the challenge of studying its demographic characteristics and how these might influence pathogen transmission. Here, we present a rapid age-grading technique that has the potential to replace traditional methods such as ovarian dissections, which are time-consuming and challenging to apply on a large scale. Using 2084 spectral data points, we trained ML models that classify the epidemiologically relevant age groups of An. funestus mosquitoes reared from wild larvae using water from the same habitats, but under laboratory conditions. The models correctly distinguished between the young An. funestus females (1–9 days old) and the older ones (10–16 days old) based on the MIR spectra indicative of the varying biochemical composition of the mosquito cuticles [57]. While this was the first demonstration of the effectiveness of this technique for predicting the age of An. funestus mosquitoes, the approach of combining infrared spectra and ML models has been widely demonstrated for predicting different indicators, including age, blood meals, infection status and insecticide resistance profiles of other Anopheles species [46,47,48]. If validated on field-collected adults, these findings could be a step towards wider applications of this approach for malaria vector surveillance in settings with different vector species.

In settings such as rural south-eastern Tanzania where An. funestus is the dominant malaria vector [17, 20], it is particularly important that vector surveillance programs are expanded to include this vector species. Indeed, the successful demonstration of this technique on An. funestus, which is one of the most efficient and also most widespread malaria vectors in Africa [58], expands the utility range of this technique for a much broader application for malaria vector surveys in different parts of Africa.

One of the key concerns regarding previous applications of MIRS-ML-based approaches for entomological assessments is that, with the exception of some cases [46], these methods have been rarely validated for wild-caught malaria vectors in field settings. Here, An. funestus mosquitoes were collected as larvae from various villages and breeding habitats, to account for genetic variation, variation in larval food sources and microbiome, and to maintain some characteristics of the natural ecosystems. The success of this analysis and the high accuracy obtained may therefore be indicative of the potential of the approach for predicting key mosquito attributes in field settings. However, it is unknown whether specific climatic factors could influence the prediction and generalizability of the MIRS-ML approach. Future studies should therefore test the generalizability of this approach across different populations of wild mosquitoes.

This study classified mosquitoes only as young (1–9 days old) or old (10–16 days old) and did not attempt to classify them at specific chronological ages because the sample size was not large enough to test it. However, the chosen age classes represent the typical epidemiological distinction relevant to the transmission of malaria parasites, which, under standard climatic conditions, requires a vector to be at least 10 days old [32]. However, it may fail to capture variations in MIR spectra or the small biochemical changes that occur within a mosquito cuticle after each ageing day (such as chronological age from 1 up to 16) [45]. Moreover, it has been demonstrated that calibrating ML models based on physiological age (which considers key life cycle processes such as blood-feeding and oviposition) may be more useful than simply relying on chronological age classifications [38, 59]. In our study, mosquitoes were all sugar-fed, and therefore physiological age was not assessed. Future efforts should assess key differences in these approaches and evaluate models trained on biological age and chronological age to determine which ones are most practical and most generalizable. An obvious next step is therefore to investigate any correlations that might exist between the machine-classified age categories and the epidemiology of malaria in human populations.

To improve the classification accuracy of our model, the XGBoost feature importance was relied upon to reduce the number of spectral features from 1665 to 100. This dimensionality reduction significantly lowered the noise and redundant features in the MIR spectra data. The important features were mostly associated with proteins, with the most influential peak (1700 cm−1) being the band associated with the amide bond from proteins. The region around 3000 cm−1, which is also related to proteins, was also found to be important in the model prediction. This implies that the model is learning from protein-based biological traits that vary depending on the age of the mosquito [46]. Moreover, when PCA was used to reduce the dimensionality of the spectra from 1665 features to eight principal components [47], the prediction accuracy matched that of the MLP model trained with the top 100 biological features as identified from the XGBoost model. This suggests that ML models may perform better when trained with fewer features that explain more variation in the data, rather than many redundant features that introduce noise into the model. Moreover, as observed previously, reducing the dimensionality of the spectra data reduces the computational resources needed to train ML models [47].

Future research should investigate the effects of rearing wild An. funestus larvae in the insectary on the predictive accuracy of the MIRS-ML approach for mosquito age-classification, as this could impact the generalizability of the findings.

Conclusions

This study demonstrates the classification of adult female An. funestus into distinct and epidemiologically relevant age categories using a MIRS-ML approach. In conjunction with prior research conducted on other Anopheles mosquitoes, this study suggests that the applicability of this approach can be extended to evaluate various entomological attributes in An. funestus. The MIRS-ML approach proves to be quick and cost-effective, and has the potential to significantly enhance An. funestus surveillance efforts, thereby contributing valuable insights to national malaria control programs, particularly in resource-constrained settings where this vector is highly prevalent. Nonetheless, further research is needed to validate the MIRS-ML approach in field conditions, using adult An. funestus populations and other vector species within malaria-endemic communities, and to examine how the machine-classified age categories correlate with the epidemiological strata of malaria in human populations.

Availability of data and materials

The mid-infrared spectral datasets generated and analysed during the current study, as well as code for the analyses, is available at [GitHub].

Abbreviations

MIRS:

Mid-infrared spectroscopy

NIRS:

Near-infrared spectroscopy

PCR:

Polymerase chain reaction

ITNs:

Insecticide-treated nets

MIRS-ML:

Mid-infrared spectroscopy and machine learning

References

  1. WHO. World malaria report. 2022.

  2. Winzeler EA, Manary MJ. Drug resistance genomics of the antimalarial drug artemisinin. Genome Biol. 2014;15:544.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Mbengue A, Bhattacharjee S, Pandharkar T, Liu H, Estiu G, Stahelin RV, et al. A molecular mechanism of artemisinin resistance in Plasmodium falciparum malaria. Nature. 2015;520:683–7.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  4. Siddiqui G, Srivastava A, Russell AS, Creek DJ. Multi-omics based identification of specific biochemical changes associated with PfKelch13-mutant artemisinin-resistant Plasmodium falciparum. J Infect Dis. 2017;215:1435–44.

    Article  CAS  PubMed  Google Scholar 

  5. Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, et al. Insecticide resistance in mosquito vectors. Nature. 2003;423:136–7.

    Article  ADS  CAS  PubMed  Google Scholar 

  6. Sokhna C, Ndiath MO, Rogier C. The changes in mosquito vector behaviour and the emerging resistance to insecticides will challenge the decline of malaria. Clin Microbiol Infect. 2013;19:902–7.

    Article  CAS  PubMed  Google Scholar 

  7. Bosco AB, Anderson K, Gresty K, Prosser C, Smith D, Nankabirwa JI, et al. Molecular surveillance reveals the presence of pfhrp2 and pfhrp3 gene deletions in Plasmodium falciparum parasite populations in Uganda, 2017–2019. Malar J. 2020;19:300.

    Article  CAS  PubMed Central  Google Scholar 

  8. Thomson R, Beshir KB, Cunningham J, Baiden F, Bharmal J, Bruxvoort KJ, et al. pfhrp2 and pfhrp3 gene deletions that affect malaria rapid diagnostic tests for Plasmodium falciparum: analysis of archived blood samples from 3 African countries. J Infect Dis. 2019;220:1444–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Funwei R, Nderu D, Nguetse CN, Thomas BN, Falade CO, Velavan TP, et al. Molecular surveillance of pfhrp2 and pfhrp3 genes deletion in Plasmodium falciparum isolates and the implications for rapid diagnostic tests in Nigeria. Acta Trop. 2019;196:121–5.

    Article  CAS  PubMed  Google Scholar 

  10. Berzosa P, González V, Taravillo L, Mayor A, Romay-Barja M, García L, et al. First evidence of the deletion in the pfhrp2 and pfhrp3 genes in Plasmodium falciparum from Equatorial Guinea. Malar J. 2020;19:99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Agaba BB, Yeka A, Nsobya S, Arinaitwe E, Nankabirwa J, Opigo J, et al. Systematic review of the status of pfhrp2 and pfhrp3 gene deletion, approaches and methods used for its estimation and reporting in Plasmodium falciparum populations in Africa: review of published studies 2010–2019. Malar J. 2019;18:355.

    Article  PubMed  PubMed Central  Google Scholar 

  12. WHO. World Malaria report 2021. 2021.

  13. Parpia AS, Ndeffo-Mbah ML, Wenzel NS, Galvani AP. Effects of response to 2014–2015 Ebola outbreak on deaths from malaria, HIV/AIDS, and tuberculosis, West Africa. Emerg Infect Dis. 2016;22:433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sherrard-Smith E, Hogan AB, Hamlet A, Watson OJ, Whittaker C, Winskill P, et al. The potential public health consequences of COVID-19 on malaria in Africa. Nat Med. 2020;26:1411–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  16. Lwetoijera DW, Harris C, Kiware SS, Dongus S, Devine GJ, McCall PJ, et al. Increasing role of Anopheles funestus and Anopheles arabiensis in malaria transmission in the Kilombero Valley, Tanzania. Malar J. 2014;13:331.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Kaindoa EW, Matowo NS, Ngowo HS, Mkandawile G, Mmbando A, Finda M, et al. Interventions that effectively target Anopheles funestus mosquitoes could significantly improve control of persistent malaria transmission in south–eastern Tanzania. PLoS ONE. 2017;12:e0177807.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Djamouko-Djonkam L, Nkahe DL, Kopya E, Talipouo A, Ngadjeu CS, Doumbe-Belisse P, et al. Implication of Anopheles funestus in malaria transmission in the city of Yaoundé, Cameroon. EDP Sciences. 2020;27:10.

  19. Ogola EO, Fillinger U, Ondiba IM, Villinger J, Masiga DK, Torto B, et al. Insights into malaria transmission among Anopheles funestus mosquitoes, Kenya. Parasit Vectors. 2018;11:577.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Mapua SA, Hape EE, Kihonda J, Bwanary H, Kifungo K, Kilalangongono M, et al. Persistently high proportions of plasmodium-infected Anopheles funestus mosquitoes in two villages in the Kilombero valley, South-Eastern Tanzania. Parasite Epidemiol Control. 2022;18:e00264.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Swai JK, Mmbando AS, Ngowo HS, Odufuwa OG, Finda MF, Mponzi W, et al. Protecting migratory farmers in rural Tanzania using eave ribbons treated with the spatial mosquito repellent, transfluthrin. Malar J. 2019;18:414.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Matowo NS, Kulkarni MA, Messenger LA, Jumanne M, Martin J, Mallya E, et al. Differential impact of dual-active ingredient long-lasting insecticidal nets on primary malaria vectors: a secondary analysis of a 3-year, single-blind, cluster-randomised controlled trial in rural Tanzania. Lancet Planet Heal. 2023;7:e370–80.

    Article  Google Scholar 

  23. Takken W, Verhulst NO. Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2013;58:433–53.

    Article  CAS  PubMed  Google Scholar 

  24. Limwagu AJ, Kaindoa EW, Ngowo HS, Hape E, Finda M, Mkandawile G, et al. Using a miniaturized double-net trap (DN-Mini) to assess relationships between indoor–outdoor biting preferences and physiological ages of two malaria vectors, Anopheles arabiensis and Anopheles funestus. Malar J. 2019;18:282.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Midega JT, Mbogo CM, Mwnambi H, Wilson MD, Ojwang G, Mwangangi JM, et al. Estimating dispersal and survival of Anopheles gambiae and Anopheles funestus along the Kenyan coast by using mark-release-recapture methods. J Med Entomol. 2007;44:923–9.

    Article  PubMed  Google Scholar 

  26. Coetzee M, Koekemoer LL. Molecular systematics and insecticide resistance in the major African malaria vector Anopheles funestus. Annu Rev Entomol. 2013;58:393–412.

    Article  CAS  PubMed  Google Scholar 

  27. Nambunga IH, Ngowo HS, Mapua SA, Hape EE, Msugupakulya BJ, Msaky DS, et al. Aquatic habitats of the malaria vector Anopheles funestus in rural south-eastern Tanzania. Malar J. 2020;19:219.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Detinova TS. Age-grouping methods in Diptera of medical importance with special reference to some vectors of malaria. Monogr Ser World Health Organ. Geneva: World Health Organization; 1962; 47:13–191.

  29. Rao V. On gonotrophic discordance among certain Indian Anopheles. Indian J Malariol. 1947;1:43–50.

    Google Scholar 

  30. Derek Charlwood J, Nenhep S, Sovannaroth S, Morgan JC, Hemingway J, Chitnis N, et al. ‘Nature or nurture’: survival rate, oviposition interval, and possible gonotrophic discordance among South East Asian anophelines. Malar J. 2016;15:356.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Charlwood JD, Tomás EVE, Andegiorgish AK, Mihreteab S, LeClair C. ‘We like it wet’: a comparison between dissection techniques for the assessment of parity in Anopheles arabiensis and determination of sac stage in mosquitoes alive or dead on collection. PeerJ. 2018;6:e5155.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Ohm JR, Baldini F, Barreaux P, Lefevre T, Lynch PA, Suh E, et al. Rethinking the extrinsic incubation period of malaria parasites. Parasit Vectors. 2018;11:178.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Chen CS, Mulla MS, March RB, Chaney JD. Cuticular hydrocarbon patterns in Culex quinquefasciatus as influenced by age, sex, and geography. Bull Soc Vector Ecol. 1990;15:129–39.

    Google Scholar 

  34. Cook PE, Hugo LE, Iturbe-Ormaetxe I, Williams CR, Chenoweth SF, Ritchie SA, et al. Predicting the age of mosquitoes using transcriptional profiles. Nat Protoc. 2007;2:2796–806.

    Article  CAS  PubMed  Google Scholar 

  35. Hugo LE, Cook PE, Johnson PH, Rapley LP, Kay BH, Ryan PA, et al. Field validation of a transcriptional assay for the prediction of age of uncaged Aedes aegypti mosquitoes in northern Australia. PLoS Negl Trop Dis. 2010;4:e608.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Wang M-H, Marinotti O, James AA, Walker E, Githure J, Yan G. Genome-wide patterns of gene expression during aging in the African malaria vector Anopheles gambiae. PLoS ONE. 2010;5:e13359.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  37. Burns DA, Ciurczak EW. Handbook of near-infrared analysis. Boca Raton: CRC Press; 2008.

    Google Scholar 

  38. Mayagaya VS, Michel K, Benedict MQ, Killeen GF, Wirtz RA, Ferguson HM, et al. Non-destructive determination of age and species of Anopheles gambiae sl using near-infrared spectroscopy. Am J Trop Med Hyg. 2009;81:622.

    Article  CAS  PubMed  Google Scholar 

  39. Ong OTW, Kho EA, Esperança PM, Freebairn C, Dowell FE, Devine GJ, et al. Ability of near-infrared spectroscopy and chemometrics to predict the age of mosquitoes reared under different conditions. Parasit Vectors. 2020;13:160.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Sikulu-Lord MT, Maia MF, Milali MP, Henry M, Mkandawile G, Kho EA, et al. Rapid and non-destructive detection and identification of two strains of Wolbachia in Aedes aegypti by near-infrared spectroscopy. PLoS Negl Trop Dis. 2016;10:e0004759.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lambert B, Sikulu-Lord MT, Mayagaya VS, Devine G, Dowell F, Churcher TS. Monitoring the age of mosquito populations using near-infrared spectroscopy. Sci Rep. 2018;8:5274.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  42. Sikulu-Lord MT, Devine GJ, Hugo LE, Dowell FE. First report on the application of near-infrared spectroscopy to predict the age of Aedes albopictus Skuse. Sci Rep. 2018;8:9590.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  43. Sikulu M, Killeen GF, Hugo LE, Ryan PA, Dowell KM, Wirtz RA, et al. Near-infrared spectroscopy as a complementary age grading and species identification tool for African malaria vectors. Parasit Vectors. 2010;3:49.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Sikulu MT, Majambere S, Khatib BO, Ali AS, Hugo LE, Dowell FE. Using a near-infrared spectrometer to estimate the age of anopheles mosquitoes exposed to pyrethroids. PLoS ONE. 2014;9:e90657.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  45. Gonzalez-Jimenez M, Babayan SA, Khazaeli P, Doyle M, Walton F, Reedy E, et al. Prediction of malaria mosquito species and population age structure using mid-infrared spectroscopy and supervised machine learning. Wellcome Open Res. 2019;4:76.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Siria DJ, Sanou R, Mitton J, Mwanga EP, Niang A, Sare I, et al. Rapid age-grading and species identification of natural mosquitoes for malaria surveillance. Nat Commun. 2022;13:1501.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  47. Mwanga EP, Siria DJ, Mitton J, Mshani IH, González-Jiménez M, Selvaraj P, et al. Using transfer learning and dimensionality reduction techniques to improve generalisability of machine-learning predictions of mosquito ages from mid-infrared spectra. BMC Bioinform. 2023;24:11.

    Article  Google Scholar 

  48. Mwanga EP, Mapua SA, Siria DJ, Ngowo HS, Nangacha F, Mgando J, et al. Using mid-infrared spectroscopy and supervised machine-learning to identify vertebrate blood meals in the malaria vector Anopheles arabiensis. Malar J. 2019;18:187.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Mwanga EP, Minja EG, Mrimi E, Jiménez MG, Swai JK, Abbasi S, et al. Detection of malaria parasites in dried human blood spots using mid-infrared spectroscopy and logistic regression analysis. Malar J. 2019;18:341.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Coetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:70.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66:804–11.

    Article  CAS  PubMed  Google Scholar 

  52. Pedregosa F, Gramfort A, Michel V, Thirion B, Grisel O, Blondel M, et al. Scikit-learn : machine learning in Python. J Mach Learn Res. 2011;12:2825–30.

    MathSciNet  Google Scholar 

  53. Abadi M, Barham P, Chen J, Chen Z, Davis A, Dean J, et al. TensorFlow: a system for large-scale machine learning. In: Proc 12th USENIX Symp Oper Syst Des Implementation, OSDI 2016. 2016.

  54. Chollet F. Keras: the Python deep learning library. Kerasio. 2015. http://keras.io.

  55. Nitish S, Geoffrey H, Alex K, Ilya S, Ruslan S. Dropout: a simple way to prevent neural networks from overfitting. J Mach Learn Res. 2014;15:1929–58.

    MathSciNet  Google Scholar 

  56. Prechelt L. Early stopping—but when? In: Montavon G, Orr GB, Müller KR, editors. Neural networks: tricks of the trade, vol. 7700. Lecture Notes in Computer Science. Berlin, Heidelberg: Springer; 2012.

    Google Scholar 

  57. Suarez E, Nguyen HP, Ortiz IP, Lee KJ, Kim SB, Krzywinski J, et al. Matrix-assisted laser desorption/ionization-mass spectrometry of cuticular lipid profiles can differentiate sex, age, and mating status of Anopheles gambiae mosquitoes. Anal Chim Acta. 2011;706:157–63.

    Article  CAS  PubMed  Google Scholar 

  58. Dia I, Guelbeogo MW, Ayala D. Advances and perspectives in the study of the malaria mosquito Anopheles funestus. London: Intech; 2013.

    Book  Google Scholar 

  59. Ntamatungiro AJ, Mayagaya VS, Rieben S, Moore SJ, Dowell FE, Maia MF. The influence of physiological status on age prediction of Anopheles arabiensis using near infra-red spectroscopy. Parasit Vectors. 2013;6:298.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors sincerely appreciate all field technicians who assisted in the collection of wild An. funestus larvae, as well as the rearing and handling of adult mosquitoes. The authors also express their gratitude to the administration team for their continuous administrative assistance. Additionally, we are grateful to the community members and local government officials in the districts of Ulanga and Kilombero for their unwavering support throughout this study.

Funding

This study was supported by a Howard Hughes Medical Institute (HHMI)-Gates International Research Scholarship (Grant No. OPP1099295) awarded to FOO and the Medical Research Council (MRC) [MR/P025501/1] awarded to FB. EPM was supported by the Wellcome Trust Masters Fellowship in Tropical Medicine & Hygiene (Grant No. 214643/Z/18/Z). FB is supported by the Academy Medical Sciences Springboard Award (ref: SBF007\100094). SAB is supported by the Bill and Melinda Gates Foundation (INV-030025) and Royal Society (ICA/R1/191238).

Author information

Authors and Affiliations

Authors

Contributions

EPM, DJS, SB, FB, and FOO conceived the study. EPM, SA, FOO, and DJS developed the study protocol. DJS collected the data, EPM carried out data analysis and ML training. EPM wrote the manuscript. EPM, DJS, SHM, IHM, MGJ, KW, SB, FB, and FOO reviewed and edited drafts of the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Emmanuel P. Mwanga.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for this study was obtained from the Ifakara Health Institute Institutional Review Board (Ref. IHI/IRB/EXT/No: 005-2018), and from the Medical Research Coordinating Committee (MRCC) at the National Institute of Medical Research (NIMR), Ref: NIMR/HQ/R.8c/Vol. II/880.

Consent for publication

All authors have read and approved the final manuscript. Permission to publish this work was obtained from the National Institute of Medical Research (NIMR), Ref: NIMR/HQ/P.12 VOL.XXXVI/48.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mwanga, E.P., Siria, D.J., Mshani, I.H. et al. Rapid classification of epidemiologically relevant age categories of the malaria vector, Anopheles funestus. Parasites Vectors 17, 143 (2024). https://doi.org/10.1186/s13071-024-06209-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-024-06209-5

Keywords