Discrimination of four Culex (Culex) species from the Neotropics based on geometric morphometrics

Wing geometric morphometrics and COI barcoding of Culex pipiens subgroup in the Republic of Korea

Two members of the Culex pipiens subgroup, Culex pallens and Culex pipiens f. molestus, are known to occur in the Republic of Korea (ROK). These species exhibit morphologically similar features and are challenging to distinguish below the species level. Therefore, this study utilized wing geometric morphometrics (GM) on the right wing of the Culex pipiens subgroup, alongside sequencing of the cytochrome c oxidase subunit I (COI) region. Mosquitoes were collected from 11 locations between June and October (2020-2022) to minimize regional and seasonal variations. Additionally, Culex pipiens f. pipiens, which is not native to the ROK, was included in the analysis. Culex tritaeniorhynchus, Aedes albopictus, and Anopheles sinensis, the primary vectors in the ROK, were used as outgroups for comparison. All three taxa in the Culex pipiens subgroup could be identified with an 82.4%-97.0% accuracy using GM. However, a comparison of the COI regions of the Culex pipiens subgroup revealed no clear differences between the taxa. These data can be used for accurate identification, contributing to effective mosquito control, in addition to providing a foundation for evolutionary and ecological studies on wing shape differences.

Geometric morphometrics to distinguish the cryptic species Anopheles minimus and An. harrisoni in malaria hot spot villages, western Thailand

Anopheles minimus Theobald 1901 and An. harrisoni Harbach & Manguin 2007 belong to the same species complex. They are morphologically similar and can exist in sympatry but have blood host preferences. The most accurate method for their identification is based on molecular techniques. Here, we measure the level of interspecific discrimination by geometric morphometry. Sixty-seven An. minimus and 22 An. harrisoni specimens were selected based on their morphological integrity and confirmed by identification polymerase chain reaction of internal transcribed spacer 2. These samples were used as reference data allowing for a morphometric identification based on geometric shape. Despite size overlap between the two species, there was a significant shape divergence allowing for differentiation of An. minimus and An. harrisoni with 90% accuracy. An intraspecific study of An. minimus showed a summer period associated to the reducing of wing size, which did not influence the shape-based differentiation of An. harrisoni. Wing venation geometry can be used to distinguish between these cryptic species mainly based on shaped divergence. This study suggests that geometric morphometrics represent a convenient low-cost method to complement morphological identification, especially concerning damaged specimens, i.e., insects having accidentally lost the anatomical features allowing a reliable morphological identification.

WingBank: A Wing Image Database of Mosquitoes

Mosquito-borne diseases affect millions of people and cause thousands of deaths yearly. Vaccines have been hitherto insufficient to mitigate them, which makes mosquito control the most viable approach. But vector control depends on correct species identification and geographical assignment, and the taxonomic characters of mosquitoes are often inconspicuous to non-taxonomists, which are restricted to a life stage and/or even damaged. Thus, geometric morphometry, a low cost and precise technique that has proven to be efficient for identifying subtle morphological dissimilarities, may contribute to the resolution of these types of problems. We have been applying this technique for more than 10 years and have accumulated thousands of wing images with their metadata. Therefore, the aims of this work were to develop a prototype of a platform for the storage of biological data related to wing morphometry, by means of a relational database and a web system named “WingBank.” In order to build the WingBank prototype, a multidisciplinary team performed a gathering of requirements, modeled and designed the relational database, and implemented a web platform. WingBank was designed to enforce data completeness, to ease data query, to leverage meta-studies, and to support applications of automatic identification of mosquitoes. Currently, the database of the WingBank contains data referring to 77 species belonging to 15 genera of Culicidae. From the 13,287 wing records currently cataloged in the database, 2,138 were already made available for use by third parties. As far as we know, this is the largest database of Culicidae wings of the world.

Comparative Phylogeography and Integrative Taxonomy of Ochlerotatus caspius (Dipera: Culicidae) and Ochlerotatus dorsalis

Given that accurately identifying pathogen vectors is vital for designing efficient mosquito control programs based on the proper surveillance of the epidemiologically important species, it has been suggested the complementary use of independently evolving genes and morphometric traits as a reliable approach for the characterization and delimitation of related species. Hence, we examined the spatial distribution of COI mtDNA and ITS2 rDNA variation from the historical perspective of Ochlerotatus caspius (Pallas, 1771) and O. dorsalis (Meigen, 1830), while simultaneously testing the utility of the two markers in integrative species delimitation when combined with phenotypic character analyses of larvae and adults. Despite the striking difference in haplotype diversity (high in COI mtDNA, low in ITS2 rDNA), no evident phylogeographic structure was apparent in the Palearctic O. caspius. The Holarctic O. dorsalis species was subdivided into two highly distinctive COI mtDNA phylogroups which corresponded to the Nearctic and Palearctic regions. Strong support for the independence of the two allopatric evolutionary lineages suggested that geographical barrier and climatic changes during Pleistocene caused vicariance of the ancestral range. COI mtDNA reliably distinguished O. caspius and O. dorsalis, while ITS2 rDNA yet again lacked the proper resolution for solving this problem. An integrative approach based on the larval and adult morphological traits have varying taxonomic applications due to their differential diagnostic values. Thus, by the implementation of an integrative taxonomic approach, we successfully detected species borders between the two epidemiologically relevant species and uncovered the presence of cryptic diversity within O. dorsalis.

Wing geometric morphometrics for identification of mosquito species (Diptera: Culicidae) of neglected epidemiological importance

Culicidae comprises more than 3500 species, some of which are responsible for the spread of various human diseases, causing millions of deaths worldwide. Correct identification of these species is essential for the development of surveillance and control strategies. The most common method of mosquito identification is based on specific traits of the external morphology of species. However, identification of mosquitoes by morphological characters can be inaccurate or even unfeasible if the specimen is damaged or there is a lack of distinguishing features, as in the case of cryptic species complexes. Wing geometric morphometrics is a reliable, affordable tool for the identification of mosquito species, including sibling species. More importantly, it can be used in addition to both traditional morphologic identification methods as well as genetic approaches. Here, wing geometric morphometrics was used to identify sixteen mosquito species from eight genera: Aedes, Coquillettidia, Culex, Limatus, Mansonia, Psorophora, Runchomyia, and Wyeomyia. The 390 specimens used here were collected in São Paulo, Brazil using CDC traps, aspiration, and Shannon traps. Allometry was assessed by multivariate regression of the Procrustes coordinates on centroid size followed by canonical variate analysis and a pairwise cross-validated reclassification test. A Neighbor-Joining tree based on Mahalanobis distances was constructed with 1,000 bootstrap replicates using MorphoJ 1.02 and Past 2.17c. The canonical variate analysis of genera resulted in distinct clusters for Culex, Limatus, and Psorophora and partial overlapping between Aedes, Coquilettidia, and Mansonia, and between Runchomyia and Wyeomyia. Pairwise cross-validated reclassification tests indicated that genera were identified with an accuracy of at least 99% and subgenera with a mean accuracy of 96% and that in 160 of the 240 possible comparisons species were identified with an accuracy of 100%. Our results show that the eight genera in the study were correctly distinguished by wing shape, as were subgenera and most species, demonstrating that wing geometric morphometrics can be used for the identification of the mosquito species studied here.

Wing Geometric Morphometrics as a Tool for the Identification of Culex Subgenus Mosquitoes of Culex (Diptera: Culicidae)

Culex is the largest subgenus within the genus Culex that includes important vectors of diseases. The correct identification of mosquitoes is critical for effective control strategies. Wing geometric morphometrics (WGM) has been used to identify mosquito species alongside traditional identification methods. Here, WGM was used for eleven Culex species from São Paulo, Brazil, and one from Esquel, Argentina. Adult mosquitoes were collected using CDC (Centers for Disease Control) traps, morphologically identified and analyzed by WGM. The canonical variate analysis (CVA) was performed and a Neighbor-joining (NJ) tree was constructed to illustrate the patterns of species segregation. A cross-validated reclassification test was also carried out. From 110 comparisons in the cross-validated reclassification test, 87 yielded values higher than 70%, with 13 comparisons yielding 100% reclassification scores. Culexquinquefasciatus yielded the highest reclassification scores among the analyzed species, corroborating with the results obtained by the CVA, in which Cx. quinquefasciatus was the most distinct species. The high values obtained at the cross-validated reclassification test and in the NJ analysis as well as the segregation observed at the CVA made it possible to distinguish among Culex species with high degrees of confidence, suggesting that WGM is a reliable tool to identify Culex species of the subgenus Culex.

Size and shape analysis of Trioza erytreae Del Guercio (Hemiptera: Triozidae), vector of citrus huanglongbing disease

BACKGROUND

The African citrus triozid (ACT) Trioza erytreae Del Guercio (Hemiptera: Triozidae) is one of the most devastating pests of citrus with a well-known role as a vector of the phloem-limited bacteria (Candidatus Liberibacter africanus) associated with huanglongbing (citrus greening disease), currently considered the world's most serious disease of citrus. Although the pest can successfully develop and reproduce on non-citrus host plants, there is no documented information on the geometric morphometry of ACT. We determine the effect of host plants on ACT morphometry under controlled laboratory conditions using traditional and geometric analysis.

RESULTS

ACT reared on C. limon and Citroncirus spp. was significantly larger than when reared on the other host plant species. ACT reared on C. anisata and C. tangelo was consistently smaller than that reared on M. koenigii and C. sinensis. Based on warped outline drawings, ACT reared on Citroncirus spp. and Murraya koenigii had narrower wings than when reared on C. anisata, C. limon, and C. sinensis with slightly broader wing patterns.

CONCLUSION

This study clearly demonstrates that host plant species affect morphometric variation in ACT, which might have a direct impact on fitness parameters of the pest as well as its potential for dispersion. Wing shape and size appear to be useful in separating populations of ACT into different groups. © 2018 Society of Chemical Industry.

Can you make morphometrics work when you know the right answer? Pick and mix approaches for apple identification

Morphological classification of living things has challenged science for several centuries and has led to a wide range of objective morphometric approaches in data gathering and analysis. In this paper we explore those methods using apple cultivars, a model biological system in which discrete groups are pre-defined but in which there is a high level of overall morphological similarity. The effectiveness of morphometric techniques in discovering the groups is evaluated using statistical learning tools. No one technique proved optimal in classification on every occasion, linear morphometric techniques slightly out-performing geometric (72.6% accuracy on test set versus 66.7%). The combined use of these techniques with post-hoc knowledge of their individual successes with particular cultivars achieves a notably higher classification accuracy (77.8%). From this we conclude that even with pre-determined discrete categories, a range of approaches is needed where those categories are intrinsically similar to each other, and we raise the question of whether in studies where potentially continuous natural variation is being categorised the level of match between categories is routinely set too high.

A potential pitfall in studies of biological shape: Does size matter?

The number of published studies using geometric morphometrics (GM) for analysing biological shape has increased steadily since the beginning of the 1990s, covering multiple research areas such as ecology, evolution, development, taxonomy and palaeontology. Unfortunately, we have observed that many published studies using GM do not evaluate the potential allometric effects of size on shape, which normally require consideration or assessment. This might lead to misinterpretations and flawed conclusions in certain cases, especially when size effects explain a large part of the shape variation. We assessed, for the first time and in a systematic manner, how often published studies that have applied GM consider the potential effects of allometry on shape. We reviewed the 300 most recent published papers that used GM for studying biological shape. We also estimated how much of the shape variation was explained by allometric effects in the reviewed papers. More than one-third (38%) of the reviewed studies did not consider the allometric component of shape variation. In studies where the allometric component was taken into account, it was significant in 88% of the cases, explaining up to 87.3% of total shape variation. We believe that one reason that may cause the observed results is a misunderstanding of the process that superimposes landmark configurations, i.e. the Generalized Procrustes Analysis, which removes isometric effects of size on shape, but not allometric effects. Allometry can be a crucial component of shape variation. We urge authors to address, and report, size effects in studies of biological shape. However, we do not propose to always remove size effects, but rather to evaluate the research question with and without the allometric component of shape variation. This approach can certainly provide a thorough understanding of how much size contributes to the observed shaped variation.

Cranial size and shape variation in isolated populations of the Olkhon mountain vole (Alticola olchonensis Litvinov, 1960)

The Olkhon mountain vole (Alticola olchonensis) is an endemic species of the Lake Baikal area with an extremely restricted range. We investigated the pattern of differentiation of cranial shape and size in five isolated insular populations of A. olchonensis from the Baikal islands (Olkhon, Hubyn, Borokchin, Ogoy, and Zamogoy). The ventral aspect of the cranium was analysed using landmark-based geometric morphometric methods While the sexes of A. olchonensis did not differ regarding cranium size and shape, multivariate statistical analyses showed that there were inter-island differences in skull morphology. Voles from the four small islands were similar in cranium size; the large skull size of voles inhabiting the Olkhon Island may be due to the more favourable habitat conditions on Olkhon. Shape differences among island populations were not associated with the geographic distances between islands. The voles from Borokchin Island differ the most in skull shape compared to other island populations. The combination of reduced gene flow among island populations (due to isolation and drift during population fluctuations) may be a suitable explanation for the significant differences in skull shape among island populations of A. olchonensis.

Molecular identification of two Culex (Culex) species of the neotropical region (Diptera: Culicidae)

Culex bidens and C. interfor, implicated in arbovirus transmission in Argentina, are sister species, only distinguishable by feature of the male genitalia; however, intermediate specimens of the species in sympatry have been found. Fourth-instar larvae and females of both species share apomorphic features, and this lack of clear distinction creates problems for specific identification. Geometric morphometric traits of these life stages also do not distinguish the species. The aim of the present study was to assess the taxonomic status of C. bidens and C. interfor using two mitochondrial genes and to determine the degree of their reproductive isolation using microsatellite loci. Sequences of the ND4 and COI genes were concatenated in a matrix of 993 nucleotides and used for phylogenetic and distance analyses. Bayesian and maximum parsimony inferences showed a well resolved and supported topology, enclosing sequences of individuals of C. bidens (0.83 BPP, 73 BSV) and C. interfor (0.98 BPP, 97 BSV) in a strong sister relationship. The mean K2P distance within C. bidens and C. interfor was 0.3% and 0.2%, respectively, and the interspecific variation was 2.3%. Bayesian clustering also showed two distinct mitochondrial lineages. All sequenced mosquitoes were successfully identified in accordance with the best close match algorithm. The low genetic distance values obtained indicate that the species diverged quite recently. Most morphologically intermediate specimens of C. bidens from Córdoba were heterozygous for the microsatellite locus GT51; the significant heterozygote excess observed suggests incomplete reproductive isolation. However, C. bidens and C. interfor should be considered good species: the ventral arm of the phallosome of the male genitalia and the ND4 and COI sequences are diagnostic characters.

Host plant affects morphometric variation of Diaphorina citri (Hemiptera: Liviidae)

The Asian citrus psyllid (ACP), Diaphorina citri Kuwayama, is one of the most serious citrus pests worldwide due to its role as vector of huanglongbing or citrus greening disease. While some optimal plant species for ACP oviposition and development have been identified, little is known of the influence of host plants on ACP size and shape. Our goal was to determine how size and shape of ACP wing and body size varies when development occurs on different host plants in a controlled rearing environment. ACP were reared on six different rutaceous species; Bergera koenigii, Citrus aurantifolia, Citrus macrophylla, Citrus maxima, Citrus taiwanica and Murraya paniculata. Adults were examined for morphometric variation using traditional and geometric analysis based on 12 traits or landmarks. ACP reared on C. taiwanica were consistently smaller than those reared on the other plant species. Wing aspect ratio also differed between C. maxima and C. taiwanica. Significant differences in shape were detected with those reared on M. paniculata having narrower wings than those reared on C. macrophylla. This study provides evidence of wing size and shape differences of ACP based on host plant species which potentially may impact dispersal. Further study is needed to determine if behavioral and physiological differences are associated with the observed phenotypic differences.

The wing venation patterns to identify single tsetse flies

This is the first study to explore the potential of various geometric morphometrics methods to help the morphological diagnostic of tsetse species, vectors of human and animal trypanosomiases in sub-Saharan Africa. We compared landmarks, semilandmarks and outlines techniques on male and female samples of species, and suggested adapted strategies according to the countries and their own Glossina fauna. We could compare up to 7 taxa belonging to the three main subgenera of the Glossina genus: Nemorhina (5 species), Glossina (1 species) and Austenina (1 species). Our sample included the major vectors of sleeping sickness: G. palpalis palpalis, G. p. gambiensis, G. fuscipes fuscipes and G. f. quanzensis, as well as two important vectors of African animal trypanosomoses: G. tachinoides and Glossina morsitans submorsitans. The average level of correct species recognition by the wing shape was satisfactory, and slightly higher for females than for males. The best scores of correct assignment, in both sexes, were obtained by the contour technique (96% of correct attribution in females, 92% in males), slightly higher than for semilandmarks (95% and 91%) or landmarks (94% and 89%) techniques. We made our images of wings freely available to be used as reference images (http://mome-clic.com), and we describe the conditions and the analytical steps to be followed to identify unknown specimens using external reference images. Under adequate conditions, such use of reference images obtained from a free access server could help species identification of new samples anywhere in Africa.

Morphometric Wing Characters as a Tool for Mosquito Identification

Mosquitoes are responsible for the transmission of important infectious diseases, causing millions of deaths every year and endangering approximately 3 billion people around the world. As such, precise identification of mosquito species is crucial for an understanding of epidemiological patterns of disease transmission. Currently, the most common method of mosquito identification relies on morphological taxonomic keys, which do not always distinguish cryptic species. However, wing geometric morphometrics is a promising tool for the identification of vector mosquitoes, sibling and cryptic species included. This study therefore sought to accurately identify mosquito species from the three most epidemiologically important mosquito genera using wing morphometrics. Twelve mosquito species from three epidemiologically important genera (Aedes, Anopheles and Culex) were collected and identified by taxonomic keys. Next, the right wing of each adult female mosquito was removed and photographed, and the coordinates of eighteen digitized landmarks at the intersections of wing veins were collected. The allometric influence was assessed, and canonical variate analysis and thin-plate splines were used for species identification. Cross-validated reclassification tests were performed for each individual, and a Neighbor Joining tree was constructed to illustrate species segregation patterns. The analyses were carried out and the graphs plotted with TpsUtil 1.29, TpsRelw 1.39, MorphoJ 1.02 and Past 2.17c. Canonical variate analysis for Aedes, Anopheles and Culex genera showed three clear clusters in morphospace, correctly distinguishing the three mosquito genera, and pairwise cross-validated reclassification resulted in at least 99% accuracy; subgenera were also identified correctly with a mean accuracy of 96%, and in 88 of the 132 possible comparisons, species were identified with 100% accuracy after the data was subjected to reclassification. Our results showed that Aedes, Culex and Anopheles were correctly distinguished by wing shape. For the lower hierarchical levels (subgenera and species), wing geometric morphometrics was also efficient, resulting in high reclassification scores.

Genetic and phenotypic variation in central and northern European populations of Aedes (Aedimorphus) vexans (Meigen, 1830) (Diptera, Culicidae)

The floodwater mosquito Aedes vexans can be a massive nuisance in the flood plain areas of mainland Europe, and is the vector of Tahyna virus and a potential vector of Dirofilaria immitis. This epidemiologically important species forms three subspecies worldwide, of which Aedes vexans arabiensis has a wide distribution in Europe and Africa. We quantified the genetic and phenotypic variation in Ae. vexans arabiensis in populations from Sweden (northern Europe), Hungary, and Serbia (central Europe). A landscape genetics approach (FST , STRUCTURE, BAPS, GENELAND) revealed significant differentiation between northern and southern populations. Similar to genetic data, wing geometric morphometrics revealed two different clusters, one made by Swedish populations, while another included Hungarian and Serbian populations. Moreover, integrated genetic and morphometric data from the spatial analysis suggested groupings of populations into three clusters, one of which was from Swedish and Hungarian populations. Data on spatial analysis regarding an intermediate status of the Hungarian population was supported by observed Isolation-by-Distance patterns. Furthermore, a low proportion of interpopulation vs intrapopulation variance revealed by AMOVA and low-to-moderate FST values on a broader geographical scale indicate a continuous between-population exchange of individuals, including considerable gene flow on the regional scale, are likely to be responsible for the maintenance of the observed population similarity in Aе. vexans. We discussed data considering population structure in the light of vector control strategies of the mosquito from public health importance.