- Open Access
Wing geometry of Culex coronator (Diptera: Culicidae) from South and Southeast Brazil
Parasites & Vectorsvolume 7, Article number: 174 (2014)
The Coronator Group encompasses Culex coronator Dyar & Knab, Culex camposi Dyar, Culex covagarciai Forattini, Culex ousqua Dyar, Culex usquatissimus Dyar, Culex usquatus Dyar and Culex yojoae Strickman. Culex coronator has the largest geographic distribution, occurring in North, Central and South America. Moreover, it is a potential vector-borne mosquito species because females have been found naturally infected with several arboviruses, i.e., Saint Louis Encephalitis Virus, Venezuelan Equine Encephalitis Virus and West Nile Virus. Considering the epidemiological importance of Cx. coronator, we investigated the wing shape diversity of Cx. coronator from South and Southeast Brazil, a method to preliminarily estimate population diversity.
Field-collected immature stages of seven populations from a large geographical area in Brazil were maintained in the laboratory to obtain both females and males linked with pupal and/or larval exuviae. For each individual female, 18 landmarks of left wings were marked and digitalized. After Procrustes superimposition, discriminant analysis of shape was employed to quantify wing shape variation among populations. The isometric estimator centroid size was calculated to assess the overall wing size and allometry.
Wing shape was polymorphic among populations of Cx. coronator. However, dissimilarities among populations were higher than those observed within each population, suggesting populational differentiation in Cx. coronator. Morphological distances between populations were not correlated to geographical distances, indicating that other factors may act on wing shape and thus, determining microevolutionary patterns in Cx. coronator. Despite the population differentiation, intrapopulational wing shape variability was equivalent among all seven populations.
The wing variability found in Cx. coronator populations brings to light a new biological problem to be investigated: the population genetics of Cx. coronator. Because of differences in the male genitalia, we also transferred Cx. yojoae to the Apicinus Subgroup.
Culex (Culex) coronator Dyar and Knab, 1906 belongs to the Coronator Group, which comprises six other species (Culex usquatus Dyar, 1918, Culex usquatissimus Dyar, 1922, Culex ousqua Dyar, 1918, Culex camposi Dyar, 1925, Culex covagarciai Forattini, 1965, and Culex yojoae Strickman, 1989) . Accurate identification of these species depends on examining traits of the dissected male genitalia because both adult females and immature stages are morphologically indistinguishable. Culex coronator was described based on morphological characteristics of larva collected in St. Joseph District in Trinidad. Adult stages of this species were previously identified as Culex secutor Theobald. In subsequent years, several species were described and their similarities to Cx. coronator were observed [2–8]. Later, only Cx. coronator, Cx. usquatus, Cx. usquatissimus, Cx. ousqua, Cx. camposi and Culex covagarciai were considered valid species because of differences in male genitalia traits and geographical distribution. Morphological similarities of immature stages and adult females suggested that Cx. coronator comprises a species complex [9, 10]. Recently, Harbach  proposed a preliminary classification for the genus Culex based on morphological characteristics and transferred Cx.yojoae to the Coronator Group.
Culex coronator has the largest distribution among species of the Coronator Group. It occurs in Central and South America (Argentina, Belize, Bolivia, Brazil, Colombia, Costa Rica, El Salvador, French Guiana, Guatemala, Honduras, Nicaragua, Panama, Paraguay, Peru, Suriname, Trinidad and Tobago, and Venezuela), Mexico and United States [11, 12]. In the latter country, this species is considered an invasive species that had dispersed eastward [13–16]. Culex coronator occurs in a large range of larval habitats that vary from natural to artificial containers, temporary ground pools, in full and partial shade or full sun, in sylvatic, rural and urban environments . In Brazil, it was found in anthropic environments . Females are predominantly nocturnal, exophilic and low anthropophilic [19, 20].
Females of this species have been found naturally infected with several arboviruses, including the Saint Louis Encephalitis Virus (SLEV) in Brazil and in Trinidad and Tobago, Venezuelan Equine Encephalitis Virus (VEEV) in Mexico  and West Nile Virus in the United States [22, 23]. In southeast Brazil, few human cases of encephalitis caused by SLEV have been reported [24, 25]. However, it is possible that the prevalence of SLEV infections is underestimated because most infections are either unnoticeable or misdiagnosed because of cross-reactivity among different flaviviruses, including Dengue Virus, which are endemic in Brazil . Moreover, VEEV Subtype IF caused encephalitis in humans in Vale do Ribeira, southeastern Atlantic Forest, Brazil. Potential determinants of the epidemics in the region of Vale do Ribeira are fully unknown, including mosquito vector species .
Dujardin , believed that phenotypic variation can improve the understanding of epidemiological importance in medically important insects, once they were submitted to various control pressures, such as landscape changes. In this sense, geometric morphometrics can be used to quantify this variation. Currently, this method has been largely employed both in distinction of closely related species [28–31] as in microevolution studies [32–36]. Considering the lack of biological and ecological studies regarding Cx. coronator, and the epidemiological importance of the species, the main objective of this study was to address variability of wing morphology across a transect to better understand the population structure and microevolutionary patterns of Cx. coronator populations in South and Southeastern Brazil. The circulation of SLVE and VEEV in Brazil and the fact of Rio Grande do Sul State is a migration route of Northern Hemisphere birds, such research may be useful to advance the knowledge about this potential vector.
Mosquito immature stages were collected in 7 localities between 2010 and 2012 from October (spring) to January (summer) in localities in South and Southeastern Brazil (Table 1, Figure 1). Larvae and pupae were maintained in the laboratory at constant temperature and humidity to obtain males and females associated with their respective exuviae. Field collected larvae were transferred to double distilled water and artificially fed once a day with a mixture of finely ground Tetramin for tropical fish, Tetramin food for Beta fish, and dried pollen. A maximum of 30 larvae were kept in each receptacle with 400 ml of distilled water and fed with the same amount of food once a day. Water from each container was replaced daily to avoid development of pathogenic microorganisms. Male genitalia traits were employed for species identification and to support female identification.
Characterization of sampling localities
Pariquera-Açu municipality (PARI population) is situated in Ribeira Valley, Southeastern São Paulo State, within the Atlantic Forest domain. Because of intense anthropic pressures in the natural environment, most of the Atlantic Forest has either been removed or is under advanced secondary stage of succession. However, it remains better preserved in small forest fragments on the east coast of Brazil, with the largest continuous forest remaining in Southeastern São Paulo and northeast Paraná States. Larvae and pupae of Cx. coronator were collected in cattle ranches, ground pools, shallow ground depressions, ground holes, in natural and artificial habitats, in the proximities of forest.
Parque Ecológico do Tietê (PET population) is embedded in an urban environment, at the border of São Paulo and Guarulhos municipalities in the Tiete River floodplain. PET is a reforested area with 14 million squares meters of native floodplain vegetation. Pupae and larvae were collected in ground pools, in a restricted area in the park.
Rio de Janeiro municipality (RJ population) has a heterogeneous landscape in the Atlantic rainforest domains, with lowland coastal areas, mangrove and coastal mountains. Embedded in the southern border of the city, in the coastal mountains, there is a secondary forest preserved area named Tijuca Forest that is part of the Serra do Mar biodiversity corridor. Culex coronator specimens were captured in natural and artificial ground pools in two separated localities. One locality was a rural area situated in the proximity of the foot of the massif of Pedra Branca, which is also part of Serra do Mar biodiversity corridor. The second collection was carried out in an urban area, in Serra da Carioca, Tijuca Forest.
Conceição do Mato Dentro municipality (MG population) is located in the eastern portion of the Serra do Cipó and Espinhaço Mountains. It is inside a large ecological corridor formed by a mountain range that runs from the central region of Minas Gerais State, to the north of Bahia State. This is a transition area between Interior Atlantic Forest and Cerrado, with pastures as the major environmental modifier. Immatures of Cx. coronator were found in animal tracks and in rock holes near waterfalls.
Tijucas do Sul and Joinville municipalities (PR and SC populations) are embedded in a transition region between dense and semidecidous rain forests, with the family agriculture as the major environmental modifier. Sampling in Paraná State was in the first plateau  along the sides of secondary, dirty roads, in altitudes above 900 m. In contrast, in Santa Catarina States, immature collections were conducted in lower altitude areas (207 meter altitude) along the Atlantic coast.
Maquiné municipality (RS population) is a northeast rural area in Rio Grande do Sul State, in the south limit of the Atlantic Forest domain. Maquiné includes forest preserved areas, for instance, the Biological Reserve of Serra Geral and Rota do Sol, and it is situated in a transition zone between the coastal lowlands and the slopes of Serra Geral mountain. The vegetation in the area is characterized by a dense rain forest in the east and semidecidous forest in the west side Serra Geral. Culex coronator larvae and pupae were collected in rural areas, in marshes and in artificial man made ground holes.
Wing preparation and data collection
Left wings of the adult females were removed and mounted on flat microscope slides under a 0.08 - 0.12 mm glass coverslip using Canada balsam. Microscope slides of wings, male genitalia, pupal and fourth-instar larval exuviae were deposited in the Coleção Entomológica de Referência, Faculdade de Saúde Pública, Universidade de São Paulo (FSP-USP), Brazil. Wing images were obtained with a Leica M205C, under 50x magnification. In each image, 18 landmarks (Figure 2) were digitalized using tpsDig software v2.0 . Culex coronator populations included in the study were obtained in Pariquera-Açu (PARI, n = 23) and Parque Ecológico do Tietê (PET, n = 27) - both in São Paulo State- and in the States of Rio de Janeiro (RJ, n = 23), Minas Gerais (MG, n = 27), Santa Catarina (SC, n = 29), Paraná (PR, n = 24) and Rio Grande do Sul (RS, n = 23). Additionally, 23 samples of Culex quinquefasciatus, from Santa Vitória do Palmar municipality , Rio Grande do Sul State, were employed as outgroup in the Neighbour-Joining analyses .
To investigate the distribution of wing shape and to address specific questions about intraspecific heterogeneities in Cx. coronator, a generalized Procrustes distance analysis was performed in COV . A Procrustes distance matrix was constructed and employed to generate a neighbor-joining (NJ) topology in PHYLIP version 3.7 . NJ topology was visualized in TreeView . Correlation between Procrustes and geographic distance was addressed using the Mantel’s Test.
The overall size of left wing was estimated by the centroid size (CS) in MOG . CS was estimated according to Bookstein . Descriptive statistic graphs were constructed in Statistica 7 (StatSoft). The ANOVA and Tukey-Kramer tests were employed to ascertain the significance of differences between averages of the CS among populations. Both tests were conducted in GraphPad InStat v.3 (GraphPad Software, San Diego California USA). To determine potential association between mean centroid size and latitudinal variation, a simple regression analysis was conducted in Statistica 7 (StatSoft). The allometric effect was verified and removed from further analyses employing the regression of Procrustes coordinates (dependent variable) against CS (independent variables).
To address variation in wing shape among groups, a discriminant analysis of principal components derived from Procrustes superimposition method was performed using Statistica 7 software. Only the set of 12 first Relative Warps were employed in discriminant analysis, since they explained 94.25% of the total variation of wing shape. Accuracy of population recognition yielded by morphometrics was evaluated by a reclassification test. For this test, the similarity of wing shape of each individual was compared with the average of wing shape of each population, using Procrustes distances in MorphoJ software . In order to find the two most divergent wing shapes, we performed a Canonical Variate Analysis (CVA) also in MorphoJ program. The most distant samples in morphospace of the CVAs were taken and, then, an analysis of spline plots in tpsSplin v3.2 program was carried out, to assess wing deformations between those specimens. Morphological diversity was estimated by the “amount of dispersion” of the individuals in the morphospace of principal components (PCs). The dispersion was estimated using the following procedure. The Positional (Cartesian) coordinates of plots in the morphospace of PCs (each one corresponding to a single mosquito individual) were digitalized using TpsDig software, similar to the wing landmarks. The centroid size vector of a set of individuals in the morphospace (a population) was estimated in TpsRelW software . Following this, the value of this vector was used to evaluate morphological diversity within each population (MD). The amount of dispersion of individuals (of a single set) in the morphospace of PCs is proportional to the morphological variability of that population set.
CS values ranged from 1.36 mm to 2.25 mm in Cx. coronator populations employed in this study. The lowest mean for the CS was found in RJ (1.45 mm) populations, while the highest mean CS was observed in PR (2.11 mm). The remaining populations had the following CVs means: 1.62 mm (MG), 1.99 mm (PARI), 1.51 mm (PET), 1.79 mm (SC) and 2.01 mm (RS). The population with the highest variation was PARI (1.79 mm to 2.2 mm) (Figure 3). Statistical significance of the differences between mean centroid sizes among populations is in Table 2. Results of the simple regression analysis showed a positive association between latitudinal gradient and mean centroid size (Figure 4), with moderate correlation (R2 = 0.50) that was not significant (P = 0.072).
The NJ tree generated using the Procrustes distance matrix (Figure 5) showed two population clusters. One cluster included populations from MG and RJ, whereas the second cluster was comprised of individuals from the remaining 5 populations.
Likewise, southern populations (PR, SC and RS) showed the highest similarity, with PR and RS morphologically more similar to each other than to SC. PARI population was found closer to southern than to southeastern populations, whereas PET was the most distinct grouping within the cluster. The most divergent samples from CVA were from RJ and SC. Thin Plate Splines analysis showed most deformations in the medial part of the wing (landmarks 13, 14, 15 and 16) and landmark 2 - Figure 6. Results of regression analysis carried out to investigate evolutionary integration in the wing shape and wing size were statistically significant, with moderate allometric effect (5.11%). After removing the allometric effect, results of the permutation test showed that differences between MG and the remaining populations were all statistically significant (Table 3). Results of cross-validated reclassification tests are in Additional file 1. Intrapopulational diversity indexes for all 7 populations were similar (see Additional file 2).
Wing morphology analysis of Cx. coronator populations collected from South to Southeastern Brazil were polymorphic showing population structured with microevolutionary patterns. Although Procrustes distances were not correlated with geographical distances, our results clearly separate Cx. coronator into two subgroups that are correlated to the macro-geographical regions where samples were collected. As such, one subgroup includes individuals from the “southeastern populations” (RJ and MG), whereas the second subgroup is composed of the “southern populations” (PARI, PET, PR, SC and RS). Morphological differentiation between both southeastern and southern Cx. coronator populations may be explained by the presence of high altitude areas of Serra da Mantiquera, along the borders of São Paulo, Rio de Janeiro and Minas Gerais States. The mountain range may preclude dispersal of Cx. coronator through those areas, causing isolation of RJ and MG from the southern populations. Thus, it is likely that a shallow, incipient speciation process is occurring in those two population subgroups.
Considering the subgroup composed of RJ and MG populations, we found that wing shape variability within MG is slightly higher than in the RJ population. It is noteworthy that representatives from MG populations were collected in the ecotone between two biodiversity hotspots, Cerrado and Atlantic Forest, within the ecological corridor of the Cerrado. In the ecotone, important ecological and evolutionary processes may occur and induce genetic diversity and wing shape variability, in response to highly heterogeneous physiognomies and landscapes present in the zone [45, 46]. In contrast, the RJ population was from localities embedded in the coastal Atlantic Forest where ecological determinants are less heterogeneous than in the transition zone between Cerrado and Atlantic Forest. In a study that used geometric morphometric methods to investigate evolutionary association between biome and female wing shape in Anopheles darlingi, a primary mosquito vector of human Plasmodium, Motoki et al. suggested that wing shape variations could be caused by distinct selective pressures present in Cerrado and interior and coastal Atlantic Forest. Based on both the results of the present study and the findings of Motoki et al. , it is possible that shallow wing shape heterogeneities between MG and RJ populations were generated by evolutionary pressures present in the Cerrado and Atlantic Forest, rather than isolation by distance. Similarly, a study with dragonfly associated heterogeneities observed in the forewing shape with selective pressures represented by distinct landscapes, with the wing becoming broader from forested to open areas . Although the effect of selection on the wing landmarks is not fully understood, wing shape characteristics have been effectively employed to address population structure associated with landscapes, habitat split and isolation in several arthropod species [32–34, 48, 49].
In the NJ topology, within the clade composed of southern populations (Figure 5), PET was separate from the remaining southern population, which formed a secondary grouping (PARI, PR, SC, and RS). Accordingly, PARI population was demonstrated to be morphologically more similar to SC than to PET and PR populations. Environmental differences between PET and remaining southern sampling localities, in contrast, are enormous. For instance, samples from PET populations were obtained in a recently reforested area of semi deciduous rain forests in the floodplain of the highly polluted Tiete River, at the border of São Paulo and Guarulhos municipalities. This area suffered and is still under strong ecological pressures caused by intensive anthropic changes in the environment and pollution. Distinctly, individuals from PARI were captured in the proximity of the borders of São Paulo and Paraná States, in a rural area embedded within the coastal Atlantic Forest (Figure 1). On the other hand, ecological exposure factors are similar for PR, SC and RS populations. Field collections in Paraná, Santa Catarina and Rio Grande do Sul States, were in areas of semi-deciduous forest in a subtropical climate. Similar ecological conditions may facilitate both the establishment and the dispersion of individuals that are better adapted to environmental characteristics present in semi-deciduous forest. Consequently, gene flow is still occurring between those populations which are undergoing similar evolutionary pressures. However, ecological factors between PARI and SC are distinct and did not explain a close morphological similarity between them. Likely, the similarity may be caused by retention of ancestral wing shape polymorphism. In spite of being able to occupy a wide range of larval habitats in heterogeneous landscapes within distinct biomes physiognomies, wing shape diversity in Cx. coronator was similar both within and in southern populations. So, it is possible that dispersal and gene flow, in both cases, are being somehow facilitated.
Our results also showed a positive association between wing size and latitude, with the average centroid size increasing in higher latitudes. According to Bergman’s rule, endothermic vertebrates in higher latitudes tend to be bigger, as a result of colder climates. Later, Ray  stated that ectothermic invertebrates such as insects also would fit this rule. However, Mousseau  suggested that the body size in ectothermic animals decreased with increasing latitude, and consequently, with decreasing temperature, and proposed the reverse Bergmann’s rule. Recently, Blanckenhorn & Demont  demonstrated that, for arthropods, both theories occurs in nature. Meanwhile, larger insects, with longer development time, follow the reverse Bergmann’s rule, while smaller insects, with faster development periods, would fit the Bergman’s rule. Here, the analyzed sample populations of Cx. coronator clearly follows the Bergman’s rule.
Currently, association between wing and body size in vectorial capacity is not completely understood. Meanwhile, body size can have a major effect on the longevity, fecundity and size of blood meal of vector-borne disease insects [53, 54], influencing the vectorial capacity of populations. Additionally, several studies have shown that heterogeneities in the wing size may be a plasticity response to environmental factors, such as relative humidity, temperature and food availability [55–57]. Nevertheless, morphological plasticity has an important role in ecology because populations that have higher adaptative responses can occupy larger ecological niches . Consequently, a better understanding of the causes and consequences of phenotypic plasticity may be important to explore mechanisms of evolution. Moreover, plasticity is a result of evolutionary integrations between environmental factors and genotypes [59–61]. Thereby, environmental determinants can facilitate the perpetuation of some populations under directional selection .
Hitherto, regional dissimilarities did not indicate Cx. coronator as a group of cryptic species. Nevertheless, we cannot discard the hypothesis that the polymorphisms observed among populations may be the substratum of an incipient speciation process. Recently, Harbach  compiled morphological data to propose a preliminary classification for the genus Culex, including the subgenus Culex. Accordingly, the Coronator Group is composed of Cx. camposi, Cx. covagarciai, Cx. ousqua, Cx. usquatissimus, Cx. usquatus and Cx. yojoae. The latter species was described by Strikman  based on female, male, fourth-instar larva and pupae from La Joya, Cortes Department, Honduras. Comparisons among characteristics of the male genitalia of Cx. yojoae and the remaining Neotropical Culex (Culex) species show that Cx. yojoae is more similar to Culex chidesteri and Culex pseudostigmatosoma than to any species of both the Coronator Group. It is noteworthy that Cx. chidesteri was included in the Apicinus Subgroup, whereas Cx. pseudostigmatosoma was considered a species of the Tarsalis Subgroup, both subgroups of the Pipiens Group, by Harbach . Presence of foliform g setae in the distal division of the subapical lobe of the gonocoxite, the shape of the ventral and dorsal divisions of the lateral plate of the phalosome and the presence of spatulate spines in the paraproct crown suggest that Cx. chidesteri, Cx. yojoae and Cx. pseudostigmatosoma may comprise a species complex. We are herein transferring both Cx. yojoae and Cx. pseudostigmatosoma to the Apicinus Subgroup. Notwithstanding, it is obvious that the subgenus Culex needs further and more detailed studies before proposing a new internal classification for it.
Results of the present study show that Cx. coronator is a polymorphic species with morphological dissimilarities in both wing shape and size. Considering that wing shape has a polygenic basis , differences among populations may be caused by polymorphisms in those genes. Moreover, it is plausible to suppose that heterogeneities in wing size are correlated to plasticity response to ecological determinants of the larval habitats. In both cases, further studies, including molecular data, will be necessary to provide a more comprehensive understanding of the variability of the wing traits of a largely dispersed mosquito vector species. While identifying the Culex specimens collected for the present study, we checked our samples with all species of the Coronator Group as defined by Harbach . In considering only adult females characteristics discussed by Strikman , he concluded that Cx. yojoae was a species of the Coronator group. However, most morphological traits used to identify species of the genus Culex are from male genitalia, thus we concluded that there is a misinterpretation in Harbach’s  revision of Culex. Consequently, we strongly suggest that Cx. yojoae belongs to Apicinus Subgroup due to the similarities between the male genitalia of this species and Cx. chidesteri.
Harbach RE: Classification within the cosmopolitan genus Culex (Diptera: Culicidae): The foundation for molecular systematics and phylogenetic research. Acta Trop. 2011, 120: 1-14. 10.1016/j.actatropica.2011.06.005.
Dyar H: The larvae of Culicidae classified as independent organisms. Insec Inscit Menst. 1906, 6: 86-111.
Dyar H: New American mosquitoes (Diptera, Culicidae). Insec Inscit Menst. 1918, 6: 86-111.
Dyar H: Note on the male genitalia of Culex coronator and allied forms. Insec Inscit Menst. 1922, 10: 18-19.
Dyar H: Some mosquitoes from Ecuador (Diptera, Culicidae). Insec Inscit Menst. 1925, 13: 27-31.
Anduze P: Estudios de Entomologia Medica en el Estado Merida –Venezuela: La Fauna Culicidiana. Boletin de Entomologia Venezolana. 1943, 2: 189-196.
Vargas L, Palácios M: Culex (Culex) coronator mooseri n. subsp. de Mexico. Rev Inst Sal Enf Trop. 1954, 14: 33-37.
Forattini OP: Entomologia Médica. 1965, São Paulo: Editora da Universidade de São Paulo Press
Bram R: Classification of Culex subgenus Culex (Diptera: Culicidae). Proc US Nat Mus. 1967, 120: 1-123.
Foratinni OP: Culicidologia Médica. 2002, São Paulo: Editora da Universidade de São Paulo Press
Walter Reed Biosystematics Unit. [http://www.wrbu.org/]
Bond JG, Casas-Martínez M, Quiroz-Martínez H, Novelo-Gutiérrez R, Marina CF, Ulloa A, Orozco-Bonilla A, Muñoz M, Williams T: Diversity of mosquitoes and the aquatic insects associated with their oviposition sites along the Pacific coast of Mexico. Parasit Vectors. 2014, 7: 41-10.1186/1756-3305-7-41.
Debboun M, KUHR DD, Rueda LM, Pecor JE: First Record of Culex (Culex) coronator in Louisiana, USA. J Am Mosq Control Assoc. 2005, 21: 455-457. 10.2987/8756-971X(2006)21[455:FROCCC]2.0.CO;2.
Smith JP, Walsh JD, Cope EH, Tennant RAJ, Kozak JA, Darsie RJ: Culex Coronator Dyar and Knab: a new Florida species record. J Am Mosq Control Assoc. 2006, 22: 330-332. 10.2987/8756-971X(2006)22[330:CCDAKA]2.0.CO;2.
Gray KM, Burkett-Cadena ND, Eubanks MD: Distribution Expansion of Culex coronator in Alabama. J Am Mosq Control Assoc. 2008, 24: 585-587. 10.2987/08-5778.1.
Moulis RA, Russell JD, Lewandowski HBJ, Thompson PS, Heusel JL: Culex coronator in Coastal Georgia and South Carolina. J Am Mosq Control Assoc. 2008, 24: 588-590. 10.2987/5766.1.
Pecor JE, Harbach RE, Peyton EL, Roberts DR, Rejmankova E, Manguin S, Palanko J: Mosquito studies in Belize, Central America: records, taxonomic notes, and a checklist of species. J Am Mosq Control Assoc. 2002, 18: 241-276.
Ribeiro AF, Urbinatti PR, Duarte AMRC, Paula MB, Pereira DM, Mucci LF, Fernandes A, Mello MHSH, Matos MOJ, Oliveira RC, Natal D, Malafronte RS: Mosquitoes in degraded and preserved areas of the Atlantic Forest and potential for vector-borne disease risk in the municipality of São Paulo, Brazil. J Vector Ecol. 2012, 37: 316-326. 10.1111/j.1948-7134.2012.00233.x.
Carpenter SJ, Lacasse WJ: Mosquitoes of North America. 1955, London: University of California Press
Consoli RAGB, Lourenço-de-Oliveira R: Principais mosquitos de importância sanitária no Brasil. 1994, Fiocruz: Rio de Janeiro: Brasil Press, 97-114.
Mackay AJ, Kramer WL, Meece JK, Brumfield RT, Foil LD: Host Feeding Patterns of Culex Mosquitoes (Diptera: Culicidae) in East Baton Rouge Parish. Louisiana J Med Entomol. 2010, 47: 238-248. 10.1603/ME09168.
Centers for Diasease Control and Prevention. [http://www.cdc.gov/ncidod/dvbid/westnile/mosquitoSpecies.htm#06, 2008]
Unlu I, Wayne LKW, Foil LD RAFe: Detection of West Nile Virus RNA in Mosquitoes and Identification of Mosquito Blood Meals Collected at Alligator Farms in Louisiana. J Med Entomol. 2010, 47: 625-633. 10.1603/ME09087.
Santos CLSSM: Genetic characterization of St. Louis encephalitis virus isolated from human in São Paulo, Brazil. Mem Inst Oswaldo Cruz. 2006, 101: 57-63.
Mondini ABR, Cardeal ILS, Santos TMILS, Eduardo L, Nunes SHP, Silva GCD, Sarkis Madrid MCFS, Rahal P, Figueiredo LT, Neto FC, Nogueira ML: Simultaneous infection by DENV-3 and SLEV in Brazil. J Clin Virol. 2007, 40: 84-86. 10.1016/j.jcv.2007.06.007.
Figueiredo LTM: Emergent arboviruses in Brazil. Rev Inst Med Trop. 2007, 40: 224-229.
Dujardin JP: Morphometrics applied to medical entomology. Infec Gen Evol. 2008, 8: 875-890. 10.1016/j.meegid.2008.07.011.
Márquez EJ-ON, Gómez-Palacio A, Dujardin JP: Morphometric and molecular differentiation of a Rhodnius robustus-like form from R. robustus Larousse, 1927 and R. prolixus Stal, 1859 (Hemiptera, Reduviidae). Acta Trop. 2011, 120: 103-109. 10.1016/j.actatropica.2011.06.009.
Vidal PO, Peruzin MC, Suesdek L: Wing diagnostic characters for Culex quinquefasciatus and Culex nigripalpus (Diptera, Culicidae). Rev Bras Entomol. 2011, 55: 134-137. 10.1590/S0085-56262011000100022.
Ruangsittichai JAC, Dujardin JP: Interspecific and sexual shape variation in the filariasis vectors Mansonia dives and Ma. bonneae. Infec Gen Evol. 2011, 11: 2089-2094. 10.1016/j.meegid.2011.10.002.
Lorenz CMT, Sallum MAM, Suesdek L: Morphometrical diagnosis of the malaria vectors Anopheles cruzii, An. homunculus and An. bellator. Parasit Vectors. 2012, 5: 257-10.1186/1756-3305-5-257.
Vicente JL, Sousa CA, Alten B, Caglar SS, Falcutá E, Latorre JM, Toty C, Barré H, Demirci B, Luca MD, Toma L, Alves R, Salgueiro P, Silva TL, Bargues MD, Mas-Coma S, Boccolini D, Romi R, Nicolescu G, Rosário VE, Ozer N, Fontenille D, Pinto J: Genetic and phenotypic variation of the malaria vector Anopheles atroparvus in southern Europe. Malaria J. 2011, 10: 5-10.1186/1475-2875-10-5.
Motoki MT, Suesdek L, Bergo ES, Sallum MAM: Wing geometry of Anopheles darlingi Root (Diptera: Culicidae) in five major Brazilian ecoregions. Infec Gen Evol. 2012, 12: 1246-1252. 10.1016/j.meegid.2012.04.002.
Demirci B, Lee Y, Lanzaro GC, Alten B: Altitudinal genetic and morphometric variation among populations of Culex theileri Theobald (Diptera: Culicidae) from northeastern Turkey. J Vector Ecol. 2012, 37: 197-209. 10.1111/j.1948-7134.2012.00217.x.
Prudhomme JGF, Rahola N, Ouanaimi F, Guernaoui S, Boumezzough A, Bañuls AL, Sereno D, Alten B: Wing size and shape variation of Phlebotomus papatasi (Diptera: Psychodidae) populations from the south and north slopes of the Atlas Mountains in Morocco. J Vector Ecol. 2012, 37: 137-147. 10.1111/j.1948-7134.2012.00210.x.
Kaba D, Ravel S, Acapovi-Yao G, Solano P, Allou K, Bosson-Vanga H, Gardes L, N’Goran EK, Schofield CJ, Kon M, Dujadin JP: Phenetic and genetic structure of tsetse fly populations (Glossina palpalis palpalis) in southern Ivory Coast. Parasit Vectors. 2012, 5: 153-10.1186/1756-3305-5-153.
Rohlf FJ: Shape statistics: procrustes superimpositions and tangent spaces. J Classif. 1999, 16: 197-223. 10.1007/s003579900054.
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: 672-676. 10.1590/S0074-02762010000500012.
Dujardin JP: Morphometrics in Medical Entomology. [http://mome-clic.com/]
Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.6. 2005, Seattle: Department of Genome Sciences, University of Washington, Distributed by the author
Page RDM: TREEVIEW: an application to display phylogenetic trees on personal computers. Comp Applic Biosci. 1996, 12: 357-358.
Bookstein FL: Morphometric tools for landmark data: Geometry and Biology. 1991, New York: Cambridge University Press
Klingenberg CP: MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Res. 2011, 11: 353-357. 10.1111/j.1755-0998.2010.02924.x.
Schilthuizen M: Ecotone: speciation-prone. Trends Ecol Evol. 2000, 15: 130-131. 10.1016/S0169-5347(00)01839-5.
Silva JMC, Bates JM: Biogeographic patterns and conservation in the South American Cerrado: a tropical savanna hotspot. BioScience. 2002, 52: 225-234. 10.1641/0006-3568(2002)052[0225:BPACIT]2.0.CO;2.
Outomuro D, Dijkstra K-D, Johansson F: Habitat variation and wing coloration affect wing shape evolution in dragonflies. J Evol Biol. 2013, 26: 1866-1874. 10.1111/jeb.12203.
Schachter-Broide JGR, Kitron U, Dujardin JP: Temporal variations of wing size and shape of Triatoma infestans (Hemiptera: Reduviidae) populations from northwestern Argentina using geometric morphometry. J Med Entomol. 2009, 46: 994-1000. 10.1603/033.046.0504.
Bouyer JRS, Dujardin JP, Meeus T, Vial L, Thévenon S, Guerrini L, Sidibé I, Solano P: Population structuring of Glossina palpalis gambiensis (Diptera: Glossinidae) according to landscape fragmentation in the Mouhoun river, Burkina Faso. J Med Entomol. 2007, 44: 788-795. 10.1603/0022-2585(2007)44[788:PSOGPG]2.0.CO;2.
Ray C: The application of Bergmann’s and Allen’s rules to poikilotherms. J Morph. 1960, 106: 85-108. 10.1002/jmor.1051060104.
Mousseau TA: Ectotherms follow the converse Bergmann’s rule. Evolution. 1997, 51: 630-632. 10.2307/2411138.
Blanckenhorn WU, Demont M: Bergmann and converse Bergmann latitudinal clines in arthropods: two ends of a continuum?. Integr Comp Biol. 2004, 44: 413-424. 10.1093/icb/44.6.413.
Silver JB, Renshaw M: Are mosquitoes from lower latitudes smaller?. Med Vet Entomol. 1999, 13: 219-223. 10.1046/j.1365-2915.1999.00157.x.
Ameneshewa B, Servide MW: The relationship between female body size and survival rate of the malaria vector Anopheles arabienses in Ethiopa. Med Vet Entomol. 1996, 10: 170-172. 10.1111/j.1365-2915.1996.tb00724.x.
Jirakanjanakit NLS, Thongrungkiat S, Apiwathnasorn C, Singhaniyom S, Bellec C, Dujardin JP: Influence of larval density or food variation on the geometry of the wing of Aedes (Stegomyia) aegypti. Trop Med Int Health. 2007, 12: 1354-1360. 10.1111/j.1365-3156.2007.01919.x.
Morales-Vargas ER, Ya-Umphan P, Phumala-Morales N, Komalamisra N, Dujardin JP: Climate associated size and shape changes in Aedes aegypti (Diptera: Culicidae) populations from Thailand. Infect Genet Evol. 2010, 10: 580-585. 10.1016/j.meegid.2010.01.004.
Ayala D, Caro-Riaño H, Dujardin JP, Rahola N, Simard F, Fontenille D: Chromosomal and environmental determinants of morphometric variation in natural populations of the malaria vector Anopheles funestus in Cameroon. Infect Genet Evol. 2011, 11: 940-947. 10.1016/j.meegid.2011.03.003.
Dujardin JP: Modern morphometrics of medically important insects. Genetics and Evolution of Infectious diseases. Edited by: Tibayrenc M. 2011, Burlington: Elsevier, Chapter 16: 473-501. pp749. ISBN: 978-0-12-384890-1
Via S, Lande R: Genotype-environment interaction and the evolution of phenotypic plasticity. Evolution. 1985, 39: 505-522. 10.2307/2408649.
Carreira VP, Imberti MA, Mensch J, Fanara JJ: Gene-by-temperature interactions and candidate plasticity genes for morphological traits in Drosophila melanogaster. PLoS ONE. 2013, 8: e70851-10.1371/journal.pone.0070851.
Price TD, Qvarnström A, Irwin DE: The role of phenotypic plasticity in driving genetic evolution. Proc R Soc Lond B. 2003, 270: 1433-1440. 10.1098/rspb.2003.2372.
Strickman D: Culex pseudostigmatosoma, Cx. yojoae, and Cx. aquarius: new Central American species in the subgenus Culex (DIPTERA: CULICIDAE). Mosq Syst. 1989, 21: 143-171.
Carreira VP, Soto IM, Mensch J, Fanara JJ: Genetic basis of wing morphogenesis in Drosophila: sexual dimorphism and non-allometric effects of shape variation. BMC Dev Biol. 2011, 11: 32-10.1186/1471-213X-11-32.
We are grateful to the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP No. 2010/02658-5 to BDS; FAPESP No. 2011/20397-7 to MAMS; FAPESP No. 2006/02622-6 to LS), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, No. 301666/2011-3 to MAMS), CAPES Research Grants# 23038.005274/2011-24 and 1275/2011. We are indebted to Denise Sant’Ana and Sandra S. Nagaki who kindly helped us with field collections.
We declare that we don’t have competing interests.
MAMS, BDS, MTM conceived the study. BDS collected material in field and morphometric data. BDS and LS performed wing morphometric analyses and interpretation of the results. All authors read and agreed with the final version of the manuscript.