Research Article | Open Access
Tanawat Chaiphongpachara, "Comparison of Landmark- and Outline-Based Geometric Morphometrics for Discriminating Mosquito Vectors in Ratchaburi Province, Thailand", BioMed Research International, vol. 2018, Article ID 6170502, 10 pages, 2018. https://doi.org/10.1155/2018/6170502
Comparison of Landmark- and Outline-Based Geometric Morphometrics for Discriminating Mosquito Vectors in Ratchaburi Province, Thailand
It is often challenging to identify mosquito vectors in the field based on morphological features due to their similar morphologies and difficulties in obtaining undamaged samples but is required for their successful control. Geometric morphometrics (GM) overcomes this issue by analyzing a suite of traits simultaneously and has the added advantages of being easy to use, low cost, and quick. Therefore, this research compared the efficiency and precision of landmark- and outline-based GM techniques for separating species of mosquitoes in Huay Nam Nak village, Ratchaburi Province, Thailand. This research collected 273 individuals belonging to seven species: Anopheles barbirostris, An. subpictus, Culex quinquefasciatus, Cx. vishnui, Cx. whitmorei, Aedes aegypti, and Ae. albopictus. Both landmark-based and outline-based GM techniques could identify malaria vectors in this area to the genus level successfully and were also very effective for identifying the malaria vectors Anopheles spp. and the dengue vectors Aedes spp. to the species level. However, they were less effective for distinguishing between species of Culex. Therefore, GM represents a valuable tool for the identification of mosquito vectors in the field, which will facilitate their successful control.
Mosquitoes (Diptera: Culicidae) are medically important small insects that transmit many diseases to humans, particularly in tropical and subtropical regions . Globally, there are 3,400 species of mosquitoes in 42 genera, many of which are important vectors of pathogens . Mosquito-borne diseases, which include dengue fever (DF), chikungunya, malaria, Japanese encephalitis, and filariasis , account for 17% of all infectious diseases worldwide and are responsible for approximately 1.4 million deaths per year , making them a serious public health issue.
Mosquito-borne diseases are one of the most important health concerns in Thailand. According to the Bureau of Epidemiology of Thailand, DF had the highest incidence rate in 2014 at 35.67 cases per 100,000 people, followed by malaria (17.48), chikungunya (0.29), filariasis (0.03), and Japanese encephalitis (0.02) . Ratchaburi Province on the western border of Thailand has one of the highest incidences of mosquito-borne diseases, particularly DF and malaria, with morbidity rates of 30.68 and 8.95 cases per 100,000 population, respectively. Therefore, it shows that Ratchaburi Province has an epidemic of mosquito-borne diseases still, making it important that control methods are found to reduce the number of patients in this area.
Mosquito control is one important strategy for controlling mosquito-borne disease epidemics [6–8]. However, since different species of mosquitoes have different characteristics, such as behaviors, breeding sites, and epidemiologies, a knowledge of which vector species are present in the area [9, 10] and the epidemiological patterns of disease transmission is required to ensure the planned control method is appropriate [11, 12]. The It is can be challenging to identify mosquito vectors in the field based on morphological features because many species are cryptic species, sibling species, or isomorphic species with similar morphologies . Furthermore, it is more difficult to obtain a mosquito sample in the field than in the laboratory because the external characteristics are often damaged during trapping and transportation . One solution to this is the use of high-efficiency molecular techniques for mosquito identification. However, these have very high associated costs [14, 15], making them unsuitable for use in the field where many samples are collected. Thus, a new technique for identifying mosquito vectors in the field is required.
Geometric morphometrics (GM) is a fast, inexpensive technique [14, 16] used to analyze the size and shape of individuals based on a suite of traits. GM has been applied widely in a number of fields, including entomology (e.g., mosquito , blow fly , bee , and eggs of Triatominae ). Furthermore, several studies have used the GM method to classify species and to examine variation among medically important mosquitoes that are morphologically similar or sibling species [14, 15, 21–24]. GM analyses can be conducted using landmark- or outline-based methods , each of which has different advantages depending on the characteristics and specificity of the sample. The landmark-based approach uses the coordinates of landmarks to analyze the morphology  and is popular in the field of medical entomology, with many studies having shown it can distinguish different species of mosquito vectors successfully, such as Anopheles spp. and Aedes spp. [14, 21, 23]. By contrast, the outline-based approach uses contour data  and has also been shown to be useful for identifying some mosquito species, such as Aedes scutellaris . However, surprisingly, the specificity with which mosquito vectors can be identified differs between the two techniques [14, 23] and, therefore, a comparison of their performance in the field is required before they can be used.
Thus, the aim of this study was to compare the efficiency of landmark- and outline-based GM techniques to separate species of mosquito vectors in Huay Nam Nak village in Ratchaburi Province, Thailand, an endemic area for mosquito-borne diseases. The result from this research will serve as a guideline for choosing the best GM technique for the identification of vectors in the field in Thailand to facilitate the control of vector-borne diseases.
2. Materials and Methods
2.1. Study Site and Mosquito Collection
Mosquitoes were collected from Huay Nam Nak village (13°22′36.0′′N, 99°16′34.9′′E) in the Suan Phueng District of Ratchaburi Province, Thailand, during June to August 2015 (Figure 1). The study site consisted mainly of wooden houses and agricultural fields, with streams flowing through and mountains and hilly forests surrounding the village.
Three Mosquito Magnet® Independence Mosquito Traps (Woodstream Corporation, USA) were placed approximately 5 meters away from houses in the village for 24 hours (6:00 am to 6:00 am of the following day) per week over the 3-month study period to collect nocturnal and diurnal mosquito vectors. Mosquito samples were collected from each trap in the morning (6.00 am) and transported to the laboratory of the College of Allied Health Science, Suan Sunandha Rajabhat University, Samut Songkhram Province for morphological identification. Female mosquitoes were then morphologically identified to species using the illustrated keys to the mosquitoes of Thailand .
2.2. Mosquito Wing Preparation
Only the right wings of female mosquitoes were analyzed by GM. The right wing of each individual was dissected under a Nikon AZ 100M stereomicroscope (Nikon Corp., Tokyo, Japan) and placed between a glass microscope slide and coverslip using Hoyer’s solution as a mounting medium. Then, all of the wing samples for each species were photographed using a Nikon DS-Ri1 SIGHT digital camera connected to a Nikon Eclipse E600 microscope (Nikon Corp., Tokyo, Japan) under 40× magnification alongside a 1-mm scale bar and were analyzed using the Collecting Landmarks for Identification and Characterization (CLIC) Program. Since the aim of this study was to compare the effectiveness of landmark- and outline-based GM for species discrimination, the same set of wing pictures was used for each species in the two analyses outlined below.
2.3. Landmark-Based GM Analysis
A total of 17 landmarks were selected based on the ease with which they could be plotted across all mosquito species and their low likelihood of being damaged, e.g., by wing vein intersection. The positions of these landmarks on the wing of each individual were digitized (Figure 2). The wing size was then computed as the centroid size (CS), which is the square root of the sum of the squared distances between each individual landmark and the center of the landmark configuration [14, 27]. The wing shape variables (partial warps [PW]) were computed as principal components (PCs) of the PW (known as relative warps [RW]) after Generalised Procrustes Analysis (GPA) (the statistical procedure of superimposition) and discriminant analysis (DA) (or Canonical Variate Analysis [CVA]) was then used to analyze the shape variables for the separation of each species.
2.4. Outline-Based GM Analysis
The outline used for analysis was constructed using coordinates along the contour of the lower part of the mosquito wing (Figure 3). This analytical contour is the most complete part of the wing in all mosquitoes and also lacks thick scales and is difficult to tear. The wing size was estimated as the length of the perimeter of the contour, which was separate from the set of shape variables used. The shape variables were constructed using the normalized elliptic Fourier coefficients (NEF).
2.5. Statistical Analysis
Prior to conducting the GM analyses, this research assessed the quality and measurement error of the digitized landmarks and digitized contour used for the landmark-based and outline-based GM analyses, respectively, by calculating the repeatability index (R). The repeatability of each coordinate was tested for 70 randomly selected wings (10 wings per species). In order to do this, measurements of the coordinates used in the original wing photographs were repeated to give a total of 140 images. The repeatability was then computed based on an ANOVA design. Prior to analysis, the error was reduced by averaging the two digitizations [23, 25].
The variation in wing CS (in mm) of mosquito samples each species was illustrated by quantile plots between P 25% and P 75%. Differences in the wing CS (for landmark-based GM) and the length of the perimeter (for outline-based GM) among species were analyzed using nonparametric permutation tests (1,000 runs) with Bonferroni correction test and a significance level of p < 0.05.
Wing shape variation among species was also visualized by superimposition the mean landmark configurations using Procrustes superimposition in the landmark-based analysis and using EFA in the outline-based analysis. The shape variables for each species were separated using DA and displayed as a factor map. The Mahalanobis distance was calculated from DA to assess the degree of similarity between population and differences in the Mahalanobis distance (i.e., wing shape) among species were computed using nonparametric permutation tests (1,000 runs) with Bonferroni correction test and a significance level of p < 0.05. Then, validated reclassification, whereby each individual was allocated to its closest group according to the Mahalanobis distance without being used to help determine the group center, was performed , and a single linkage hierarchical classification tree was created based on the bootstrap technique on input data within species, according to Morales Vargas et.al.  and visualized using XYOM.
Finally, since size and shape are not independent attributes, this research examined the relationship between the two (the allometric effect) also using linear correlation  after the GM analysis of the size variables.
Both the landmark- and outline-based GM analyses were performed using the CLIC package version 97, which was developed by Professor Jean-Pierre Dujardin [13, 29] and is freely available at http://xyom-clic.eu/. Various modules of the CLIC package were used, including the COO module to digitize the landmarks or pseudolandmarks. The TET module was used to transform the data for analysis. The MOG and FOG modules were used to construct the size and shape variables, perform principle component analyses (PCA) and DA, compute Procrustes distances, and create quantile plots in the landmark- and outline-based analyses, respectively. The VAR module was used to analyze allometry and the statistical significance of differences in the size variables among species. The PAD module was used to analyze the statistical significance of differences in the shape variables among species. The single linkage hierarchical classification tree as built by the recent online morphometric package, XYOM (https://xyom.io), was compared to the UPGMA tree as computed by the R software (https://cran.r-project.org/).
This research collected 273 mosquitoes belonging to seven species within three genera (Table 1).
The two sets of repeated measurements from the same images used in the landmark- and outline-based GM analyses showed good repeatability scores for both size and shape. In the landmark-based analysis, the repeatability of the CS was 0.98, while in the outline-based analysis, the repeatability of the perimeter length was 0.97.
This research found a weak relationship between size and shape, with the allometric residuals of the discriminant factors explaining 8% of the variation in the landmark-based GM and 5% of the variation in the outline-based GM. However, since both of the components are important for species identification, neither was removed from the GM analyses.
3.3. Wing Size
In the landmark-based analysis, An. barbirostris had the largest wings (mean = 3.54 mm) and Ae. aegypti had the smallest (mean = 2.19 mm) (Table 2 and Figure 4). In addition, An. barbirostris had significant intraspecific variation in mean CS that was the highest among all species.
Species with different superscripts letters had significantly different wing sizes at p < 0.05. Min, minimum; Max, maximum; S.D., standard deviation.
In the outline-based analysis, An. barbirostris again had the largest wings (mean = 4.51 mm) while Cx. whitmorei had the smallest (mean = 3.39 mm). Both An. barbirostris and Ae. albopictus exhibited significant intraspecific variation in the mean perimeter length of the contour.
3.4. Wing Shape
Superimposition of the mean wing shapes of each mosquito species using Procrustes superimposition showed the positions of landmarks 1, 7, 8, 12, 13, and 17 varied among species in the landmark-based analysis (Figure 5), and the outlines of Ae. aegypti and Ae. albopictus were distinct from the other species in the outline-based analysis (Figure 6).
Comparison of the factor maps derived from DA showed both the landmark- and outline-based GM analyses gave very similar results, particularly at the genus level (Figure 7). Based on the Mahalanobis distances, wing shape was significantly different among all species in both the landmark- and outline-based GM analyses (nonparametric permutation test, 1,000 cycles, p < 0.05; Table 3).The greatest Mahalanobis distances were between Ae. aegypti and An. subpictus (8.86) in the landmark-based analysis and between Ae. albopictus and An. barbirostris (22.67) in the outline-based analysis.
indicates statistical significance at p< 0.05.
The validated reclassification scores were high (>80% accuracy) for five species using the landmark-based approach (An. barbirostris [96%], An. subpictus [81%], Cx. whitmorei [91%], Ae. aegypti [88%], and Ae. albopictus [90%]) and four species using the outline-based approach (An. barbirostris [93%], An. subpictus [86%], Cx. vishnui [85%], Ae. aegypti [86%], and Ae. albopictus [95%]). The validated reclassification scores were lowest for Cx. vishnui with the landmark-based approach (65%) and for Cx. quinquefasciatus with the outline-based approach (42%) (Table 4). Single Linkage Hierarchical classification trees separated each species of mosquito (Figure 8). The trees that were produced using the two different approaches were very similar.
This research analyzed 273 wings across seven species of mosquitoes using two different GM approaches. Mosquito wings are almost bidimensional and relatively rigid, which reduces mistakes when digitizing them for GM analysis , as reflected by the good repeatability scores for both size and shape in the present study.
4.1. Wing Size
This result found few differences between the landmark- and outline-based analyses in terms of the wing CS, with both methods showing An. barbirostris and Aedes spp. had significantly larger and smaller wings, respectively, than the other species. However, these species also exhibited significant intraspecific differences in size, likely due to variation in factors such as temperature, humidity, and food availability . Because of such variation, it has been suggested previously shape is more appropriate than size for distinguishing among morphologically similar species, and is also more informative in terms of the genetics and evolution of organisms [15, 16, 23]. However, wing size is also useful in the initial identification of species, particularly since some species of mosquitoes are larger than many other species; for example, An. barbirostris was clearly larger than Culex species in this study. Similarly, Stanford et al.  found wing beat frequencies, which mediate assortative mating, are related to size and consequently are unique for each species.
4.2. Wing Shape
The shape and venation of mosquito wings are unique characteristics that can be used to separate different genera and species . This research found little difference between the landmark- and outline-based methods when considering the wing shape of mosquitoes. For both methods, factor maps derived from DA showed there was no overlap between genera but some overlap between species (Figure 7). However, wing shape was found to be significantly different among all seven species using both methods. Similarly, Wilke et al.  found GM was a good tool (100% accuracy) for species identification at the genus level for Aedes, Anopheles, and Culex, but was also efficient at the subgenus and species levels. In both GM analyses, the reclassification scores were higher for Anopheles spp. and Aedes spp. (>80% accuracy in all species) than for Culex spp., with the exception of Cx. whitmorei for the landmark-based analysis (91%) and Cx. vishnui for the outline-based analysis (85%). Consistent with this, the single linkage hierarchical classification trees clearly separated the three genera in the landmark- and outline-based analysis (Figure 8(a)). Superimposition of the mean landmark configurations of each mosquito species showed landmarks 2, 3, 5, 6, 7, 9, and 10 were in similar places for Culex spp. (particularly Cx. quinquefasciatus and Cx. vishnui) and the other species (Figure 5).
These findings show GM can be used successfully to classify Anopheles spp., which are often implicated in malaria transmission; An. barbirostris is a suspected vector of Plasmodium spp. in Thailand , while An. subpictus is not considered a malaria vector in Thailand but is a vector in Sri Lanka . GM is also useful for identifying Aedes spp., including Ae. aegypti and Ae. albopictus, which are dengue and chikungunya vectors, respectively . However, GM had low classification rates for Culex spp., including Cx. quinquefasciatus, Cx. vishnui, and Cx. whitmorei, which have been incriminated as Japanese encephalitis and filariasis vectors .
Our findings demonstrate both GM approaches can be used to identify mosquitoes in the study area, particularly to the genus level. Currently, the landmark-based GM approach is the most popular for helping to identify species and investigating variation among vectors [14, 21, 23] as it is less time-consuming, requiring the definition of only a few analytical points for analysis, and is a powerful approach. By contrast, few studies use the outline-based GM approach because it requires a lot of time and samples than the landmark-based approach. However, this approach has the advantage of not requiring the use of specific locations . Furthermore, this research found that the outline-based GM approach was better and more effective at discriminating some species, such as Cx. vishnui.
Both the landmark- and outline-based GM methods are practical and effective for discriminating between species of mosquitoes in the study area. However, since each species of mosquito has a unique wing identity, the best method needs to be selected to suit the species occurring in a particular vector control area. GM has the advantages of being easy to use, low cost, and quick, as well as not requiring advanced entomological skills, and these advantages make it particularly attractive for use in the field to facilitate the control of mosquito vectors.
The data supporting the conclusions of this article are provided within the article. The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The author declares that there are no conflicts of interest.
The author would like to thank College of Allied Health Science, Suan Sunandha Rajabhat University, Thailand, for their kind support of our research. This work was supported by Suan Sunandha Rajabhat University, Bangkok, Thailand.
- R. Killick-Kendrick, Medical Entomology for Students, 1996.
- M. Service, Medical Entomology for Students, 4th edition, 2008.
- M. A. Tolle, “Mosquito-borne diseases,” Current Problems in Pediatric and Adolescent Health Care, vol. 39, no. 4, pp. 97–140, 2009.
- World Health Organization, Mosquito Born Diseases, WHO, 2016.
- Ministry of Public Health (MOPH), “Annual Report,” Bureau of Epidemiology, MOPH, Thailand, 2014.
- T. Chaiphongpachara, P. Bunyuen, and K. K. Chansukh, “Development of a more effective mosquito trapping box for vector control,” The Scientific World Journal, vol. 2018, Article ID 6241703, 8 pages, 2018.
- T. Chaiphongpachara, S. Laojun, and C. Kunphichayadecha, “Effect of the CDC light trap on control of nocturnal mosquitoes in coastal Samut Songkhram Province, Thailand,” Biodiversitas Journal of Biological Diversity, vol. 19, no. 5, pp. 1750–1754, 2018.
- T. Chaiphongpachara, O. Padidpoo, K. K. Chansukh, and S. Sumruayphol, “Efficacies of five edible mushroom extracts as odor baits for resting boxes to attract mosquito vectors: A field study in Samut Songkhram Province, Thailand,” Tropical Biomedicine, vol. 35, no. 3, pp. 653–663, 2018.
- C. Tananchai, R. Tisgratog, W. Juntarajumnong et al., “Species diversity and biting activity of Anopheles dirus and Anopheles baimaii (Diptera: Culicidae) in a malaria prone area of western Thailand,” Parasites & Vectors, vol. 5, no. 1, article 211, 2012.
- T. Chaiphongpachara and S. Sumruayphol, “Species diversity and distribution of mosquito vectors in coastal habitats of Samut Songkhram province, Thailand,” Tropical Biomedicine, vol. 34, no. 3, pp. 524–532, 2017.
- K. Tainchum, M. Kongmee, S. Manguin, M. J. Bangs, and T. Chareonviriyaphap, “Anopheles species diversity and distribution of the malaria vectors of Thailand,” Trends in Parasitology, vol. 31, no. 3, pp. 109–119, 2015.
- T. Chaiphongpachara, S. Pimsuka, W. S. Na Ayudhaya, and W. Wassanasompong, “The application of geographic information system in dengue haemorrhagic fever risk assessment in samut Songkhram province, Thailand,” International Journal of GEOMATE, vol. 12, no. 30, pp. 53–60, 2017.
- J.-P. Dujardin, “Morphometrics applied to medical entomology,” Infection, Genetics and Evolution, vol. 8, no. 6, pp. 875–890, 2008.
- S. Sumruayphol, C. Apiwathnasorn, J. Ruangsittichai et al., “DNA barcoding and wing morphometrics to distinguish three Aedes vectors in Thailand,” Acta Tropica, vol. 159, pp. 1–10, 2016.
- J. Ruangsittichai, C. Apiwathnasorn, and J.-P. Dujardin, “Interspecific and sexual shape variation in the filariasis vectors Mansonia dives and Ma. bonneae,” Infection, Genetics and Evolution, vol. 11, no. 8, pp. 2089–2094, 2011.
- C. Garros and J. P. Dujardin, “Genetic and Phenetic Approaches to Anopheles Systematics,” in Anopheles Mosquitoes - New Insights into Malaria Vectors, pp. 81–105, 2013.
- C. Lorenz, F. Almeida, F. Almeida-Lopes et al., “Geometric morphometrics in mosquitoes: What has been measured?” Infection, Genetics and Evolution, vol. 54, pp. 205–215, 2017.
- N. Sontigun, K. L. Sukontason, B. K. Zajac et al., “Wing morphometrics as a tool in species identification of forensically important blow flies of Thailand,” Parasites & Vectors, vol. 10, no. 1, 2017.
- D. A. De Souza, Y. Wang, O. Kaftanoglu et al., “Morphometric identification of queens, workers, intermediates in in vitro reared honey bees (Apis mellifera),” PLoS ONE, vol. 10, no. 4, 2015.
- S. Santillán-Guayasamín, A. G. Villacís, M. J. Grijalva, and J.-P. Dujardin, “The modern morphometric approach to identify eggs of Triatominae,” Parasites & Vectors, vol. 10, no. 1, pp. 1–10, 2017.
- G. F. Gómez, E. J. Márquez, L. A. Gutiérrez, J. E. Conn, and M. M. Correa, “Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation,” Acta Tropica, vol. 135, no. 1, pp. 75–85, 2014.
- M. T. Motoki, L. Suesdek, E. S. Bergo, and M. A. M. Sallum, “Wing geometry of Anopheles darlingi Root (Diptera: Culicidae) in five major Brazilian ecoregions,” Infection, Genetics and Evolution, vol. 12, no. 6, pp. 1246–1252, 2012.
- J.-P. Dujardin, “Modern morphometrics of medically important insects,” in Genetics and Evolution of Infectious Diseases, pp. 473–501, 2011.
- N. Jaramillo-O, J.-P. Dujardin, D. Calle-Londoño, and I. Fonseca-González, “Geometric morphometrics for the taxonomy of 11 species of Anopheles (Nyssorhynchus) mosquitoes,” Medical and Veterinary Entomology, vol. 29, no. 1, pp. 26–36, 2015.
- J.-P. Dujardin, D. Kaba, P. Solano, M. Dupraz, K. D. McCoy, and N. Jaramillo-O, “Outline-based morphometrics, an overlooked method in arthropod studies?” Infection, Genetics and Evolution, vol. 28, pp. 704–714, 2014.
- R. Rattanarithikul, B. A. Harrison, P. Panthusiri, and R. E. Coleman, “Illustrated keys to the mosquitoes of Thailand. I. Background; geographic distribution; lists of genera, subgenera, and species; and a key to the genera,” Southeast Asian Journal of Tropical Medicine and Public Health, vol. 36, no. 1, pp. 1–80, 2005.
- J.-P. Dujardin et al., “The exchangeability of shape,” BMC Research Notes, vol. 3, no. 1, 2010.
- R. E. Morales Vargas, N. Phumala-Morales, T. Tsunoda, C. Apiwathnasorn, and J.-P. Dujardin, “The phenetic structure of Aedes albopictus,” Infection, Genetics and Evolution, vol. 13, no. 1, pp. 242–251, 2013.
- M. R. Sanford, B. Demirci, C. D. Marsden, Y. Lee, A. J. Cornel, and G. C. Lanzaro, “Morphological differentiation may mediate mate-choice between incipient species of Anopheles gambiae s.s.,” PLoS ONE, vol. 6, no. 11, 2011.
- A. B. B. Wilke, R. De Oliveira Christe, L. C. Multini et al., “Morphometric wing characters as a tool for mosquito identification,” PLoS ONE, vol. 11, no. 8, 2016.
- P. Sriwichai, Y. Samung, S. Sumruayphol et al., “Natural human Plasmodium infections in major Anopheles mosquitoes in western Thailand,” Parasites & Vectors, vol. 9, no. 1, article 17, 2016.
- S. N. Surendran, P. J. Jude, and R. Ramasamy, “Variations in salinity tolerance of malaria vectors of the Anopheles subpictus complex in Sri Lanka and the implications for malaria transmission,” Parasites & Vectors, vol. 4, no. 1, article no. 117, 2011.
Copyright © 2018 Tanawat Chaiphongpachara. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.