Fang evolution in venomous snakes: Adaptation of 3D tooth shape to the biomechanical properties of their prey

How venom pore placement may influence puncture performance in snake fangs

When designing experimental studies, it is important to understand the biological context of the question being asked. For example, many biological puncture experiments embed the puncture tool to a standardized depth based on a percentage of the total tool length, to compare the performance between tools. However, this may not always be biologically relevant to the question being asked. To understand how definitions of penetration depth may influence comparative results, we performed puncture experiments on a series of venomous snake fangs using the venom pore location as a functionally relevant depth standard. After exploring variation in pore placement across snake phylogeny, we compared the work expended during puncture experiments across a set of snake fangs using various depth standards: puncture initiation, penetration to a series of depths defined by the venom pore and penetration to 15% of fang length. Contrary to our hypothesis, we found almost no pattern in pore placement between clades, dietary groups or venom toxicity. Rank correlation statistics of our experimental energetics results showed no difference in the broad comparison of fangs when different puncture depth standards were used. However, pairwise comparisons between fangs showed major shifts in significance patterns between the different depth standards used. These results imply that the interpretation of experimental puncture data will heavily depend upon which depth standard is used during the experiments. Our results illustrate the importance of understanding the biological context of the question being addressed when designing comparative experiments.

Bayesian Prediction of Multivariate Ecology from Phenotypic Data Yields New Insights into the Diets of Extant and Extinct Taxa

AbstractMorphology often reflects ecology, enabling the prediction of ecological roles for taxa that lack direct observations, such as fossils. In comparative analyses, ecological traits, like diet, are often treated as categorical, which may aid prediction and simplify analyses but ignores the multivariate nature of ecological niches. Furthermore, methods for quantifying and predicting multivariate ecology remain rare. Here, we ranked the relative importance of 13 food items for a sample of 88 extant carnivoran mammals and then used Bayesian multilevel modeling to assess whether those rankings could be predicted from dental morphology and body size. Traditional diet categories fail to capture the true multivariate nature of carnivoran diets, but Bayesian regression models derived from living taxa have good predictive accuracy for importance ranks. Using our models to predict the importance of individual food items, the multivariate dietary niche, and the nearest extant analogs for a set of data-deficient extant and extinct carnivoran species confirms long-standing ideas for some taxa but yields new insights into the fundamental dietary niches of others. Our approach provides a promising alternative to traditional dietary classifications. Importantly, this approach need not be limited to diet but serves as a general framework for predicting multivariate ecology from phenotypic traits.

Armed to the teeth: The underestimated diversity in tooth shape in snakes and its relation to feeding behavior and diet

The structure, composition, and shape of teeth have been related to dietary specialization in many vertebrate species, but comparative studies on snakes' teeth are lacking. Yet, snakes have diverse dietary habits that may impact the shape of their teeth. We hypothesize that prey properties, such as hardness and shape, as well as feeding behavior, such as aquatic or arboreal predation, or holding vigorous prey, impose constraints on the evolution of tooth shape in snakes. We compared the morphology of the dentary teeth of 63 species that cover the phylogenetic and dietary diversity of snakes, using 3D geometric morphometrics and linear measurements. Our results show that prey hardness, foraging substrate, and the main feeding mechanical challenge are important drivers of tooth shape, size, and curvature. Overall, long, slender, curved teeth with a thin layer of hard tissue are observed in species that need to maintain a grip on their prey. Short, stout, less curved teeth are associated with species that undergo high or repeated loads. Our study demonstrates the diversity of tooth morphology in snakes and the need to investigate its underlying functional implications to better understand the evolution of teeth in vertebrates.

The effect of structural curvature on the load-bearing characteristics of biomechanical elements

Miniature, sharped-edge, curved-shape biomechanical elements appear in various biological systems and grant them diverse functional capabilities, such as mechanical defense, venom injection, and frictional support. While these biomechanical elements demonstrate diverse curved shapes that span from slightly curved needle-like elements (e.g., stingers), through moderately curved anchor-like elements (e.g., claws), to highly curved hook-like elements (e.g., fangs)-the curvature effect on the load-bearing capabilities of these biomechanical elements are yet mostly unknown. Here, we employ structural-mechanical modeling to explore the relationships between the curved shapes of biomechanical elements on their local deformation mechanisms, overall elastic stiffness, and reaction forces on a target surface. We found that the curvature of the biomechanical element is a prime modulator of its load-bearing characteristics that substantially affect its functional capabilities. Slightly curved elements are preferable for penetration states with optimal load-bearing capabilities parallel to their tips but possess high directional sensitivity and degraded capabilities for scratching states; contrary, highly curved elements are suitable for combined penetration-scratching states with mild directional sensitivity and optimal load-bearing capabilities in specialized angular orientation to their tips. These structural-mechanical principles are tightly linked to the intrinsic functional roles of biomechanical elements in diverse natural systems, and their synthetic realizations may promote new engineering designs of advanced biomedical injections, functional surfaces, and micromechanical devices.

Morphological quantification of the maxillary canine tooth in the domestic dog (Canis lupus familiaris)

Canine tooth shape is known to vary with diet and killing behavior in wild animals and the relationship between form and function is driven in part by selective pressure. However, comparative investigation of the domestic dog (Canis lupus familiaris) is of interest. How do they compare to their wild counterparts? This study sought to quantify and characterize the morphology of the canine tooth in the domestic dog, and to provide a preliminary investigation into the variance in canine tooth morphology across individual dogs of varying breeds. Three-dimensional (3D) models generated from micro-computed tomography (µ-CT) studies of 10 mature maxillary canine teeth from the domesticated dog (Canis lupus familiaris) were used to quantify key morphological features and evaluate variance among dogs. Results show that, utilizing modern imaging and model building software, the morphology of the canine tooth can be comprehensively characterized and quantified. Morphological variables such as second moment of area and section modulus (geometrical parameters related to resistance to bending), as well as aspect ratio, ridge sharpness, cusp sharpness and enamel thickness are optimized in biomechanically critical areas of the tooth crown to balance form and function. Tooth diameter, second moment of area, section modulus, cross sectional area, tooth volume and length as well as enamel thickness are highly correlated with body weight. In addition, we found preliminary evidence of morphological variance across individual dogs. Quantification of these features provide insight into the balance of form and function of the canine tooth in wild and domesticated canids. In addition, results suggest that variance between dogs exist in some morphological features and most morphological features are highly correlated with body weight.

Shedding Light on the Dentition and Venom Delivery System of the Rear-Fanged Snake, Galvarinus chilensis chilensis (Serpentes: Dipsadidae: Tachymenini) from Chile

Although the rear-fanged snake Galvarinus chilensis chilensis (formerly named Tachymenis ch. chilensis) causes ophidian accidents with clinical importance in Chile, the anatomical and histological characterizations of the venom delivery system (venom gland and fang) of this species still remain unknown. This study describes the dentition and characteristics of fangs and their ontogenetic variations in G. ch. chilensis. Moreover, histological and histochemistry analyses of the venom glands of this species are presented. Using micro-computed tomography and scanning electron microscopy, the dentitions of neonates, juveniles, and adults were analyzed, and no ontogenetic variations in teeth length and number present in the dentary and maxilla were observed. Moreover, we found three types of basic teeth, with distributional patterns conserved in all ontogenetic categories. The fangs exhibited a groove from the base to the middle. At the end of the groove, prominent ridges are formed. The fang and groove lengths were significantly distinct between ontogenetic categories. No differences between females and males were observed. Histologically, we found that the venom gland is close to the fangs and has a seromucous composition. Our results describe, for the first time, the distributional pattern and characteristics of the dentition and venom delivery system of the poorly studied snake G. ch. chilensis.

Taking a stab at modelling canine tooth biomechanics in mammalian carnivores with beam theory and finite-element analysis

Canine teeth are vital to carnivore feeding ecology, facilitating behaviours related to prey capture and consumption. Forms vary with specific feeding ecologies; however, the biomechanics that drive these relationships have not been comprehensively investigated. Using a combination of beam theory analysis (BTA) and finite-element analysis (FEA) we assessed how aspects of canine shape impact tooth stress, relating this to feeding ecology. The degree of tooth lateral compression influenced tolerance of multidirectional loads, whereby canines with more circular cross-sections experienced similar maximum stresses under pulling and shaking loads, while more ellipsoid canines experienced higher stresses under shaking loads. Robusticity impacted a tooth's ability to tolerate stress and appears to be related to prey materials. Robust canines experience lower stresses and are found in carnivores regularly encountering hard foods. Slender canines experience higher stresses and are associated with carnivores biting into muscle and flesh. Curvature did not correlate with tooth stress; however, it did impact bending during biting. Our simulations help identify scenarios where canine forms are likely to break and pinpoint areas where this breakage may occur. These patterns demonstrate how canine shape relates to tolerating the stresses experienced when killing and feeding, revealing some of the form-function relationships that underpin mammalian carnivore ecologies.

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.

Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins

Background

The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes.

Results

Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae – Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae – Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms.

Conclusions

We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator–prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01208-9.

Fang shape varies with ontogeny and sex in the venomous elapid snake Pseudonaja affinis

A predator's preferred prey often changes over the course of its life as it grows from an inexperienced juvenile through to a sexually mature adult. For species with highly specialised feeding strategies, this may require its anatomy to change over the course of its life. The dugite (Pseudonaja affinis, Günther 1872) is a venomous snake from Australia that displays such a diet shift, with juveniles feeding on small reptiles, while adults mainly target mammals. We examined the morphology of fangs across both sexes and throughout ontogeny using geometric morphometrics and cross-sectional sharpness measurements of key functional regions on these teeth. This highlighted key differences in shape that likely relate to the varied properties of their adult and juvenile diet. We found that juveniles display a more robust and blunter fang, which likely relates to feeding on scaly lizard prey, whereas adults have slender fangs with sharper tips, which reflects their diet of softer mammalian prey. There were also differences between males and females, with male snakes having significantly more slender fangs than females, which might be an indication of niche partitioning between the sexes. Using snout-vent length as a proxy for age, we found that the ontogenetic shift in fang shape occurs when P. affinis is around 60 cm long, corresponding with previous studies that found this size to be the moment where these snakes switch from their juvenile to adult diet.

Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake (Crotalus horridus)

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)–hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.

The killer’s toolkit: remarkable adaptations in the canine teeth of mammalian carnivores

Often the first point of contact between predator and prey, mammalian canine teeth are essential for killing, dismembering and consuming prey. Yet despite their importance, few associations among shape, function and phylogeny are established. We undertook the first comprehensive analysis of canine tooth shape across predatory mammals (Carnivora, Didelphimorphia and Dasyuromorphia), integrating shape analysis with function of this fundamental feature. Shape was quantified using three-dimensional geometric morphometrics and cross-sectional sharpness. Canines vary in three main ways (sharpness, robustness and curvature) which vary with diet, killing behaviour and phylogeny. Slender, sharp canines are associated with carnivores such as felids that target the neck of their prey and primarily consume the ‘softer’ parts of a carcass. Robust, blunt canines are found in mustelids and dasyurids that typically consume ‘harder’ materials, such as bone, or bite into skulls. Differences in the killing behaviours of felids and canids probably result in more curved canines in the latter, which act as hooks to hold prey. We find functional specialization in the upper and lower canines of individuals and across the major mammalian clades. These patterns demonstrate how canine teeth are adapted to suit diverse diets and hunting styles, enabling mammals to become some of nature's most successful predators.