EFFECT OF BODY POSITION ON ECHOCARDIOGRAPHIC PARAMETERS IN PRAIRIE RATTLESNAKES (CROTALUS VIRIDIS).

Echocardiography is a noninvasive diagnostic tool that can provide instantaneous information about cardiac function, but it is uncommonly used by veterinarians to assess reptilian patients. Echocardiograms were performed on 14 clinically healthy, adult prairie rattlesnakes (Crotalus viridis), and cardiac measurements were taken in a horizontal and vertical position. Cardiac parameters including ventricular volume in systole and diastole, as well as the diameter of the left atrium, pulmonary artery, and paired left and right aortic arches were obtained. No evidence of cardiac disease was noted in any of the study animals. Males had a greater percentage of ventricular volume change (VVC) than females in the vertical position (P = 0.043). The percentage of ventricular volume change was significantly lower in the horizontal compared with the vertical position (P = 0.032) and was not different by sex. For the short-axis views, the right atrial diameter and short-axis ventricular area in diastole and systole were significantly smaller in snakes in the vertical compared with the horizontal position. This study is the first to obtain echocardiographic measurements in North American vipers and adds to the understanding of techniques used to evaluate the cardiac function of these species.

Intraspecific sequence and gene expression variation contribute little to venom diversity in sidewinder rattlesnakes ( Crotalus cerastes).

Traits can evolve rapidly through changes in gene expression or protein-coding sequences. However, these forms of genetic variation can be correlated and changes to one can influence the other. As a result, we might expect traits lacking differential expression to preferentially evolve through changes in protein sequences or morphological adaptation. Given the lack of differential expression across the distribution of sidewinder rattlesnakes ( Crotalus cerastes), we tested this hypothesis by comparing the coding regions of genes expressed in the venom gland transcriptomes and fang morphology. We calculated Tajima's D and FST across four populations comparing toxin and nontoxin loci. Overall, we found little evidence of directional selection or differentiation between populations, suggesting that changes to protein sequences do not underlie the evolution of sidewinder venom or that toxins are under extremely variant selection pressures. Although low-expression toxins do not have higher sequence divergence between populations, they do have more standing variation on which selection can act. Additionally, we found significant differences in fang length among populations. The lack of differential expression and sequence divergence suggests sidewinders-given their generalist diet, moderate gene flow and environmental variation-are under stabilizing selection which functions to maintain a generalist phenotype. Overall, we demonstrate the importance of examining the relationship between gene expression and protein-coding changes to understand the evolution of complex traits.

Interfang Distances of Rattlesnakes: Sexual, Interspecific, and Body Size-related Variation, and Implications for Snakebite Research and Management.

OBJECTIVES: Snakebite severity corresponds to size of snake because the amount of venom a snake injects is positively associated with snake size. Because fang marks are often present on snakebite patients, we tested whether the relationship between snake length and distance between fang puncture wounds can be generalized for rattlesnakes of genus Crotalus. METHODS: We measured 2 interfang distances from 79 rattlesnakes of both sexes, 5 species, and varying body length: 1) distance between fang bases in anesthetized snakes, and 2) distance between fang punctures in a membrane-covered beaker bitten defensively. RESULTS: Statistical analyses supported our 2 hypotheses, that 1) body size-related fang divergence during fang protraction (ie, anterolateral movement during fang erection), and 2) the relationship between snake length and interfang distance are similar between the sexes and among different rattlesnake species. We therefore derived a general equation to estimate snake length based on distance between fang marks, and recommended 5 snake size categories: very small (<10 mm), small (10-15 mm), medium (15-20 mm), large (20-25 mm), and very large (>25 mm). CONCLUSIONS: The distance between fang marks on a snakebite patient may be used to estimate the size or size category of the offending snake, which in some cases may have predictive value for overall clinical severity of a given envenomation. Assessing interfang distance from puncture wounds can improve snakebite research and anticipation of snakebite severity.

Molecular basis of infrared detection by snakes.

Snakes possess a unique sensory system for detecting infrared radiation, enabling them to generate a 'thermal image' of predators or prey. Infrared signals are initially received by the pit organ, a highly specialized facial structure that is innervated by nerve fibres of the somatosensory system. How this organ detects and transduces infrared signals into nerve impulses is not known. Here we use an unbiased transcriptional profiling approach to identify TRPA1 channels as infrared receptors on sensory nerve fibres that innervate the pit organ. TRPA1 orthologues from pit-bearing snakes (vipers, pythons and boas) are the most heat-sensitive vertebrate ion channels thus far identified, consistent with their role as primary transducers of infrared stimuli. Thus, snakes detect infrared signals through a mechanism involving radiant heating of the pit organ, rather than photochemical transduction. These findings illustrate the broad evolutionary tuning of transient receptor potential (TRP) channels as thermosensors in the vertebrate nervous system.

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus).

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age-specific selection is constrained to act on the same genome by altering age-specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype-phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.

Evaporative respiratory cooling augments pit organ thermal detection in rattlesnakes.

Rattlesnakes use their facial pit organs to sense external thermal fluctuations. A temperature decrease in the heat-sensing membrane of the pit organ has the potential to enhance heat flux between their endothermic prey and the thermal sensors, affect the optimal functioning of thermal sensors in the pit membrane and reduce the formation of thermal "afterimages", improving thermal detection. We examined the potential for respiratory cooling to improve strike behaviour, capture, and consumption of endothermic prey in the South American rattlesnake, as behavioural indicators of thermal detection. Snakes with a higher degree of rostral cooling were more accurate during the strike, attacking warmer regions of their prey, and relocated and consumed their prey faster. These findings reveal that by cooling their pit organs, rattlesnakes increase their ability to detect endothermic prey; disabling the pit organs caused these differences to disappear. Rattlesnakes also modify the degree of rostral cooling by altering their breathing pattern in response to biologically relevant stimuli, such as a mouse odour. Our findings reveal that low humidity increases their ability to detect endothermic prey, suggesting that habitat and ambush site selection in the wild may be influenced by external humidity levels as well as temperature.

Directional sensitivity in the thermal response of the facial pit in western diamondback rattlesnakes (Crotalus atrox).

Recent work published in the accompanying paper used a combination of 3D morphological reconstruction to define optical spread functions and heat transfer physics to study how external heat energy would reach the sensory membrane within the facial pit of pitvipers. The results from all of the species examined indicated asymmetric directional sensitivity, e.g. the pit would preferentially respond to stimuli located below and behind the snake. The present study was intended as a test of these findings through a quantitative neurophysiological analysis of directional sensitivity in the facial pit of the western diamondback rattlesnake, Crotalus atrox. An infrared emitter was positioned through a coordinate system (with varying angular orientations and distances) and the response it evoked measured through neurophysiological recordings of a trigeminal nerve branch composed of the afferents from the sensory membrane of the facial pit. Significant differences were found in the strength of the membrane's neural response to a constant stimulus presented at different orientations (relative to the facial pit opening) and over different distances. The peak sensitivity (at 12 deg above and 20 deg in front of the facial pit opening) was in good agreement with the predicted directional sensitivities based on optical spread functions and 3D topography. These findings support the hypothesis that the topography, and functional performance, of the facial pit has undergone an adaptive radiation within the pit vipers, and that differences in the behavioral ecology of the pit vipers (i.e. terrestrial versus arboreal) are reflected within the facial pits.

Experimentally altered navigational demands induce changes in the cortical forebrain of free-ranging northern pacific rattlesnakes (Crotalus o. oreganus).

The hippocampus of birds and mammals plays a crucial role in spatial memory and navigation. The hippocampus exhibits plasticity in adulthood in response to diverse environmental factors associated with spatial demands placed on an animal. The medial and dorsal cortices of the telencephalon of squamate reptiles have been implicated as functional homologues to the hippocampus. This study sought to experimentally manipulate the navigational demands placed on free-ranging northern Pacific rattlesnakes (Crotalus o. oreganus) to provide direct evidence of the relationship between spatial demands and neuroplasticity in the cortical telencephalon of the squamate brain. Adult male rattlesnakes were radio-tracked for 2 months, during which time 1 of 3 treatments was imposed weekly, namely 225-meter translocation in a random direction, 225-meter walk and release at that day's capture site (handling control) or undisturbed (control). Snakes were then sacrificed and the brains were removed and processed for histological analysis of cortical features. The activity range was larger in the translocated (Tr) group compared to the handled (Hd) and undisturbed control (Cn) groups when measured via 95% minimum convex polygon (MCP). At the 100% MCP level, Tr snakes had larger activity ranges than the Cn snakes only. The volume of the medial cortex (MC) was larger in the Tr group compared to the Cn group. The MC of Hd snakes was not significantly different from that of either of the other groups. No differences in dorsal cortex (DC) or lateral cortex volumes were detected among the groups. Numbers of 5-bromo-2'-deoxyuridine (BrdU)-labeled cells in the MC and DC 3 weeks after BrdU injection were not affected by treatment. This study establishes a causal relationship between navigational demands and greater MC volume in a free-ranging reptile.

Crotalus Durissus Ruruima: Current Knowledge on Natural History, Medical Importance, and Clinical Toxinology.

Crotalus durissus ruruima is a rattlesnake subspecies mainly found in Roraima, the northernmost state of Brazil. Envenomings caused by this subspecies lead to severe clinical manifestations (e.g. respiratory muscle paralysis, rhabdomyolysis, and acute renal failure) that can lead to the victim's death. In this review, we comprehensively describe C. d. ruruima biology and the challenges this subspecies poses for human health, including morphology, distribution, epidemiology, venom cocktail, clinical envenoming, and the current and future specific treatment of envenomings by this snake. Moreover, this review presents maps of the distribution of the snake subspecies and evidence that this species is responsible for some of the most severe envenomings in the country and causes the highest lethality rates. Finally, we also discuss the efficacy of the Brazilian horse-derived antivenoms to treat C. d. ruruima envenomings in Roraima state.

Large snake size suggests increased snakebite severity in patients bitten by rattlesnakes in Southern california.

OBJECTIVE: To correlate rattlesnake size and other characteristics of envenomation with the severity of envenomation. METHODS: We retrospectively reviewed 145 charts of patients bitten by rattlesnakes in Southern California between 1995 and 2004, measuring Snakebite Severity Scores (SSS) and characteristics of envenomation that might be correlated with snakebite severity, including rattlesnake size, rattlesnake species, patient size, and anatomic location of the bite. Outcomes measured included SSS, complications of envenomation, number of vials of antivenom used, and length of hospital stay. RESULTS: Of the patients bitten by rattlesnakes, 81% were men, and 79% of bites were on the upper extremities. Fifty-five percent of bites were provoked by the patient, and 44% were unprovoked. Neither location of snakebite nor provocation of snakebite affected the SSS. Only 1 patient had a snakebite without envenomation, and only 1 patient died from envenomation. Rattlesnake size was positively correlated with SSS, and SSS was positively correlated with the number of vials of antivenom used and with the length of hospital stay. Rattlesnake species and patient mass did not affect SSS. CONCLUSIONS: Larger rattlesnakes cause more severe envenomations, which contradicts popular belief.

Side-impact collision: mechanics of obstacle negotiation in sidewinding snakes.

Snakes excel at moving through cluttered environments, and heterogeneities can be used as propulsive contacts for snakes performing lateral undulation. However, sidewinding, which is often associated with sandy deserts, cuts a broad path through its environment that may increase its vulnerability to obstacles. Our prior work demonstrated that sidewinding can be represented as a pair of orthogonal body waves (vertical and horizontal) that can be independently modulated to achieve high maneuverability and incline ascent, suggesting that sidewinders may also use template modulations to negotiate obstacles. To test this hypothesis, we recorded overhead video of four sidewinder rattlesnakes (Crotalus cerastes) crossing a line of vertical pegs placed in the substrate. Snakes used three methods to traverse the obstacles: a Propagate Through behavior in which the lifted moving portion of the snake was deformed around the peg and dragged through as the snake continued sidewinding (115/160 runs), Reversal turns that reorient the snake entirely (35/160), or switching to Concertina locomotion (10/160). The Propagate Through response was only used if the anterior-most region of static contact would propagate along a path anterior to the peg, or if a new region of static contact could be formed near the head to satisfy this condition; otherwise, snakes could only use Reversal turns or switch to Concertina locomotion. Reversal turns allowed the snake to re-orient and either escape without further peg contact or re-orient into a posture amenable to using the Propagate Through response. We developed an algorithm to reproduce the Propagate Through behavior in a robophysical model using a modulation of the two-wave template. This range of behavioral strategies provides sidewinders with a versatile range of options for effectively negotiating obstacles in their natural habitat, as well as provide insights into the design and control of robotic systems dealing with heterogeneous habitats.

Evolutionary morphology of the rattlesnake style.

BACKGROUND: The rattlesnake rattling system is an evolutionary novelty that includes anatomical, behavioral, and physiological modifications of the generalized pitviper tail. One such modification, the formation of a bony clublike style at the terminal region of the caudal vertebrae, has not previously been examined in a phylogenetic context. Here we used skeletal material, cleared and stained preparations, and radiographs of whole preserved specimens to examine interspecific variation in style morphology among 34 rattlesnake species. RESULTS: Evolutionary Principal Components Analysis revealed an inverse relationship between caudal segmental counts and style size, supporting the hypothesis that bone from caudal vertebral elements was reallocated to style formation during the evolution of this structure. Most of the basal rattlesnake species have small styles consisting of few compacted vertebral elements; however, early in the rattlesnake radiation there appears to have been two independent transitions to relatively large, pronged styles consisting of multiple coalesced vertebrae (once in Sistrurus catenatus, and once in Crotalus following the divergence of the Mexican long-tailed rattlesnakes). In terms of style shape, the two most divergent species, C. catalinensis and C. ericsmithi, provide insight into the possible relationship between style and rattle matrix morphology and lineage-specific evolutionary strategies for retaining rattle segments. CONCLUSION: The considerable interspecific variation in rattle morphology appears to correspond to variation in the bony style. We hypothesize that style morphology evolves indirectly as an integrated module responding to adaptive evolution on matrix morphology.

Evidence for divergent patterns of local selection driving venom variation in Mojave Rattlesnakes (Crotalus scutulatus).

Snake venoms represent an enriched system for investigating the evolutionary processes that lead to complex and dynamic trophic adaptations. It has long been hypothesized that natural selection may drive geographic variation in venom composition, yet previous studies have lacked the population genetic context to examine these patterns. We leverage range-wide sampling of Mojave Rattlesnakes (Crotalus scutulatus) and use a combination of venom, morphological, phylogenetic, population genetic, and environmental data to characterize the striking dichotomy of neurotoxic (Type A) and hemorrhagic (Type B) venoms throughout the range of this species. We find that three of the four previously identified major lineages within C. scutulatus possess a combination of Type A, Type B, and a 'mixed' Type A + B venom phenotypes, and that fixation of the two main venom phenotypes occurs on a more fine geographic scale than previously appreciated. We also find that Type A + B individuals occur in regions of inferred introgression, and that this mixed phenotype is comparatively rare. Our results support strong directional local selection leading to fixation of alternative venom phenotypes on a fine geographic scale, and are inconsistent with balancing selection to maintain both phenotypes within a single population. Our comparisons to biotic and abiotic factors further indicate that venom phenotype correlates with fang morphology and climatic variables. We hypothesize that links to fang morphology may be indicative of co-evolution of venom and other trophic adaptations, and that climatic variables may be linked to prey distributions and/or physiology, which in turn impose selection pressures on snake venoms.

The venom gland transcriptome of the Desert Massasauga rattlesnake (Sistrurus catenatus edwardsii): towards an understanding of venom composition among advanced snakes (Superfamily Colubroidea).

BACKGROUND: Snake venoms are complex mixtures of pharmacologically active proteins and peptides which belong to a small number of superfamilies. Global cataloguing of the venom transcriptome facilitates the identification of new families of toxins as well as helps in understanding the evolution of venom proteomes. RESULTS: We have constructed a cDNA library of the venom gland of a threatened rattlesnake (a pitviper), Sistrurus catenatus edwardsii (Desert Massasauga), and sequenced 576 ESTs. Our results demonstrate a high abundance of serine proteinase and metalloproteinase transcripts, indicating that the disruption of hemostasis is a principle mechanism of action of the venom. In addition to the transcripts encoding common venom proteins, we detected two varieties of low abundance unique transcripts in the library; these encode for three-finger toxins and a novel toxin possibly generated from the fusion of two genes. We also observed polyadenylated ribosomal RNAs in the venom gland library, an interesting preliminary obsevation of this unusual phenomenon in a reptilian system. CONCLUSION: The three-finger toxins are characteristic of most elapid venoms but are rare in viperid venoms. We detected several ESTs encoding this group of toxins in this study. We also observed the presence of a transcript encoding a fused protein of two well-characterized toxins (Kunitz/BPTI and Waprins), and this is the first report of this kind of fusion in a snake toxin transcriptome. We propose that these new venom proteins may have ancillary functions for envenomation. The presence of a fused toxin indicates that in addition to gene duplication and accelerated evolution, exon shuffling or transcriptional splicing may also contribute to generating the diversity of toxins and toxin isoforms observed among snake venoms. The detection of low abundance toxins, as observed in this and other studies, indicates a greater compositional similarity of venoms (though potency will differ) among advanced snakes than has been previously recognized.

Soft epidermis of a scaleless snake lacks beta-keratin.

Beta-keratins are responsible for the mechanical resistance of scales in reptiles. In a scaleless crotalus snake (Crotalus atrox), large areas of the skin are completely devoid of scales, and the skin appears delicate and wrinkled. The epidermis of this snake has been assessed for the presence of beta-keratin by immunocytochemistry and immunoblotting using an antibody against chicken scale beta-keratin. This antibody recognizes beta-keratins in normal snake scales with molecular weights of 15-18 kDa and isoelectric points at 6.8, 7.5, 8.3 and 9.4. This indicates that beta-keratins of the stratum corneum are mainly basic proteins, so may interact with cytokeratins of the epidermis, most of which appear acidic (isoelectric points 4.5-5.5). A beta-layer and beta-keratin immunoreactivity are completely absent in moults of the scaleless mutant, and the corneous layer comprises a multi-layered alpha-layer covered by a flat oberhautchen. In conclusion, the present study shows that a lack of beta-keratins is correlated with the loss of scales and mechanical protection in the skin of this mutant snake.

Deconstructing a Species-Complex: Geometric Morphometric and Molecular Analyses Define Species in the Western Rattlesnake (Crotalus viridis).

Morphological data are a conduit for the recognition and description of species, and their acquisition has recently been broadened by geometric morphometric (GM) approaches that co-join the collection of digital data with exploratory 'big data' analytics. We employed this approach to dissect the Western Rattlesnake (Crotalus viridis) species-complex in North America, currently partitioned by mitochondrial (mt)DNA analyses into eastern and western lineages (two and seven subspecies, respectively). The GM data (i.e., 33 dorsal and 50 lateral head landmarks) were gleaned from 2,824 individuals located in 10 museum collections. We also downloaded and concatenated sequences for six mtDNA genes from the NCBI GenBank database. GM analyses revealed significant head shape differences attributable to size and subspecies-designation (but not their interactions). Pairwise shape distances among subspecies were significantly greater than those derived from ancestral character states via squared-change parsimony, with the greatest differences separating those most closely related. This, in turn, suggests the potential for historic character displacement as a diversifying force in the complex. All subspecies, save one, were significantly differentiated in a Bayesian discriminant function analysis (DFA), regardless of whether our priors were uniform or informative (i.e., mtDNA data). Finally, shape differences among sister-clades were significantly greater than expected by chance alone under a Brownian model of evolution, promoting the hypothesis that selection rather than drift was the driving force in the evolution of the complex. Lastly, we combine head shape and mtDNA data so as to derived an integrative taxonomy that produced robust boundaries for six OTUs (operational taxonomic units) of the C. viridis complex. We suggest these boundaries are concomitant with species-status and subsequently provide a relevant nomenclature for its recognition and representation.

Body size variation among mainland populations of the western rattlesnake (Crotalus viridis).

In general, squamate reptiles follow the converse to Bergmann's rule, attaining smaller sizes in cooler environments, whereas other vertebrate groups follow Bergmann's rule, attaining larger sizes in cooler areas. Intensive studies of body size evolution for species of squamates are necessary to understand the processes responsible for this trend. Here I present data on body size variation among mainland populations of the western rattlesnake, Crotalus viridis. This species consists of two well-differentiated phylogenetic clades, therefore all analyses were performed for the C. viridis group as a whole and separately for each of the two clades within the C. viridis group. Although both phylogenetic and nonphylogenetic analyses were performed, the data did not show phylogenetic conservatism, and therefore the nonphylogenetic results are preferred. I found no significant relationships between mean adult female snout-vent length and any of the physical and climatic variables that were examined for the C. viridis group using simple linear regression analysis. Examined separately, I found that individuals of the western clade, C. oreganus, were smaller in cooler and more seasonal environments, whereas individuals of the eastern clade. C. viridis sensu stricto, were larger in cooler and more seasonal areas. Thus, the observed size trends were in opposite directions for the two clades. Multiple regression analysis revealed that seasonality was a stronger predictor of body size variation than was temperature for both clades. The differences in body size trends between these clades may be due to differences in mortality rates among populations.

Fang tip spread, puncture distance, and suction for snake bite.

We measured the distance between fang tip punctures in defensive bites by western diamondback rattlesnakes (Crotalus atrox) and the distance between their retracted fangs. Because the fang tips at penetration average 112% further apart than their bases at rest, The Extractor, a device widely marketed in the United States for snake bite first aid, will not simultaneously cover both punctures of most adult New World pitvipers.

Primary culture of venom gland cells from the South American rattlesnake (Crotalus durissus terrificus).

Primary cultures of venom gland cells from the South American rattlesnake (Crotalus durissus terrificus) were attempted. At first, six different cell types were obtained including potentially secreting epithelial-like cells. Nonepithelial cell cultures were later invaded by fibroblast-like cells. Cultures of epithelial-like gland cells were successfully maintained, after testing different culture conditions by varying the media, incubation temperature, use of dissociating agents and adhesion substrates. The best results were achieved using plates precoated with rattlesnake skin collagen and incubation in CMRL 1415 modified for snake gland cells plus 10% fetal calf serum at 30 degrees C. The presence of venom could be demonstrated in the supernatant of five out of six epithelial-like gland cell cultures tested by ELISA, in the very first passages. After the third passage, however, venom amounts dropped to undetectable values. A total of 23 venom gland cell lines were obtained and are kept frozen in the laboratory; among them, five epithelial-like gland cell lines with up to 12 passages, that were continuously cultured for more than 30 weeks. The methodology described here was successfully applied to C. d. terrificus kidney cells culturing, developed to be used as negative control.

Mating systems, reproductive success, and sexual selection in secretive species: a case study of the western diamond-backed rattlesnake, Crotalus atrox.

Long-term studies of individual animals in nature contribute disproportionately to our understanding of the principles of ecology and evolution. Such field studies can benefit greatly from integrating the methods of molecular genetics with traditional approaches. Even though molecular genetic tools are particularly valuable for species that are difficult to observe directly, they have not been widely adopted. Here, we used molecular genetic techniques in a 10-year radio-telemetric investigation of the western diamond-backed rattlesnake (Crotalus atrox) for an analysis of its mating system and to measure sexual selection. Specifically, we used microsatellite markers to genotype 299 individuals, including neonates from litters of focal females to ascertain parentage using full-pedigree likelihood methods. We detected high levels of multiple paternity within litters, yet found little concordance between paternity and observations of courtship and mating behavior. Larger males did not father significantly more offspring, but we found evidence for size-specific male-mating strategies, with larger males guarding females for longer periods in the mating seasons. Moreover, the spatial proximity of males to mothers was significantly associated with reproductive success. Overall, our field observations alone would have been insufficient to quantitatively measure the mating system of this population of C. atrox, and we thus urge more widespread adoption of molecular tools by field researchers studying the mating systems and sexual selection of snakes and other secretive taxa.