Wing sexual dimorphism of pathogen-vector culicids
© Virginio et al.; licensee BioMed Central. 2015
Received: 5 August 2014
Accepted: 25 February 2015
Published: 14 March 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.
KeywordsCulex Aedes Anopheles Ochlerotatus Culicidae Mosquitoes Geometric morphometrics
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
Data for Culicidae species collected
No. of specimens
Type of breeding
Rio Pinheiros, SP
-23.595138 S; -46.694258 W
17 F / 17 M
Parque Ecológico do Tietê, SP
-23.480094 S; -46.509274 W
16 F / 20 M
-22.954916 S; -45.543534 W
24 F / 25 M
São José do Rio Preto, SP
-20.810039 S; -49.368546 W
25 F / 25 M
-22.905906 S; -47.069657 W
25 F / 25 M
An. albitarsis l.s.
-20.030129 S; -49.021425 W
15 F / 09 M
-24.695063 S; -47.870972 W
32 F / 24 M
An. triannulatus l.s.
-20.030129 S; -49.021425 W
22 F / 17 M
Adult and Immature
An. strodei l.s.
-20.030129 S; -49.021425 W
22 F / 16 M
-25.012376 S, -47.935381 W
40 F / 33 M
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.
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 .
Phenetic differentiation between sexes in each of the species
Number of individuals
Cross validation score
Without allometry (%)
With allometry (%)
An. albitarsis l.s.
An. strodei l.s.
An. triannulatus l.s.
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.
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.
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).
- 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-110Google Scholar
- 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.0001Google Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Darwin CR. The descent of man, and selection in relation to sex. 2ªth ed. London: John Murray, Albermarle Street; 1874.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Wilkinson GS, Dodson GN. Function and evolution of antlers and eye stalks in flies. Cambridge, UK: Cambridge University Press; 1997.View ArticleGoogle Scholar
- Forattini OP. Culicidologia médica, vol. 2. 2002Google Scholar
- Devicari M, Lopes AR, Suesdek L. Wing sexual dimorphism in Aedes scapularis (Diptera: Culicidae). Biota Neotrop. 2011;11:165–9.View ArticleGoogle Scholar
- 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.Google Scholar
- Alexander RD, Brown Jr W. Mating behavior and the origin of insect wings. Ann Arbor, Michigan: University of Michigan Museum of Zoology; 1963.Google Scholar
- 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.View ArticleGoogle Scholar
- Klingenberg CP, Debat V, Roff DA. Quantitative genetics of shape in cricket wings: developmental integration in a functional structure. Evolution. 2010;64:2935–51.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Dujardin JP. Morphometrics applied to medical entomology. Infect Genet Evol. 2008;8:875–90.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Klingenberg CP. Evolution and development of shape: integrating quantitative approaches. Nat Rev Genet. 2010;11:623–35.PubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Monteiro LR, Reis SF. Princípios de morfometria geométrica. Holos: Ribeirão Preto; 1999.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- Rohlf FJ. Relative warp analysis and example of its application to mosquito wing. Madrid: Museo Nacional de Ciencias Naturales; 1993.Google Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Vidal PO, Carvalho E, Suesdek L. Temporal variation of wing geometry in Aedes albopictus. Mem Inst Oswaldo Cruz. 2012;107:1030–4.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Vidal PO, Peruzin MC, Suesdek L. Wing diagnostic characters for Culex quinquefasciatus and Culex nigripalpus (Diptera, Culicidae). Revista Brasileira Entomol. 2011;55:134–7.View ArticleGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Armell JH. A review of the Scapularis group of Aedes (Ochlerotatus), vol. 13. Contri Amer Ent. 1976.Google Scholar
- 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.Google Scholar
- Lane J. Neotropical culicidae, vol. 2. Sao Paulo: University of Sao Paulo; 1953.Google Scholar
- Zavortink TJ. Mosquito studies (Diptera: Culicidae) XXIX. A review of the subgenus Kerteszia of Anopheles. Contri Amer Ent. 1973;9:1–54.Google Scholar
- Rohlf FJ. Shape statistics: procrustes superimpositions and tangent spaces. J Classif. 1999;16:197–223.View ArticleGoogle Scholar
- Klingenberg CP. MORPHOJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011;11:353–7.View ArticlePubMedGoogle Scholar
- Felsenstein J. PHYLIP (Phylogeny Inference Package). 3.6 edn. University of Washington, Seattle: Distributed by the author. Department of Genome Sciences; 2005Google Scholar
- Dujardin JP. COO, MOG and COV for windows. Marseille: Institut de Recherches pourle Développement; 2010.Google Scholar
- 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.Google Scholar
- Klingenberg CP, Zaklan SD. Morphological integration between developmental compartments in the Drosophila wing. Evolution. 2000;54(4):1273–85.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Robertson SP, Caprio MA, Faver MK. Heritability of wing-beat frequency in Anopheles quadrimaculatus. J Am Mosq Control Assoc. 2002;18(4):316–20.PubMedGoogle Scholar
- 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.Google Scholar
- Abbott JK, Bedhomme S, Chippindale AK. Sexual conflict in wing size and shape in Drosophila melanogaster. J Evol Biol. 2010;23:1989–97.View ArticlePubMedGoogle Scholar
- Reeve JP, Fairbairn DJ. Change in sexual size dimorphism as a correlated response to selection on fecundity. Heredity. 1999;83:697–706.View ArticlePubMedGoogle Scholar
- Moller AP, Zamora-Muñoz C. Antennal asymmetry and sexual selection in a cerambycid beetle. Anim Behav. 1997;54:1509–15.View ArticlePubMedGoogle Scholar
- Forrest TG. Insect size tactics and developmental strategies. Oecologia. 1987;73:178–84.View ArticleGoogle Scholar
- Allen CE, Zwaan BJ, Brakefield PM. Evolution of sexual dimorphism in the Lepidoptera. Annu Rev Entomol. 2011;56:445–64.View ArticlePubMedGoogle Scholar
- McAlpine DK. Observations on sexual behavior in some Australian Platystomatidae (Diptera: Schizophora). Rec Aust Mus. 1973;29(1):1–10.View ArticleGoogle Scholar
- Sivinski J. Ornaments in the Diptera. Fla Entomol. 1997;80(2):142–64.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.