Wing sexual dimorphism of pathogen-vector culicids
Parasites & Vectorsvolume 8, Article number: 159 (2015)
Sexual dimorphism in animals has been studied from different perspectives for decades. In 1874 Darwin hypothesized that it was related to sexual selection, and even after nearly 140 years, when additional empirical data has become available and the subject has been investigated from a contemporary viewpoint, this idea is still supported. Although mosquito (Culicidae) wings are of great importance as they play a sex-specific role, little is known about wing sexual dimorphism in these pathogen-vector insects. Detection and characterization of wing sexual dimorphism in culicids may indirectly enhance our knowledge of their epidemiology or reveal sex-linked genes, aspects that have been discussed by vector control initiatives and developers of genetically modified mosquitoes.
Using geometric morphometrics, we carried out a comparative assessment of wing sexual dimorphism in ten culicid species of medical/veterinary importance from genera Culex, Aedes, Anopheles and Ochlerotatus collected in Brazil.
Discriminant analysis revealed significant sexual dimorphism in all the species studied, indicating that phenotypic expression of wing shape in mosquitoes is indeed sex-specific. A cross-validated test performed to reclassify the sexes with and without allometry yielded very similar results. Mahalanobis distances among the ten species showed that the species had different patterns of shape sexual dimorphism and that females are larger than males in some species.
Wing morphology differed significantly between species. The finding of sexual dimorphism in all the species would suggest that the wing geometry of Culicidae is canalized. Although sexual dimorphism is prevalent, species-specific patterns occur. Allometry was not the main determinant of sexual dimorphism, which suggests that sexual selection or other evolutionary mechanisms underlie wing sexual dimorphism in these insects.
Typically, morphological sexual dimorphism (SD) is dichotomically classified into size and shape dimorphism; however, there is also known to be interaction between body size and shape (allometric SD). Although in some cases SD appears to be the result of allometry , examples of SD that can be attributed purely to size  or shape  have also been reported. For decades, SD in animals has been investigated from either an evolutionary or ecological point of view. About 140 years ago, Darwin  hypothesized that SD might be related to sexual selection, i.e., the process whereby the maintenance of sex-specific traits is driven by the preference of the opposite sex. Contemporary views of SD take into account empirical data and new constraints and concepts but preserve the underlying idea behind Darwin’s hypothesis .
Insects have been used as models for investigating SD as sex-related differences can occur in various organs, such as the eyes , legs [7,8] and head capsule . Other structures and features, such as hairs on the antennae, body size, mouthparts, genitalia  and wings [11,12] can also be observed and help with sexing. Insect wings may have evolved through selective action associated with mating behavior [13,14]. Sexual behavior, such as courtship songs, can depend on sex-specific wing geometry .
Few studies have been published to date on wing SD in Culicidae, although this family includes insects which are vectors of etiologic agents of serious human diseases. Wing SD in this family has been described occasionally in taxonomic keys, but its specific patterns and variability are not well known. Wing sexual size dimorphism (SSD) and sexual shape dimorphism (SShD) have already been reported in Oc. scapularis , but a comparative perspective has yet to be explored. Investigation of wing SD in culicids may indirectly enrich our knowledge of their epidemiology given that both the epidemiological relevance of and role played by wings are particular to each sex. Only females are hematophagous and pathogen-competent and use their wings to ensure an accurate approach to other animals and suck their blood, while males can copulate with several mates and use wing beats to attract the opposite sex during courtship. Wing SD merits study because it may underlie sex-linked genetic markers, which are especially useful for the development of genetically-modified mosquitoes [16,17].
An initial approach to the study of wing SD in Culicidae would be to detect and quantify wing SD in species that are representative of the main taxonomic subgroups. This would yield information about the frequency of SD and subsequently provide insights into the role of phylogenetic constraints and/or species-specific adaptation in the evolution of SD. Nowadays, investigation of wing SD in many samples is more feasible than in the past because of the availability of geometric morphometrics, a cheap and highly accurate technique that has become increasingly popular [11,18-33].
We decided to investigate whether male and female culicids have different wing characteristics, a hypothesis that was proposed following a case study of Oc. scapularis . We also hypothesized that different species of Culicidae have different degrees of SD, as has been reported for Drosophilidae . To this end we investigated ten Culicidae species of medical/veterinary importance from the two main subfamilies of Culicidae using geometric morphometric analysis of wing shape and wing size separately.
Collection of the biological samples
Mosquitoes of the genus Aedes (Aedes aegypti, Aedes albopictus), Anopheles (Anopheles albitarsis l.s., Anopheles cruzii, Anopheles homunculus, Anopheles strodei l.s., Anopheles triannulatus l.s.), Culex (Culex quinquefasciatus, Culex nigripalpus) and Ochlerotatus (Ochlerotatus scapularis) were collected between 2007 and 2013 in the States of São Paulo and Minas Gerais (Table 1). The immature stages (eggs, larvae and pupae) were collected from artificial and natural breeding sites and maintained in the laboratory under standard temperature and humidity conditions (25 ± 1°C; 80 ± 10%) with a photoperiod of 12:12 (light:dark) until the emergence of adult mosquitoes. Adult mosquitoes were collected with entomological aspirators and white Shannon-type traps. All individuals were identified at the species level [10,34-37] and stored in 70% ethanol.
Wings were detached from the thorax of each individual (both males and females) and mounted between a slide and coverslip as described by Lorenz et al. . Images of the wings were captured by a Leica DFC320 digital camera coupled to a Leica S6 stereoscope at 40X magnification.
The coordinates of eighteen landmarks at the vein intersections on these images were digitized using TpsDig v. 1.4 software . This set of landmarks has proven to be sufficiently sensitive to describe SD  (Figure 1).
Most of the analyses were performed with right wings; left wings were only used when the corresponding right wing was damaged.
To analyze the contribution of size and shape to SD separately, we removed the allometric effect in all the analyses. Wing shape and size SD was determined in the different species. Morphological similarity among the ten species sampled was also examined. This was done separately for males and females.
Variations in wing shape (partial warps) were determined by Procrustes superimpositions through generalized least squares, eliminating the differences in orientation, position and isometric size . Multivariate regression of Procrustes coordinates versus centroid size (CS) was used to remove allometric effects, and a permutation test with 10,000 randomizations was applied using MorphoJ software  to test the significance of the allometry. Differences in wing shape SD were determined by canonical variate analysis using MorphoJ software. Wing shape dissimilarity between females and males of each species was estimated by Mahalanobis distances (MD) and compared in a permutation test with 10,000 randomizations using MorphoJ software.
Morphological divergence among the samples was illustrated by UPGMA phenograms constructed using PHYLIP (Phylogeny Inference Package) v.3.6 distributed by J. Felsenstein , Department of Genome Sciences - University of Washington, Seattle/WA). To test for dissimilarity between males and females, individuals were reclassified according to their similarity to each group using the MD.
The global size of the wing was estimated using CS, the square root of the sum of the squared distances of all landmarks from the centroid. GraphPad InStat v.3.01 (San Diego, CA) was used to compare CS scores between samples. The unpaired T-test was used in populations that had a Gaussian distribution and the Mann-Whitney test in populations that did not. The CS of all the species were represented graphically with the MOG v.82 program .
We first removed the allometric effect from all the analyses as, although the study was a macroevolutionary one, we wanted to analyze shape and size separately. Canonical variate analysis revealed differences in wing shape between all the species analyzed for both males and females, thus we could see that all species have clearly sexually dimorphic wings (Figure 2). MD (without allometry) ranged from 7.34 to 34.06 (p < 0.0001), indicating a great variability in SD between species. Ae. albopictus and Cx. quinquefasciatus were the least and most sexually dimorphic species, respectively.
Table 2 shows the phenetic differentiation between sexes in each species, indicates that allometry was low (<10%) or not significant and did not contribute significantly to SShD (see the MD columns).
Comparison of the extremes of differentiation for males and females revealed major displacement of the landmarks between the sexes (Figure 3). As in Drosophila [1,42-44], the most variable landmarks were in the proximal and distal regions of the radial and medium veins. Females were slightly wider and significantly shorter than males.
The phenograms of the pairwise distances for the ten species (sexes separated) corroborated the SD found in the other analyses. The “Anophelinae - males” cluster was symmetrical to the “Anophelinae - females” cluster, but clusters comprising Culicinae representatives were not congruent between the sexes (Figure 4).
A cross-validated reclassification test based on the MD was performed with an accuracy of 53.3 to 100%. The same analysis was carried out taking into account the allometric effect to identify whether allometry contributed significantly to SD in the species studied. The score in this analysis was 58.8 to 100%. Both analyses (with and without the effects of allometry) yielded very similar results except for the species Cx. quinquefasciatus, An. albitarsis l.s. and Ae aegypti, which had high scores when allometric effects were included. Analysis of CS indicated SSD in all species except An. homunculus (Figure 5).
Wing CS for each species lay in the following ranges (in mm): Ae. albopictus (F: 1.85- 2.58, M: 1.46-1.91), Cx. quinquefasciatus (F: 3.40-3.85, M: 2.64-3.02), Ae. aegypti (F: 1.72-2.92, M: 1.46-1.91), Cx. nigripalpus (F: 2.69-3.40, M: 2.39-2.97), An. triannulatus l.s. (F: 2.30-2.83, M: 2.01-2.71), An. strodei l.s. (F: 2.33-3.23, M: 2.07-3.01), An. albitarsis l.s. (F:2.57-3.14, M:2.30-3.14), Oc. scapularis (F: 2.37-3.15, M: 2.43-2.94), An. homunculus (F: 2.37-2.87, M: 2.45-2.84), An. cruzii (F: 1.51-2.45, M: 1.49-2.57). The species with the lowest wing CS in males were Ae. albopictus and Ae. aegypti, while the corresponding species for females was An. cruzii. The species with the largest wing CS in males and females were An. albitarsis l.s. and Cx. quinquefasciatus, respectively.
The graphical representation of the ratio of mean female CS to mean male CS (Figure 6) shows that females were larger than males in most species. The exceptions were An. homunculus (ratio = 1), for which the analysis failed to reveal SSD, and An. cruzii (ratio = 0.78), for which females were smaller than males.
Our results indicated wing SD in the ten mosquito species analyzed, suggesting that the phenotypic expression of wing shape is sex-specific. This is in accordance with the findings reported by Devicari et al.  in Oc. scapularis.
There was significant SShD in all the taxonomic groups analyzed, showing a very marked pattern of SD in the family Culicidae. There may be some evolutionary canalization of wing shape that keeps wings sexually dimorphic in so many species, including even those that are phylogenetically more distant in the family. Such canalization does not appear to be driven by allometry because the allometric effects in all the species were low (and sometimes not significant) and only exerted a marginal influence on SShD. Because of the low allometric values observed (as previously detected in An. cruzii by Lorenz et al. ), we cannot assert that SD in culicids is only an allometric effect, as some authors have conjectured . It is possible that the ubiquity of SShD is maintained by the different ways in which wings are used by each sex, i.e., the overall wing shape may be adapted to each sex-specific function.
Although, unlike for crickets , we do not know how wing shape determines wing beat dynamics in culicids, this interpretation has a parallel in the evidence of sex-specific selection of wing beats in some Anopheles species [45,46]. We cannot discard the possible existence of some sex-specific non-adaptive developmental constraint behind the supposed canalization.
The incongruence between Culicinae clusters in the phenograms suggests that species of this subfamily have different patterns of sex dimorphism. Accordingly, Culicinae is the most heterogeneous cluster in terms of the degree of SShD and SSD. As Culicinae is also the most macroevolutionarily derived taxon, the phylogeny of SD merits further investigation.
In some species of insects, males can be larger than females but have less body mass . In our study, SSD was common to most species, although it manifested differently. In most of the species, the females were larger, as reported for other groups of insects such as Ophion intricatus , Drosophila melanogaster [48,49], Cx. quinquefasciatus , Ae. albopictus , Stenurella melanura , Scapteriscus acletus and Scapteriscus vicinus . The fact that SSD did not follow the female > male pattern in An. cruzii and was not present in An. homunculus may be the result of evolutionary and environmental factors particular to these species, which are evolutionarily and ecologically closely related . SSD is probably more variable and less canalized than SShD, i.e., size is evolutionarily less stable than shape. This corroborates the findings of Devicari et al.  and confirms the theory put forward by Dujardin  according to which wing size in insects may be plastic and more influenced by non-genetic factors. Curiously, An. cruzii also had low SShD (as well as “inverted” SSD), making it a particularly unusual species compared with the others. We do not yet know the causes of this variability. Although some studies suggest that SSD and SShD can be adaptive [1-3,39], a clear explanation for these types of SD has yet to be proposed.
Although we detected several features of wing SD in Culicidae species here, an explanation for the observed patterns has yet to be formulated. Considering that SD is prevalent and not primarily due to allometry, we believe that complex evolutionary mechanisms are responsible for the maintenance of SD in mosquitoes. The interspecific divergence of SD has also yet to be explained. It has been suggested that sexual selection and mating system are primary forces that direct the evolution of SD in insects , as occurs in Drosophilidae, the best known sister group of Culicidae. Although it has been suggested that ecological SD occurs, it does not appear to be the most important form of SD . However, although many researchers support the idea that the evolution of SD is directed by sexual selection [9,53,54], others believe it may be an effect of natural selection . Chenoweth et al. , in experimental studies of Drosophila serrata, found that SD tends to increase under the pressure of sexual selection and decrease under the pressure of natural selection. Nevertheless, as we have described here many peculiarities of wing SD in culicids, we do not know to what extent the considerations regarding other insects are generalizable to Culicidae.
Establishing answers to the following questions could help clarify whether SShD is driven by sexual selection: Does wing geometry influence wing beat and other sexually selectable wing traits? Which genes are involved in determining wing shape in culicids?
What is the inheritance and expressivity of these genes? Until recently, little was known about wing SD in Culicidae. It is now clear that this subject is an open field for further research that will eventually enrich our knowledge of the biology, evolution and epidemiology of these mosquitoes.
We believe that wing characters provide useful information for the study of SD and are a practical tool for identifying the sex of culicids. Our findings of apparent canalization of SShD in all the species studied and high variability of SSD appear to reflect complex underlying evolutionary factors. Further, more detailed studies of the determinants of SD may provide essential information for an understanding of the biology of mosquitoes and subsequently help improve vector control methods.
Gidaszewski NA, Baylac M, Klingenberg CP. Evolution of sexual dimorphism of wing shape in the Drosophila melanogaster subgroup. BMC Evol Biol. 2009;9: doi:10.1186/1471-2148-9-110
Fairbairn DJ, Blanckenhorn WU, Székely T. Sex, size and gender roles: evolutionary studies of sexual size dimorphism. Oxford Scholarship Online. 2007;doi:10.1093/acprof:oso/9780199208784.001.0001
Benítez HA, Bravi R, Parra LE, Sanzana MJ, Sepulveda-Zuniga E. Allometric and non-allometric patterns in sexual dimorphism discrimination of wing shape in Ophion intricatus: might two male morphotypes coexist. J Insect Sci. 2013;13:143.
Darwin CR. The descent of man, and selection in relation to sex. 2ªth ed. London: John Murray, Albermarle Street; 1874.
Fairbairn DJ. Allometry for sexual size dimorphism: pattern and process in the coevolution of body size in males and females. Annu Rev Ecol Syst. 1997;28:659–87.
Meyer-Rochow VB, Reid WA. Male and female eyes of the Antarctic midge Belgica antarctica (Diptera, Hironomidae) - a scanning electron microscope study. Appl Entomol Zool. 1994;29:439–42.
Adler PH, Adler CRL. Mating behavior and the evolutionary significance of mate guarding in 3 species of crane flies (Diptera: Tipulidae). J Insect Behav. 1991;4(5):619–32.
Eberhard WG. Physical restraint or stimulation? The function(s) of the modified front legs of male Archisepsis diversiformis (Diptera, Sepsidae). J Insect Behav. 2002;15(6):831–50.
Wilkinson GS, Dodson GN. Function and evolution of antlers and eye stalks in flies. Cambridge, UK: Cambridge University Press; 1997.
Forattini OP. Culicidologia médica, vol. 2. 2002
Devicari M, Lopes AR, Suesdek L. Wing sexual dimorphism in Aedes scapularis (Diptera: Culicidae). Biota Neotrop. 2011;11:165–9.
Benítez HA, Parra LE, Sepulveda E, Sanzana MJ. Geometric perspectives of sexual dimorphism in the wing shape of Lepidoptera: the case of Synneuria sp. (Lepidoptera: Geometridae). J Entomol Res Soc. 2011;13:53–60.
Alexander RD, Brown Jr W. Mating behavior and the origin of insect wings. Ann Arbor, Michigan: University of Michigan Museum of Zoology; 1963.
Shevtsova E, Hansson C, Janzen DH, Kjaerandsen J. Stable structural color patterns displayed on transparent insect wings. PNAS Early Edition. 2011;108(2):668–73.
Klingenberg CP, Debat V, Roff DA. Quantitative genetics of shape in cricket wings: developmental integration in a functional structure. Evolution. 2010;64:2935–51.
Sperança MA, Capurro ML. Perspectives in the control of infectious diseases by transgenic mosquitoes in the post-genomic era – a review. Mem Inst Oswaldo Cruz. 2007;102(4):425–33.
Scali C, Catteruccia F, Li Q, Crisanti A. Identification of sex-specific transcripts of the Anopheles gambiae double sex gene. J Exp Biol. 2005;208:3701–9.
Dujardin JP. Morphometrics applied to medical entomology. Infect Genet Evol. 2008;8:875–90.
Henry A, Thongsripong P, Fonseca-Gonzalez I, Jaramillo-Ocampo N, Dujardin JP. Wing shape of dengue vectors from around the world. Infect Genet Evol. 2010;10(2):207–14.
N. J, Leemingsawat S, Thongrungkiat S, Apiwathnasorn C, Singhaniyom S, Bellec C, et al. Influence of larval density or food variation on the geometry of the wing of Aedes (Stegomyia) aegypti. Trop Med Int Health. 2007;12:1354–60.
Kaba D, Ravel S, Acapovi-Yao G, Solano P, Allou K, Bosson-Vanga H, et al. Phenetic and genetic structure of tsetse fly populations (Glossina palpalis palpalis) in southern Ivory Coast. Parasit Vectors. 2012;5:153.
Klingenberg CP. Evolution and development of shape: integrating quantitative approaches. Nat Rev Genet. 2010;11:623–35.
Lorenz C, Marques TC, Sallum MA, Suesdek L. Morphometrical diagnosis of the malaria vectors Anopheles cruzii. An. homunculus and An. bellator. Parasit Vectors. 2012;5:257.
Monteiro LR, Reis SF. Princípios de morfometria geométrica. Holos: Ribeirão Preto; 1999.
Motoki MT, Suesdek L, Bergo ES, Sallum MA. Wing geometry of Anopheles darlingi Root (Diptera: Culicidae) in five major Brazilian ecoregions. Infect Genet Evol. 2012;12(6):1246–52.
Rohlf FJ. Relative warp analysis and example of its application to mosquito wing. Madrid: Museo Nacional de Ciencias Naturales; 1993.
Vicente JL, Sousa CA, Alten B, Caglar SS, Falcuta E, Latorre JM, et al. Genetic and phenotypic variation of the malaria vector Anopheles atroparvus in southern Europe. Malar J. 2011;10:5.
Vidal PO, Carvalho E, Suesdek L. Temporal variation of wing geometry in Aedes albopictus. Mem Inst Oswaldo Cruz. 2012;107:1030–4.
Vidal PO, Suesdek L. Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti. Infect Genet Evol. 2012;12:591–6.
Vidal PO, Peruzin MC, Suesdek L. Wing diagnostic characters for Culex quinquefasciatus and Culex nigripalpus (Diptera, Culicidae). Revista Brasileira Entomol. 2011;55:134–7.
Demari-Silva B, Suesdek L, Sallum MAM, Marrelli MT. Wing geometry of Culex coronator (Diptera: Culicidae) from South and Southeast Brazil. Parasit Vectors. 2014;7:174.
Lorenz C, Marques TC, Sallum MAM, Suesdek L. Altitudinal population structure and microevolution of the malaria vector Anopheles cruzii (Diptera: Culicidae). Parasit Vectors. 2014;7:581.
Morais SA, Moratore C, Suesdek L, Marrelli MT. Genetic-morphometric variation in Culex quinquefasciatus from Brazil and La Plata. Argentina. Mem Inst Oswaldo Cruz. 2010;105(5):672–6.
Armell JH. A review of the Scapularis group of Aedes (Ochlerotatus), vol. 13. Contri Amer Ent. 1976.
Consoli RAGB, Lourenço-de-Oliveira R. Principais mosquitos de importância sanitária no Brasil. [online]. Rio de Janeiro: Editora FIOCRUZ; 1994. p. 228. ISBN 85-85676-03-5. Available from SciELO Books.
Lane J. Neotropical culicidae, vol. 2. Sao Paulo: University of Sao Paulo; 1953.
Zavortink TJ. Mosquito studies (Diptera: Culicidae) XXIX. A review of the subgenus Kerteszia of Anopheles. Contri Amer Ent. 1973;9:1–54.
Rohlf FJ. Shape statistics: procrustes superimpositions and tangent spaces. J Classif. 1999;16:197–223.
Klingenberg CP. MORPHOJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011;11:353–7.
Felsenstein J. PHYLIP (Phylogeny Inference Package). 3.6 edn. University of Washington, Seattle: Distributed by the author. Department of Genome Sciences; 2005
Dujardin JP. COO, MOG and COV for windows. Marseille: Institut de Recherches pourle Développement; 2010.
Dworkin I, Gibson G. Epidermal growth factor receptor and transforming growth factor-b signaling contributes to variation for wing shape in Drosophila melanogaster. G3 (Bethesda). 2006;173:1417–31.
Klingenberg CP, Zaklan SD. Morphological integration between developmental compartments in the Drosophila wing. Evolution. 2000;54(4):1273–85.
Gilchrist AS, Azevedo RBR, Partridge L, O’Higginsc P. Adaptation and constraint in the evolution of Drosophila melanogaster wing shape. Evol Dev. 2000;2(2):114–24.
Caprio MA, Huang JX, Faver MK, Moore A. Characterization of male and female wing-beat frequencies in the Anopheles quadrimaculatus complex in Mississippi. J Am Mosq Control Assoc. 2001;17:186–9.
Robertson SP, Caprio MA, Faver MK. Heritability of wing-beat frequency in Anopheles quadrimaculatus. J Am Mosq Control Assoc. 2002;18(4):316–20.
Cepeda-Pizarro J, Vásquez H, Veas H, Colon G. Relaciones entre tamaño corporal y biomasa en adultos de Tenebrionidae (Coleoptera) de la estepa costera del margen meridional del desierto chileno. Rev Chil Hist Nat. 1996;69:67–76.
Abbott JK, Bedhomme S, Chippindale AK. Sexual conflict in wing size and shape in Drosophila melanogaster. J Evol Biol. 2010;23:1989–97.
Reeve JP, Fairbairn DJ. Change in sexual size dimorphism as a correlated response to selection on fecundity. Heredity. 1999;83:697–706.
Moller AP, Zamora-Muñoz C. Antennal asymmetry and sexual selection in a cerambycid beetle. Anim Behav. 1997;54:1509–15.
Forrest TG. Insect size tactics and developmental strategies. Oecologia. 1987;73:178–84.
Allen CE, Zwaan BJ, Brakefield PM. Evolution of sexual dimorphism in the Lepidoptera. Annu Rev Entomol. 2011;56:445–64.
McAlpine DK. Observations on sexual behavior in some Australian Platystomatidae (Diptera: Schizophora). Rec Aust Mus. 1973;29(1):1–10.
Sivinski J. Ornaments in the Diptera. Fla Entomol. 1997;80(2):142–64.
Chenoweth SF, Rundle HD, Blows MW. Genetic constraints and the evolution of display trait sexual dimorphism by natural and sexual selection. Am Nat. 2008;171(1):22–34.
We would like to thank Aristides Fernandes for taxonomic identification of the specimens and Fernanda Silva Almeida for assembling the wings. We are also grateful to Nelly A. Gidaszewski for reviewing our manuscript and for the invaluable suggestions she made. Financial support was provided by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP grants no. 2010/15039-1 and no. 2006/02622-5) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES grant no. 23038.005274/2011-24).
The authors declare that they have no competing interests.
FV, POV and LS conceived the study, carried out the data analysis, interpreted the results and wrote the manuscript. All the authors approved the final version of the manuscript.