Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus).

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Independent Recruitment of Different Types of Phospholipases A2 to the Venoms of Caenophidian Snakes: The Rise of PLA2-IIE within Pseudoboini (Dipsadidae).

Snake venoms harbor a wide and diverse array of enzymatic and nonenzymatic toxic components, allowing them to exert myriad effects on their prey. However, they appear to trend toward a few optimal compositional scaffolds, dominated by four major toxin classes: SVMPs, SVSPs, 3FTxs, and PLA2s. Nevertheless, the latter appears to be restricted to vipers and elapids, as it has never been reported as a major venom component in rear-fanged species. Here, by investigating the original transcriptomes from 19 species distributed in eight genera from the Pseudoboini tribe (Dipsadidae: Xenodontinae) and screening among seven additional tribes of Dipsadidae and three additional families of advanced snakes, we discovered that a novel type of venom PLA2, resembling a PLA2-IIE, has been recruited to the venom of some species of the Pseudoboini tribe, where it is a major component. Proteomic and functional analyses of these venoms further indicate that these PLA2s play a relevant role in the venoms from this tribe. Moreover, we reconstructed the phylogeny of PLA2s across different snake groups and show that different types of these toxins have been recruited in at least five independent events in caenophidian snakes. Additionally, we present the first compositional profiling of Pseudoboini venoms. Our results demonstrate how relevant phenotypic traits are convergently recruited by different means and from homologous and nonhomologous genes in phylogenetically and ecologically divergent snake groups, possibly optimizing venom composition to overcome diverse adaptative landscapes.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers.

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

MITGARD: an automated pipeline for mitochondrial genome assembly in eukaryotic species using RNA-seq data.

MOTIVATION: Over the past decade, the field of next-generation sequencing (NGS) has seen dramatic advances in methods and a decrease in costs. Consequently, a large expansion of data has been generated by NGS, most of which have originated from RNA-sequencing (RNA-seq) experiments. Because mitochondrial genes are expressed in most eukaryotic cells, mitochondrial mRNA sequences are usually co-sequenced within the target transcriptome, generating data that are commonly underused or discarded. Here, we present MITGARD, an automated pipeline that reliably recovers the mitochondrial genome from RNA-seq data from various sources. The pipeline identifies mitochondrial sequence reads based on a phylogenetically related reference, assembles them into contigs, and extracts a complete mtDNA for the target species. RESULTS: We demonstrate that MITGARD can reconstruct the mitochondrial genomes of several species throughout the tree of life. We noticed that MITGARD can recover the mitogenomes in different sequencing schemes and even in a scenario of low-sequencing depth. Moreover, we showed that the use of references from congeneric species diverging up to 30 million years ago (MYA) from the target species is sufficient to recover the entire mitogenome, whereas the use of species diverging between 30 and 60 MYA allows the recovery of most mitochondrial genes. Additionally, we provide a case study with original data in which we estimate a phylogenetic tree of snakes from the genus Bothrops, further demonstrating that MITGARD is suitable for use on biodiversity projects. MITGARD is then a valuable tool to obtain high-quality information for studies focusing on the phylogenetic and evolutionary aspects of eukaryotes and provides data for easily identifying a sample using barcoding, and to check for cross-contamination using third-party tools.

Total-evidence phylogeny and evolutionary morphology of New World pitvipers (Serpentes: Viperidae: Crotalinae).

Crotalines (pitvipers) in the Americas are distributed from southern Canada to southern Argentina, and are represented by 13 genera and 163 species that constitute a monophyletic group. Their phylogenetic relationships have been assessed mostly based on DNA sequences, while morphological data have scarcely been used for phylogenetic inquiry. We present a total-evidence phylogeny of New World pitvipers, the most taxon/character comprehensive phylogeny to date. Our analysis includes all genera, morphological data from external morphology, cranial osteology and hemipenial morphology, and DNA sequences from mitochondrial and nuclear genes. We performed analyses with parsimony as an optimality criterion, using different schemes for character weighting. We evaluated the contribution of the different sources of characters to the phylogeny through analyses of reduced datasets and calculation of weighted homoplasy and retention indexes. We performed a morphological character analysis to identify synapomorphies for the main clades. In terms of biogeography, our results support a single colonization event of the Americas by pitvipers, and a cladogenetic event into a Neotropical clade and a North American/Neotropical clade. The results also shed light on the previously unstable position of some taxa, although they could not sufficiently resolve the position of Bothrops lojanus, which may lead to the paraphyly of either Bothrops or Bothrocophias. The morphological character analyses demonstrated that an important phylogenetic signal is contained in characters related to head scalation, the jaws and the dorsum of the skull, and allowed us to detect morphological convergences in external morphology associated with arboreality.

Evolutionary allometry and ecological correlates of fang length evolution in vipers.

Traits for prey acquisition form the phenotypic interface of predator-prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems.

The phylogenetic position of ridley's worm lizard reveals the complex biogeographic history of New World insular amphisbaenids.

The archipelago of Fernando de Noronha (FN) is located in the Equatorial South Atlantic Ocean, at 375 km off the northeastern coast of Brazil. Its endemic vertebrate land fauna is restricted to only six species, and three main hypotheses have been proposed to explain their presence in the archipelago. These hypotheses suggest FN had alternative biogeographic connections with: 1) the West Indies; 2) the South America mainland; or 3) Africa. Here, we evaluate for the first time the phylogenetic position of Amphisbaena ridleyi within the diversity of Amphisbaenia, and we infer the biogeographic processes that explain its presence in FN and its relationship with amphisbaenids from the West Indies. We analyzed a comprehensive multilocus dataset for Amphisbaenidae using maximum likelihood and time-calibrated Bayesian phylogenetic approaches. Based on our time-calibrated tree, we tested different biogeographic scenarios through historical biogeographic analyses. Our phylogenetic results for the high-level relationships of Amphisbaenia can be parenthetically summarized as (Rhineuridae, (Blanidae, (Bipedidae, (Cadeidae, (Trogonophidae, (Amphisbaenidae))))). Nine highly supported groups of species were recovered among the mainland South American amphisbaenids (SAA), whereas two phylogenetically distant groups of species were inferred for the West Indies: 1) WIC01, an Oligocene lineage present in Cuba and Hispaniola, which is the sister group of all other SAA groups; and 2) WIC02, a Miocene lineage that is restricted to southern Hispaniola and is closely related to Am. ridleyi. We estimated two events of transatlantic dispersal of amphisbaenians from Africa to West Indies: the dispersal of Cadeidae during the Middle Eocene, and the dispersal of the ancestor of Amphisbaena during the transition Eocene/Oligocene. These events were likely affected by the North Equatorial and South Equatorial currents, respectively, which have been flowing westwards since the Paleocene. The ancestral cladogenesis of Amphisbaena during the Late Oligocene is likely related to overwater dispersal events, or alternatively can be associated with the fragmentation of GAARlandia, when WIC01 was isolated in the West Indies, while the remaining groups of Amphisbaena diversified throughout the South American continent. During the Late Miocene, the ancestor of WIC02 dispersed from northern South America to the West Indies, while Am. ridleyi dispersed from the same region to FN. The overwater dispersal of WIC02 was driven by the North Brazilian Current and the dispersal of Am. ridleyi was likely influenced by the periodic shifts in direction and strength estimated for the North Equatorial Counter-Current during the Late Miocene.

When the tail shakes the snake: phylogenetic affinities and morphology of Atractus badius (Serpentes: Dipsadidae), reveals some current pitfalls on the snake's genomic age.

Atractus badius has a long and controversial nomenclatural history due to both its antiquity and the brevity of its original description. This species was described based on two syntypes from Java. Later, a lectotype was designated and the distribution range restricted to the Guiana Shield. Although this species has been repeatedly recorded throughout Amazonia and the Andes, these records have been erroneously assigned to A. badius because of a considerable level of confusion in the literature. We found 13 additional specimens of this poorly known snake, expanding our knowledge on its morphological variability, phylogenetic relationships and distribution. In this paper, we report new localities and data on meristic, morphometric, coloration, scales micro-ornamentation, osteology and hemipenial morphology to A. badius. We also determine the phylogenetic position of A. badius as nested in a composite Guiana Shield clade. In addition, we discuss recent advances on the systematics of Atractus and comment on some 'new arrangements' with respect to previously recognized species groups and available morphological evidence. Finally, we highlight the importance of accurate voucher identification before promoting taxonomic changes or implementing nomenclatural acts derived from new phylogenetic hypotheses, and notably propose the synonymy of the recently described A. pyroni with A. roulei.

Replacement and Parallel Simplification of Nonhomologous Proteinases Maintain Venom Phenotypes in Rear-Fanged Snakes.

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components.

Phylogenomics, Biogeography, and Morphometrics Reveal Rapid Phenotypic Evolution in Pythons After Crossing Wallace's Line.

Ecological opportunities can be provided to organisms that cross stringent biogeographic barriers towards environments with new ecological niches. Wallace's and Lyddeker's lines are arguably the most famous biogeographic barriers, separating the Asian and Australo-Papuan biotas. One of the most ecomorphologically diverse groups of reptiles, the pythons, is distributed across these lines, and are remarkably more diverse in phenotype and ecology east of Lydekker's line in Australo-Papua. We used an anchored hybrid enrichment approach, with near complete taxon sampling, to extract mitochondrial genomes and 376 nuclear loci to resolve and date their phylogenetic history. Biogeographic reconstruction demonstrates that they originated in Asia around 38-45 Ma and then invaded Australo-Papua around 23 Ma. Australo-Papuan pythons display a sizeable expansion in morphological space, with shifts towards numerous new adaptive optima in head and body shape, coupled with the evolution of new micro-habitat preferences. We provide an updated taxonomy of pythons and our study also demonstrates how ecological opportunity following colonization of novel environments can promote morphological diversification in a formerly ecomorphologically conservative group. [Adaptive radiation; anchored hybrid enrichment; biogeography; morphometrics; snakes.].

Interrogating Genomic-Scale Data for Squamata (Lizards, Snakes, and Amphisbaenians) Shows no Support for Key Traditional Morphological Relationships.

Genomics is narrowing uncertainty in the phylogenetic structure for many amniote groups. For one of the most diverse and species-rich groups, the squamate reptiles (lizards, snakes, and amphisbaenians), an inverse correlation between the number of taxa and loci sampled still persists across all publications using DNA sequence data and reaching a consensus on the relationships among them has been highly problematic. In this study, we use high-throughput sequence data from 289 samples covering 75 families of squamates to address phylogenetic affinities, estimate divergence times, and characterize residual topological uncertainty in the presence of genome-scale data. Importantly, we address genomic support for the traditional taxonomic groupings Scleroglossa and Macrostomata using novel machine-learning techniques. We interrogate genes using various metrics inherent to these loci, including parsimony-informative sites (PIS), phylogenetic informativeness, length, gaps, number of substitutions, and site concordance to understand why certain loci fail to find previously well-supported molecular clades and how they fail to support species-tree estimates. We show that both incomplete lineage sorting and poor gene-tree estimation (due to a few undesirable gene properties, such as an insufficient number of PIS), may account for most gene and species-tree discordance. We find overwhelming signal for Toxicofera, and also show that none of the loci included in this study supports Scleroglossa or Macrostomata. We comment on the origins and diversification of Squamata throughout the Mesozoic and underscore remaining uncertainties that persist in both deeper parts of the tree (e.g., relationships between Dibamia, Gekkota, and remaining squamates; among the three toxicoferan clades Iguania, Serpentes, and Anguiformes) and within specific clades (e.g., affinities among gekkotan, pleurodont iguanians, and colubroid families).

Diversification in vipers: Phylogenetic relationships, time of divergence and shifts in speciation rates.

Snakes of the cosmopolitan family Viperidae comprise around 329 venomous species showing a striking heterogeneity in species richness among lineages. While the subfamily Azemiopinae comprises only two species, 70% of all viper species are arranged in the subfamily Crotalinae or the "pit vipers". The radiation of the pit vipers was marked by the evolution of the heat-sensing pits, which has been suggested to be a key innovation for the successful diversification of the group. Additionally, only crotalines were able to successfully colonize the New World. Here, we present the most complete molecular phylogeny for the family to date that comprises sequences from nuclear and mitochondrial genes representing 79% of all living vipers. We also investigated the time of divergence between lineages, using six fossils to calibrate the tree, and explored the hypothesis that crotalines have undergone an explosive radiation. Our phylogenetic analyses retrieved high support values for the monophyly of the family Viperidae, subfamilies Viperinae and Crotalinae, and 22 out of 27 genera, as well as well-supported intergeneric relationships throughout the family. We were able to recover a strongly supported sister clade to the New World pit vipers that comprises Gloydius, Ovophis, Protobothrops and Trimeresurus gracilis. Our results agree in many aspects with other studies focusing on the phylogenetics of vipers, but we recover new relationships as well. Despite the addition of new sequences we were not able to resolve some of the poor supported relationships previously suggested. Time of divergence estimates suggested that vipers started to radiate around the late Paleocene to middle Eocene with subfamilies most likely dating back to the Eocene. The invasion of the New World might have taken place sometime close to the Oligocene/Miocene boundary. Diversification analyses suggested a shift in speciation rates during the radiation of a sub-clade of pit vipers where speciation rates rapidly increased but slowed down toward the present. Thus, the evolution of the loreal pits alone does not seem to explain their explosive speciation rates. We suggest that climatic and geological changes in Asia and the invasion of the New World may have also contributed to the speciation shift found in vipers.

Consuming viscous prey: a novel protein-secreting delivery system in neotropical snail-eating snakes.

BACKGROUND: Efficient venom delivery systems are known to occur only in varanoid lizards and advanced colubroidean snakes among squamate reptiles. Although components of these venomous systems might have been present in a common ancestor, the two lineages independently evolved strikingly different venom gland systems. In snakes, venom is produced exclusively by serous glands in the upper jaw. Within the colubroidean radiation, lower jaw seromucous infralabial glands are known only in two distinct lineages-the basal pareatids and the more advanced Neotropical dipsadines known as "goo-eating snakes". Goo-eaters are a highly diversified, ecologically specialized clade that feeds exclusively on invertebrates (e.g., gastropod molluscs and annelids). Their evolutionary success has been attributed to their peculiar feeding strategies, which remain surprisingly poorly understood. More specifically, it has long been thought that the more derived Dipsadini genera Dipsas and Sibynomorphus use glandular toxins secreted by their infralabial glands to extract snails from their shells. RESULTS: Here, we report the presence in the tribe Dipsadini of a novel lower jaw protein-secreting delivery system effected by a gland that is not functionally related to adjacent teeth, but rather opens loosely on the oral epithelium near the tip of the mandible, suggesting that its secretion is not injected into the prey as a form of envenomation but rather helps control the mucus and assists in the ingestion of their highly viscous preys. A similar protein-secreting system is also present in the goo-eating genus Geophis and may share the same adaptive purpose as that hypothesized for Dipsadini. Our phylogenetic hypothesis suggests that the acquisition of a seromucous infralabial gland represents a uniquely derived trait of the goo-eating clade that evolved independently twice within the group as a functionally complex protein-secreting delivery system. CONCLUSIONS: The acquisition by snail-eating snakes of such a complex protein-secreting system suggests that the secretion from the hypertrophied infralabial glands of goo-eating snakes may have a fundamental role in mucus control and prey transport rather than envenomation of prey. Evolution of a functional secretory system that combines a solution for mucus control and transport of viscous preys is here thought to underlie the successful radiation of goo-eating snakes.

Highly conserved and extremely variable: The paradoxical pattern of toxin expression revealed by comparative venom-gland transcriptomics of Phalotris (Serpentes: Dipsadidae).

Although non-front fanged snakes account for almost two-thirds of snake diversity, most studies on venom composition and evolution focus exclusively on front-fanged species, which comprise most of the clinically relevant accidents. Comprehensive reports on venom composition of non-front fanged snakes are still scarce for several groups. In this study, we address such shortage of knowledge by providing new insights about the venom composition among species of Phalotris, a poorly studied Neotropical dipsadid genus. Phalotris are known for their specialized venom delivery system and toxic venoms, which can cause life-threatening accidents in humans. We evaluate the venom-gland transcriptome of Phalotris, comparing the following three South American species: P. reticulatus for the Araucaria Pine forests, P. lemniscatus for the Pampa grasslands, and P. mertensi for the Brazilian Cerrado. Our results indicate similar venom profiles, in which they share a high expression level of Kunitz-type inhibitors (KUNZ). On the other hand, comparative analyses revealed substantial differences in the expression levels of C-type lectins (CTL) and snake venom metalloproteinases (SVMP). The diverse set of SVMP and CTL isoforms shows signals of positive selection, and we also identified truncated forms of type III SVMPs, which resemble type II and type I SVMPs of viperids. Additionally, we identified a CNP precursor hosting a proline-rich region containing a BPP motif resembling those commonly detected in viperid venoms with hypotensive activity. Altogether, our results suggest an evolutionary history favoring high expression levels of few KUNZ isoforms in Phalotris venoms, contrasting with a highly diverse set of SVMP and CTL isoforms. Such diversity can be comparable with the venom variability observed in some viperids. Our findings highlight the extreme phenotypic diversity of non-front fanged snakes and the importance to allocate greater effort to study neglected groups of Colubroidea.

Venom Composition of Neglected Bothropoid Snakes from the Amazon Rainforest: Ecological and Toxinological Implications.

Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora, Bothrops taeniatus, Bothrops bilineatus smaragdinus, Bothrops brazili, and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.

Molecular systematics and historical biogeography of tree boas (Corallus spp.).

Inferring the evolutionary and biogeographic history of taxa occurring in a particular region is one way to determine the processes by which the biodiversity of that region originated. Tree boas of the genus Corallus are an ancient clade and occur throughout Central and South America and the Lesser Antilles, making it an excellent group for investigating Neotropical biogeography. Using sequenced portions of two mitochondrial and three nuclear loci for individuals of all recognized species of Corallus, we infer phylogenetic relationships, present the first molecular analysis of the phylogenetic placement of the enigmatic C. cropanii, develop a time-calibrated phylogeny, and explore the biogeographic history of the genus. We found that Corallus diversified within mainland South America, via over-water dispersals to the Lesser Antilles and Central America, and via the traditionally recognized Panamanian land bridge. Divergence time estimates reject the South American Caribbean-Track as a general biogeographic model for Corallus and implicate a role for events during the Oligocene and Miocene in diversification such as marine incursions and the uplift of the Andes. Our findings also suggest that recognition of the island endemic species, C. grenadensis and C. cookii, is questionable as they are nested within the widely distributed species, C. hortulanus. Our results highlight the importance of using widespread taxa when forming and testing biogeographic hypotheses in complex regions and further illustrate the difficulty of forming broadly applicable hypotheses regarding patterns of diversification in the Neotropical region.

A new species of lizard genus Potamites from Ecuador (Squamata, Gymnophthalmidae).

Potamitesflavogularis sp. nov. is described from the Napo and Tungurahua Provinces around 1800 m elevation in eastern Ecuador. The new species is closely related, sibling, and sympatric to Potamites cochranae to which it has been previously confused. It is characterized by the absence of isolated basal flounces of spines and presence of calcareous spinules on flounces of the hemipenis, a short (1,30-1,41 times SVL) and slightly compressed tail without tubercles, tympanum slightly recessed, subimbricate ventral scales, lateral body scales lacking conspicuous enlarged tubercles, four longitudinal rows of dorsal tubercles, 6 transverse series of ventral scales, absence of intercalated scales along sides of tail, and absence of tubercles on sides of neck and gular regions. Like their congeners, the new species was found close to vegetation surrounding streams in primary and secondary forests.

Patterns of morphological variation and ecological correlates in the skull of vipers (Serpentes: Viperidae).

The skull of vipers is a highly kinetic anatomical structure involved in envenomating and consuming of prey. Morphological knowledge about the viperid skull is based on studies on some groups of species, but information on its variation within the whole family and its functional morphology is still scarce. In this study, we aimed to explore variation in skull morphology among species of the three subfamilies of Viperidae, and test whether that variation correlates with macrohabitat and diet. We performed quantitative analyses of the viperid skull based on broad taxonomic sampling and two methodological approaches: linear and geometric morphometrics. The results of both approaches showed that much of the variation lies in differences of shape and relative size of the premaxilla, the nasals, the frontals, and the parietals. The results indicated that phylogeny and size influence the shape of the skull, but we also found evidence of morphological differentiation between arboreal and terrestrial species and in species with mammal specialist diet. Our findings imply that, besides evolutionary allometry and phylogenetic signal, demands of particular diets coupled with use of certain habitats have in part shaped morphological evolution of the viperid skull.

Comparing morphological and secretory aspects of cephalic glands among the New World coral snakes brings novel insights on their biological roles.

Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers.

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.