Deep evolutionary lineages in a Western Mediterranean snake (Vipera latastei/monticola group) and high genetic structuring in Southern Iberian populations

Venom variation among the three subspecies of the North African mountain viper Vipera monticola (Saint-Girons 1954)

The North African mountain viper (Vipera monticola) is a medically relevant venomous snake distributed in Morocco, Algeria, and Tunisia. Three subspecies of V. monticola, exhibiting differences in morphotypes and dietary regimes, are currently recognised: V. m. monticola, V. m. atlantica, and V. m. saintgironsi. Through the application of snake venomics, we analysed the venoms of specimens of Moroccan origin belonging to each of the three subspecies. Snake venom metalloproteinase (svMP), snake venom serine protease (svSP), C-type lectin and C-type lectin-related proteins (CTL), and phospholipase A2 (PLA2) were predominant, with PLA2 being the most abundant toxin family overall. Disintegrins (DI) and cysteine-rich secretory proteins (CRISP) were exclusive to V. m. monticola and V. m. atlantica, while l-amino-acid oxidases (LAAO) were only found in V. m. saintgironsi. The differences detected in the venom profiles, as well as in presence/absence and relative abundances of toxin families, indicate the occurrence of intraspecific venom variation within V. monticola. The identified patterns of venom similarity between subspecies seem to align more with their phylogenetic relationships than with the reported differences in their feeding habits.

Isolation and characterization of polymorphic microsatellite loci for the three Iberian vipers, Vipera aspis, V. Latastei and V. seoanei by Illumina MiSeq sequencing

BACKGROUND

European vipers (genus Vipera) are a well-studied taxonomic group, but the low resolution of nuclear sanger-sequenced regions has precluded thorough studies at systematic, ecological, evolutionary and conservation levels. In this study, we developed novel microsatellite markers for the three Iberian vipers, Vipera aspis, V. latastei and V. seoanei, and assessed their polymorphism in north-central Iberian populations.

METHODS AND RESULTS

Genomic libraries were developed for each species using an Illumina Miseq sequencing approach. From the 70 primer pairs initially tested, 48 amplified reliably and were polymorphic within species. Cross-species transferability was achieved for 31 microsatellites loci in the three target species and four additional loci that were transferable to one species only. The 48 loci amplified in average seven alleles, and detected average expected and observed heterozygosities of 0.7 and 0.55, in the three genotyped populations/species (26 V. aspis, 20 V. latastei and 10 V. seoanei).

CONCLUSIONS

Our study provides a selection of 48 polymorphic microsatellite markers that will contribute significantly to current knowledge on genetic diversity, gene flow, population structure, demographic dynamics, systematics, reproduction and heritability in these species, and potentially in other congeneric taxa.

Interpopulational variation and ontogenetic shift in the venom composition of Lataste's viper (Vipera latastei, Boscá 1878) from northern Portugal

Lataste's viper (Vipera latastei) is a venomous European viper endemic to the Iberian Peninsula, recognised as medically important by the World Health Organization. To date, no comprehensive characterisation of this species' venom has been reported. Here, we analysed the venoms of juvenile and adult specimens of V. latastei from two environmentally different populations from northern Portugal. Using bottom-up venomics, we produced six venom proteomes (three per population) from vipers belonging to both age classes (i.e., two juveniles and four adults), and RP-HPLC profiles of 54 venoms collected from wild specimens. Venoms from juveniles and adults differed in their chromatographic profiles and relative abundances of their toxins, suggesting the occurrence of ontogenetic changes in venom composition. Specifically, snake venom metalloproteinase (SVMP) was the most abundant toxin family in juvenile venoms, while snake venom serine proteinases (SVSPs), phospholipases A₂ (PLA₂s), and C-type lectin-like (CTLs) proteins were the main toxins comprising adult venoms. The RP-HPLC venom profiles were found to vary significantly between the two sampled localities, indicating geographic variability. Furthermore, the presence/absence of certain peaks in the venom chromatographic profiles appeared to be significantly correlated also to factors like body size and sex of the vipers. Our findings show that V. latastei venom is a variable phenotype. The intraspecific differences we detected in its composition likely mirror changes in the feeding ecology of this species, taking place during different life stages and under different environmental pressures. SIGNIFICANCE: Lataste's viper (Vipera latastei) is a medically important viper endemic to the Iberian Peninsula, inhabiting different habitats and undergoing a marked ontogenetic dietary shift. In the current study, we report the first proteomic analysis of V. latastei venom from two environmentally different localities in northern Portugal. Our bottom-up venomic analyses show that snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins are the major components of adult V. latastei venom. The comparative analysis of young and adult venoms suggests the occurrence of ontogenetic shift in toxin abundances, with snake venom metalloproteinases (SVMPs) being the predominant toxins in juvenile venoms. Moreover, geographic venom variation between the two studied populations is also detected, with our statistical analyses suggesting that factors like body size and sex of the vipers are possibly at play in its determination. Our work represents the first assessment of the composition of V. latastei venom, and the first step towards a better understanding of the drivers behind its variability.

Origin, extinction and ancient DNA of a new fossil insular viper: molecular clues of overseas immigration

Viperinae is a subfamily of viperid snakes whose fossil record in the Mediterranean islands is, until now, restricted to 12 palaeontological deposits on seven islands. Revision of the material excavated 30 years ago from the Middle/Late Pleistocene–Holocene deposit of Es Pouàs [Eivissa (= Ibiza), Balearic Islands, western Mediterranean] revealed about 6000 bones of a small-sized viper across different stratigraphic levels. Its morphological characteristics are different enough to known species of Vipera to warrant the description of a new species, but the nearly complete mitochondrial genome obtained from this snake based on a sample dated to 16 130 ± 45 bp, suggested it belonged to a new insular population of Lataste’s viper (Vipera latastei), Vipera latastei ebusitana subsp. nov. Phylogenetic analysis indicates that the dispersal of the ancestors of V. l. ebusitana to Eivissa, most probably from a north-east Iberian population, occurred via overwater colonization < 1.5 Mya, well after the Messinian Salinity Crisis (5.97–5.32 Mya) when land bridges allowed terrestrial colonization of the Balearic Islands by mainland faunas. The morphological differences between V. l. ebusitana and the Iberian populations suggest that it is a new dwarf taxon resulting from insular evolutionary processes, becoming extinct shortly after the first human arrival to this island about 4000 years ago.

Genes on the edge: A framework to detect genetic diversity imperiled by climate change

Ongoing global warming is disrupting several ecological and evolutionary processes, spanning different levels of biological organization. Species are expected to shift their ranges as a response to climate change, with relevant implications to peripheral populations at the trailing and leading edges. Several studies have analyzed the exposure of species to climate change but few have explored exposure at the intraspecific level. We introduce a framework to forecast exposure to climate change at the intraspecific level. We build on existing methods by combining correlative species distribution models, a model of species range dynamics, and a model of phylogeographic interpolation. We demonstrate the framework by applying it to 20 Iberian amphibian and reptile species. Our aims were to: (a) identify which species and intraspecific lineages will be most exposed to future climate change; (b) test if nucleotide diversity at the edges of species ranges are significantly higher or lower than on the overall range; and (c) analyze if areas of higher species gain, loss, and turnover coincide with those predicted for lineages richness and nucleotide diversity. We found that about 80% of the studied species are predicted to contract their range. Within each species, some lineages were predicted to contract their range, while others were predicted to maintain or expand it. Therefore, estimating the impacts of climate change at the species level only can underestimate losses at the intraspecific level. Some species had significant high amount of nucleotide at the trailing or leading edge, or both, but we did not find a consistent pattern across species. Spatial patterns of species richness, gain, loss, and turnover were fairly concurrent with lineages richness and nucleotide diversity. Our results support the need for increased attention to intraspecific diversity regarding monitoring and conservation strategies under climate change.

Chasing the phantom: biogeography and conservation of Vipera latastei-monticola in the Maghreb (North Africa)

The Maghreb region (North Africa) constitutes a major component of the Mediterranean Basin biodiversity hotspot. During the last centuries, a consistent human population growth has led to an unprecedented rate of habitat transformation and loss in the region and thus, threatening its biodiversity. The Western Mediterranean viper Vipera latastei-monticola inhabits humid and subhumid areas in the main mountain ranges of the Maghreb, facing such threatening factors; however, its elusive character and rarity hindered data collection for distinct biological purposes. Here, we study the biogeographical patterns and conservation status of the Maghrebian V. latastei-monticola resulting from recent sampling campaigns in Morocco and Tunisia. We update species distribution, and integrate phylogeographic and ecological niche modelling analyses at both species and lineage level to identify suitable areas, and to evaluate the impact of anthropogenic transformation and level of protection of their suitable space. We identified four highly divergent mitochondrial lineages, including a new lineage endemic to the Western High Atlas, with allopatric distributions and restricted to mountain ranges, supporting the role of mountains as past climatic refugia. Despite the remoteness of suitable areas, we report widespread habitat degradation and identify the low effectiveness of the current protected areas system in preserving the species and lineages range. Our study shows the urgent need to apply management actions for the long-term conservation of this vulnerable species and suggests a revaluation of the specific status of V. monticola, as these populations likely represent an ecotype of V. latastei.

Integrative inference of population history in the Ibero-Maghrebian endemic Pleurodeles waltl (Salamandridae)

Inference of population histories from the molecular signatures of past demographic processes is challenging, but recent methodological advances in species distribution models and their integration in time-calibrated phylogeographic studies allow detailed reconstruction of complex biogeographic scenarios. We apply an integrative approach to infer the evolutionary history of the Iberian ribbed newt (Pleurodeles waltl), an Ibero-Maghrebian endemic with populations north and south of the Strait of Gibraltar. We analyzed an extensive multilocus dataset (mitochondrial and nuclear DNA sequences and ten polymorphic microsatellite loci) and found a deep east-west phylogeographic break in Iberian populations dating back to the Plio-Pleistocene. This break is inferred to result from vicariance associated with the formation of the Guadalquivir river basin. In contrast with previous studies, North African populations showed exclusive mtDNA haplotypes, and formed a monophyletic clade within the Eastern Iberian lineage in the mtDNA genealogy. On the other hand, microsatellites failed to recover Moroccan populations as a differentiated genetic cluster. This is interpreted to result from post-divergence gene flow based on the results of IMA2 and Migrate analyses. Thus, Moroccan populations would have originated after overseas dispersal from the Iberian Peninsula in the Pleistocene, with subsequent gene flow in more recent times, implying at least two trans-marine dispersal events. We modeled the distribution of the species and of each lineage, and projected these models back in time to infer climatically favourable areas during the mid-Holocene, the last glacial maximum (LGM) and the last interglacial (LIG), to reconstruct more recent population dynamics. We found minor differences in climatic favourability across lineages, suggesting intraspecific niche conservatism. Genetic diversity was significantly correlated with the intersection of environmental favourability in the LIG and LGM, indicating that populations of P. waltl are genetically more diverse in regions that have remained environmentally favourable through the last glacial cycle, particularly southern Iberia and northern Morocco. This study provides novel insights into the relative roles of geology and climate on the biogeography of a biodiversity hotspot.

Patterns of Genetic Variability in Island Populations of the Cane Toad (Rhinella marina) from the Mouth of the Amazon

The Amazonian coast has several unique geological characteristics resulting from the interaction between drainage pattern of the Amazon River and the Atlantic Ocean. It is one of the most extensive and sedimentologically dynamic regions of the world, with a large number of continental islands mostly formed less than 10,000 years ago. The natural distribution of the cane toad (Rhinella marina), one of the world's most successful invasive species, in this complex Amazonian system provides an intriguing model for the investigation of the effects of isolation or the combined effects of isolation and habitat dynamic changes on patterns of genetic variability and population differentiation. We used nine fast-evolving microsatellite loci to contrast patterns of genetic variability in six coastal (three mainlands and three islands) populations of the cane toad near the mouth of the Amazon River. Results from Bayesian multilocus clustering approach and Discriminant Analyses of Principal Component were congruent in showing that each island population was genetically differentiated from the mainland populations. All FST values obtained from all pairwise comparisons were significant, ranging from 0.048 to 0.186. Estimates of both recent and historical gene flow were not significantly different from zero across all population pairs, except the two mainland populations inhabiting continuous habitats. Patterns of population differentiation, with a high level of population substructure and absence/restricted gene flow, suggested that island populations of R. marina are likely isolated since the Holocene sea-level rise. However, considering the similar levels of genetic variability found in both island and mainland populations, it is reliable to assume that they were also isolated for longer periods. Given the genetic uniqueness of each cane toad population, together with the high natural vulnerability of the coastal regions and intense human pressures, we suggest that these populations should be treated as discrete units for conservation management purposes.

PSMC analysis of effective population sizes in molecular ecology and its application to black-and-white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (Ne ) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in Ne over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black-and-white flycatchers. Ne curves of Atlas, collared, pied and semicollared flycatcher converged 1-2 million years ago (Ma) at an Ne of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate Ne trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100-200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in Ne differed among species, and for pied flycatcher, the temporal dynamics of Ne differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per-site filter of ≥10 reads and no more than 25% of missing data.

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species' evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.

Out of Africa: did Emys orbicularis occidentalis cross the Strait of Gibraltar twice?

The narrow Strait of Gibraltar has separated the African and European continents since the Miocene (5.3 Mya), with a different degree of permeability for Mediterranean taxa. Southern and northern regions of the Iberian Peninsula and Morocco, respectively, are key areas to evaluate the colonization dynamics and biogeographic history of taxa occurring at both sides of this strait. The Ibero-Maghrebian subspecies of the European pond turtle, Emys orbicularis occidentalis, is patchily distributed and threatened throughout most of the Iberian Peninsula and northern Morocco and its origin is thought to be in North Africa. Here we expand the geographic sampling across the Iberian Peninsula and Morocco, with special emphasis in the southern tip of the peninsula and northern Morocco, and analyze mtDNA sequences of 183 E. o. occidentalis to better understand the complex biogeographic history of this subspecies. We provide for the first time evidence for shared haplotypes of Iberian and North African pond turtles, with an additional haplotype in the southern Iberian Peninsula derived from Moroccan haplotypes. This supports the hypothesis that the Strait of Gibraltar constitutes no significant biogeographic barrier for E. orbicularis. However, the newly discovered shared, or extremely similar, haplotypes of European pond turtles from the southern Iberian Peninsula and Morocco suggest either that at least two independent natural colonization waves from Morocco have reached the Iberian Peninsula or that Moroccan turtles were accidentally or deliberately introduced there.

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three-way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.