The wing venation patterns to identify single tsetse flies

Entomological study on vectorial density, temporal variation of tsetse fly and other biting flies in intervention and non-intervention areas of South Omo Zone, Ethiopia

A repeated cross-sectional entomological survey was conducted to estimate Glossina (tsetse) and other biting flies density, their seasonal variation and associated risk factors in intervention and non- intervention areas of South Omo Zone, Southwest Ethiopia from January 2019-November 2019. In both dry and wet seasons, a total of 96 NGU traps (64 traps in tsetse intervention districts and 32 traps in tsetse non- intervention districts) were deployed at an interval of about 100-200 m in purposively selected and suspected tsetse habitats. Thus, Glossina pallidipes was found to be the only cyclical vector along with mechanical vectors of Tabanus, Stomoxys and Haematopota. In tsetse intervention areas, G. pallidipes apparent density of 2.64 F/T/D and 0.42 F/T/D was recorded in dry and wet season respectively. Mechanical vectors (dry; wet) of Tabanus (205; 155), Stomoxys (34; 54) and Haematopota (50; 33) were also recorded in tsetse intervened areas. Whereas, in non- intervention areas, apparent density of G. pallidipes was 2.03 F/T/D and 0.56 F/T/D, respectively in dry and wet season. Similarly, Tabanus (22; 56), Stomoxys (10; 8) and Haematopota (5; 7) respectively in dry and wet (dry; wet) season were recorded in tsetse non- intervention areas. According to Negative Binomial Regression (NBR), season was the only variable significantly affecting (P < 0.05) the Glossina count in the current study area. Accordingly, the incidence G. pallidipes during wet season was decreased by the factor of 0.21 (CI; 0.097-0.47) when compared to its incidence in dry season by holding other variables constant. In conclusion, cyclical vectors were playing vital role in transmission of trypanosomosis in South Omo Zone along with numerous mechanical vectors even though there have been vector intervention activities in the areas. Therefore, strong, sustainable, environmentally friend and community participating vector control strategies should be followed to tackle the vector distribution in the area.

Geometry of posterior larval spiracles to identify medically and forensically important calliphorids in Thailand

Fly identification is the primary step of analysis in forensic entomology. Although morphology and molecular techniques are considered satisfactory methods, some constraints may arise from a financial or even human point of view. Over the past decade, the geometric morphometric approach has been increasingly advocated for the classification and identification of arthropods. This study explored the method for species identification of 800 third-instar larvae of eight blow fly species of medical and forensic importance: Chrysomya chani Kurahashi, Chrysomya megacephala (Fabricius), Chrysomya (Ceylonomyia) nigripes Aubertin, Chrysomya pinguis (Walker), Chrysomya (Achoetandrus) rufifacies (Macquart), Hemipyrellia ligurriens (Wiedemann), Lucilia cuprina (Wiedemann), and Lucilia porphyrina (Walker). Based on the posterior spiracles geometry, the cross-validation revealed a relatively high percentage of correct classification in most species, ranking from 86% to 100%. The results of this study confirmed that the geometric morphometric (GM) analysis of posterior spiracles might be utilized as a larva identification tool. Therefore, this GM method represents one way of overcoming difficulties with the identification of blow fly larvae and can support further studies of these flies.

Application of biomolecular techniques on tsetse fly puparia for species identification at larvipostion sites

Puparia are commonly found in tsetse fly larviposition sites during studies on larval ecology. This chitinous shell is representative of past or ongoing exploitation of these sites by tsetse flies. The morphological characteristics of the puparium are not sufficiently distinctive to allow identification of the species. This study explores the applicability of biomolecular techniques on empty puparia for tsetse fly species identification. Five techniques were compared for DNA extraction from tsetse fly puparia, 1/Chelex® 100 Resin, 2/CTAB, 3/Livak's protocol, 4/DEB + proteinase K and 5/QIAamp® DNA Mini kit, using two homogenisation methods (manual and automated). Using a combination of two primer pairs, Chelex, CTAB, and DEB + K proved the most efficient on fresh puparia with 90, 85, and 70% samples identified, respectively. Shifting from fresh to one- to nine-month-old puparia, the Chelex method gave the best result allowing species identification on puparia up to seven months old. The subsequent testing of the Chelex extraction protocol identified 152 (60%) of 252 field-collected puparia samples at species level. The results show that reliable genetic identification of tsetse flies species can be performed from empty puparia, what can prove of great interest for future ecological studies on larviposition sites. The Chelex technique was the most efficient for DNA extraction, though the age-limit of the samples stood at seven months, beyond which DNA degradation probably compromises the genetic analysis.

Wing morphometrics of biting midges (Diptera: Culicoides) of veterinary importance in Madagascar

Biting midges are vectors of arboviruses such as bluetongue virus, bovine ephemeral fever virus, Akabane virus, African horse sickness virus, epizootic haemorrhagic disease virus and Schmallenberg virus. Fast and accurate identification of biting midges is crucial in the study of Culicoides-borne diseases. Morphological identification of biting midges has revealed the presence of cryptic species. A total of 20 species are reported in Madagascar. In this study, we assessed wing morphometric analysis for identification of seven species namely C. dubitatus Kremer, Rebholtz-Hirtzel and Delécolle, C. enderleini Cornet and Brunhes, C. kibatiensis Goetghebuer, C. miombo Meiswinkel, C. moreli Clastrier, C. nevilli Cornet and Brunhes, and C. zuluensis de Meillon. Culicoides enderleini, C. miombo, C. moreli, C. nevilli and C. zuluensis are vectors diseases. A molecular approach, based on the cytochrome oxidase I gene (Cox1), was used for species delimitation. The molecular analysis presented seven different clades grouped two-by-two according to morphological characters. A total of 179 wing images were digitised. We found morphometric variation among seven species based on 11 landmarks and two outlines. Wing shape variation plots showed that species overlapped with species belonging to the same group. The cross-validation revealed a relatively high percentage of correct classification in most species, ranging from 91.3% to 100% for landmarks; 60% to 82.6% for outlines-1 and 77.1% to 91.3% for outlines-2. Our study suggests that wing geometric morphometric analysis is a robust tool for reliable "Moka Fohy" identification in Madagascar. This inexpensive and simple method is a precise supplement to morphological identification, with reaches the accuracy of Cox1 barcoding.

Wing geometric morphometrics to distinguish and identify Haematobosca flies (Diptera: Muscidae) from Thailand

The hematophagous flies of the genus Haematobosca Bezzi, 1907 (Diptera: Muscidae) are important ectoparasites in domestic animals and wildlife. Two species of this genus have been recorded in Thailand, viz., Haematobosca sanguinolenta (Austen, 1909) and Haematobosca aberrans (Pont, Duvallet & Changbunjong, 2020). They have a similar morphology and coexist in the same habitat. The correct species identification of these flies is extremely important for understanding disease epidemiology and developing effective control measures. Geometric morphometrics (GM) has been confirmed to be a useful tool for differentiating and identifying morphologically similar insect species. Therefore, GM was used to distinguish and identify H. sanguinolenta and H. aberrans in Thailand. Adult flies of both sexes were collected using Nzi traps, morphologically identified, and analyzed by landmark-based GM of the wing. Results showed that GM was highly effective in distinguishing the two Haematobosca species based on their wing shape, with an overall accuracy score of 99.3%. We also revealed that our study material could be used as reference data to identify new field specimens collected from other geographic locations. We propose that wing GM can be used as a supplement to conventional morphology identification, particularly for Haematobosca specimen that has been damaged or has lost its diagnostic characteristics due to specimen collection and processing in the field.

Deep learning approaches to landmark detection in tsetse wing images

Morphometric analysis of wings has been suggested for identifying and controlling isolated populations of tsetse (Glossina spp), vectors of human and animal trypanosomiasis in Africa. Single-wing images were captured from an extensive data set of field-collected tsetse wings of species Glossina pallidipes and G. m. morsitans. Morphometric analysis required locating 11 anatomical landmarks on each wing. The manual location of landmarks is time-consuming, prone to error, and infeasible for large data sets. We developed a two-tier method using deep learning architectures to classify images and make accurate landmark predictions. The first tier used a classification convolutional neural network to remove most wings that were missing landmarks. The second tier provided landmark coordinates for the remaining wings. We compared direct coordinate regression using a convolutional neural network and segmentation using a fully convolutional network for the second tier. For the resulting landmark predictions, we evaluate shape bias using Procrustes analysis. We pay particular attention to consistent labelling to improve model performance. For an image size of 1024 × 1280, data augmentation reduced the mean pixel distance error from 8.3 (95% confidence interval [4.4,10.3]) to 5.34 (95% confidence interval [3.0,7.0]) for the regression model. For the segmentation model, data augmentation did not alter the mean pixel distance error of 3.43 (95% confidence interval [1.9,4.4]). Segmentation had a higher computational complexity and some large outliers. Both models showed minimal shape bias. We deployed the regression model on the complete unannotated data consisting of 14,354 pairs of wing images since this model had a lower computational cost and more stable predictions than the segmentation model. The resulting landmark data set was provided for future morphometric analysis. The methods we have developed could provide a starting point to studying the wings of other insect species. All the code used in this study has been written in Python and open sourced.

Phenotypic divergence in an island bee population: Applying geometric morphometrics to discriminate population-level variation in wing venation

Phenotypic divergence is an important consequence of restricted gene flow in insular populations. This divergence can be challenging to detect when it occurs through subtle shifts in morphological traits, particularly in traits with complex geometries, like insect wing venation. Here, we employed geometric morphometrics to assess the extent of variation in wing venation patterns across reproductively isolated populations of the social sweat bee, Halictus tripartitus. We examined wing morphology of specimens sampled from a reproductively isolated population of H. tripartitus on Santa Cruz Island (Channel Islands, Southern California). Our analysis revealed significant differentiation in wing venation in this island population relative to conspecific mainland populations. We additionally found that this population-level variation was less pronounced than the species-level variation in wing venation among three sympatric congeners native to the region, Halictus tripartitus, Halictus ligatus, and Halictus farinosus. Together, these results provide evidence for subtle phenotypic divergence in an island bee population. More broadly, these results emphasize the utility and potential of wing morphometrics for large-scale assessment of insect population structure.

Species Discrimination of Stomoxys Flies S. bengalensis, S. calcitrans, and S. sitiens (Diptera: Muscidae) Using Wing Geometric Morphometrics

The flies of the genus Stomoxys Geoffroy, 1762 (Diptera: Muscidae), are regarded as pests of veterinary and medical importance. In Thailand, Stomoxys calcitrans (Linnaeus, 1758) is the most abundant species and is widely distributed throughout the country. This Stomoxys species can coexist with two other morphologically similar species: Stomoxys bengalensis Picard, 1908, and Stomoxys sitiens Rondani, 1873. Hence, discriminating using morphological characteristics is difficult, especially if the specimen is damaged or loses its diagnostic characteristics. In this study, we evaluated the effectiveness of the landmark-based geometric morphometric (GM) approach to discriminate among the three Stomoxys spp.: S. bengalensis, S. calcitrans, and S. sitiens. Left-wing images of S. bengalensis (n = 120), S. calcitrans (n = 150), and S. sitiens (n = 155) were used for the GM analyses. The results of the wing shape analyses revealed that the GM approach was highly effective for discriminating three Stomoxys, with high accuracy scores ranging from 93.75% to 100%. This study adds to the evidence that landmark-based GM is an excellent alternative approach for discriminating Stomoxys species.

Wing morphology variations in Culicoides circumscriptus from France

The biting midge Culicoides circumscriptus Kieffer, 1918 is a European widespread vector of avian malaria throughout the continent and is a possible vector of Akabane virus and Bluetongue virus. This species populates a wide range of environments in contrasting ecological settings often exposed to strong seasonal fluctuations. The main goals of this study were to investigate C. circumscriptus phenotypic variation at three departments in France (Corsica Island, Moselle and Var) and to determine if its phenotypes vary with the environment. Culicoides circumscriptus wing phenotypes were analyzed using a geometric morphometric approach based on anatomical landmarks and outlines of the wing. Dendogram trees based on landmarks and the outlines-2 set (cell m4) showed similar topologies and separated populations of C. circumscriptus. In contrast, another set of outlines-1 (covering the r-m cross vein, M, radiale and arculus) presented a different hierarchical clustering tree. The phenotypic variation observed in C. circumscriptus indicated that these populations are exposed to environmental and ecological pressures. Our results suggest the presence of phenotypic plasticity in this species.

Evaluating the efficacy of various traps in catching tsetse flies at Nech Sar and Maze National Parks, Southwestern Ethiopia: An Implication for Trypanosoma Vector Control

Tsetse flies are the vector of protozoan parasite of the genus Trypanosoma, the causative agent of human African sleeping sickness and animal trypanosomiasis. Traps such as Nguruman (NGU), biconical and sticky traps are in use for tsetse flies sampling and monitoring. However, there is no evidence regarding their comparative efficiency in catching flies using olfactory cues. Therefore, the present study aimed to evaluate the efficiency of different types of traps in catching tsetse flies at Nech Sar and Maze National Parks, Southwestern Ethiopia. The study was done for six consecutive months from February to July 2019. Briefly, a 3×4 Latin square design was performed, and tsetse flies were collected for three days each month in four different vegetation types, including wood grassland, bush land, forest, and riverine forest. To avoid trapping position bias, rotation of traps has been done every day. Almost all (99.5%) of the flies were Glossina pallidipes and the remaining were G. fuscipes. The latter were present only at Maze national park. NGU traps were the most efficient type with 12.1 flies/trap/day at Nech Sar National Park and it was 2.2 flies/trap/day at Maze National Park followed by biconical and sticky traps. The number of tsetse flies collected by biconical trap was three-fold lower than NGU trap, and it was four-fold lower in sticky trap than NGU trap in both Nech Sar and Maze National Parks. A substantial number (41%) of G. pallidipes were collected from woody grassland (WGL). In conclusion, G. pallidipes monitoring and evaluation activities could consider NGU trap model as it performed better in most vegetation types in the region.

Wing Interferential Patterns (WIPs) and machine learning, a step toward automatized tsetse (Glossina spp.) identification

A simple method for accurately identifying Glossina spp in the field is a challenge to sustain the future elimination of Human African Trypanosomiasis (HAT) as a public health scourge, as well as for the sustainable management of African Animal Trypanosomiasis (AAT). Current methods for Glossina species identification heavily rely on a few well-trained experts. Methodologies that rely on molecular methodologies like DNA barcoding or mass spectrometry protein profiling (MALDI TOFF) haven't been thoroughly investigated for Glossina sp. Nevertheless, because they are destructive, costly, time-consuming, and expensive in infrastructure and materials, they might not be well adapted for the survey of arthropod vectors involved in the transmission of pathogens responsible for Neglected Tropical Diseases, like HAT. This study demonstrates a new type of methodology to classify Glossina species. In conjunction with a deep learning architecture, a database of Wing Interference Patterns (WIPs) representative of the Glossina species involved in the transmission of HAT and AAT was used. This database has 1766 pictures representing 23 Glossina species. This cost-effective methodology, which requires mounting wings on slides and using a commercially available microscope, demonstrates that WIPs are an excellent medium to automatically recognize Glossina species with very high accuracy.

Intraspecific variation in wing geometry among Tabanus rubidus (Diptera: Tabanidae) populations in Thailand

Tabanus rubidus (Wiedemann, 1821) (Diptera: Tabanidae) is a hematophagous insect of veterinary and medical importance and is the predominant Tabanus spp. in Thailand. It is a potential mechanical vector of Trypanosoma evansi, which causes surra in domestic and wild animals. Wing geometric morphometrics is widely used as morphological markers for species identification and to assess the insect population structure. Herein, we investigated the intraspecific variation in wing geometry among T. rubidus populations in Thailand using landmark-based geometric morphometric analysis. Tabanus rubidus females were collected from five populations in four geographical regions in Thailand. The left wings of 240 specimens were removed and digitized using 22 landmarks for analysis. While wing size variations were found between some populations, wing shape variations were detected in all. These intraspecific variations in T. rubidus populations indicate an adaptive response to the local environmental conditions.

Landmark Data to Distinguish and Identify Morphologically Close Tabanus spp. (Diptera: Tabanidae)

Tabanus spp., also known as horse flies (Diptera: Tabanidae), are important vectors of several animal pathogens. Adult females of Tabanus&nbsp;megalops and Tabanus striatus, which are members of the T. striatus complex, are morphologically similar and hence difficult to distinguish using morphological characteristics. In addition, molecular identification by DNA barcoding is also unable to distinguish these species. These two species can occur sympatrically with Tabanus rubidus, which is morphologically similar to T. megalops and T. striatus. Wing geometric morphometrics has been widely used in various insects to distinguish morphologically similar species. This study explored the effectiveness of landmark-based geometrics at distinguishing and identifying T. megalops, T. rubidus, and T. striatus in Thailand. Specimens were collected from different geographical regions of Thailand, and only unambiguously identified specimens were used for geometric morphometric analyses. Left wings of females of T. megalops (n = 160), T. rubidus (n = 165), and T. striatus (n = 85) were photographed, and 22 wing landmarks were used for the analysis. Wing shape was able to distinguish among species with high accuracy scores, ranging from 94.38% to 99.39%. We showed that morphologically very close species of Tabanus can be reliably distinguished by the geometry of their wing venation, and we showed how our experimental material could be used as a reference to tentatively identify new field collected specimens.

Molecular mechanisms underlying simplification of venation patterns in holometabolous insects

How mechanisms of pattern formation evolve has remained a central research theme in the field of evolutionary and developmental biology. The mechanism of wing vein differentiation in Drosophila is a classic text-book example of pattern formation using a system of positional-information, yet very little is known about how species with a different number of veins pattern their wings, and how insect venation patterns evolved. Here, we examine the expression pattern of genes previously implicated in vein differentiation in Drosophila in two butterfly species with more complex venation Bicyclus anynana and Pieris canidia. We also test the function of some of these genes in B. anynana. We identify both conserved as well as new domains of decapentaplegic, engrailed, invected, spalt, optix, wingless, armadillo, blistered, and rhomboid gene expression in butterflies, and propose how the simplified venation in Drosophila might have evolved via loss of decapentaplegic, spalt and optix gene expression domains, silencing of vein inducing programs at Spalt-expression boundaries, and changes in gene expression of vein maintenance genes.

Effectiveness of landmark- and semi-landmark-based geometric morphometric to identify four species of Culex mosquitoes in Thailand

Objective:

This research aims to study the effectiveness of landmark- and semi-landmark-based geometric morphometric (LMSL) in the identification of four Culex vectors, namely, C. quinquefasciatus, C. visnui, C. sitiens, and C. whitmorei in Thailand and also compared the potential between LMSL and the landmark-based geometric morphometric method (LM).

Materials and Methods:

The overall size of the wing sample was estimated by the centroid size. Shape variables were computed as principal components of the “partial warp” calculated after generalized procrustes analysis of raw coordinates. Discriminant analysis of the canonical variables performed to explore the shape dissimilarity between Culex species has been shown as a factor map and to calculate the Mahalanobis distance. Size and shape differences based on pairwise Mahalanobis distances were tested using non-parametric methods (1,000 cycles) with Bonferroni correction at a p-value of <0.05.

Results:

A total of 120 individuals were used that were divided into 30 individuals per Culex species. The mean CS of C. sitiens had the largest wings followed by C. visnui, C. quinquefasciatus, and C. whitmorei in LM and LMSM. The patterns of statistical difference in CS of both methods were similar and wing shapes among Culex species were different based on a comparison of pairwise Mahalanobis distances (p < 0.05) in both methods. For the cross-validated reclassification test, LM provided Culex species separation ranging from 54% to 84% and 51% to 93% for LMSM.

Conclusion:

Thus, LMSM is another option to use for the identification in mosquito vectors that have a curved line on the wing specific to the species.

Morphometric diagnosis of Glossina palpalis (Diptera: Glossinidae) population structure in Ghana

OBJECTIVE

This study aimed to identify isolated population(s) of Glossina palpalis in Ghana using geometric morphometrics to evaluate variations in wing-shape and size between populations of the fly from three regions.

RESULTS

Wing shape of G. palpalis tsetse flies from the Northern, Western and Eastern Regions varied significantly between each other. Populations from the Northern and Western Regions varied the most (Mahalanobis Distance = 54.20). The least variation was noticed between populations from the Western and Eastern Regions (MD = 1.99). On morphospace, the Northern population clearly separated from the Eastern and Western populations both of which overlapped. Wing centroid size also significantly varied among populations. Reclassification scores were satisfactory reaching 100% for the Northern population. The Northern population of G. palpalis is possibly isolated from the Western and Eastern Region populations. Meanwhile, a panmictic relationship could be on-going between the Western and Eastern populations. We speculate that geographical distance and subspecific difference between populations are among factors responsible for observed pattern of wing shape variations among the studied populations. The implications of results regarding choice of control strategy and limitations of the study are discussed.

Identification of Muscidae (Diptera) of medico-legal importance by means of wing measurements

Cadavers attract numerous species and genera of Muscidae, both regular elements of carrion insect assemblages, and accidental visitors. Identification of adult Muscidae may be considered difficult, particularly by non-experts. Since species identification is a vital first step in the analysis of entomological material in any forensic entomology orientated experiment and real cases, various alternative methods of species identification have been proposed. We investigated possibility of semiautomated identification by means of wing measurements as an alternative for classic morphology and DNA-based approaches. We examined genus-level identification success for 790 specimens representing 13 genera of the most common European cadavers visiting Muscidae. We found 99.8% of examined specimens correctly identified to the genus-level. Without error, the following were identified: Azelia, Eudasyphora, Graphomya, Hydrotaea, Musca, Muscina, Mydaea, Neomyia, Polietes, Stomoxys and Thricops. Genus-level misidentifications were found only in Helina and Phaonia. Discrimination of examined material on the species level within Hydrotaea (318 specimens representing eight species) and Muscina (163 specimens representing four species) showed lower, yet still high average identification success, 97.2 and 98.8%, respectively. Our results revealed relatively high success in both genus and species identification of Muscidae of medico-legal importance. Semiautomated identification by means of wing measurements can be used by non-experts and does not require sophisticated equipment. This method will facilitate the identification of forensically relevant muscids in comparison to more difficult and more time-consuming identification approaches based on taxonomic keys or DNA-based methods. However, for unambiguous identification of some taxa, we recommend complementary use of identification keys.