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.     

Geographical patterns of morphological variation and environmental correlates in contact zones: a multi-scale approach using two Mediterranean vipers (Serpentes)

Adaptation to similar selective pressures can explain morphological convergence between closely related species in contact zones. Geostatistics and Geographical Information System were used to identify multi-scale patterns of morphological variability and test the hypothesis of morphological convergence due to local environmental pressures in the contact areas between Vipera aspis and Vipera latastei in the Iberian Peninsula. Nine morphological traits from 630 and 362 vipers at regional and local scale, respectively, were interpolated by Kriging to generate surfaces of morphological variation. Kriging is a geostatistical algorithm that allows investigating the spatial structure of data with statistical models. At both scales, a convergent north–south pattern in morphological variability was observed and the contact areas were identified as integration zones where intermediate vipers are found. Significant correlations were found between surfaces of univariate and multivariate traits, with precipitation and temperature seasonality. Thus, several morphological traits were apparently under local environmental selection. Nevertheless, the influence of biotic pressures and gene flow on morphological convergence of vipers in contact zones deserves further study.     

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.     

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.     

Experimental studies of adaptation in Clarkia xantiana. I. Sources of trait variation across a subspecies border

Abstract Both genetic differentiation and phenotypic plasticity might be expected to affect the location of geographic range limits. Co‐gradient variation (CoGV), plasticity that is congruent with genetic differentiation, may enhance performance at range margins, whereas its opposite, counter‐gradient variation (CnGV) may hinder performance. Here we report findings of reciprocal transplant experiments intended to tease apart the roles of differentiation and plasticity in producing phenotypic variation across a geographic border between two plant subspecies. Clarkia xantiana ssp. xantiana and C. xantiana ssp. parviflora are California‐endemic annuals that replace each other along a west‐east gradient of declining precipitation. We analyzed variation in floral traits, phenological traits, and vegetative morphological and developmental traits by sowing seeds of 18 populations (six of ssp. xantiana and 12 of ssp. parviflora) at three sites (one in each subspecies' exclusive range and one in the subspecies' contact zone), in two growing seasons (an exceptionally wet El Niño winter and a much drier La Niña winter). Significant genetic differences between subspecies appeared in 11 of 12 traits, and differences were of the same sign as in nature. These findings are consistent with the hypothesis that selection is responsible for subspecies differences. Geographic variation within subspecies over part of the spatial gradient mirrored between‐subspecies differences present at a larger scale. All traits showed significant plasticity in response to spatial and temporal environmental variation. Plasticity patterns ranged from spatial and temporal CoGV (e.g., in node of first flower), to spatial CnGV (e.g., in flowering time), to patterns that were neither CoGV nor CnGV (the majority of traits). Instances of CoGV may reflect adaptive plasticity and may serve to increase performance under year‐to‐year environmental variation and at sites near the subspecies border. However, the presence of spatial CnGV in some critical traits suggests that subspecies ranges may also be constrained by patterns of plasticity.     

An agonist of CXCR4 induces a rapid recovery from the neurotoxic effects of Vipera ammodytes and Vipera aspis venoms

People bitten by Alpine vipers are usually treated with antivenom antisera to prevent the noxious consequences caused by the injected venom. However, this treatment suffers from a number of drawbacks and additional therapies are necessary. The venoms of Vipera ammodytes and of Vipera aspis are neurotoxic and cause muscle paralysis by inducing neurodegeneration of motor axon terminals because they contain a presynaptic acting sPLA₂ neurotoxin. We have recently found that any type of damage to motor axons is followed by the expression and activation of the intercellular signaling axis consisting of the CXCR4 receptor present on the membrane of the axon stump and of its ligand, the chemokine CXCL12 released by activated terminal Schwann cells. We show here that also V. ammodytes and V. aspis venoms cause the expression of the CXCL12-CXCR4 axis. We also show that a small molecule agonist of CXCR4, dubbed NUCC-390, induces a rapid regeneration of the motor axon terminal with functional recovery of the neuromuscular junction. These findings qualify NUCC-390 as a promising novel therapeutics capable of improving the recovery from the paralysis caused by the snakebite of the two neurotoxic Alpine vipers.

     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ between Iris hexagona and Iris fulva, and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, with I. hexagona alleles tending to be favoured over those of I. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.     

Movement of a Heliconius hybrid zone over 30 years: A Bayesian approach

Hybrid zones have long been of interest to biologists as natural laboratories where we can gain insight into the processes of adaptation and speciation. Repeated sampling of individual hybrid zones has been particularly useful in elucidating the dynamic balance between selection and dispersal that maintains most hybrid zones. Here, we revisit a hybrid zone between Heliconius erato butterflies in Panamá for a third time over more than 30 years. We combine a novel Bayesian extension of stepped‐cline hybrid zone models with environmental data to understand the genetic and environmental causes of cline dynamics in this species. The cline has continued to move west, likely due to dominance drive, but has slowed and broadened. Environmental analyses suggest that widespread deforestation in Panamá could be leading to decreased avian predation and relaxed selection, causing the observed changes in cline dynamics.     

Population structure in the barn swallow,Hirundo rustica: a comparison between neutral DNA markers and quantitative traits

Local adaptation to variable environments can generate clinal variation in morphological traits. Alternatively, similar patterns of clinal variation may be generated simply as a result of genetic drift/migration balance. Teasing apart these different processes is a continuing focus in evolutionary ecology. We compare genetic differentiation at molecular loci and quantitative traits to analyse the effect of these different processes in a morphological latitudinal cline of the barn swallow, Hirundo rustica, breeding across Europe. The results obtained show no structuring at neutral microsatellite loci, which contrasts with positive structuring at five quantitative morphometric traits. This supports the hypothesis that the observed morphometric cline in barn swallows is the result of selection acting in a spatially heterogeneous environment. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 306–314.     

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.     

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.     

Sex chromosome inversions enforce reproductive isolation across an avian hybrid zone

Across hybrid zones, the sex chromosomes are often more strongly differentiated than the autosomes. This is regularly attributed to the greater frequency of reproductive incompatibilities accumulating on sex chromosomes and their exposure in the heterogametic sex. Working within an avian hybrid zone, we explore the possibility that chromosome inversions differentially accumulate on the Z chromosome compared to the autosomes and thereby contribute to Z chromosome differentiation. We analyse the northern Australian hybrid zone between two subspecies of the long‐tailed finch (Poephila acuticauda), first described based on differences in bill colour, using reduced‐representation genomic sequencing for 293 individuals over a 1,530‐km transect. Autosomal differentiation between subspecies is minimal. In contrast, 75% of the Z chromosome is highly differentiated and shows a steep genomic cline, which is displaced 350 km to the west of the cline in bill colour. Differentiation is associated with two or more putative chromosomal inversions, each predominating in one subspecies. If inversions reduce recombination between hybrid incompatibilities, they are selectively favoured and should therefore accumulate in hybrid zones. We argue that this predisposes inversions to differentially accumulate on the Z chromosome. One genomic region affecting bill colour is on the Z, but the main candidates are on chromosome 8. This and the displacement of the bill colour and Z chromosome cline centres suggest that bill colour has not strongly contributed to inversion accumulation. Based on cline width, however, the Z chromosome and bill colour both contribute to reproductive isolation established between this pair of subspecies.     

A hybrid zone between Chorthippus parallelus parallelus and Chorthippus parallelus erythropus (Orthoptera: Acrididae): morphological and electrophoretic characters

Two subspecies of the grasshopper Chorthippus parullelus have distributions that abut in the Pyrenees. We have studied morphological and electrophoretic variation in II populations near Mont Louis, Pyrénées-Orientales, France. The two subspecies differ in several morphological characters and at one enzyme locus, esterase-2, and in this area they hybridize forming a cline in morphology less than 5 km wide. Examination of museum specimens suggests that this hybrid zone extends for most of the length of the Pyrenees, possibly breaking down in the extreme west. As the two subspecies are known to differ in mating behaviour this hybrid zone is well suited to the study of reinforcement.     

小叶锦鸡儿天然居群叶形态性状变异研究

为揭示小叶锦鸡儿(Caragana microphylla)天然居群叶形态性状的变异规律及其生态适应性特征,该研究以10个小叶锦鸡儿天然居群为对象,通过多重比较、巢式方差分析、相关性分析、聚类分析和主成分分析等方法,对7个叶形态性状进行分析。结果表明：(1)小叶锦鸡儿叶形态性状在居群内和居群间均存在极显著差异(P < 0.01),平均变异系数为10.13%,不同性状的变异幅度为6.23%~12.78%;平均叶形态性状的表型分化系数为43.62%,居群内变异(30.09%)大于居群间变异(24.91%),说明居群内是其叶形态性状变异的主要来源。(2)相关性分析表明,环境因子对小叶锦鸡儿的叶形态性状变异有很大的影响,在地理空间上主要呈现出沿海拔梯度的变异模式;主成分分析的结果显示,小叶宽、叶柄宽和叶柄长对小叶锦鸡儿叶形态变异起主导作用;利用欧式距离对小叶锦鸡儿居群进行UPGMA聚类分析结果显示,基于叶形态性状和环境因子可分别将小叶锦鸡儿10个居群分为3类和2类,Mantel检验结果表明,小叶锦鸡儿的叶形态性状变异不存在地理连续性。研究结果为小叶锦鸡儿的适应性进化和开发利用提供了理论依据。     

Natural hybridization in eastern-Mediterranean firs: The case of Abies borisii-regis

The genus Abies is represented in southern Balkans by A. alba, A. cephalonica, and A. borisii-regis. To infer the status of southern-Balkans firs, as well as the extent and patterns of introgression within this taxonomical complex, we analyzed genetic variation patterns of 29 indigenous fir populations in Bulgaria, Macedonia, Greece, and Calabria using a combination of maternally and biparentally inherited markers. Three mitochondrial lineages were observed, one comprising Calabrian populations and two distributed in the Balkans, coinciding with A. alba and A. cephalonica. The boundary between lineages is sharp; only two populations containing a mixture of haplotypes were found. Bayesian analysis of population structure based on seven nuclear microsatellite (nSSR) loci revealed the existence of two clusters whose proportions exhibited a latitudinal cline with a width of 2.3° ( ≈ 255 km). Populations in the center of the latitudinal cline exhibit the most symmetrical, the flattest, and the broadest distributions of cluster proportions within individual tree genomes. A neighbor-net network reflects the cline resulting from the Bayesian analysis. The observed variation patterns are not consistent with the hypothesis of A. borisii-regis as a monophyletic taxon or a stabilized hybridogenous species resulted from ancient hybridization; the taxon rather is a product of recent introgression.     

Taxonomy of the bean goose-pink-footed goose

The bean goose Anser fabalis and the pink-footed goose A. brachyrhynchus breed in the tundra and taiga zones of Eurasia and eastern Greenland, and the taxonomy of the group based on morphology has been controversial. We investigated the phylogenetic relationships within the bean goose–the pink-footed goose complex using mitochondrial control region sequences of 199 individuals collected from the breeding areas in the Palaearctic and Eastern Nearctic. We found three mitochondrial clades geographically distributed to (1) Greenland, Iceland and Svalbard (A. brachyrhynchus), (2) the eastern taiga zone (former subspecies A. fabalis middendorffii), and (3) the western taiga and the tundra zone (subspecies A. fabalis rossicus, serrirostris and fabalis). MtDNA phylogeny suggests that morphological affinities between the taxa, e.g. in the bill structure, result from convergent evolution due to adaptation to similar habitats. Although a latitudinal cline in morphology was observed, clear phylogenetic discontinuities exist in the taiga and tundra zones supporting a species status for brachyrhynchus and middendorffii.