Abstract
Background: Several mosquito population suppression strategies based on the rearing and release of sterile males have provided promising results. However, the lack of an efficient male selection method has hampered the expansion of these approaches into large-scale operational programmes. Currently, most of these programmes targeting Aedes mosquitoes rely on sorting methods based on the sexual size dimorphism (SSD) at the pupal stage. The currently available sorting methods have not been developed based on biometric analysis, and there is therefore potential for improvement. We applied an automated pupal size estimator developed by Grupo Tragsa with laboratory samples of Anopheles arabiensis, Aedes albopictus, Ae. polynesiensis, and three strains of Ae. aegypti. The frequency distribution of the pupal size was analyzed. We propose a general model for the analysis of the frequency distribution of mosquito pupae in the context of SSD-sorting methods, which is based on a Gaussian mixture distribution functions, thus making possible the analysis of performance (% males recovery) and purity (% males on the sorted sample). Results: For the three Aedes species, the distribution of the pupae size can be modeled by a mixture of two Gaussian distribution functions and the proposed model fitted the experimental data. For a given population, each size threshold is linked to a specific outcome of male recovery. Two dimensionless parameters that measure the suitability for SSD-based sorting of a specific batch of pupae are provided. The optimal sorting results are predicted for the highest values of SSD and lowest values of intra-batch variance. Rearing conditions have a strong influence in the performance of the SSD-sorting methods and non-standard rearing can lead to increase pupae size heterogeneity. Conclusions: Sex sorting of pupae based on size dimorphism can be achieved with a high performance (% males recovery) and a reasonably high purity (% males on the sorted sample) for the different Aedes species and strains. The purity and performance of a sex sorting operation in the tested Aedes species are linked parameters whose relation can be modeled. The conclusions of this analysis are applicable to all the existing SSD-sorting methods. The efficiency of the SSDsorting methods can be improved by reducing the heterogeneity of pupae size within rearing containers. The heterogeneity between batches does not strongly affect the quality of the sex sorting, as long as a specific separation threshold is not pre-set before the sorting process. For new developments, we recommend using adaptive and precise threshold selection methods applied individually to each batch or to a mix of batches. Adaptive and precise thresholds will allow the sex-sorting of mixed batches in operational conditions maintaining the target purity at the cost of a reduction in performance. We also recommend a strategy whereby an acceptable level of purity is pre-selected and remains constant across the different batches of pupae while the performance varies from batch to batch to fit with the desired purity.