A user-friendly software to easily count Anopheles egg batches
© Mollahosseini et al.; licensee BioMed Central Ltd. 2012
Received: 14 March 2012
Accepted: 10 May 2012
Published: 19 June 2012
Studies on malaria vector ecology and development/evaluation of vector control strategies often require measures of mosquito life history traits. Assessing the fecundity of malaria vectors can be carried out by counting eggs laid by Anopheles females. However, manually counting the eggs is time consuming, tedious, and error prone.
In this paper we present a newly developed software for high precision automatic egg counting. The software written in the Java programming language proposes a user-friendly interface and a complete online manual. It allows the inspection of results by the operator and includes proper tools for manual corrections. The user can in fact correct any details on the acquired results by a mouse click. Time saving is significant and errors due to loss of concentration are avoided.
The software was tested over 16 randomly chosen images from 2 different experiments. The results show that the proposed automatic method produces results that are close to the ground truth.
The proposed approaches demonstrated a very high level of robustness. The adoption of the proposed software package will save many hours of labor to the bench scientist. The software needs no particular configuration and is freely available for download on: http://w3.ualg.pt/∼hshah/eggcounter/.
Understanding the ecology and evolution of malaria vector species and populations is a key factor in controlling the disease they carry . In several contexts, the estimation of the fitness of Anopheles mosquitoes is required. For instance, this may help to understand adaptation of different populations to given environmental conditions [2–4] or the opposite [5–7], interactions between the vectors and the parasites [8–12] or the vertebrate hosts . More directly related to vector control, fitness measures can help deciphering the effect of insecticides , including new classes of insect growth regulators , effect of insecticide resistance [16, 17], genetic manipulation of vector populations to make them resistant to parasites [18, 19], or to impede their reproductive success. The fecundity is a key variable to assess the fitness of populations. In mosquitoes, it can be measured at several stages, but the number of eggs laid by females is often considered as a good estimate. Depending on the case study, the experiments may involve hundreds of females, possibly over many gonotrophic cycles and generations. As female Anopheles can lay over one hundred eggs every 2-3 days, one experiment only, often involves counting thousands of eggs [20, 21]. Such a task is obviously tedious, time-consuming and prone to errors. These facts make the identification and counting of eggs a natural task for automation.
As software capable of correctly counting Anopheles eggs laid on filter paper was not available, we developed a fully automatic software together with its own user-friendly interface. Egg-Counter v1.0 has the capacity of counting individual eggs as well as eggs laid in piles. At first, it counts the eggs laid on the paper, excluding the debris. Then, the user is able to modify any details of the presented result by means of a provided user-friendly interface. The software was tested by comparing the automatic and manually corrected results with the number of manually counted eggs, and a precision of 98% was obtained. Egg-Counter is fully written in Java and, therefore, completely platform-independent; it can be downloaded, installed and used within a couple of minutes.
The rest of this paper is organized as follows: In Methods, the proposed method for processing the images is elaborated. A brief description of the software provided and its features are introduced in Software. Experimental results on 16 different images, taken in different experiments are presented in Results and Discussion.
Where T auto is the value of T that makes C B (T) + C F (T) maximum, p(g) is the probability of intensity level g and p(T) is the cumulative probability function. Unlike the original method , here controls the maximum value of threshold in order to remove high exposure in the images caused by flash where empirically . Figure 2(c) shows the resulting binary map of the sample image. As it is shown eggs and debris are properly extracted, after the automatic threshold is found.
In this step eggs and debris are detected by extracting eight-connected neighborhood components. Putting it simply, connected component labeling is the process of assigning labels to the foreground elements in such a way that adjacent foreground elements are assigned to the same label . Here, adjacent means 8-adjacent pixels neighborhood.
Egg Size Estimation
where is the centroid of the object and represents whole set of pixel coordinates inside the object. The value of the eccentricity falls in the interval [0,1]. A value of one indicates a perfect circle and zero indicates a line segment. In this step, the objects with high eccentricity, e 1<e<e 2, are selected as eggs. Values of e 1e 2can be set from 0.95 up to 0.99. Experimental results show that more than 80% of the components containing one egg are in this range. Finally, the mode of the size of the candidate objects is chosen as the estimated size for one egg.
In this step, the objects that have at least X pixels with C> C E in their radii r are considered as debris. Parameter C E is 0.15 by default, but the application allows the user to tune it, resulting in a change of the detection sensitivity. Because eggs have a dark color and the paper is white, the rate of true negatives (considering non-debris areas as debris) is negligible and it only happens for the eggs that are inside the yellow brownish area.
Egg Count and Noise Removal
where A(i) is the area of the object, S is the estimated size of one egg, ∥.∥ is the round function and C is the accumulative ratio which should be less than one. The accumulative ratio is necessary due to the fact that regions in the intersection of the eggs inside a pile are shadowed by the eggs and are usually considered as a connected object in the binary map. Therefore, multiplying by this ratio allows better estimation of the number of eggs.
It should be mentioned that thresholding may produce noise, especially in the grid lines of the paper. There are many methods for noise reduction such as low-pass filters and other smoothing operators. However, these methods highly depend on the size of the filter. Here components smaller than half of the estimated size of one egg, are considered as noise in the Equation (10), and are removed.
Results and Discussion
Results of 16 different images
702 / 29
711 / 9
342 / 20
366 / 18
608 / 40#
399 / 620
620 / 35
391 / 21
277 / 400
400 / 35
372 / 21
314 / 363
363 / 25
219 / 20
219 / 18
460 / 27
455 / 20
086 / 15
89 / 18
114 / 20
119 / 18
280 / 21
272 / 20
352 / 25
360 / 20
480 / 25
489 / 20
022 / 06
22 / 15
297 / 20
297 / 115
331 / 21
331 / 110
311 / 20
311 / 110
As it is shown in the error rate column, the precision of the software after a first correction is on average more than 98%, whereas the time saved is enormous. For image number 3, the automatic processing was very poor and the user had to manually set the sensitivity parameters. It should be mentioned that a proper sensitivity can be easily found by just adjusting a sliding bar. Manually counting all eggs of an image containing more than 500 eggs can take up to 30 minutes but, the same image, can be processed by Egg-Counter v1.0 in only a few seconds.
The results show that the proposed automatic method produces results that are close to the ground truth, and by making a few adjustments, mainly in the regions of debris and the paper’s grid lines, these results get even closer to the ground truth. Multiplying the time gained by the number of images and experiments, results in economizing an extremely large amount of time. The sensitivity value of the test results for assessing the eggs is around 0.945 in the fully automatic use of the software and is higher after a few mouse clicks. The specificity value, in this case, makes no sense since the true negatives are not actually known. Elimination of thousands of noise pixels during different stages of the analysis does not allow a reliable calculation of the specificity of the software. Positive predictive value of the test is set around 0.921 for the full automatic version.
The proposed algorithm works well, if eggs are clearly darker than the background paper, even if they are brown or another dark color. To introduce better and higher contrast and a wider color difference between the paper and laid eggs, we suggest using white or semi-white paper. In addition, it is highly recommended to take images with a camera of at least 5M pixels, using the lowest ISO and maximum shutter speed and avoiding re-sizing the images afterwards. Depending on your Laboratory lightning conditions, you might need to use a flash.
In this article, a promising new software package for automatically counting Anopheles eggs was presented. A very intuitive graphical user interface is available, that allows correction of the detection results, if needed. Additionally, in the parameters section of the interface it is possible to control the sensitivity to egg detection. Nevertheless, the proposed approaches involved in the processing of such images demonstrated a very high level of robustness, as it was shown by the results. The adoption of the proposed software package will save many hours of labor for the bench scientist.
Ali Mollahosseini, António dos Anjos and Hamid Reza Shahbazkia contributed with the proposed methods and software and Marie Rossignol, Cédric Pennetier, Anna Cohuet and Fabrice Chandre for testing and providing the necessary material. All authors read and approved the final version of the manuscript.
We would like to acknowledge IRD MIVEGEC Montpellier and DEEI UALG for having provided the necessary backup. We would like to thank Prof. Chris Arme, the editor of Parasites & Vectors, and the anonymous reviewers for their insightful and constructive comments.
- Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, Gimnig JE, Fish D, Killeen GF: Ecology: A prerequisite for malaria elimination and eradication. PLoS Med. 2010, 7 (8): e1000303-10.1371/journal.pmed.1000303.PubMed CentralView ArticlePubMedGoogle Scholar
- Yaro AS, Dao A, Adamou A, Crawford JE, Traoré SF, Touré AM, Gwadz R, Lehmann T: Reproductive output of female Anopheles gambiae (Diptera: Culicidae): comparison of molecular forms. J Med Entomol. 2006, 43 (5): 833-839. 10.1603/0022-2585(2006)43[833:ROOFAG]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Sharma A, Parasher H, Singh OP, Adak T: Species B of Anopheles culicifacies (Diptera: Culicidae) is reproductively less fit than species A and C of the complex. Acta Tropica. 2009, 112 (3): 316-319. 10.1016/j.actatropica.2009.08.006.View ArticlePubMedGoogle Scholar
- Hogg JC, Thomson MC, Hurd H: Comparative fecundity and associated factors for two sibling species of the Anopheles gambiae complex occurring sympatrically in The Gambia. Med Vet Entomol. 1996, 10 (4): 385-391. 10.1111/j.1365-2915.1996.tb00761.x.View ArticlePubMedGoogle Scholar
- Manda H, Gouagna LC, Foster WA, Jackson RR, Beier JC, Githure JI, Hassanali A: Effect of discriminative plant-sugar feeding on the survival and fecundity of Anopheles gambiae. Malaria Journal. 2007, 6: 113-124. 10.1186/1475-2875-6-113.PubMed CentralView ArticlePubMedGoogle Scholar
- Gary RE, Foster WA: Effects of available sugar on the reproductive fitness and vectorial capacity of the malaria vector Anopheles gambiae (Diptera: Culicidae). J Med Entomol. 2001, 38: 22-28. 10.1603/0022-2585-38.1.22.View ArticlePubMedGoogle Scholar
- Pfaehler O, Oulo DO, Gouagna LC, Githure J, Guerin PM: Influence of soil quality in the larval habitat on development of Anopheles gambiae Giles. J Vector Ecol. 2006, 31 (2): 400-405. 10.3376/1081-1710(2006)31[400:IOSQIT]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Yaro AS, Touré AM, Guindo A, Coulibaly MB, Dao A, Diallo M, Traoré SF: Reproductive success in Anopheles arabiensis and the M and S molecular forms of Anopheles gambiae: Do natural sporozoite infection and body size matter?. Acta Tropica. 2012, 122: 87-93. 10.1016/j.actatropica.2011.12.005.View ArticlePubMedGoogle Scholar
- Hogg JC, Hurd H: The effects of natural Plasmodium falciparum infection on the fecundity and mortality of Anopheles gambiae s. l. in north east Tanzania. Parasitology. 1997, 114 (Pt4): 325-331.View ArticlePubMedGoogle Scholar
- Hogg JC, Hurd H: Malaria-induced reduction of fecundity during the first gonotrophic cycle of Anopheles stephensi mosquitoes. Med Vet Entomol. 1995, 9 (2): 176-180. 10.1111/j.1365-2915.1995.tb00175.x.View ArticlePubMedGoogle Scholar
- Ferguson HM, Read AF: Why is the effect of malaria parasites on mosquito survival still unresolved?. Trends in Parasitology. 2002, 18 (6): 256-261. 10.1016/S1471-4922(02)02281-X.View ArticlePubMedGoogle Scholar
- Ahmed AM, Hurd H: Immune stimulation and malaria infection impose reproductive costs in Anopheles gambiae via follicular apoptosis. Microbes and Infection. 2006, 8 (2): 308-315. 10.1016/j.micinf.2005.06.026.View ArticlePubMedGoogle Scholar
- Lyimo IN, Keegan SP, Ranford-Cartwright LC, Ferguson HM: The impact of uniform and mixed species blood meals on the fitness of the mosquito vector Anopheles gambiae s.s: does a specialist pay for diversifying its host species diet?. Journal of Evolutionary Biology. 2012, 25 (3): 452-460. 10.1111/j.1420-9101.2011.02442.x.View ArticlePubMedGoogle Scholar
- Mosqueira B, Duchon S, Chandre F, Hougard JM, Carnevale P, Mas-Coma S: Efficacy of an insecticide paint against insecticide-susceptible and resistant mosquitoes - part 1: laboratory evaluation. Malaria Journal. 2010, 9: 340-346. 10.1186/1475-2875-9-340.PubMed CentralView ArticlePubMedGoogle Scholar
- Dhadialla TS, Carlson GR, Le DP: New insecticides with ecdysteroidal and juvenile hormone activity. Annu Rev Entomol. 1998, 43: 545-569. 10.1146/annurev.ento.43.1.545.View ArticlePubMedGoogle Scholar
- Okoye PN, Brooke BD, Hunt RH, Coetzee M: Relative developmental and reproductive fitness associated with pyrethroid resistance in the major southern African malaria vector, Anopheles funestus. Bulletin of Entomological Research. 2007, 97 (6): 599-605.View ArticlePubMedGoogle Scholar
- Djogbénou L, Labbé P, Chandre F, Pasteur N, Weill M: Ace-1 duplication in Anopheles gambiae: a challenge for malaria control. Malaria Journal. 2009, 8: 70-76. 10.1186/1475-2875-8-70.PubMed CentralView ArticlePubMedGoogle Scholar
- Li C, Marrelli MT, Yan G, Jacobs-Lorena M: Fitness of transgenic Anopheles stephensi mosquitoes expressing the SM1 peptide under the control of a vitellogenin promoter. J Hered. 2008, 99 (3): 275-282. 10.1093/jhered/esn004.PubMed CentralView ArticlePubMedGoogle Scholar
- Lambrechts L, Koella JC, Boëte C: Can transgenic mosquitoes afford the fitness cost?. Trends in Parasitology. 2008, 24: 4-7. 10.1016/j.pt.2007.09.009.View ArticlePubMedGoogle Scholar
- Munhenga G, Brooke BD, Chirwa TF, Hunt RH, Coetzee M, Govender D, Koekemoer LL: Evaluating the potential of the sterile insect technique for malaria control: relative fitness and mating compatibility between laboratory colonized and a wild population of Anopheles arabiensis from the Kruger National Park, South Africa. Parasites & Vectors. 2011, 4: 208-11. 10.1186/1756-3305-4-208.View ArticleGoogle Scholar
- Kweka EJ, Owino EA, Mwang’onde BJ, Mahande AM, Nyindo M, Mosha F: The role of cow urine in the oviposition site preference of culicine and Anopheles mosquitoes. Parasites & Vectors. 2011, 4: 184-190. 10.1186/1756-3305-4-184.View ArticleGoogle Scholar
- Yen JC, Chang FJ, Chang S: A new criterion for automatic multilevel thresholding. Image Processing, IEEE Transactions on. 1995, 4 (3): 370-378. 10.1109/83.366472.View ArticleGoogle Scholar
- González R, Woods R: 2001, Addison-Wesley Longman Publishing Co., Inc., Boston, MA, USAGoogle Scholar
- Burger W, Burge MJ: Digital Image Processing: an algorithmic introduction using Java. 2008, Springer, New YorkView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.