Mosquitoes of the Maculipennis complex in Northern Italy

Characterization of environmental drivers influencing the abundance of Anopheles maculipennis complex in Northern Italy

BACKGROUND

In Italy, malaria was endemic until the 1970s, when it was declared eradicated by WHO. Nowadays, with the persistence of competent mosquito populations, the effect of climate change, and increased possibility of importing malaria parasites from endemic counties due to growing migration, a malaria resurgence in Italy has become more likely. Hence, enhancing the understanding of the current distribution of the Anopheles maculipennis complex and the factors that influence the presence of this malaria vector is crucial, especially in Northern Italy, characterised by a high density of both human population and livestock.

METHODS

To assess the presence and abundance of malaria vectors, a 4-year field survey in the plain areas of Lombardy and Emilia-Romagna region in Italy was conducted. Every sampling point was characterised in space by the land use in a 500-m radius and in time considering meteorological data collected in the short and long time periods before sampling. We combined the results of a linear regression model with a random forest analysis to understand the relative importance of the investigated niche dimensions in determining Anopheles mosquito presence and abundance.

RESULTS

The estimated normalised variable importance indicates that rice fields were the most important land use class explaining the presence of Anopheles, followed by transitional woodlands and shrubland. Farm buildings were the third variable in terms of importance, likely because of the presence of animal shelters, followed by urbanised land. The two most important meteorological variables influencing the abundance of Anopheles in our study area were mean temperature in the 24 h before the sampling date and the sum of degree-days with temperature between 18 °C and 30 °C in the 14 days before the sampling date.

CONCLUSIONS

The results obtained in this study could be helpful in predicting the risk of autochthonous malaria transmission, based on local information on land cover classes that might facilitate the presence of malaria vectors and presence of short- and medium-term meteorological conditions favourable to mosquito development and activity. The results can support the design of vector control measures through environmental management.

Trends in mosquito species distribution modeling: insights for vector surveillance and disease control

Species distribution modeling (SDM) has become an increasingly common approach to explore questions about ecology, geography, outbreak risk, and global change as they relate to infectious disease vectors. Here, we conducted a systematic review of the scientific literature, screening 563 abstracts and identifying 204 studies that used SDMs to produce distribution estimates for mosquito species. While the number of studies employing SDM methods has increased markedly over the past decade, the overwhelming majority used a single method (maximum entropy modeling; MaxEnt) and focused on human infectious disease vectors or their close relatives. The majority of regional models were developed for areas in Africa and Asia, while more localized modeling efforts were most common for North America and Europe. Findings from this study highlight gaps in taxonomic, geographic, and methodological foci of current SDM literature for mosquitoes that can guide future efforts to study the geography of mosquito-borne disease risk.

Ecological Niche Modelling Approaches: Challenges and Applications in Vector-Borne Diseases

Vector-borne diseases (VBDs) pose a major threat to human and animal health, with more than 80% of the global population being at risk of acquiring at least one major VBD. Being profoundly affected by the ongoing climate change and anthropogenic disturbances, modelling approaches become an essential tool to assess and compare multiple scenarios (past, present and future), and further the geographic risk of transmission of VBDs. Ecological niche modelling (ENM) is rapidly becoming the gold-standard method for this task. The purpose of this overview is to provide an insight of the use of ENM to assess the geographic risk of transmission of VBDs. We have summarised some fundamental concepts and common approaches to ENM of VBDS, and then focused with a critical view on a number of crucial issues which are often disregarded when modelling the niches of VBDs. Furthermore, we have briefly presented what we consider the most relevant uses of ENM when dealing with VBDs. Niche modelling of VBDs is far from being simple, and there is still a long way to improve. Therefore, this overview is expected to be a useful benchmark for niche modelling of VBDs in future research.

Unsupervised machine learning and geometric morphometrics as tools for the identification of inter and intraspecific variations in the Anopheles Maculipennis complex

Geometric morphometric analysis was combined with two different unsupervised machine learning algorithms, UMAP and HDBSCAN, to visualize morphological differences in wing shape among and within four Anopheles sibling species (An. atroparvus, An. melanoon, An. maculipennis s.s. and An. daciae sp. inq.) of the Maculipennis complex in Northern Italy. Specifically, we evaluated: 1) wing shape variation among and within species; 2) the consistencies between groups of An. maculipennis s.s. and An. daciae sp. inq. identified based on COI sequences and wing shape variability; and 3) the spatial and temporal distribution of different morphotypes. UMAP detected at least 13 main patterns of variation in wing shape among the four analyzed species and mapped intraspecific morphological variations. The relationship between the most abundant COI haplotypes of An. daciae sp. inq. and shape ordination/variation was not significant. However, morphological variation within haplotypes was reported. HDBSCAN also recognized different clusters of morphotypes within An. daciae sp. inq. (12) and An. maculipennis s.s. (4). All morphotypes shared a similar pattern of variation in the subcostal vein, in the anal vein and in the radio-medial cross-vein of the wing. On the contrary, the marginal part of the wings remained unchanged in all clusters of both species. Any spatial-temporal significant difference was observed in the frequency of the identified morphotypes.  Our study demonstrated that machine learning algorithms are a useful tool combined with geometric morphometrics and suggest to deepen the analysis of inter and intra specific shape variability to evaluate evolutionary constrains related to wing functionality.

Comparative efficacy of BG-Sentinel 2 and CDC-like mosquito traps for monitoring potential malaria vectors in Europe

Background

Different trapping devices and attractants are used in the mosquito surveillance programs currently running in Europe. Most of these devices target vector species belonging to the genera Culex or Aedes, and no studies have yet evaluated the effectiveness of different trapping devices for the specific targeting of Anopheles mosquito species, which are potential vectors of malaria in Europe. This study aims to fill this gap in knowledge by comparing the performance of trapping methods that are commonly used in European mosquito surveillance programs for Culex and Aedes for the specific collection of adults of species of the Anopheles maculipennis complex.

Methods

The following combinations of traps and attractants were used: (i) BG-Sentinel 2 (BG trap) baited with a BG-Lure cartridge (BG + lure), (ii) BG trap baited with a BG-Lure cartridge and CO₂ (BG + lure + CO₂), (iii) Centers for Disease Control and Prevention-like trap (CDC trap) baited with CO₂ (CDC + CO₂), (iv) CDC trap used with light and baited with BG-Lure and CO₂ (CDC light + lure + CO₂). These combinations were compared in the field using a 4 × 4 Latin square study design. The trial was conducted in two sites in northeastern Italy in 2019. Anopheles species were identified morphologically and a sub-sample of An. maculipennis complex specimens were identified to species level by molecular analysis.

Results

Forty-eight collections were performed on 12 different trapping days at each site, and a total of 1721 An. maculipennis complex specimens were captured. The molecular analysis of a sub-sample comprising 254 specimens identified both Anopheles messeae/Anopheles daciae (n = 103) and Anopheles maculipennis sensu stricto (n = 8) at site 1, while at site 2 only An. messeae/An. daciae (n = 143) was found. The four trapping devices differed with respect to the number of An. messeae/An. daciae captured. More mosquitoes were caught by the BG trap when it was used with additional lures (i.e. BG + lure + CO₂) than without the attractant, CO₂ [ratio_{BG+lure vs BG+lure+CO2} = 0.206, 95% confidence interval (CI) 0.101–0.420, P < 0.0001], while no significant differences were observed between CDC + CO₂ and CDC light + lure + CO₂ (P = 0.321). The addition of CO₂ to BG + lure increased the ability of this combination to capture An. messeae/An. daciae by a factor of 4.85, and it also trapped more mosquitoes of other, non-target species (Culex pipiens, ratio_{BG+lure vs BG+lure+CO2} = 0.119, 95% CI 0.056–0.250, P < 0.0001; Ochlerotatus caspius, ratio_{BG+lure vs BG+lure+CO2} = 0.035, 95% CI 0.015–0.080, P < 0.0001).

Conclusions

Our results show that both the BG-Sentinel and CDC trap can be used to effectively sample An. messeae/An. daciae, but that the combination of the BG-Sentinel trap with the BG-Lure and CO₂ was the most effective means of achieving this. BG + lure + CO₂ is considered the best combination for the routine monitoring of host-seeking An. maculipennis complex species such as An. messeae/An. daciae. The BG-Sentinel and CDC traps have value as alternative methods to human landing catches and manual aspiration for the standardized monitoring of Anopheles species in Europe.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05285-9.

Updated occurrence and bionomics of potential malaria vectors in Europe: a systematic review (2000-2021)

Despite the eradication of malaria across most European countries in the 1960s and 1970s, the anopheline vectors are still present. Most of the malaria cases that have been reported in Europe up to the present time have been infections acquired in endemic areas by travelers. However, the possibility of acquiring malaria by locally infected mosquitoes has been poorly investigated in Europe, despite autochthonous malaria cases having been occasionally reported in several European countries. Here we present an update on the occurrence of potential malaria vector species in Europe. Adopting a systematic review approach, we selected 288 papers published between 2000 and 2021 for inclusion in the review based on retrieval of accurate information on the following Anopheles species: An. atroparvus, An. hyrcanus sensu lato (s.l.), An. labranchiae, An. maculipennis sensu stricto (s.s.), An. messeae/daciae, An. sacharovi, An. superpictus and An. plumbeus. The distribution of these potential vector species across Europe is critically reviewed in relation to areas of major presence and principal bionomic features, including vector competence to Plasmodium. Additional information, such as geographical details, sampling approaches and species identification methods, are also reported. We compare the information on each species extracted from the most recent studies to comparable information reported from studies published in the early 2000s, with particular reference to the role of each species in malaria transmission before eradication. The picture that emerges from this review is that potential vector species are still widespread in Europe, with the largest diversity in the Mediterranean area, Italy in particular. Despite information on their vectorial capacity being fragmentary, the information retrieved suggests a re-definition of the relative importance of potential vector species, indicating An. hyrcanus s.l., An. labranchiae, An. plumbeus and An. sacharovi as potential vectors of higher importance, while An. messeae/daciae and An. maculipennis s.s. can be considered to be moderately important species. In contrast, An. atroparvus and An. superpictus should be considered as vectors of lower importance, particularly in relation to their low anthropophily. The presence of gaps in current knowledge of vectorial systems in Europe becomes evident in this review, not only in terms of vector competence but also in the definition of sampling approaches, highlighting the need for further research to adopt the appropriate surveillance system for each species.

Diversity and distribution of larval habitats of mosquitoes (Diptera: Culicidae) in northern Spain: from urban to natural areas

Studies of the biodiversity of mosquito larval habitats are important for vector-borne disease control programs and help to improve vector distribution maps. This study was designed to investigate the geographical distribution of mosquito species and their larval habitats in urban, rural, and natural areas in northern Spain. Pre-imaginal stages were collected over two sampling periods (spring and summer) in 2019. In the laboratory, immature specimens were reared to the adult stage for species identification by morphological taxonomy and/or molecular methods. In total, 2,182 specimens belonging to 13 different native mosquito species of five genera were collected from 135 sampling points of which 59.2% harbored larvae. Culex pipiens s.l. was the most abundant species (45.1%), followed by Culex torrentium (12.3%), Anopheles maculipennis s.l. (10.2%), Culex hortensis (9.5%), and nine other species with lower frequencies that accounted for less than 25%. By molecular identification, An. maculipennis s.s. was identified as the only species within the An. maculipennis species complex and Cx. pipiens pipiens as the predominant subspecies of the Cx. pipiens species complex. Margins in large sunlit water bodies were the most suitable sites for An. maculipennis s.l., whereas Cx. pipiens s.l. was present in both natural and artificial habitats. The larval site index, species richness, and relative abundance of the mosquitoes were determined based on the characteristics of the sites where they were collected. The public health importance and ecology of some identified mosquitoes is also discussed.

Entomological Surveillance in Former Malaria-endemic Areas of Southern Italy

Malaria still represents a potential public health issue in Italy, and the presence of former Anopheles vectors and cases imported annually merit continuous surveillance. In areas no longer endemic, the concurrent presence of gametocyte carriers and competent vectors makes re-emergence of local transmission possible, as recently reported in Greece. In October 2017, due to the occurrence of four suspected introduced malaria cases in the province of Taranto (Apulia region), entomological investigations were performed to verify the involvement of local anopheline species. In 2019–2020 entomological surveys were extended to other areas historically prone to malaria between the provinces of Taranto and Matera and the province of Foggia (Gargano Promontory). Resting mosquitoes were collected in animal shelters and human dwellings, larvae were sampled in natural and artificial breeding sites, and specimens were both morphologically and molecularly identified. A total of 2228 mosquitoes were collected, 54.3% of which were anophelines. In all the investigated areas, Anopheles labranchiae was the most widespread species, while Anopheles algeriensis was predominant at the Gargano sites, and Anopheles superpictus and Anopheles plumbeus were recorded in the province of Matera. Our findings showed a potentially high receptivity in the surveyed areas, where the abundance of the two former malaria vectors, An. labranchiae and An. superpictus, is related to environmental and climatic parameters and to anthropic activities.

New Cytogenetic Photomap and Molecular Diagnostics for the Cryptic Species of the Malaria Mosquitoes Anopheles messeae and Anopheles daciae from Eurasia

The Eurasian malaria vector Anopheles messeae is a widely spread and genetically diverse species. Five widespread polymorphic chromosomal inversions were found in natural populations of this mosquito. A cryptic species, Anopheles daciae, was differentiated from An. messeae by the presence of several nucleotide substitutions in the Internal Transcribed Spacer 2 (ITS2) region of ribosomal DNA. However, because of the absence of a high-quality reference cytogenetic map, the inversion polymorphisms in An. daciae and An. messeae remain poorly understood. Moreover, a recently determined heterogeneity in ITS2 in An. daciae questioned the accuracy of the previously used Restriction Fragment Length Polymorphism (RFLP) assay for species diagnostics. In this study, a standard-universal cytogenetic map was constructed based on orcein stained images of chromosomes from salivary glands for population studies of the chromosomal inversions that can be used for both An. messeae and An. daciae. In addition, a new ITS2-RFLP approach for species diagnostics was developed. Both methods were applied to characterize inversion polymorphism in populations of An. messeae and An. daciae from a single location in Western Siberia in Russia. The analysis demonstrates that cryptic species are remarkably different in their frequencies of chromosomal inversion variants. Our study supports previous observations that An. messeae has higher inversion polymorphism in all autosomes than the cryptic species An. daciae.

Geometric morphometrics and machine learning as tools for the identification of sibling mosquito species of the Maculipennis complex (Anopheles)

Geometric morphometrics allows researchers to use the specific software to quantify and to visualize morphological differences between taxa from insect wings. Our objective was to assess wing geometry to distinguish four Anopheles sibling species of the Maculipennis complex, An. maculipennis s. s., An. daciae sp. inq., An. atroparvus and An. melanoon, found in Northern Italy. We combined the geometric morphometric approach with different machine learning alghorithms: support vector machine (SVM), random forest (RF), artificial neural network (ANN) and an ensemble model (EN). Centroid size was smaller in An. atroparvus than in An. maculipennis s. s. and An. daciae sp. inq. Principal component analysis (PCA) explained only 33% of the total variance and appeared not very useful to discriminate among species, and in particular between An. maculipennis s. s. and An. daciae sp. inq. The performance of four different machine learning alghorithms using procrustes coordinates of wing shape as predictors was evaluated. All models showed ROC-AUC and PRC-AUC values that were higher than the random classifier but the SVM algorithm maximized the most metrics on the test set. The SVM algorithm with radial basis function allowed the correct classification of 83% of An. maculipennis s. s. and 79% of An. daciae sp. inq. ROC-AUC analysis showed that three landmarks, 11, 16 and 15, were the most important procrustes coordinates in mean wing shape comparison between An. maculipennis s. s. and An. daciae sp. inq. The pattern in the three-dimensional space of the most important procrustes coordinates showed a clearer differentiation between the two species than the PCA. Our study demonstrated that machine learning algorithms could be a useful tool combined with the wing geometric morphometric approach.