Phylogenomic insights into the diversity and evolution of Palearctic vipers.

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of ∼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

Contribution of host species and pathogen clade to snake fungal disease hotspots in Europe.

Infectious diseases are influenced by interactions between host and pathogen, and the number of infected hosts is rarely homogenous across the landscape. Areas with elevated pathogen prevalence can maintain a high force of infection and may indicate areas with disease impacts on host populations. However, isolating the ecological processes that result in increases in infection prevalence and intensity remains a challenge. Here we elucidate the contribution of pathogen clade and host species in disease hotspots caused by Ophidiomyces ophidiicola, the pathogen responsible for snake fungal disease, in 21 species of snakes infected with multiple pathogen strains across 10 countries in Europe. We found isolated areas of disease hotspots in a landscape where infections were otherwise low. O. ophidiicola clade had important effects on transmission, and areas with multiple pathogen clades had higher host infection prevalence. Snake species further influenced infection, with most positive detections coming from species within the Natrix genus. Our results suggest that both host and pathogen identity are essential components contributing to increased pathogen prevalence.

Ophiodimyces ophiodiicola, Etiologic Agent of Snake Fungal Disease, in Europe since Late 1950s.

The fungus Ophiodimyces ophiodiicola is the etiologic agent of snake fungal disease. Recent findings date US occurrence at least as far back as 1945. We analyzed 22 free-ranging snakes with gross lesions consistent with snake fungal disease from museum collections from Europe. We found 5 positive samples, the oldest collected in 1959.

Phenological and elevational shifts of plants, animals and fungi under climate change in the European Alps.

Mountain areas are biodiversity hotspots and provide a multitude of ecosystem services of irreplaceable socio-economic value. In the European Alps, air temperature has increased at a rate of about 0.36°C decade-1 since 1970, leading to glacier retreat and significant snowpack reduction. Due to these rapid environmental changes, this mountainous region is undergoing marked changes in spring phenology and elevational distribution of animals, plants and fungi. Long-term monitoring in the European Alps offers an excellent natural laboratory to synthetize climate-related changes in spring phenology and elevational distribution for a large array of taxonomic groups. This review assesses the climatic changes that have occurred across the European Alps during recent decades, spring phenological changes and upslope shifts of plants, animals and fungi from evidence in published papers and previously unpublished data. Our review provides evidence that spring phenology has been shifting earlier during the past four decades and distribution ranges show an upwards trend for most of the taxonomic groups for which there are sufficient data. The first observed activity of reptiles and terrestrial insects (e.g. butterflies) in spring has shifted significantly earlier, at an average rate of -5.7 and -6.0 days decade-1 , respectively. By contrast, the first observed spring activity of semi-aquatic insects (e.g. dragonflies and damselflies) and amphibians, as well as the singing activity or laying dates of resident birds, show smaller non-significant trends ranging from -1.0 to +1.3 days decade-1 . Leaf-out and flowering of woody and herbaceous plants showed intermediate trends with mean values of -2.4 and -2.8 days decade-1 , respectively. Regarding species distribution, plants, animals and fungi (N = 2133 species) shifted the elevation of maximum abundance (optimum elevation) upslope at a similar pace (on average between +18 and +25 m decade-1 ) but with substantial differences among taxa. For example, the optimum elevation shifted upward by +36.2 m decade-1 for terrestrial insects and +32.7 m decade-1 for woody plants, whereas it was estimated to range between -1.0 and +11 m decade-1 for semi-aquatic insects, ferns, birds and wood-decaying fungi. The upper range limit (leading edge) of most species also shifted upslope with a rate clearly higher for animals (from +47 to +91 m decade-1 ) than for plants (from +17 to +40 m decade-1 ), except for semi-aquatic insects (-4.7 m decade-1 ). Although regional land-use changes could partly explain some trends, the consistent upward shift found in almost all taxa all over the Alps is likely reflecting the strong warming and the receding of snow cover that has taken place across the European Alps over recent decades. However, with the possible exception of terrestrial insects, the upward shift of organisms seems currently too slow to track the pace of isotherm shifts induced by climate warming, estimated at about +62 to +71 m decade-1 since 1970. In the light of these results, species interactions are likely to change over multiple trophic levels through phenological and spatial mismatches. This nascent research field deserves greater attention to allow us to anticipate structural and functional changes better at the ecosystem level.

A powerful long metabarcoding method for the determination of complex diets from faecal analysis of the European pond turtle (Emys orbicularis, L. 1758).

High-throughput sequencing has become an accurate method for the identification of species present in soil, water, faeces, gut or stomach contents. However, information at the species level is limited due to the choice of short barcodes and based on the idea that DNA is too degraded to allow longer sequences to be amplified. We have therefore developed a long DNA metabarcoding method based on the sequencing of short reads followed by de novo assembly, which can precisely identify the taxonomic groups of organisms associated with complex diets, such as omnivorous individuals. The procedure includes 11 different primer pairs targeting the COI gene, the large subunit of the ribulose-1,5-bisphosphate carboxylase gene, the maturase K gene, the 28S rRNA and the trnL-trnF chloroplastic region. We validated this approach using 32 faeces samples from an omnivorous reptile, the European pond turtle (Emys orbicularis, L. 1758). This metabarcoding approach was assessed using controlled experiments including mock communities and faecal samples from captive feeding trials. The method allowed us to accurately identify prey DNA present in the diet of the European pond turtles to the species level in most of the cases (82.4%), based on the amplicon lengths of multiple markers (168-1,379 bp, average 546 bp), and produced by de novo assembly. The proposed approach can be adapted to analyse various diets, in numerous conservation and ecological applications. It is consequently appropriate for detecting fine dietary variations among individuals, populations and species as well as for the identification of rare food items.

Greater topoclimatic control of above- versus below-ground communities.

Assessing the degree to which climate explains the spatial distributions of different taxonomic and functional groups is essential for anticipating the effects of climate change on ecosystems. Most effort so far has focused on above-ground organisms, which offer only a partial view on the response of biodiversity to environmental gradients. Here including both above- and below-ground organisms, we quantified the degree of topoclimatic control on the occurrence patterns of >1,500 taxa and phylotypes along a c. 3,000 m elevation gradient, by fitting species distribution models. Higher model performances for animals and plants than for soil microbes (fungi, bacteria and protists) suggest that the direct influence of topoclimate is stronger on above-ground species than on below-ground microorganisms. Accordingly, direct climate change effects are predicted to be stronger for above-ground than for below-ground taxa, whereas factors expressing local soil microclimate and geochemistry are likely more important to explain and forecast the occurrence patterns of soil microbiota. Detailed mapping and future scenarios of soil microclimate and microhabitats, together with comparative studies of interacting and ecologically dependent above- and below-ground biota, are thus needed to understand and realistically forecast the future distribution of ecosystems.

Population Genetic Analyses Using 10 New Polymorphic Microsatellite Loci Confirms Genetic Subdivision within the Olm, Proteus anguinus.

We provide a comparative population genetic study of the elusive amphibian, Proteus anguinus, by comparing the genetic diversity and divergence among 4 cave populations (96 individuals) sampled in the Dinaric Karst of Croatia. We developed 10 variable microsatellite markers using pyrosequencing and applied them to the 4 selected populations belonging to 4 different cave systems. The results showed strong genetic differentiation between the 4 caves corroborating with previous findings suggesting that Proteus might comprise several unrecognized taxa. Our results confirmed that gene flow should be high within the caves, whereas it is low between hydrographic systems since geological periods. Finally, we conclude that the high genetic subdivision suggests the necessity of treating the 4 studied Proteus populations as evolutionary significant units.

Hybridization patterns in two contact zones of grass snakes reveal a new Central European snake species.

Recent studies found major conflicts between traditional taxonomy and genetic differentiation of grass snakes and identified previously unknown secondary contact zones. Until now, little is known about gene flow across these contact zones. Using two mitochondrial markers and 13 microsatellite loci, we examined two contact zones. One, largely corresponding to the Rhine region, involves the western subspecies Natrix natrix helvetica and the eastern subspecies N. n. natrix, whereas in the other, more easterly, contact zone two lineages meet that are currently identified with N. n. natrix and N. n. persa. This second contact zone runs across Central Europe to the southern Balkans. Our analyses reveal that the western contact zone is narrow, with parapatrically distributed mitochondrial lineages and limited, largely unidirectional nuclear gene flow. In contrast, the eastern contact zone is very wide, with massive nuclear admixture and broadly overlapping mitochondrial lineages. In combination with additional lines of evidence (morphology, phylogeny, divergence times), we conclude that these differences reflect different stages in the speciation process and that Natrix helvetica should be regarded as a distinct species. We suggest a nomenclatural framework for presently recognized grass snake taxa and highlight the need for reconciling the conflicts between genetics and taxonomy.

Breeding system, shell size and age at sexual maturity affect sperm length in stylommatophoran gastropods.

BACKGROUND: Sperm size and quality are key factors for fertilization success. There is increasing empirical evidence demonstrating that sperm form and function are influenced by selective pressures. Theoretical models predict that sperm competition could favour the evolution of longer sperm. In hermaphrodites, self-fertilizing species are expected to have shorter sperm than cross-fertilizing species, which use sperm stored from several mating partners for the fertilization of their eggs and thus are exposed to intense sperm competition. We tested this hypothesis by comparing original data on sperm length in 57 species of simultaneously hermaphroditic stylommatophoran gastropods from Europe and South America with respect to the species' breeding system. We used 28S rRNA nuclear and COI mitochondrial sequence data to construct a molecular phylogeny. Phylogenetic generalized linear models were applied to examine the potential influence of morphological and life-history characters. RESULTS: The best-fit model revealed that the breeding system and age at sexual maturity influence sperm length in gastropods. In general, species with predominant cross-fertilization had longer sperm than species with predominant self-fertilization or a mixed breeding system. Across species with shells (snails), sperm length also increased with shell size. CONCLUSIONS: Our study provides evidence that sperm length in stylommatophoran gastropods is influenced by the risk of sperm competition, as well as by age at sexual maturity and shell size. This finding extends present knowledge of sperm evolution to a group of so far poorly studied simultaneous hermaphrodites.

DNA Quantity and Quality in Remnants of Traffic-Killed Specimens of an Endangered Longhorn Beetle: A Comparison of Different Methods.

The sampling of living insects should be avoided in highly endangered species when the sampling would further increase the risk of population extinction. Nonlethal sampling (wing clips or leg removals) can be an alternative to obtain DNA of individuals for population genetic studies. However, nonlethal sampling may not be possible for all insect species. We examined whether remnants of traffic-killed specimens of the endangered and protected flightless longhorn beetle Iberodorcadion fuliginator (L., 1758) can be used as a resource for population genetic analyses. Using insect fragments of traffic-killed specimens collected over 15 yr, we determined the most efficient DNA extraction method in relation to the state of the specimens (crushed, fragment, or intact), preservation (dried, airtight, or in ethanol), storage duration, and weight of the sample by assessing the quantity and quality of genomic DNA. A modified cetyltrimethyl ammonium bromide method provided the highest recovery rate of genomic DNA and the largest yield and highest quality of DNA. We further used traffic-killed specimens to evaluate two DNA amplification techniques (quantitative polymerase chain reaction [qPCR] and microsatellites). Both qPCR and microsatellites revealed successful DNA amplification in all degraded specimens or beetle fragments examined. However, relative qPCR concentration and peak height of microsatellites were affected by the state of specimen and storage duration but not by specimen weight. Our investigation demonstrates that degraded remnants of traffic-killed beetle specimens can serve as a source of high-quality genomic DNA, which allows to address conservation genetic issues.

Diversifying selection and color-biased dispersal in the asp viper.

BACKGROUND: The presence of intraspecific color polymorphism can have multiple impacts on the ecology of a species; as a consequence, particular color morphs may be strongly selected for in a given habitat type. For example, the asp viper (Vipera aspis) shows a high level of color polymorphism. A blotched morph (cryptic) is common throughout its range (central and western Europe), while a melanistic morph is frequently found in montane populations, presumably for thermoregulatory reasons. Besides, rare atypical uniformly colored individuals are known here and there. Nevertheless, we found in a restricted treeless area of the French Alps, a population containing a high proportion (>50%) of such specimens. The aim of the study is to bring insight into the presence and function of this color morph by (i) studying the genetic structure of these populations using nine microsatellite markers, and testing for (ii) a potential local diversifying selection and (iii) differences in dispersal capacity between blotched and non-blotched vipers. RESULTS: Our genetic analyses support the occurrence of local diversifying selection for the non-blotched phenotype. In addition, we found significant color-biased dispersal, blotched individuals dispersing more than atypical individuals. CONCLUSION: We hypothesize that, in this population, the non-blotched phenotype possess an advantage over the typical one, a phenomenon possibly due to a better background matching ability in a more open habitat. In addition, color-biased dispersal might be partly associated with the observed local diversifying selection, as it can affect the genetic structure of populations, and hence the distribution of color morphs.

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central-marginal hypothesis.

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central-marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold-adapted species likely used two isolated glacial refugia in southern France, in permafrost-free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.

Influence of climate on the presence of colour polymorphism in two montane reptile species.

The coloration of ectotherms plays an important role in thermoregulation processes. Dark individuals should heat up faster and be able to reach a higher body temperature than light individuals and should therefore have benefits in cool areas. In central Europe, montane local populations of adder (Vipera berus) and asp viper (Vipera aspis) exhibit a varying proportion of melanistic individuals. We tested whether the presence of melanistic V. aspis and V. berus could be explained by climatic conditions. We measured the climatic niche position and breadth of monomorphic (including strictly patterned individuals) and polymorphic local populations, calculated their niche overlap and tested for niche equivalency and similarity. In accordance with expectations, niche overlap between polymorphic local populations of both species is high, and even higher than that of polymorphic versus monomorphic montane local populations of V. aspis, suggesting a predominant role of melanism in determining the niche of ectothermic vertebrates. However, unexpectedly, the niche of polymorphic local populations of both species is narrower than that of monomorphic ones, indicating that colour polymorphism does not always enable the exploitation of a greater variability of resources, at least at the intraspecific level. Overall, our results suggest that melanism might be present only when the thermoregulatory benefit is higher than the cost of predation.

High genetic differentiation among French populations of the Orsini's viper (Vipera ursinii ursinii) based on mitochondrial and microsatellite data: implications for conservation management.

The Orsini's viper (Vipera ursinii) is one of the most threatened snakes in Europe due to its highly fragmented distribution and specific open environment (steppic habitat) requirement. French populations are isolated on top of mountain massifs of the southern Prealps/Alps. Mitochondrial sequences (cytochrome b) and 6 microsatellite loci have been used to estimate the levels of genetic diversity and isolation within and among 11 French fragmented populations (a total of 157 individuals). Eleven cytochrome b haplotypes with a limited divergence were observed (mean divergence between haplotypes: 0.31%). However, we detected considerable genetic differentiation among populations (global F(ST) = 0.76 and 0.26 for mitochondrial and nuclear DNA, respectively). Results indicate that 3 populations possibly went through a bottleneck and 1 population showed low genetic diversity compared with the others. Although a significant isolation by distance was detected for both markers, strong differentiation was also observed between geographically close populations, probably due to the ragged landscape that constitutes a serious barrier to gene flow owing to the limited dispersal capability of the viper. Despite some discrepancies between the 2 markers, 8 Management Units have been identified and should be considered for future management projects.

Male reproductive success and multiple paternity in wild, low-density populations of the adder (Vipera berus).

We studied for the first time the occurrence of multiple paternity, male reproductive success, and neonate survival in wild, low-density adder (Vipera berus) populations using 13 microsatellite loci. Paternity was assigned for 15 clutches, collected during 3 years. Our data demonstrated that multiple paternity can occur at a high level (69%) in natural populations of V. berus, even if the density of adults is low. The high proportion of multiple sired clutches was comparable to the proportion observed in captive populations. Male reproductive success significantly increased with body length, and only the largest males successfully sired entire clutches. Finally, no relationship was detected between the number of fathers per clutch and neonate survival. These results suggest that multiple matings could be beneficial in populations with high level of inbreeding or low male fecundity.