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.

Variations in age structure and growth in congeners Lacerta viridis and Lacerta media.

Determining the age of any species allows it to be analyzed from the ontogenetic, demographic, and ecological perspectives. In the present study, we tested the hypothesis that the age structure of congener species (Lacerta media and Lacerta viridis) with the same ecological niche may vary in different areas. In this context, we applied skeletochronology method to reveal various demographic parameters, such as age structure, longevity, age at sexual maturity, growth rate, survival rate, adult life expectancy, and the relationship between age and body size in the green lizard, L. viridis, and the medium lizard, L. media. In L. media and L. viridis, the maximum lifespan was 10 and 8 years, respectively. The mean age and body size of females were significantly greater than those of males in L. media. However, in the examined L. viridis population, no appreciable variation in mean age or body size was found to exist between the sexes. It was estimated that the green lizards reach maturity at the age of 2 or 3 years. However, the L. media reached sexual maturity approximately 1 year later than the congener. The body size markedly increased with age in males for both studied populations. However, in females, body size positively increased with age only in L. media. The approach of skeletochronology that we utilized in this study to assess age structure makes it simple to gather a variety of time-dependent ecological data for such ectothermic species.

Extended Snake Venomics by Top-Down In-Source Decay: Investigating the Newly Discovered Anatolian Meadow Viper Subspecies, Vipera anatolica senliki.

Herein, we report on the venom proteome of Vipera anatolica senliki, a recently discovered and hitherto unexplored subspecies of the critically endangered Anatolian meadow viper endemic to the Antalya Province of Turkey. Integrative venomics, including venom gland transcriptomics as well as complementary bottom-up and top-down proteomics analyses, were applied to fully characterize the venom of V. a. senliki. Furthermore, the classical top-down venomics approach was extended to elucidate the venom proteome by an alternative in-source decay (ISD) proteomics workflow using the reducing matrix 1,5-diaminonaphthalene. Top-down ISD proteomics allows for disulfide bond counting and effective de novo sequencing-based identification of high-molecular-weight venom constituents, both of which are difficult to achieve by commonly established top-down approaches. Venom gland transcriptome analysis identified 96 toxin transcript annotations from 18 toxin families. Relative quantitative snake venomics revealed snake venom metalloproteinases (42.9%) as the most abundant protein family, followed by several less dominant toxin families. Online mass profiling and top-down venomics provide a detailed insight into the venom proteome of V. a. senliki and facilitate a comparative analysis of venom variability for the closely related subspecies, Vipera anatolica anatolica.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders.

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus' distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

Environmental temperatures shape thermal physiology as well as diversification and genome-wide substitution rates in lizards.

Climatic conditions changing over time and space shape the evolution of organisms at multiple levels, including temperate lizards in the family Lacertidae. Here we reconstruct a dated phylogenetic tree of 262 lacertid species based on a supermatrix relying on novel phylogenomic datasets and fossil calibrations. Diversification of lacertids was accompanied by an increasing disparity among occupied bioclimatic niches, especially in the last 10 Ma, during a period of progressive global cooling. Temperate species also underwent a genome-wide slowdown in molecular substitution rates compared to tropical and desert-adapted lacertids. Evaporative water loss and preferred temperature are correlated with bioclimatic parameters, indicating physiological adaptations to climate. Tropical, but also some populations of cool-adapted species experience maximum temperatures close to their preferred temperatures. We hypothesize these species-specific physiological preferences may constitute a handicap to prevail under rapid global warming, and contribute to explaining local lizard extinctions in cool and humid climates.

Intact protein mass spectrometry reveals intraspecies variations in venom composition of a local population of Vipera kaznakovi in Northeastern Turkey.

We report on the variable venom composition of a population of the Caucasus viper (Vipera kaznakovi) in Northeastern Turkey. We applied a combination of venom gland transcriptomics, de-complexing bottom-up and top-down venomics. In contrast to sole bottom-up venomics approaches and gel or chromatography based venom comparison, our combined approach enables a faster and more detailed comparison of venom proteomes from multiple individuals. In total, we identified peptides and proteins from 15 toxin families, including snake venom metalloproteinases (svMP; 37.8%), phospholipases A2 (PLA2; 19.0%), snake venom serine proteinases (svSP; 11.5%), C-type lectins (CTL; 6.9%) and cysteine-rich secretory proteins (CRISP; 5.0%), in addition to several low abundant toxin families. Furthermore, we identified intraspecies variations of the venom composition of V. kaznakovi, and find these were mainly driven by the age of the animals, with lower svSP abundance detected in juveniles. On the proteoform level, several small molecular weight toxins between 5 and 8 kDa in size, as well as PLA2s, drove the differences observed between juvenile and adult individuals. This study provides novel insights into the venom variability of V. kaznakovi and highlights the utility of intact mass profiling for fast and detailed comparison of snake venom. BIOLOGICAL SIGNIFICANCE: Population level and ontogenetic venom variation (e.g. diet, habitat, sex or age) can result in a loss of antivenom efficacy against snakebites from wide ranging snake populations. The current state of the art for the analysis of snake venoms are de-complexing bottom-up proteomics approaches. While useful, these have the significant drawback of being time-consuming and following costly protocols, and consequently are often applied to pooled venom samples. To overcome these shortcomings and to enable rapid and detailed profiling of large numbers of individual venom samples, we integrated an intact protein analysis workflow into a transcriptomics-guided bottom-up approach. The application of this workflow to snake individuals of a local population of V. kaznakovi revealed intraspecies variations in venom composition, which are primarily explained by the age of the animals, and highlighted svSP abundance to be one of the molecular drivers for the compositional differences observed.

Cytotoxicity of Silica Nanoparticles with Transcaucasian Nose-Horned Viper, Vipera ammodytes transcaucasiana, Venom on U87MG and SHSY5Y Neuronal Cancer Cells.

Highly bioactive compounds of the snake venom make them particular sources for anticancer agent development. They contain very rich peptide-protein structures. Therefore, they are very susceptible to environmental conditions such as temperature, pH, and light. In this study, Vipera ammodytes transcaucasiana venom was encapsulated in PAMAM-G4 dendrimer by sol-gel method in order to prevent degradation of venom contents from the environmental conditions. For this purpose, nanoparticles were prepared by sol-gel methodology and SEM analyses were performed. U87MG and SHSY5Y neuronal cancer cell lines were treated with different concentrations of venom-containing nanoparticles and cytotoxicity was determined by MTT assay. IC50 values of nanoparticles with snake venom were calculated as 37.24 and 44.64 µg/ml for U87MG and SHSY5Y cells, respectively. The IC50 values of nanoparticles with snake venom were calculated as 10.07 and 7.9 µg/ml for U87MG and SHSY5Y cells, respectively. As a result, nanoparticles with V. a. transcaucasiana venom showed remarkably high cytotoxicity. Encapsulation efficiency of nanoparticles with 1 mg/ml snake venom was determined as %67 via BCA™ protein analysis. In conclusion, this method is found to be convenient and useful for encapsulating snake venom as well as being suitable for drug delivery systems.

Comparative Venomics of the Vipera ammodytes transcaucasiana and Vipera ammodytes montandoni from Turkey Provides Insights into Kinship.

The Nose-horned Viper (Vipera ammodytes) is one of the most widespread and venomous snakes in Europe, which causes high frequent snakebite accidents. The first comprehensive venom characterization of the regional endemic Transcaucasian Nose-horned Viper (Vipera ammodytes transcaucasiana) and the Transdanubian Sand Viper (Vipera ammodytes montandoni) is reported employing a combination of intact mass profiling and bottom-up proteomics. The bottom-up analysis of both subspecies identified the major snake protein families of viper venoms. Furthermore, intact mass profiling revealed the presence of two tripeptidic metalloprotease inhibitors and their precursors. While previous reports applied multivariate analysis techniques to clarify the taxonomic status of the subspecies, an accurate classification of Vipera ammodytestranscaucasiana is still part of the ongoing research. The comparative analysis of the viper venoms on the proteome level reveals a close relationship between the Vipera ammodytes subspecies, which could be considered to clarify the classification of the Transcaucasian Nose-horned Viper. However, the slightly different ratio of some venom components could be indicating interspecific variations of the two studied subspecies or intraspecies alternations based on small sample size. Additionally, we performed a bioactivity screening with the crude venoms against several human cancerous and non-cancerous cell lines, which showed interesting results against a human breast adenocarcinoma epithelial cell line. Several fractions of Vipera a. transcaucasiana demonstrated a strong cytotoxic effect on triple negative MDA MB 231 breast cancer cells.

Biological activities of skin and parotoid gland secretions of bufonid toads (Bufo bufo, Bufo verrucosissimus and Bufotes variabilis) from Turkey.

Toad glandular secretions and skin extractions contain numerous natural agents which may provide unique resources for novel drug development. Especially the skin-parotoid gland secretions of toads from genus Bufo contain as many as 86 different types of active compounds, each with the potential of becoming a potent drug. In the present study, crude skin-parotoid gland secretions from Bufo bufo, Bufo verrucosissimus and Bufotes variabilis from Turkey were screened against various cancer cells together with normal cells using MTT assay. Furthermore, the antimicrobial properties of skin secretions were tested on selected bacterial and fungal species for assessing the possible medical applications. Antimicrobial activity of skin secretions was studied by determining minimal inhibitory concentration (MIC) in broth dilution method. Hemolytic activity of each skin-secretion was also estimated for evaluating pharmaceutical potential. Both skin-parotoid gland secretions showed high cytotoxic effect on all cancerous and non-cancerous cell lines with IC50 values varying between <0.1µg/ml and 6.02µg/ml. MIC results of antimicrobial activity tests were found to be between 3.9µg/ml and 250µg/ml. No hemolytic activities on rabbit red blood cells at concentrations between 0.5µg/ml and 50µg/ml were observed. In conclusion, skin-parotoid secretions of bufonid toads might be remarkable candidates for anti-cancer and antimicrobial agents without hemolytic activities.