Inferring evolutionary scenarios with geostatistics and geographical information systems for the viperid snakes Vipera latastei and Vipera monticola

The present study aimed to infer evolutionary scenarios for Vipera latastei and Vipera monticola in the Iberian Peninsula and the Maghreb through the identification of spatial patterns in morphological character variation and biogeographic patterns in morphological variability distribution. Ten morphological traits from 630 vipers were analysed with geostatistic and ecological niche modelling in a geographical information system. Interpolation by Kriging was used to generate surfaces of morphological variation, which were combined with spatial principal components analysis (SPCA). Putative morphological differentiated groups generated by SPCA maps were tested with discriminant function analysis (DFA). Maximum entropy modelling and nine environmental variables were used to identify factors limiting the distribution of groups and areas for their potential occurrence. Groups supported by DFA were: Western Iberia, Eastern Iberia, Rif plus Middle Atlas, Algeria, and High Atlas. Their distribution is influenced by common environmental factors such as precipitation. Areas of probable sympatry between Iberian groups matched the morphological clines observed by geostatistics tools. Geographic variation patterns in V. latastei‐monticola are probably due to vicariant separation of Iberian and African populations during the opening of the Strait of Gibraltar, and population refugia during the Quaternary glaciations with secondary contact. The taxonomic status of northern Morocco and Algerian groups should be further investigated. We conclude that geostatistics and niche‐modelling tools are adequate to infer morphological variability across wide geographic ranges of species. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 790–806.     

The use of geostatistics and GIS for evolutionary history studies: the case of the nose‐horned viper (Vipera ammodytes) in the Balkan Peninsula

Geostatistics and geographical information system (GIS) procedures are novel techniques helpful for the identification of environmental correlates sustaining contact zones among subspecies or closely related species. In this paper, we tried to infer evolutionary scenarios for Vipera ammodytes across the European part of its distribution area using geostatistics and ecological niche‐based models, hence trying to solve several biogeographical questions that remained unclear after the application of classical morphological tools and genetic analyses. Eleven morphological traits from 871 vipers were analysed with geostatistics and ecological niche‐based modelling. Interpolation by kriging was used to generate surfaces of morphological variation, which were combined with spatial principal components analysis (SPCA). SPCA maps were used to test putative morphological differentiated groups with discriminant function analysis (DFA). Maximum entropy modelling and seven environmental variables were used to identify factors limiting the distribution of groups and areas for the potential occurrence of such groups. Three patterns of morphological variation were observed: a north‐west/south‐west cline, transition zones with steep clines of variation in a west–east arc, and particular character traits that disturbed the general cline. SPCA identified between three and nine putative population groups, of which three were supported by DFA. Areas of potential occurrence of these groups were coherent with the range of the three subspecies of V. ammodytes currently recognized. The distribution of all subspecies was mostly related to precipitation in the driest month. Areas of probable sympatry between subspecies are generally small and restricted. The main patterns of geographic variation of morphological characters for V. ammodytes were similar to the patterns obtained for Vipera latastei and Vipera monticola; the same environmental factors limit the distribution of differentiated groups of vipers in the Balkans and the Iberian Peninsula. The influence of humidity on the variation of morphological traits in spatially separated viper taxa from the two European peninsulas coincides with their phylogenetic relatedness. Geostatistics and GIS procedures were successful in the identification of environmental correlates sustaining contact zones among V. ammodytes subspecies in the Balkans. The same techniques should be applied for studying other parapatric forms and refugia regions. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 651–666.     

GIS-based niche models identify environmental correlates sustaining a contact zone between three species of European vipers

The current range of European vipers is mostly parapatric but local‐scale allopatric distribution is common and few cases of sympatry are known. In the High Course of Ebro River, northern Spain, there is a contact zone between Vipera aspis, V. latastei, and V. seoanei. Sympatry was detected between aspis and latastei and also specimens with intermediate morphological traits. Presence‐data at a local scale (1 × 1 km) and ecological niche‐based models manipulated in a GIS were used to (1) identify how environmental factors correlate with the distribution of the three vipers and with the location of the sympatry area, and (2) identify potential areas for viper occurrence and sympatry. Ensemble for casting with 10 Maximum Entropy models identified a mixture of topographical (altitude, slope), climatic (precipitation, evapotranspiration, and minimum and maximum temperature), and habitat factors (land cover) as predictors for viper occurrence. Similar predicted probabilities according to the variation of some environmental factors (indicating probable sympatry) were observed only for aspis‐latastei and aspis‐seoanei. In fact, areas of probable occurrence of vipers were generally allopatric but probable sympatry between vipers was identified for aspis‐latastei in 76 UTM 1 × 1 km squares, for aspis‐seoanei in 23 squares, and latastei‐seoanei in two squares. Environmental factors correlate with the location of this contact zone by shaping the species range: some enhance spatial exclusion and constrain distribution to spatially non‐overlapping ranges, while others allow contact between species. The distribution in the contact zone apparently results from the balance between the pressures exerted by the different environmental factors and in the sympatry area probably by interspecific competition. Further ecological and genetical data are needed to evaluate the dynamics of the probable hybrid zone. GIS and niche‐modelling tools proved to be powerful tools to identify environmental factors sustaining the location of contact zones.     

Integrating classical and spatial multivariate analyses for assessing morphological variability in the endemic Iberian viper Vipera seoanei

Historical and ecological processes have deeply affected biogeographic patterns of animals. Studying morphological variability of species, using classical and spatial analyses, can elucidate these patterns and give insights on both processes. Morphological variability of the endemic Iberian viper Vipera seoanei is examined to identify morphological coherent groups, biogeographic patterns and the putative role of abiotic pressures in the geographic variation of morphological variation. Results from classic and spatial multivariate analyses over 27 morphometric traits for 468 specimens from the global range of the species were integrated. Classic analyses reported large morphological variability and confirmed the differentiation of two coherent groups, which are representatives of current subspecies. Spatial analyses reported a geographic gradient pattern from western Cantabrian Mountains to the rest of the study area. Areas of high morphological variability were found, and two spatial coherent groups with an integration zone were recognized. Significant spatial correlations and trends suggest that some traits could be under selection and may display adaptations to local environments. Although observed patterns can be attributed to Pleistocene climatic cycles, an adaptive diversification of the species is supported. The combination of classical and spatial multivariate analyses is a useful methodology to identify morphological patterns and infer underlying factors.     

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.

     

A mitochondrial DNA phylogeny of the endangered vipers of the Vipera ursinii complex

The last two populations of the Hungarian meadow viper Vipera ursinii rakosiensis were thought to persist in the steppe fragments of Hungary until meadow vipers were discovered in central Romania (Transylvania), suggesting a possible existence of remnant populations elsewhere. We assessed the phylogenetic position of the Transylvanian vipers using 2030 bp of mitochondrial DNA sequence. We showed that they were closely related to the Hungarian vipers, while those from northeastern Romania (Moldavia) and Danube Delta belonged to the subspecies Vipera ursinii moldavica. Montane subspecies from Europe (Vipera ursinii ursinii and Vipera ursinii macrops) formed a sister clade to the two lowland subspecies. Vipera renardi formed a sister clade to V. ursinii, with populations from the Greater Caucasus (Vipera renardi lotievi) and Tien Shan (Vipera renardi tienshanica) as the sister group to Vipera renardi renardi, and Vipera renardi eriwanensis from the Lesser Caucasus as the most basal taxon in the species. Our results illustrate that the divergence between the lowland and montane populations occurred separately in each species and several times in V. renardi. We demonstrated that the recently discovered Transylvanian population is the third surviving population of V. u. rakosiensis and the only known population outside of Hungary.     

Differentiation of <Emphasis Type="Italic">Harmonia axyridis</Emphasis> Pall. according to polymorphic morphological traits and variability of the mitochondrial COI gene

The geographic variability of Harmonia axyridis Pall. has been studied with regard to three morphological features. A study of the cytochrome C oxidase subunit I (COI) polymorphism was begun. Two geographically separated zones inside the Russian part of the geographic range of the species are recognized. Interpopulational differences, as judged using polymorphic traits, are very small within each zone. In contrast, the differences between the populations of the western and eastern zones reach the subspecies level of. A zone of clinal variability of morphological traits is noted in Transbaikalia. It obviously resulted from secondary contact between the western and eastern groups of populations, which were separated during the last glaciation and then rejoined.     

Molecular phylogeny of Vipera Laurenti, 1768 and the related genera Macrovipera (Reuss, 1927) and Daboia (Gray, 1842), with comments about neurotoxic Vipera aspis aspis populations

We used mtDNA sequences (cytochrome b and NADH dehydrogenase subunit 2) to reconstruct molecular phylogenies of Vipera sensu lato, Vipera sensu stricto, and Vipera aspis. Three major clades were identified within the Vipera s.l. group: (1) the European vipers, (2) the oriental vipers, consisting of Montivipera (Vipera 2) plus Macrovipera lebetina, and (3) a group of Asian and North African vipers consisting of Daboia russelii, V. palaestinae, and Macrovipera mauritanica. We also distinguished three clades within the monophyletic European Vipera group: V. ammodytes, V. aspis, and V. latastei, and Pelias with monophyly of Vipera 1 uncertain. Within V. aspis, the specimens collected in France formed the sister group of an Italian clade. The "neurotoxic" French population of V. aspis, which has a specific venom profile, separated from other French V. aspis early in the history of this group.     

Are the opposing selection pressures on exposed and protected shores sufficient to maintain genetic differentiation between gastropod populations with high intermigration rates?

Throughout the world intertidal gastropods living on exposed rocky shores differ strikingly in a number of morphological and life history traits from those on protected shores. Where surf is heavy gastropods tend to be smaller and to have thinner and smoother shells with larger apertures than do those from sheltered areas where crab predation is more intense. These morphological differences can occur within a species and there is evidence that they can be partially genetic and partially environmental. In addition the convergence of shell features in each habitat suggests that there are consistent differences between the selective pressures on exposed shores and the selective pressures on protected shores. I constructed a simulation model for a polygenic trait that experiences different selective pressures on exposed and sheltered shores. The results show that genetic differences can be maintained between the two populations despite high intermigration rates. Replacement of a portion of the random environmental variance with adaptive environmental variance reduces the effect of selection and thus the size of the difference maintained between the two populations. Genetic differentiation between exposed and protected populations can persist for significant periods of time and may have sometimes been the first step in speciation.     

Contrasting post-settlement selection results in many-to-one mapping of high performance phenotypes in the Hawaiian waterfall-climbing goby <Emphasis Type="Italic">Sicyopterus stimpsoni</Emphasis>

Natural selection drives adaptive evolution, but contrasting environmental pressures may lead to trade-offs between phenotypes that confer different performances. Such trade-offs may weaken the strength of selection and/or generate complex fitness surfaces with multiple local optima that correspond to different selection regimes. We evaluated how differences in patterns of phenotypic selection might promote morphological differences between subpopulations of the amphidromous Hawaiian waterfall-climbing goby, Sicyopterus stimpsoni. We conducted laboratory experiments on fish from the islands of Kaua‘i and Hawai‘i (the “Big Island”) to compare patterns of linear and nonlinear selection, and the opportunity for selection, that result from two contrasting pressures, predator evasion and waterfall climbing, which vary in intensity between islands. We found directional and nonlinear selection were strongest when individuals were exposed to their primary selective pressures (predator evasion on Kaua‘i, waterfall climbing on the Big Island). However, the opportunity for selection was greater for the non-primary pressure: climbing on Kaua‘i, predator evasion on the Big Island. Canonical rotation of the nonlinear gamma matrix demonstrated that individuals from Kaua‘i and the Big Island occupy regions near their local fitness peaks for some traits. Therefore, selection for predator evasion on Kaua‘i and climbing on the Big Island may be less effective in promoting morphological changes in this species, because variation of functionally important traits in their respective environments may have been reduced by directional or stabilizing selection. These results demonstrate that despite constraints on the opportunities for selection, population differences in phenotypic traits can arise due to differences in selective regimes. For S. stimpsoni, sufficient variation exists in other locomotor traits, allowing for necessary levels of performance in the contrasting selective regime (i.e., climbing on Kaua‘i and predator evasion on the Big Island) through many-to-one-mapping, which may be essential for the survival of local populations in an evanescent island environment.     

Assessing concurrent patterns of environmental niche and morphological evolution among species of horned lizards (Phrynosoma)

The prediction that variation in species morphology is related to environmental features has long been of interest to ecologists and evolutionary biologists. Many studies have demonstrated strong associations between morphological traits and local habitat characteristics, but few have considered the extent to which morphological traits may be associated with environmental features across broad geographic areas. Here, we use morphological, environmental and phylogenetic data compiled from Phrynosoma species to examine morphological and climatic variation across the geographic ranges of these species in an evolutionary context. We find significant phylogenetic signal in species’ environmental niches, but not in morphological traits. Furthermore, we demonstrate a significant correlation between species’ environmental niches and morphological traits when phylogenetic history is accounted for in the analysis. Our results suggest the importance of climatic variables in influencing morphological variation among species, and have implications for understanding how species distributions are constrained by environmental variation.     

Phylogeny of the vipers of the Caucasus (Reptilia, Viperidae)

Phylogenetic relationships for five taxa of Palearctic vipers (genus Vipera ) from the Caucasian region were revealed by cladistic analyses of separate and combined morphological and biochemical characters. The different data sets yielded largely congruent cladograms. Vipera berus from Sweden was included as an ingroup and V. aspis was used for outgroup comparison. For V. kaznakovi and V. dinniki , three and four different sub-populations, respectively, were treated as independent terminal taxa in the analyses. The most parsimonious cladograms confirmed the systematic positions of these populations, discussed in a recent study, and support the hypothesis that the montane populations of the western main Caucasus comprise one polymorphie species: V. dinniki.
Analyses of combined biochemical and morphological data generated two equally parsimonious cladograms (for all ingroups compared), but yielded only one fully resolved topology when ingroups were condensed to the species level: ( berus ((renardi ('ursinii'-eriwanensis )))( dinnikikaznakovi )).     

When do meadow vipers (Vipera ursinii) become sexually dimorphic? – ontogenetic patterns of sexual size dimorphisms

Contrary to an increasing number of papers that document sexual dimorphism in size (and/or shape) in adults, studies dealing with sex differences in newborn and juvenile snakes are surprisingly scarce. Data about ontogenetic shifts in sexual dimorphism are generally lacking and hence, it is unclear whether sex differences are set at birth or arise post‐natally. In this study, we analyzed patterns of sexual dimorphism in body size, head dimensions and tail length (TL) among newborn, subadult and adult meadow vipers (Vipera ursinii) from the Bjelasica Mt. in Montenegro. Patterns of sexual size dimorphisms differed among traits. There was no significant difference in head dimension of males and females, but adult snakes were sexually dimorphic in body size. Sexual differences in TL were evident since birth but changed in degree throughout ontogeny. Neonate meadow vipers presented highly significant inter‐litter variation in the sexual dimorphism of all traits we have measured. Such family effects may have an important influence on extent of inter‐sexual differences in snakes and should be included in analyses of sexual dimorphism.     

Aposematism and crypsis are not enough to explain dorsal polymorphism in the Iberian adder

Aposematic organisms can show phenotypic variability across their distributional ranges. The ecological advantages of this variability have been scarcely studied in vipers. We explored this issue in Vipera seoanei, a species that exhibits five geographically structured dorsal colour phenotypes across Northern Iberia: two zigzag patterned (Classic and Cantabrica), one dorsal-strip patterned (Bilineata), one even grey (Uniform), and one melanistic (Melanistic). We compared predation rates (raptors and mammals) on plasticine models resembling each colour phenotype in three localities. Visual modelling techniques were used to infer detectability (i.e. conspicuousness) of each model type for visually guided predators (i.e. diurnal raptors). We hypothesize that predation rates will be lower for the two zigzag models (aposematism hypothesis) and that models with higher detectability would show higher predation rates (detectability hypothesis). Classic and Bilineata models were the most conspicuous, while Cantabrica and Uniform were the less. Melanistic presented an intermediate conspicuousness. Predation rate was low (3.24% of models) although there was variation in attack frequency among models. Zigzag models were scarcely predated supporting the aposematic role of the zigzag pattern in European vipers to reduce predation (aposematism hypothesis). From the non-zigzag models, high predation occurred on Bilineata and Melanistic models, and low on Uniform models, partially supporting our detectability hypothesis. These results suggest particular evolutionary advantages for non-zigzag phenotypes such as better performance of Melanistic phenotypes in cold environments or better crypsis of Uniform phenotypes. Polymorphism in V. seoanei may respond to a complex number of forces acting differentially across an environmental gradient.     

A standardized assessment of forest mammal communities reveals consistent functional composition and vulnerability across the tropics

The understanding of global diversity patterns has benefitted from a focus on functional traits and how they relate to variation in environmental conditions among assemblages. Distant communities in similar environments often share characteristics, and for tropical forest mammals, this functional trait convergence has been demonstrated at coarse scales (110–200 km resolution), but less is known about how these patterns manifest at fine scales, where local processes (e.g. habitat features and anthropogenic activities) and biotic interactions occur. Here, we used standardized camera trapping data and a novel analytical method that accounts for imperfect detection to assess how the functional composition of terrestrial mammal communities for two traits – trophic guild and body mass – varies across 16 protected areas in tropical forests and three continents, in relation to the extent of protected habitat and anthropogenic pressures. We found that despite their taxonomic differences, communities generally have a consistent trophic guild composition, and respond similarly to these factors. Insectivores were found to be sensitive to the size of protected habitat and surrounding human population density. Body mass distribution varied little among communities both in terms of central tendency and spread, and interestingly, community average body mass declined with proximity to human settlements. Results indicate predicted trait convergence among assemblages at the coarse scale reflects consistent functional composition among communities at the local scale, suggesting that broadly similar habitats and selective pressures shaped communities with similar trophic strategies and responses to drivers of change. These similarities provide a foundation for assessing assemblages under anthropogenic threats and sharing conservation measures.     

Spatio-temporal dynamics of multi-trophic communities reveal ecosystem-wide functional reorganization

Large-scale alterations in marine ecosystems as a response to environmental and anthropogenic pressures have been documented worldwide. Yet, these are primarily investigated by assessing abundance fluctuations of a few dominant species, which inadequately reflect ecosystem-wide changes. In addition, it is increasingly recognized that it is not species identity per se, but their traits that determine environmental responses, biological interactions and ecosystem functioning. In this study, we investigated long-term, spatio-temporal variability in trait composition across multiple organism groups to assess whether functional changes occur in a similar way across trophic levels and whether shifts in trait composition link to environmental change. We combined extensive trait datasets with long-term surveys (30–40 yr) of four organism groups (phytoplankton, zooplankton, benthic invertebrates and fish) in three environmentally distinct areas of a large marine ecosystem. We found similar temporal trajectories in the community weighted mean trait time-series of the different trophic groups, revealing ecosystem-wide functional changes. The traits involved and their dynamics differed between areas, concurrent with climate-driven changes in temperature and salinity, as well as more local dynamics in nutrients and oxygen. This finding highlights the importance of considering both global climate, as well as local external drivers when studying ecosystem changes. Using a multi-trophic trait-based approach, our study demonstrates the importance of integrating community functional dynamics across multiple trophic levels to capture ecosystem-wide responses which could, ultimately, help moving towards a holistic understanding, assessment and management of marine ecosystems.     

Patterns of morphological and genetic variability in UK populations of the shore crab, Carcinus maenas Linnaeus, 1758 (Crustacea: Decapoda: Brachyura)

Previous research has identified extensive inter-population variability in the morphology of the shore crab (Carcinus maenas L.). To determine the source of this variation (genetic or environmental), morphological and genetic data were analysed from crabs collected from eight sites around the coast of the UK. Ten morphometric traits were measured from over 800 crabs and the degree of morphological similarity among sites was calculated using multivariate techniques. Allozyme electrophoresis was used to investigate patterns of genetic similarity. Extensive morphological variability was detected: eight out of the ten morphometric traits analysed were useful when discriminating between crabs from each site. Discriminant function analysis revealed that over 35% of individuals could be classified to their site of origin on the basis of their morphology. In contrast, the allozyme analysis revealed low levels of genetic variability, both within the meta-population and among the crab population at each site. Pairwise comparisons revealed a moderate correlation between the degree of morphological and genetic similarity of crabs at each site, which suggests that the observed phenotypic variability has a genetic component. However, only around 20% of the phenotypic variability detected was associated with the patterns of genetic similarity. This means that patterns of morphological variability in this species are largely determined by the local environmental conditions: local factors could have a within-generation selective influence on mean trait values or C. maenas may exhibit phenotypic plasticity.     

Morphological variability and developmental instability in subpopulations of the Eurasian badger (Meles meles) in Denmark

Aim Local populations from different geographical regions may differ in the selection regimes to which they are exposed. Differences in environmental factors and population density may affect the relative importance of different selective forces (e.g. natural vs. sexual selection). We suggest a direction of investigation concerned with the developmental instability of morphological traits. The goal is to disclose putative small‐scale geographical differences in the evolutionary forces, which may be hard to detect. Location Craniometrical investigations were carried out on ninety‐eight skulls and teeth of the Eurasian badger (Meles meles) collected during the period 1995–97 from three different populations in Denmark. One of these thrives at low population density, whereas the two others are characterized by high local density. Methods The skulls were investigated for developmental instability (DI) using fluctuating asymmetry (FA) as its estimator. FA was measured on canines, molars, premolar teeth and other skull and mandible traits. For the statistical analyses, we applied nonparametric permutation tests. Results Evidence was found suggesting differentiation among populations in mean degree of FA, and the FA values measured on canines were higher in the high‐density populations. FA of the canines was significantly higher in males than females, in contrast to FA of the other traits. Evidence of a negative relationship between canine size and their FA was found, whereas no significant correlations were found between the molar and premolar teeth measures and their FA. Main conclusions Our results suggest that canines could be under directional selection stemming from intrasexual competition, which may be stronger in high‐density zones. The other teeth investigated seem to be under a stabilizing regime hence their FA is mainly affected by environmental stresses. The negative relationship between canine size and FA found in males suggests the capacity of badgers to respond in an evolutionary way to environmental changes, despite the low genetic variability previously found at the molecular level.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

Geographic variation in genital morphology of Ciulfina praying mantids

Geographic variation in morphological traits is widespread and important to our current understanding of evolutionary processes. Although male genitalia are perhaps the most divergent morphological traits in animals, geographic variation in genital traits has received little attention and the mechanism driving such variation is unclear. The species isolation hypothesis of genital evolution makes explicit predictions about geographic variation in genitalia predicting patterns of genital divergence that reflect the risk of mating with related but incompatible species. The sexual selection and pleiotropy hypotheses, however, predict general levels of geographic variation that reflect divergent sexual selection pressures or genetic drift. To test these predictions, we investigated geographic variation in genital morphology in the praying mantid genus Ciulfina (Mantodea: Liturgusidae) using elliptic Fourier analysis. We found significant levels of geographic variation in the genital morphology of four Ciulfina species irrespective of the relative proximity of different populations to contact zones with other species. These results reject the species isolation hypothesis, and instead support either the sexual selection or pleiotropy hypotheses to explain patterns of genital evolution in this genus.