Next Article in Journal
Attributes of Yellow Traps Affecting Attraction of Diaphorina citri (Hemiptera: Liviidae)
Previous Article in Journal
Blood Meal Sources of Anopheles spp. in Malaria Endemic Areas of Honduras
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 11
Issue 7
10.3390/insects11070451
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae)

by
Tanasak Changbunjong
1,2,*,
Jiraporn Ruangsittichai
3,
Gerard Duvallet
4 and
Adrian C. Pont
5
1
Department of Pre-clinic and Applied Animal Science, Faculty of Veterinary Science, Mahidol University, Nakhon Pathom 73170, Thailand
2
The Monitoring and Surveillance Center for Zoonotic Diseases in Wildlife and Exotic Animals (MoZWE), Faculty of Veterinary Science, Mahidol University, Nakhon Pathom 73170, Thailand
3
Department of Medical Entomology, Faculty of Tropical Medicine, Mahidol University, Bangkok 10400, Thailand
4
UMR5175 CEFE, Centre d’Ecologie Fonctionnelle et Evolutive, Université Paul-Valéry, 34090 Montpellier, France
5
Oxford University Museum of Natural History, Parks Road, Oxford OX1 3PW, UK
*
Author to whom correspondence should be addressed.
Insects 2020, 11(7), 451; https://doi.org/10.3390/insects11070451
Submission received: 15 June 2020 / Revised: 8 July 2020 / Accepted: 13 July 2020 / Published: 16 July 2020
Browse Figures
Versions Notes

Abstract

:
The genus Haematobosca Bezzi, 1907 (Diptera: Muscidae) contains haematophagous flies of veterinary importance. A new fly species of this genus was recognised from northern Thailand based on morphological characters and described as Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020. In the present study, the mitochondrial cytochrome c oxidase I (COI) gene was used to confirm the morphological identification of H. aberrans. In addition, landmark-based geometric morphometrics was used to determine sexual dimorphism. The molecular analysis was conducted with 10 COI sequences. The results showed that all sequences were 100% identical. The sequence was not highly similar to reference sequences from GenBank and did not match any identified species from Barcode of Life Data Systems (BOLD). Phylogenetic analysis clearly differentiated this species from other species within the subfamily Stomoxyinae. For geometric morphometric analysis, a total of 16 wing pictures were analysed using the landmark-based approach. The results showed significant differences in wing shape between males and females, with a cross-validated classification score of 100%. The allometric analysis showed that wing shape has no correlation with size. Therefore, the COI gene is effective in species identification of H. aberrans, and geometric morphometrics is also effective in determining sexual dimorphism.

Graphical Abstract

1. Introduction

Haematobosca Bezzi, 1907 (Diptera: Muscidae) is a genus of haematophagous flies in the subfamily Stomoxyinae. The genus comprises 15 species: 11 species in the Afrotropical region, one in the Palaearctic region, two in the Holarctic region and one in the Oriental and Australasian regions [1,2,3,4]. All the species in this genus are of considerable veterinary importance, affecting both domestic animals and wildlife [1,5]. The adult flies are medium-sized flies, 3.5 to 9.0 mm in length and yellow to black in colour. The arista has dorsal and ventral hairs. The palpi are about as long as the proboscis and are grooved internally. Anterior and posterior katepisternal setae (or sternopleural bristles) are present [1]. Recently, a new species of Haematobosca was described from Thailand as Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 [6]. It can be distinguished from other Haematobosca by the absence of the anterior katepisternal seta [6]. However, one female that probably belongs to this species was previously recorded incorrectly as Stygeromyia Austen, 1907 by Tumrasvin and Shinonaga (1978) [7].
The DNA-based identification system is a useful tool for species identification, biodiversity monitoring and molecular phylogenetic investigation [8,9]. The mitochondrial cytochrome c oxidase subunit I (COI) gene with a length of about 650 base pairs (bps) is commonly used as the DNA barcode of animal species. The COI gene has been found to be effective for the identification of various insects of medical and veterinary importance such as black flies [10,11], mosquitoes [12,13,14,15], sand flies [16,17,18,19] and tabanids [20,21,22,23]. The efficiency of the COI gene for the identification of stomoxyine flies has been demonstrated by Changbunjong et al. (2016) [24].
In parallel to DNA-based identification, geometric morphometrics is also increasingly applied to distinguish morphologically similar species of insects of medical and veterinary importance [25,26,27,28,29,30]. This technique is a fast, low-cost and easy to use tool [25]. Geometric morphometric analysis can be conducted using various methods such as landmark and outline-based methods, depending on the characteristics and specifics of the specimens [25,31]. The landmark-based method uses the coordinates of landmarks to analyse the variation in size and shape, while the outline-based method uses the outline or contour data in the analysis [25,31]. Many studies have used geometric morphometric methods for describing sexual dimorphism in the wing of various insects such as blow flies [28,32], mosquitoes [33,34], moths [35,36], fruit flies [37] and stomoxyine flies [26]. Sexual dimorphism in morphological characters is the most interesting source of phenotypic variation in various taxa and has attracted considerable interest in evolutionary biology [37]. Therefore, in our study, the COI gene was used as the molecular marker for the identification of H. aberrans, collected in northern Thailand. In addition, geometric morphometric analysis was used to determine sexual dimorphism.

2. Materials and Methods

2.1. Ethics Statement

Our protocol for specimen collection was approved by the Faculty of Veterinary Science, Mahidol University Animal Care and Use Committee (Ref. MUVS-2018-10-58).

2.2. Specimen Collection and Identification

Specimens of Haematobosca aberrans were collected in the Mae Wang district, Chiang Mai Province of northern Thailand (18°38′35″ N, 98°31′27″ E) in January 2019 using Nzi traps [38] (Figure 1). All flies were collected inside a buffalo farm close to a forest area. They were euthanised in a freezer (−10 °C) and placed in individual 1.5 mL microcentrifuge tubes. They were identified using the description and taxonomic key in Pont et al. (2020) [6]. Specimens were stored at −20 °C for further laboratory investigation.

2.3. Molecular Analysis

2.3.1. DNA Extraction and COI Amplification

Genomic DNA was extracted from the legs of individual specimens using DNeasy® Blood & Tissue Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol. PCR amplification of the COI gene with a 658 base pair (bp) fragment was performed using a primer pair consisting of LepF1 (5′-ATTCAACCAATCATAAAGATATTGG-3′) and LepR1 (5′-TAAACTTCTGGATGTCCAAAAAATCA-3′) [39]. The PCR cycle was performed in a C1000TM Thermal Cycler (Bio-Rad, Hercules, CA, USA) with a total reaction volume of 50 μL containing: 5 μL of 10× PCR buffer, 2.50 μL of 50 mM MgCl2, 36.01 μL of distilled water, 0.25 μL of 10 mM dNTPs, 0.24 μL of 1.2 unit of platinum Taq polymerase (Invitrogen, Carlsbad, CA, USA), 0.50 μL of 10 mM forward primer, 0.50 μL of 10 mM reverse primer, and 5 μL of DNA template. The PCR cycling conditions consisted of initial denaturation at 94 °C for 1 min; 5 cycles of 94 °C for 30 s, annealing at 45 °C for 40 s, and extension at 72 °C for 1 min; 30–35 cycles of 94 °C for 30 s, 55 °C for 40 s, and 72 °C for 1 min; final extension at 72 °C for 10 min; followed by indefinite hold at 4 °C. Amplified PCR products were separated by gel electrophoresis on 2% agarose gel, with a GeneRuler™ 100 bp DNA ladder (Thermo Scientific, Vilnius, Lithuania), which was used to estimate the size of DNA fragments on the gel. The voucher specimens are deposited in the collection of the Faculty of Veterinary Science, Mahidol University (VSMU), Thailand.

2.3.2. DNA Sequencing and Data Analysis

DNA sequencing of PCR products was performed using an ABI 3730 XL DNA analyzer at Bioneer (Daejeon, Korea) with LepF1 and LepR1 primers. The COI sequences were edited and analysed using MEGA 6.0 (Molecular Evolutionary Genetics Analysis) software [40]. Multiple alignments of all nucleotide sequences were done using ClustalW, available within MEGA 6.0. All sequences obtained in this study were 100% identical. A single representative sequence was deposited in GenBank under the accession number MN883828. To confirm the species identification, the representative sequence was blasted against the GenBank database using the nucleotide BLAST (available at https://blast.ncbi.nlm.nih.gov/Blast.cgi) and compared with the Barcode of Life Data Systems (BOLD) using the BOLD Identification System (IDS) (available at http://v4.boldsystems.org/index.php/IDS_OpenIdEngine). A phylogenetic tree based on neighbour joining (NJ) was constructed with the Kimura two-parameter (K2P) model and bootstrapping (1000 replicates) in MEGA 6.0. The sequences of stomoxyine flies available in GenBank were included in the analysis (Table 1). Musca domestica Linnaeus, 1758 (Diptera: Muscidae) KM571920 was selected as an outgroup.

2.4. Geometric Morphometric Analysis

2.4.1. Wing Preparation

The left wing of males and females of H. aberrans was dissected from the body and mounted with Hoyer’s medium on microscope slides. All slides were photographed using a digital camera connected to a stereomicroscope (Nikon AZ 100, Nikon Corp, Tokyo, Japan). A total of 16 wing pictures (9 males and 7 females) were analysed using the landmark-based method.

2.4.2. Landmark-Based Method

The coordinates of ten wing landmarks (Figure 2) based on the study of Changbunjong et al. (2016) [26] were digitised for geometric morphometric analysis. The wing size was estimated using the centroid size (CS) derived from coordinates of all landmarks. The centroid size is defined as the square root of the sum of the squared distances between the centre of the configuration of landmarks and each separate landmark [47]; it was computed after Generalised Procrustes Analysis (GPA), also known as Procrustes superimposition [48]. The centroid size difference was compared between sexes by non-parametric tests based on 1000 permutations, with Bonferroni correction for a test of significance at p-value < 0.05. The wing shape variables were computed as principal components of the partial warp (PW) scores or relative warps (RWs); they were also computed after Generalised Procrustes Analysis. The Procrustes superimposition provided configurations for visual comparison of the mean anatomical landmarks between the sexes. The relative warps were used to explore the morphospace and used as an input for the discriminant analysis (DA) or canonical variate analysis (CVA), which was illustrated by the factor map. Mahalanobis distance between the sexes was calculated from DA analysis. The statistical significance of shape differences between sexes was obtained by non-parametric analyses based on 1000 permutations, with Bonferroni correction for a test of significance at p-value < 0.05. A cross-validated classification was used to test the accuracy of classification between sexes based on Mahalanobis distances. The allometric effect, the effect of size on shape variation, was performed by linear regression of the Procrustes shape coordinates on centroid size and then estimated by the determination coefficient and statistical significance of the model at p-value < 0.05.

2.4.3. Morphometric Software

Geometric and multivariate analyses were performed using the CLIC package version 99 [25], freely available at https://xyom-clic.eu.

3. Results

A total of 16 specimens of H. aberrans, comprising 9 males (56.25%) and 7 females (43.75%), were collected in this study. For genetic analysis, only 10 specimens (5 males, 5 females) were used in the current study. A standardised 658 bp fragment of the mitochondrial COI gene was successfully amplified from all specimens. No stop codons, insertions or deletions were observed in these sequences. All sequences were 100% identical. The sequence had high A + T content (69.3%), with nucleotide composition of A = 30.2%, T = 39.1%, C = 14.4% and G = 16.3%. The sequence showed 92.52% similarity to Neivamyia flavicornis (Malloch, 1928) in GenBank but did not match any identified species in BOLD. The sequence of H. aberrans (MN883828) was used in the phylogenetic analysis. The phylogenetic relationship of the COI sequence of H. aberrans and other related species was established using an NJ tree (Figure 3). Haematobosca aberrans was clearly separated from other species within the subfamily Stomoxyinae.
For geometric morphometric analysis, a total of 16 wing pictures (9 males and 7 females) were analysed using the landmark-based approach. The wing size (centroid size) variation between sexes of H. aberrans is illustrated in Figure 4. The mean wing CS of males and females was 2.98 ± 0.13 mm and 3.05 ± 0.22 mm, respectively. There was no statistically significant difference between the sexes. The visual comparison of the mean anatomical landmark positions between the sexes showed most visible landmark displacements in the upper and lower parts of the wing (landmarks 1, 7, 8, 9, 10) (Figure 5). Morphospace based on wing shape variables from principal components and the factor map of the discriminant factors showed a clear difference between the sexes (Figure 6 and Figure 7). Based on the Mahalanobis distance, the wing shape was significantly different between the sexes of H. aberrans (Mahalanobis distance = 7.17). The accuracy score after a cross-validated classification test was 100% in both males and females. The wing shape has no allometric effect. The determination coefficient was 10% but was not statistically significant.

4. Discussion

Our study has revealed updated data about stomoxyine flies in Thailand. The new species Haematobosca aberrans was described recently in terms of its morphological characteristics [6] and was found in the Chiang Mai Province of northern Thailand. However, this species was apparently previously recorded in the Kanchanaburi Province of western Thailand but under the name Stygeromyia sp. [7]. Buffaloes may be important hosts of this fly in the studied area. However, buffaloes have also been reported to be hosts that can be attacked by the species H. stimulans [1], while other hosts of Haematobosca spp. include cattle, horses, wildlife and humans [1,2,3,5,49]. In the present study, the habitat of the buffalo farm consisted of a human settlement and of forest, which may provide suitable breeding sites for these flies. In the genus Haematobosca, only H. sanguinolenta had been reported in Thailand. This species was found in livestock farms, wildlife conservation areas and national parks [5]. In our present study, we found that H. sanguinolenta and H. aberrans could be found together in the same habitat of the livestock farm. These results suggest that livestock farms are sites with the greatest diversity of stomoxyine species [5].
The results obtained from this study showed that the alignment of COI sequences of H. aberrans was straightforward and these sequences lacking stop codons, insertions or deletions indicated that all sequences were functional mitochondrial COI sequences. The nucleotide sequences showed a high A + T content (69.3%), which is consistent with previous studies on stomoxyine flies (69%) [24]. One hundred percent similarity could be detected between the different H. aberrans sequences in the present study. This may be due to the low number of samples and/or to the fact that all the flies were collected from the one population. Comparison between our sequence and reference sequences in GenBank and BOLD databases can be used to confirm a new species of H. aberrans as previously described by Pont et al. (2020) [4]. The phylogenetics based on the NJ tree indicated that the COI gene could distinguish among stomoxyine species.
There is a correlation between the morphological characters already described [6] and the molecular results presented in our study. H. aberrans could be separated from other species by morphological characters. However, H. aberrans and N. flavicornis were clustered together in the phylogenetic tree but they received low bootstrap support. Based on the morphological characters, H. aberrans and N. flavicornis are easily distinguished by their strap-like palpi, which are not grooved internally in N. flavicornis but are more dilated terminally and grooved internally in H. aberrans. In addition, the anterior katepisternal seta is absent in H. aberrans but present in N. flavicornis [1,3]. Within the genus Haematobosca, the main morphological character distinguishing H. aberrans from other species is the absence of the anterior katepisternal seta [6]. According to the key by Zumpt (1973) [1], the species H. aberrans resembles species of the genus Stygeromyia Austen, 1907 because the anterior katepisternal seta is absent. However, H. aberrans has long hairs on the upper and lower sides of the arista, as in other species of Haematobosca, and this differs from Stygeromyia where there are long hairs only on the upper side. In addition, the meron and prosternum are bare in H. aberrans and other Haematobosca but are setulose in Stygeromyia [6]. So far as the genetics are concerned, there is little information on Haematobosca and other stomoxyine species, especially Stygeromyia spp., in the GenBank database that is available for comparison with our new sequences.
The geometric morphometric study was carried out to compare the wing size and wing shape between sexes of H. aberrans. Landmark-based geometric morphometrics has been mostly used in sexual dimorphism studies [32,33,36]. In the present study, we revealed a sexual dimorphism in wing shape, indicating that the phenotypic expression of wing shape is a sex-specific difference. This is consistent with a previous report on stomoxyine flies: Stomoxys indicus, S. pullus and S. uruma [26]. The sexual dimorphism of wing shape has been reported in other insect vectors such as mosquitoes of the genus Aedes (Ae. aegypti (Linnaeus, 1762), Ae. albopictus (Skuse, 1894)), Anopheles (An. albitarsis Lynch Arribalzaga, 1878, An. cruzii Dyar & Knab, 1908, An. homunculus Komp, 1937, An. strodei Root, 1926, An. triannulatus Komp, 1937), Culex (Cx. quinquefasciatus Say, 1823, Cx. nigripalpus Theobald, 1901) and Ochlerotatus (Oc. scapularis (Rondani, 1848)) [34]. Furthermore, sexual dimorphism of both wing size and wing shape has also been found in some mosquito species: Mansonia bonneae Edwards, 1930 and M. dives (Schiner, 1868) [33]. The sexual size dimorphism may be mainly affected by environmental factors such as food quantity/quality or temperature [50], while sexual shape dimorphism is often a passive consequence of sexual size dimorphism (allometry) [33]. The sexual shape dimorphism in our study is not correlated with size (non-allometric effect). However, non-allometric effects were probably playing a role in the sexual shape dimorphism of H. aberrans. According to Gidaszewski et al. (2009) [37], flight behaviour and mating systems may be related to variations in wing shape.

5. Conclusions

This study has provided genetic and phenotypic information on H. aberrans in Thailand. The COI gene can be used to differentiate Haematobosca spp. and other related species within the subfamily Stomoxyinae. It has also confirmed that H. aberrans is distinct and different from the other known species. Landmark-based geometric morphometrics is an effective way of separating between the sexes of this species. The sexual shape and size dimorphism could be a basis for further studies on the influence of sex-specific divergence in response to environmental conditions and the species identification of Haematobosca spp. Due to the small number of specimens in this study, further studies including more specimens are recommended to increase the reliability of wing shape for sexual discrimination. Also, behavioral and ecological studies are required to confirm the existence of sexual dimorphism. Additionally, a more comprehensive phylogeny with various gene markers, different populations of the same species, further species from each genus, and other related species will provide additional information on genetic relationships within the Stomoxyinae.

Author Contributions

Conceptualisation, T.C.; methodology, T.C.; validation, J.R., G.D., A.C.P.; investigation, T.C.; resources, T.C.; data curation, T.C.; writing—original draft preparation, T.C.; writing—review and editing, T.C., J.R., G.D., A.C.P.; project administration, T.C.; funding acquisition, T.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the Faculty of Veterinary Science, Mahidol University (0517.131/00192) and partially supported by the Thailand Research Fund and the Office of the Higher Education Commission (grant no. MRG6180193).

Acknowledgments

We thank the Monitoring and Surveillance Center for Zoonotic Diseases in Wildlife and Exotic Animals, Faculty of Veterinary Science, Mahidol University, for supplying specimens.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zumpt, F. The Stomoxyine Biting Flies of the World; Gustav Fisher Verlag: Stuttgart, Germany, 1973; pp. 44–71. [Google Scholar]
  2. Burger, J.F.; Anderson, J.R. Taxonomy and life history of the Moose Fly, Haematobosca alcis, and its association with the moose, Alces alces shirasi in Yellowstone National Park. Ann. Ent. Soc. Am. 1974, 67, 204–214. [Google Scholar] [CrossRef]
  3. Pont, A.C.; Mihok, S. A new species of Haematobosca Bezzi from Kenya (Diptera, Muscidae). Stud. Dipterol. 2000, 7, 25–32. [Google Scholar]
  4. Pont, A.C.; Dsouli, N. A new species of Haematobosca Bezzi from Gabon (Diptera, Muscidae). Stud. Dipterol 2009, 15, 259–266. [Google Scholar]
  5. Changbunjong, T.; Weluwanarak, T.; Ratanakorn, P.; Maneeon, P.; Ganpanakngan, M.; Apiwathnasorn, C.; Sungvornyothin, S.; Sriwichai, P.; Sumruayphol, S.; Ruangsittichai, J. Distribution and abundance of Stomoxyini flies (Diptera: Muscidae) in Thailand. Southeast Asian J. Trop. Med. Public Health 2012, 43, 1400–1410. [Google Scholar] [PubMed]
  6. Pont, A.C.; Duvallet, G.; Changbunjong, T. A new species of Haematobosca Bezzi (Diptera: Muscidae) from Thailand. Zootaxa 2020, 4763, 538–544. [Google Scholar] [CrossRef]
  7. Tumrasvin, W.; Shinonaga, S. Studies on medically important flies in Thailand. V. On 32 species belonging to the subfamilies Muscinae and Stomoxyini including the taxonomic keys (Diptera: Muscidae). Bull. Tokyo Med. Dent. Univ. 1978, 25, 201–227. [Google Scholar]
  8. Hebert, P.D.N.; Gregory, T.R. The Promise of DNA Barcoding for Taxonomy. Syst. Biol. 2005, 54, 852–859. [Google Scholar] [CrossRef]
  9. Meusnier, I.; Singer, G.A.; Landry, J.; Hickey, D.A.; Hebert, P.D.N.; Hajibabaei, M. A 683 universal DNA mini-barcode for biodiversity analysis. BMC Genom. 2008, 9, 214. [Google Scholar] [CrossRef] [Green Version]
  10. Pramual, P.; Adler, P.H. DNA barcoding of tropical black flies (Diptera: Simuliidae) of Thailand. Mol. Ecol. Resour. 2014, 14, 262–271. [Google Scholar] [CrossRef]
  11. Onder, Z.; Yildirim, A.; Duzlu, O.; Arslan, M.O.; Sari, B.; Tasci, G.T.; Ciloglu, A.; Aydin, N.P.; Inci, A.; Adler, P.H. Molecular characterization of black flies (Diptera: Simuliidae) in areas with pest outbreaks and simuliotoxicosis in Northeast Anatolia Region, Turkey. Acta Trop. 2019, 199, 105149. [Google Scholar] [CrossRef]
  12. Cywinska, A.; Hunter, F.F.; Hebert, P.D. Identifying Canadian mosquito species through DNA barcodes. Med. Vet. Entomol. 2006, 20, 413–424. [Google Scholar] [CrossRef] [PubMed]
  13. Kumar, N.P.; Rajavel, A.R.; Natarajan, R.; Jambulingam, P. DNA barcodes can distinguish species of Indian mosquitoes (Diptera: Culicidae). J. Med. Entomol. 2007, 44, 1–7. [Google Scholar] [CrossRef] [PubMed]
  14. Sumruayphol, S.; Apiwathnasorn, C.; Ruangsittichai, J.; Sriwichai, P.; Attrapadung, S.; Samung, Y.; Dujardin, J.P. DNA barcoding and wing morphometrics to distinguish three Aedes vectors in Thailand. Acta Trop. 2016, 159, 1–10. [Google Scholar] [CrossRef] [PubMed]
  15. Vadivalagan, C.; Karthika, P.; Murugan, K.; Panneerselvam, C.; Del Serrone, P.; Benelli, G. Exploring genetic variation in haplotypes of the filariasis vector Culex quinquefasciatus (Diptera: Culicidae) through DNA barcoding. Acta Trop. 2017, 169, 43–50. [Google Scholar] [CrossRef]
  16. Contreras Gutiérrez, M.A.; Vivero, R.J.; Vélez, I.D.; Porter, C.H.; Uribe, S. DNA barcoding for the identification of sand fly species (Diptera, Psychodidae, Phlebotomine) in Colombia. PLoS ONE 2014, 9, e85496. [Google Scholar] [CrossRef]
  17. Nzelu, C.O.; Cáceres, A.G.; Arrunátegui-Jiménez, M.J.; Lañas-Rosas, M.F.; Yañez-Trujillano, H.H.; Luna-Caipo, D.V.; Holguín-Mauricci, C.E.; Katakura, K.; Hashiguchi, Y.; Kato, H. DNA barcoding for identification of sand fly species (Diptera: Psychodidae) from leishmaniasis-endemic areas of Peru. Acta Trop. 2015, 145, 45–51. [Google Scholar] [CrossRef]
  18. Romero-Ricardo, L.; Lastre-Meza, N.; Pérez-Doria, A.; Bejarano, E.E. DNA barcoding to identify species of phlebotomine sand fly (Diptera: Psychodidae) in the mixed leishmaniasis focus of the Colombian Caribbean. Acta Trop. 2016, 159, 125–131. [Google Scholar] [CrossRef]
  19. Lozano-Sardaneta, Y.N.; Paternina, L.E.; Sánchez-Montes, S.; Quintero, A.; Ibáñez-Bernal, S.; Sánchez-Cordero, V.; Bejarano, E.E.; Becker, I. DNA barcoding and fauna of phlebotomine sand flies (Diptera: Psychodidae: Phlebotominae) from Los Tuxtlas, Veracruz, Mexico. Acta Trop. 2020, 201, 105220. [Google Scholar] [CrossRef]
  20. Morita, S.I.; Bayless, K.M.; Yeates, D.K.; Wiegmann, B.M. Molecular phylogeny of the horse flies: A framework for renewing tabanid taxonomy. Syst. Entomol. 2016, 41, 56–72. [Google Scholar] [CrossRef]
  21. Nitiyamatawat, E.; Changbunjong, T.; Sumruayphol, S.; Sriwichai, P.; Thaenkham, U.; Ruangsittichai, J. DNA barcoding for identification of cleg flies, Haematopota spp. (Diptera: Tabanidae) in Thailand. Thai J. Vet. Med. Suppl. 2017, 47, 303–305. [Google Scholar]
  22. Changbunjong, T.; Bhusri, B.; Sedwisai, P.; Weluwanarak, T.; Nitiyamatawat, E.; Chareonviriyaphap, T.; Ruangsittichai, J. Species identification of horse flies (Diptera: Tabanidae) in Thailand using DNA barcoding. Vet. Parasitol. 2018, 259, 35–43. [Google Scholar] [CrossRef]
  23. Changbunjong, T.; Weluwanarak, T.; Sedwisai, P.; Ruangsittichai, J.; Duvallet, G.; Chareonviriyaphap, T. New records and DNA barcoding of deer flies, Chrysops (Diptera: Tabanidae) in Thailand. Acta Trop. 2020, 210, 105532. [Google Scholar] [CrossRef] [PubMed]
  24. Changbunjong, T.; Weluwanarak, T.; Samung, Y.; Ruangsittichai, J. Molecular identification and genetic variation of stomoxyine flies (Diptera: Muscidae) in Thailand based on cox1 barcode. J. Asia Pac. Entomol. 2016, 19, 1117–1123. [Google Scholar] [CrossRef]
  25. Dujardin, J.P. Morphometrics applied to medical entomology. Infect. Genet. Evol. 2008, 8, 875–890. [Google Scholar] [CrossRef]
  26. Changbunjong, T.; Sumruayphol, S.; Weluwanarak, T.; Ruangsittichai, J.; Dujardin, J.P. Landmark and outline-based geometric morphometrics analysis of three Stomoxys flies (Diptera: Muscidae). Folia Parasitol. (Praha) 2016, 63, 37. [Google Scholar]
  27. Lorenz, C.; Almeida, F.; Almeida-Lopes, F.; Louise, C.; Pereira, S.N.; Petersen, V.; Vidal, P.O.; Virginio, F.; Suesdek, L. Geometric morphometrics in mosquitoes: What has been measured? Infect. Genet. Evol. 2017, 54, 205–215. [Google Scholar] [CrossRef]
  28. Sontigun, N.; Sukontason, K.L.; Zajac, B.K.; Zehner, R.; Sukontason, K.; Wannasan, A.; Amendt, J. Wing morphometrics as a tool in species identification of forensically important blow flies of Thailand. Parasit. Vectors 2017, 10, 229. [Google Scholar] [PubMed] [Green Version]
  29. Chaiphongpachara, T. Comparison of Landmark- and Outline-Based Geometric Morphometrics for Discriminating Mosquito Vectors in Ratchaburi Province, Thailand. Biomed. Res. Int. 2018, 2018, 6170502. [Google Scholar] [CrossRef] [PubMed]
  30. Chaiphongpachara, T.; Sriwichai, P.; Samung, Y.; Ruangsittichai, J.; Morales Vargas, R.E.; Cui, L.; Sattabongkot, J.; Dujardin, J.P.; Sumruayphol, S. Geometric morphometrics approach towards discrimination of three member species of Maculatus group in Thailand. Acta Trop. 2019, 192, 66–74. [Google Scholar] [CrossRef] [PubMed]
  31. Dujardin, J.P.; Kaba, D.; Solano, P.; Dupraz, M.; McCoy, K.D.; Jaramillo-O, N. Outline-based morphometrics, an overlooked method in arthropod studies? Infect. Genet. Evol. 2014, 28, 704–714. [Google Scholar]
  32. Chaiphongpachara, T.; Laojun, S. Using the modern morphometric approach to determine sexual dimorphism of three medically important flies (Order: Diptera) in Thailand. Biodiversitas 2019, 20, 1482–1486. [Google Scholar] [CrossRef]
  33. Ruangsittichai, J.; Apiwathnasorn, C.; Dujardin, J.P. Interspecific and sexual shape variation in the filariasis vectors Mansonia dives and Ma. bonneae. Infect. Genet. Evol. 2011, 11, 2089–2094. [Google Scholar] [CrossRef] [PubMed]
  34. Virginio, F.; Oliveira Vidal, P.; Suesdek, L. Wing sexual dimorphism of pathogen-vector culicids. Parasit. Vectors 2015, 8, 159. [Google Scholar] [CrossRef] [PubMed]
  35. de Camargo, W.R.; de Camargo, N.F.; Corrêa, D.; de Camargo, A.J.; Diniz, I.R. Sexual Dimorphism and Allometric Effects Associated with the Wing Shape of Seven Moth Species of Sphingidae (Lepidoptera: Bombycoidea). J. Insect Sci. 2015, 15, 107. [Google Scholar] [CrossRef] [PubMed]
  36. Benítez, H.A.; Vargas, H.A. Sexual dimorphism and population differentiation in the Chilean Neotropical moth Macaria mirthae (Lepidoptera, Geometridae): A wing geometric morphometric example. Rev. Bras. Entomol. 2017, 61, 365–369. [Google Scholar] [CrossRef]
  37. Gidaszewski, N.A.; Baylac, M.; Klingenberg, C.P. Evolution of sexual dimorphism of wing shape in the Drosophila melanogaster subgroup. BMC Evol. Biol. 2009, 9, 110. [Google Scholar] [CrossRef]
  38. Mihok, S. The development of a multipurpose trap (the Nzi) for tsetse and other biting flies. Bull. Entomol. Res. 2002, 92, 385–403. [Google Scholar] [CrossRef]
  39. Hebert, P.D.N.; Penton, E.H.; Burns, J.M.; Janzen, D.H.; Hallwachs, W. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc. Natl. Acad. Sci. USA 2004, 101, 14812–14817. [Google Scholar] [CrossRef]
  40. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  41. Iwasa, M.; Ishiguro, N. Genetic and morphological differences of Haematobia irritans and H. exigua and molecular phylogeny of Japanese Stomoxyini flies (Diptera, Muscidae). Med. Entomol. Zool. 2010, 61, 335–344. [Google Scholar] [CrossRef]
  42. Renaud, A.K.; Savage, J.; Adamowicz, S.J. DNA barcoding of Northern Nearctic Muscidae (Diptera) reveals high correspondence between morphological and molecular species limits. BMC Ecol. 2012, 12, 24. [Google Scholar] [CrossRef] [PubMed]
  43. Hebert, P.D.; Ratnasingham, S.; Zakharov, E.V.; Telfer, A.C.; Levesque-Beaudin, V.; Milton, M.A.; Pedersen, S.; Jannetta, P.; deWaard, J.R. Counting animal species with DNA barcodes: Canadian insects. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2016, 371, 20150333. [Google Scholar] [CrossRef] [Green Version]
  44. Sikes, D.S.; Bowser, M.; Morton, J.M.; Bickford, C.; Meierotto, S.; Hildebrandt, K. Building a DNA barcode library of Alaska’s non-marine arthropods. Genome 2017, 60, 248–259. [Google Scholar] [CrossRef] [PubMed]
  45. Kutty, S.N.; Pape, T.; Pont, A.; Wiegmann, B.M.; Meier, R. The Muscoidea (Diptera: Calyptratae) are paraphyletic: Evidence from four mitochondrial and four nuclear genes. Mol. Phylogenet. Evol. 2008, 49, 639–652. [Google Scholar] [CrossRef] [PubMed]
  46. Haseyama, K.L.; Wiegmann, B.M.; Almeida, E.A.; de Carvalho, C.J. Say goodbye to tribes in the new house fly classification: A new molecular phylogenetic analysis and an updated biogeographical narrative for the Muscidae (Diptera). Mol. Phylogenet. Evol. 2015, 89, 1–12. [Google Scholar] [CrossRef]
  47. Bookstein, F.L. Morphometric Tools for Landmark Data: Geometry and Biology, 1st ed.; Cambridge University Press: New York, NY, USA, 1991; 435p. [Google Scholar]
  48. Rohlf, F.J. Rotational fit (Procrustes) methods. In Proceedings of the Michigan Morphometrics Workshop; Rohlf, F.J., Bookstein, F.L., Eds.; The University of Michigan Museum of Zoology: Ann Arbor, MI, USA, 1990; pp. 227–236. [Google Scholar]
  49. Changbunjong, T.; Weluwanarak, T.; Sedwisai, P.; Chamsai, T. Stomoxyini fly fauna of the Khao Yai National Park, Thailand. Asian Pac. J. Trop. Dis. 2013, 3, 348–351. [Google Scholar] [CrossRef] [Green Version]
  50. Stillwell, R.C.; Blanckenhorn, W.U.; Teder, T.; Davidowitz, G.; Fox, C.W. Sex differences in phenotypic plasticity affect variation in sexual size dimorphism in insects: From physiology to evolution. Annu. Rev. Entomol. 2010, 55, 227–245. [Google Scholar] [CrossRef]
Insects 11 00451 g001 550
Figure 1. Map showing the collection site of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 in Thailand.
Figure 1. Map showing the collection site of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 in Thailand.
Insects 11 00451 g001
Insects 11 00451 g002 550
Figure 2. Ten landmarks digitised on wings of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020.
Figure 2. Ten landmarks digitised on wings of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020.
Insects 11 00451 g002
Insects 11 00451 g003 550
Figure 3. Neighbour-joining (NJ) tree based on Kimura two-parameter (K2P) model for cytochrome c oxidase I (COI) sequence of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 collected from northern Thailand and sequences obtained from GenBank. Only bootstrap values ≥50% are shown. The scale bar represents 0.01% divergence.
Figure 3. Neighbour-joining (NJ) tree based on Kimura two-parameter (K2P) model for cytochrome c oxidase I (COI) sequence of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 collected from northern Thailand and sequences obtained from GenBank. Only bootstrap values ≥50% are shown. The scale bar represents 0.01% divergence.
Insects 11 00451 g003
Insects 11 00451 g004 550
Figure 4. Centroid size variation of the wings between sexes of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020, shown as quartile boxes. Each box shows the group median separating the 25th and 75th quartiles. Vertical bars under the boxes represent the wing (units as mm).
Figure 4. Centroid size variation of the wings between sexes of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020, shown as quartile boxes. Each box shows the group median separating the 25th and 75th quartiles. Vertical bars under the boxes represent the wing (units as mm).
Insects 11 00451 g004
Insects 11 00451 g005 550
Figure 5. Configuration of the mean anatomical landmark positions connected by a straight line after Procrustes superimposition of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 in both sexes.
Figure 5. Configuration of the mean anatomical landmark positions connected by a straight line after Procrustes superimposition of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 in both sexes.
Insects 11 00451 g005
Insects 11 00451 g006 550
Figure 6. Morphospace of males and females of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 based on wing shape variables. The horizontal axis is the first relative warp (RW1), and the vertical axis is the second relative warp (RW2).
Figure 6. Morphospace of males and females of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020 based on wing shape variables. The horizontal axis is the first relative warp (RW1), and the vertical axis is the second relative warp (RW2).
Insects 11 00451 g006
Insects 11 00451 g007 550
Figure 7. Factor map of the discriminant factor separating males (black bar) and females (grey bar) of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020.
Figure 7. Factor map of the discriminant factor separating males (black bar) and females (grey bar) of Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020.
Insects 11 00451 g007
Table
Table 1. Stomoxyine flies from the GenBank database used for constructing the phylogenetic tree.
Table 1. Stomoxyine flies from the GenBank database used for constructing the phylogenetic tree.
SpeciesOriginAccession NumberReference
Haematobia irritans irritans (Linnaeus, 1758)JapanAB479523[41]
Haematobia irritans exigua de Meijere, 1906ThailandKC960695[24]
Haematobosca alcis (Snow, 1891)CanadaHM389222, KR589429 [42,43]
USAKU496750[44]
Haematobosca sanguinolenta (Austen, 1909)ThailandKC960696, KC960703[24]
VietnamAB479536[41]
Haematobosca stimulans (Meigen, 1824)NAFJ025615[45]
Haematostoma austeni Malloch, 1932ThailandKC960710[24]
Neivamyia flavicornis (Malloch, 1928)BrazilKP161712[46]
Rhinomusca dutoiti Zumpt, 1950South AfricaKP161721[46]
Stomoxys calcitrans (Linnaeus, 1758)ThailandKC960724[24]
Stomoxys indicus Picard, 1908ThailandKC960735[24]
Stomoxys pullus Austen, 1909ThailandKC960747[24]
Stomoxys sitiens Rondani, 1873ThailandKC960751[24]
Stomoxys uruma Shinonaga & Kano, 1966ThailandKC960753[24]
NA, no available data.

Share and Cite

MDPI and ACS Style

Changbunjong, T.; Ruangsittichai, J.; Duvallet, G.; Pont, A.C. Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae). Insects 2020, 11, 451. https://doi.org/10.3390/insects11070451

AMA Style

Changbunjong T, Ruangsittichai J, Duvallet G, Pont AC. Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae). Insects. 2020; 11(7):451. https://doi.org/10.3390/insects11070451

Chicago/Turabian Style

Changbunjong, Tanasak, Jiraporn Ruangsittichai, Gerard Duvallet, and Adrian C. Pont. 2020. "Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae)" Insects 11, no. 7: 451. https://doi.org/10.3390/insects11070451

APA Style

Changbunjong, T., Ruangsittichai, J., Duvallet, G., & Pont, A. C. (2020). Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae). Insects, 11(7), 451. https://doi.org/10.3390/insects11070451

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop