A genome assembly for the southern Pacific rattlesnake, Crotalus oreganus helleri, in the western rattlesnake species complex.

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake (Crotalus oreganus) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus, in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general.

A Promising Biomolecule Able to Degrade Neutrophil Extracellular Traps: CdcPDE, a Rattlesnake Phosphodiesterase.

Neutrophil extracellular traps (NETs) are an important mechanism for defense against pathogens. Their overproduction can be harmful since excessive NET formation promotes inflammation and tissue damage in several diseases. Nucleases are capable to degrade NET on basis of their DNA hydrolysis activity, including the CdcPDE, a nuclease isolated from Crotalus durissus collilineatus snake venom. Here, we report a new finding about CdcPDE activity, demonstrating its efficiency in degrading cell-free DNA from NETs, being a potential candidate to assist in therapies targeting inflammatory diseases.

Montane Rattlesnakes in México: Venoms of Crotalus tancitarensis and Related Species within the Crotalus intermedius Group.

The Crotalus intermedius group is a clade of rattlesnakes consisting of several species adapted to a high elevation habitat, primarily in México. Crotalus tancitarensis was previously classified as C. intermedius, until individuals occurring on Cerro Tancítaro in Michoacán, México, were reevaluated and classified as a new species (C. tancitarensis) based on scale pattern and geographic location. This study aimed to characterize the venom of C. tancitarensis and compare the venom profile to those of other species within the Crotalus intermedius group using gel electrophoresis, biochemical assays, reverse-phase high performance liquid chromatography, mass spectrometry, and lethal toxicity (LD50) assays. Results show that the venom profiles of species within the Crotalus intermedius group are similar, but with distinct differences in phospholipase A2 (PLA2), metalloproteinase PI (SVMP PI), and kallikrein-like serine proteinase (SVSP) activity and relative abundance. Proteomic analysis indicated that the highland forms produce venoms with 50-60 protein isoforms and a composition typical of type I rattlesnake venoms (abundant SVMPs, lack of presynaptic PLA2-based neurotoxins), as well as a diversity of typical Crotalus venom components such as serine proteinases, PLA2s, C-type lectins, and less abundant toxins (LAAOs, CRiSPs, etc.). The overall venom profile of C. tancitarensis appears most similar to C. transversus, which is consistent with a previous mitochondrial DNA analysis of the Crotalus intermedius group. These rattlesnakes of the Mexican highlands represent a radiation of high elevation specialists, and in spite of divergence of species in these Sky Island habitats, venom composition of species analyzed here has remained relatively conserved. The majority of protein family isoforms are conserved in all members of the clade, and as seen in other more broadly distributed rattlesnake species, differences in their venoms are largely due to relative concentrations of specific components.

Insights on the inhibition properties of Jatromollistatin (a cyclic heptapeptide) against Crotalus adamanteus metalloendopeptidase using molecular docking analysis.

Jatropha mollissima is endemic to Brazil and is used for traditional medicinal purposes, including the treatment of snakebite. In this study, latex obtained from this plant was fractioned using reversed-phase chromatography, and the fractions were then screened for peptides. A 755 g/mol peptide was obtained, and MS/MS analyses indicated it had a cyclic sequence (Pro-Leu-Gly-Val-Leu-Leu-Tyr). This peptide sequence was present in the Jatropha genome database, and an identity value of 90.71%, an E-value of 0.0, and a score of 883 with NO-associated protein 1/chloroplastic/mitochondria of Jatropha curcas were obtained from the NCBI nonredundant protein sequence (nr) database. Molecular docking analyses performed with the peptide against a metalloendopeptidase belonging to Crotalus adamanteus snake venom suggested the cyclic peptide establishes favorable interactions with the catalytic site of the enzyme. Therefore, it could inhibit enzyme catalysis. This belief was corroborated by the formation of 6 hydrogen bonds with the linear form of the peptide. Tighter complexation of the cyclic form (41 kcal/mol more energetic) revealed better spatial blocking. The linear form outperformed the cyclic form in complexing the required energy, recruiting more catalytic residues (6/2), and in establishing more hydrogen bonds (6/3). However, cyclic folding provided a more significant spatial block within the catalytic site. The set of results suggests that the cycle peptide, here called Jatromollistatin, which was previously described as jatrophidin and pohlianin A in two other species of Jatropha, is a promising candidate to inhibit venom proteases. This belief is corroborated by the topical use of the latex for initial treatment of snakebites.

SEQUENCE SLIDER: integration of structural and genetic data to characterize isoforms from natural sources.

Proteins isolated from natural sources can be composed of a mixture of isoforms with similar physicochemical properties that coexist in the final steps of purification. Yet, even where unverified, the assumed sequence is enforced throughout the structural studies. Herein, we propose a novel perspective to address the usually neglected sequence heterogeneity of natural products by integrating biophysical, genetic and structural data in our program SEQUENCE SLIDER. The aim is to assess the evidence supporting chemical composition in structure determination. Locally, we interrogate the experimental map to establish which side chains are supported by the structural data, and the genetic information relating sequence conservation is integrated into this statistic. Hence, we build a constrained peptide database, containing most probable sequences to interpret mass spectrometry data (MS). In parallel, we perform MS de novo sequencing with genomic-based algorithms to detect point mutations. We calibrated SLIDER with Gallus gallus lysozyme, whose sequence is unequivocally established and numerous natural isoforms are reported. We used SLIDER to characterize a metalloproteinase and a phospholipase A2-like protein from the venom of Bothrops moojeni and a crotoxin from Crotalus durissus collilineatus. This integrated approach offers a more realistic structural descriptor to characterize macromolecules isolated from natural sources.

A Clot Twist: Extreme Variation in Coagulotoxicity Mechanisms in Mexican Neotropical Rattlesnake Venoms.

Rattlesnakes are a diverse clade of pit vipers (snake family Viperidae, subfamily Crotalinae) that consists of numerous medically significant species. We used validated in vitro assays measuring venom-induced clotting time and strength of any clots formed in human plasma and fibrinogen to assess the coagulotoxic activity of the four medically relevant Mexican rattlesnake species Crotalus culminatus, C. mictlantecuhtli, C. molossus, and C. tzabcan. We report the first evidence of true procoagulant activity by Neotropical rattlesnake venom in Crotalus culminatus. This species presented a strong ontogenetic coagulotoxicity dichotomy: neonates were strongly procoagulant via Factor X activation, whereas adults were pseudo-procoagulant in that they converted fibrinogen into weak, unstable fibrin clots that rapidly broke down, thereby likely contributing to net anticoagulation through fibrinogen depletion. The other species did not activate clotting factors or display an ontogenetic dichotomy, but depleted fibrinogen levels by cleaving fibrinogen either in a destructive (non-clotting) manner or via a pseudo-procoagulant mechanism. We also assessed the neutralization of these venoms by available antivenom and enzyme-inhibitors to provide knowledge for the design of evidence-based treatment strategies for envenomated patients. One of the most frequently used Mexican antivenoms (Bioclon Antivipmyn®) failed to neutralize the potent procoagulant toxic action of neonate C. culminatus venom, highlighting limitations in snakebite treatment for this species. However, the metalloprotease inhibitor Prinomastat substantially thwarted the procoagulant venom activity, while 2,3-dimercapto-1-propanesulfonic acid (DMPS) was much less effective. These results confirm that venom-induced Factor X activation (a procoagulant action) is driven by metalloproteases, while also suggesting Prinomastat as a more promising potential adjunct treatment than DMPS for this species (with the caveat that in vivo studies are necessary to confirm this potential clinical use). Conversely, the serine protease inhibitor 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (AEBSF) inhibited the direct fibrinogen cleaving actions of C. mictlantecuhtli venom, thereby revealing that the pseudo-procoagulant action is driven by kallikrein-type serine proteases. Thus, this differential ontogenetic variation in coagulotoxicity patterns poses intriguing questions. Our results underscore the need for further research into Mexican rattlesnake venom activity, and also highlights potential limitations of current antivenom treatments.

Crotamine in Crotalus durissus: distribution according to subspecies and geographic origin, in captivity or nature

Background:Crotalus durissus is considered one of the most important species of venomous snakes in Brazil, due to the high mortality of its snakebites. The venom of Crotalus durissus contains four main toxins: crotoxin, convulxin, gyroxin and crotamine. Venoms can vary in their crotamine content, being crotamine-negative or -positive. This heterogeneity is of great importance for producing antivenom, due to their different mechanisms of action. The possibility that antivenom produced by Butantan Institute might have a different immunorecognition capacity between crotamine-negative and crotamine-positive C. durissus venoms instigated us to investigate the differences between these two venom groups. Methods:The presence of crotamine was analyzed by SDS-PAGE, western blotting and ELISA, whereas comparison between the two types of venoms was carried out through HPLC, mass spectrometry analysis as well as assessment of antivenom lethality and efficacy. Results:The results showed a variation in the presence of crotamine among the subspecies and the geographic origin of snakes from nature, but not in captive snakes. Regarding differences between crotamine-positive and -negative venoms, some exclusive proteins are found in each pool and the crotamine-negative pool presented more phospholipase A2 than crotamine-positive pool. This variation could affect the time to death, but the lethal and effective dose were not affected. Conclusion:These differences between venom pools indicate the importance of using both, crotamine-positive and crotamine-negative venoms, to produce the antivenom.(AU)

Bioaffinity Fishing Procedure Using Secretory Phospholipase A2 for Screening for Bioactive Components: Modulation of Pharmacological Effect Induced by sPLA2 from Crotalus durissus terrificus by Hispidulin from Moquiniastrum floribundum.

Bioaffinity capturing of molecules allows the discovery of bioactive compounds and decreases the need for various stages in the natural compound isolation process. Despite the high selectivity of this technique, the screening and identification methodology depends on the presence of a protein to capture potential ligands. However, some proteins, such as snake secretory phospholipase A2 (sPLA2), have never been investigated using this approach. The purpose of this study was to evaluate the use of a new method for screening natural compounds using a bioaffinity-guided ultrafiltration method on Crotalus durissus terrificus sPLA2 followed by HPLC-MS to identify the compounds, and this method could be used to discover new anti-inflammatory compounds from the various organisms originating from biodiversity. Different extracts were selected to evaluate their ability to inhibit sPLA2 activity. The extracts were incubated with sPLA2 and the resulting mixture was ultrafiltrated to elute unbound components. The resulting compounds were identified by HPLC-MS. We identified hispidulin as one of the components present in the Moquiniastrum floribundum leaf and evaluated the ability of this isolated compound to neutralize the inflammatory activity of sPLA2 from Crotalus durissus terrificus.

Cryptic diversity in the Mexican highlands: Thousands of UCE loci help illuminate phylogenetic relationships, species limits and divergence times of montane rattlesnakes (Viperidae: Crotalus).

With the continued adoption of genome-scale data in evolutionary biology comes the challenge of adequately harnessing the information to make accurate phylogenetic inferences. Coalescent-based methods of species tree inference have become common, and concatenation has been shown in simulation to perform well, particularly when levels of incomplete lineage sorting are low. However, simulation conditions are often overly simplistic, leaving empiricists with uncertainty regarding analytical tools. We use a large ultraconserved element data set (>3,000 loci) from rattlesnakes of the Crotalus triseriatus group to delimit lineages and estimate species trees using concatenation and several coalescent-based methods. Unpartitioned and partitioned maximum likelihood and Bayesian analysis of the concatenated matrix yield a topology identical to coalescent analysis of a subset of the data in bpp. ASTRAL analysis on a subset of the more variable loci also results in a tree consistent with concatenation and bpp, whereas the SVDquartets phylogeny differs at additional nodes. The size of the concatenated matrix has a strong effect on species tree inference using SVDquartets, warranting additional investigation on optimal data characteristics for this method. Species delimitation analyses suggest up to 16 unique lineages may be present within the C. triseriatus group, with divergences occurring during the Neogene and Quaternary. Network analyses suggest hybridization within the group is relatively rare. Altogether, our results reaffirm the Mexican highlands as a biodiversity hotspot and suggest that coalescent-based species tree inference on data subsets can provide a strongly supported species tree consistent with concatenation of all loci with a large amount of missing data.

Genetic diversity and structure of Crotalus triseriatus, a rattlesnake of central Mexico.

The isolated and fragmented populations are highly susceptible to stochastic events, increasing the extinction risk because of the decline in putative adaptive potential and individual fitness. The population has high heterozygosity values and a moderate allelic diversity, the heterozygosity values are higher than in most other Crotalus species and snake studies. Possibly these high levels of genetic diversity can be related to a large founder size, high effective population size, multiple paternity and overlapping generations. We did not find the genetic structuring but the effective number of alleles (Ne) was 138.1. We found evidence of bottlenecks and the majority of rattlesnakeswere unrelated, despite the small sample size, endemic status, the isolated and fragmented habitat. The genetic information provided in this study can be useful as a first approach to try to make informed conservation efforts for this species and also, important to preserve the habitat of this species; the endangered Abies-Pinus forest of the Nevado the Toluca Volcano.

An acidic model pro-peptide affects the secondary structure, membrane interactions and antimicrobial activity of a crotalicidin fragment.

In order to study how acidic pro-peptides inhibit the antimicrobial activity of antimicrobial peptides, we introduce a simple model system, consisting of a 19 amino-acid long antimicrobial peptide, and an N-terminally attached, 10 amino-acid long acidic model pro-peptide. The antimicrobial peptide is a fragment of the crotalicidin peptide, a member of the cathelidin family, from rattlesnake venom. The model pro-peptide is a deca (glutamic acid). Attachment of the model pro-peptide only leads to a moderately large reduction in the binding to- and induced leakage of model liposomes, while the antimicrobial activity of the crotalicidin fragment is completely inhibited by attaching the model pro-peptide. Attaching the pro-peptide induces a conformational change to a more helical conformation, while there are no signs of intra- or intermolecular peptide complexation. We conclude that inhibition of antimicrobial activity by the model pro-peptide might be related to a conformational change induced by the pro-peptide domain, and that additional effects beyond induced changes in membrane activity must also be involved.

Cryptic genetic diversity, population structure, and gene flow in the Mojave rattlesnake (Crotalus scutulatus).

The Mojave rattlesnake (Crotalus scutulatus) inhabits deserts and arid grasslands of the western United States and Mexico. Despite considerable interest in its highly toxic venom and the recognition of two subspecies, no molecular studies have characterized range-wide genetic diversity and population structure or tested species limits within C. scutulatus. We used mitochondrial DNA and thousands of nuclear loci from double-digest restriction site associated DNA sequencing to infer population genetic structure throughout the range of C. scutulatus, and to evaluate divergence times and gene flow between populations. We find strong support for several divergent mitochondrial and nuclear clades of C. scutulatus, including splits coincident with two major phylogeographic barriers: the Continental Divide and the elevational increase associated with the Central Mexican Plateau. We apply Bayesian clustering, phylogenetic inference, and coalescent-based species delimitation to our nuclear genetic data to test hypotheses of population structure. We also performed demographic analyses to test hypotheses relating to population divergence and gene flow. Collectively, our results support the existence of four distinct lineages within C. scutulatus, and genetically defined populations do not correspond with currently recognized subspecies ranges. Finally, we use approximate Bayesian computation to test hypotheses of divergence among multiple rattlesnake species groups distributed across the Continental Divide, and find evidence for co-divergence at this boundary during the mid-Pleistocene.

Integrated Venomics and Venom Gland Transcriptome Analysis of Juvenile and Adult Mexican Rattlesnakes Crotalus simus, C. tzabcan, and C. culminatus Revealed miRNA-modulated Ontogenetic Shifts.

Adult rattlesnakes within genus Crotalus express one of two distinct venom phenotypes, type I (hemorrhagic) and type II (neurotoxic). In Costa Rican Central American rattlesnake, ontogenetic changes in the concentration of miRNAs modulate venom type II to type I transition. Venomics and venom gland transcriptome analyses showed that adult C. simus and C. tzabcan expressed intermediate patterns between type II and type I venoms, whereas C. culminatus had a canonical type I venom. Neonate/juvenile and adult Mexican rattlesnakes showed notable inter- and intraspecific variability in the number, type, abundance and ontogenetic shifts of the transcriptional and translational venom gland activities. These results support a role for miRNAs in the ontogenetic venom compositional changes in the three congeneric Mexican rattlesnakes. It is worth noting the finding of dual-action miRNAs, which silence the translation of neurotoxic heterodimeric PLA2 crotoxin and acidic PLA2 mRNAs while simultaneously up-regulating SVMP-targeting mRNAs. Dual transcriptional regulation potentially explains the existence of mutually exclusive crotoxin-rich (type-II) and SVMP-rich (type-I) venom phenotypic dichotomy among rattlesnakes. Our results support the hypothesis that alterations of the distribution of miRNAs, modulating the translational activity of venom gland toxin-encoding mRNAs in response to an external cue, may contribute to the mechanism generating adaptive venom variability.

Is Hybridization a Source of Adaptive Venom Variation in Rattlesnakes? A Test, Using a Crotalus scutulatus × viridis Hybrid Zone in Southwestern New Mexico.

Venomous snakes often display extensive variation in venom composition both between and within species. However, the mechanisms underlying the distribution of different toxins and venom types among populations and taxa remain insufficiently known. Rattlesnakes (Crotalus, Sistrurus) display extreme inter- and intraspecific variation in venom composition, centered particularly on the presence or absence of presynaptically neurotoxic phospholipases A2 such as Mojave toxin (MTX). Interspecific hybridization has been invoked as a mechanism to explain the distribution of these toxins across rattlesnakes, with the implicit assumption that they are adaptively advantageous. Here, we test the potential of adaptive hybridization as a mechanism for venom evolution by assessing the distribution of genes encoding the acidic and basic subunits of Mojave toxin across a hybrid zone between MTX-positive Crotalus scutulatus and MTX-negative C. viridis in southwestern New Mexico, USA. Analyses of morphology, mitochondrial and single copy-nuclear genes document extensive admixture within a narrow hybrid zone. The genes encoding the two MTX subunits are strictly linked, and found in most hybrids and backcrossed individuals, but not in C. viridis away from the hybrid zone. Presence of the genes is invariably associated with presence of the corresponding toxin in the venom. We conclude that introgression of highly lethal neurotoxins through hybridization is not necessarily favored by natural selection in rattlesnakes, and that even extensive hybridization may not lead to introgression of these genes into another species.

Structural and evolutionary insights into endogenous alpha-phospholipase A2 inhibitors of Latin American pit vipers.

Phospholipases A2 are major components of snake venoms (svPLA2s) and are able to induce multiple local and systemic deleterious effects upon envenomation. Several snake species are provided with svPLA2 inhibitors (sbPLIs) in their circulating blood, which confer a natural resistance against the toxic components of homologous and heterologous venoms. The sbPLIs belong to any of three structural classes named α, ß and γ. In the present study, we identified, characterized and performed structural and evolutionary analyses of sbαPLIs transcripts and the encoded proteins, in the most common Latin American pit vipers belonging to Crotalus, Bothrops and Lachesis genera. Mutation data indicated that sbαPLIs from Latin American snakes might have evolved in an accelerated manner, similarly to that reported for sbαPLIs from Asian snakes, and possibly co-evoluted with svPLA2s in response to the diversity of target enzymes. The importance of sbαPLI trimerization for the effective binding and inhibition of acidic svPLA2s is discussed and conserved cationic residues located at the central pore of the inhibitor trimer are suggested to be a significant part of the binding site of sbαPLIs to acidic svPLA2s. Our data contribute to the current body of knowledge on the structural and evolutionary characteristics of sbPLIs, in general, and may assist in the future development of selective inhibitors for secretory PLA2 from several sources.

Diversity-dependent cladogenesis throughout western Mexico: Evolutionary biogeography of rattlesnakes (Viperidae: Crotalinae: Crotalus and Sistrurus).

Rattlesnakes (Crotalus and Sistrurus) represent a radiation of approximately 42 species distributed throughout the New World from southern Canada to Argentina. Interest in this enigmatic group of snakes continues to accrue due, in part, to their ecomorphological diversity, contributions to global envenomations, and potential medicinal importance. Although the group has garnered substantial attention from systematists and evolutionary biologists for decades, little is still known regarding patterns of lineage diversification. In addition, few studies have statistically quantified broad-scale biogeographic patterns in rattlesnakes to ascertain how dispersal occurred throughout the New World, particularly among the different major biomes of the Americas. To examine diversification and biogeographic patterns in this group of snakes we assemble a multilocus data set consisting of over 6700bp encompassing three nuclear loci (NT-3, RAG-1, C-mos) and seven mitochondrial genes (12S, 16S, ATPase6, ATPase8, ND4, ND5, cytb). Fossil-calibrated phylogenetic and subsequent diversification rate analyses are implemented using maximum likelihood and Bayesian inference, to examine their evolutionary history and temporal dynamics of diversity. Based on ancestral area reconstructions we explore dispersal patterns throughout the New World. Cladogenesis occurred predominantly during the Miocene and Pliocene with only two divergences during the Pleistocene. Two different diversification rate models, advocating diversity-dependence, are strongly supported. These models indicate an early rapid radiation followed by a recent speciation rate decline. Biogeographic analyses suggest that the high elevation pine-oak forests of western Mexico served as a major speciation pump for the majority of lineages, with the desert biome of western North America colonized independently at least twice. All together, these results provide evidence for rapid diversification of rattlesnakes throughout the Mexican highlands during the Neogene, likely in response to continual orogenesis of Mexico's major mountain systems, followed by more recent dispersal into desert and tropical biomes.

Genetic diversity and genetic structure of an endemic Mexican Dusky Rattlesnake (Crotalus triseriatus) in a highly modified agricultural landscape: implications for conservation.

It is necessary to determine genetic diversity of fragmented populations in highly modified landscapes to understand how populations respond to land-use change. This information will help guide future conservation and management strategies. We conducted a population genetic study on an endemic Mexican Dusky Rattlesnake (Crotalus triseriatus) in a highly modified landscape near the Toluca metropolitan area, in order to provide crucial information for the conservation of this species. There was medium levels of genetic diversity, with a few alleles and genotypes. We identified three genetically differentiated clusters, likely as a result of different habitat cover type. We also found evidence of an ancestral genetic bottleneck and medium values of effective population size. Inbreeding coefficients were low and there was a moderate gene flow. Our results can be used as a basis for future research and C. triseriatus conservation efforts, particularly considering that the Trans-Mexican Volcanic Belt is heavily impacted by destructive land-use practices.

Incipient speciation with biased gene flow between two lineages of the Western Diamondback Rattlesnake (Crotalus atrox).

We used mitochondrial DNA sequence data from 151 individuals to estimate population genetic structure across the range of the Western Diamondback Rattlesnake (Crotalus atrox), a widely distributed North American pit viper. We also tested hypotheses of population structure using double-digest restriction site associated DNA (ddRADseq) data, incorporating thousands of nuclear genome-wide SNPs from 42 individuals. We found strong mitochondrial support for a deep divergence between eastern and western C. atrox populations, and subsequent intermixing of these populations in the Inter-Pecos region of the United States and Mexico. Our nuclear RADseq data also identify these two distinct lineages of C. atrox, and provide evidence for nuclear admixture of eastern and western alleles across a broad geographic region. We identified contrasting patterns of mitochondrial and nuclear genetic variation across this genetic fusion zone that indicate partially restricted patterns of gene flow, which may be due to either pre- or post-zygotic isolating mechanisms. The failure of these two lineages to maintain complete genetic isolation, and evidence for partially-restricted gene flow, imply that these lineages were in the early stages of speciation prior to secondary contact.

Phylogenetic relationships of the enigmatic longtailed rattlesnakes (Crotalus ericsmithi, C. lannomi, and C. stejnegeri).

The longtailed rattlesnakes of western Mexico represent an enigmatic group of poorly known venomous snake species: Crotalus ericsmithi, C. lannomi, and C. stejnegeri. In the 120 years since their discovery, fewer than twenty individuals have been deposited in natural history collections worldwide. These three species share similar morphological traits, including a particularly long tail that has been interpreted as either an ancestral condition among rattlesnakes or as derived within the longtailed group. An understanding of the phylogenetic distinctiveness and relationships among the longtailed rattlesnakes, and their relationships to other rattlesnake groups, has previously been hampered by a dearth of comparative material and tissues for collection of DNA sequence data. Facilitated by the recent availability of tissue samples from multiple individuals of each species, we estimate the phylogenetic relationships among the longtailed rattlesnakes and their placement among other rattlesnake groups, using DNA sequence data from three mitochondrial and three nuclear gene fragments. We explore phylogenetic signal in our data using Bayesian and maximum likelihood methods, species tree analyses and hypothesis testing. Our results strongly support the monophyly of longtailed rattlesnakes and suggest the three species diverged from each other during the mid to late Pliocene or early Pleistocene (~1.5-5.6 mya). Contrary to prevailing hypotheses, we find no evidence for an early or basal divergence of the longtailed clade within the rattlesnake tree, and instead estimate that it diverged relatively recently (~6.8 mya) from its sister lineage, composed of the diamondback rattlesnakes (C. atrox group) and the prairie rattlesnakes (C. viridis group). With our added sampling of lineages and identification of previously used problematic sequences, we provide a revised hypothesis for relationships among Crotalus species, yet underscore the need for future studies and new data to deliver a well-supported robust estimate of rattlesnake relationships.

Phylogenetic analysis of ß-defensin-like genes of Bothrops, Crotalus and Lachesis snakes.

Defensins are components of the vertebrate innate immune system; they comprise a diverse group of small cationic antimicrobial peptides. Among them, ß-defensins have a characteristic ß-sheet-rich fold plus six conserved cysteines with particular spacing and intramolecular bonds. They have been fully studied in mammals, but there is little information about them in snakes. Using a PCR approach, we described 13 ß-defensin-like sequences in Bothrops and Lachesis snakes. The genes are organized in three exons and two introns, with exception of B.atrox_defensinB_01 which has only two exons. They show high similarities in exon 1, intron 1 and intron 2, but exons 2 and 3 have undergone accelerated evolution. The theoretical translated sequences encode a pre-ß-defensin-like molecule with a conserved signal peptide and a mature peptide. The signal peptides are leucine-rich and the mature ß-defensin-like molecules have a size around 4.5 kDa, a net charge from +2 to +11, and the conserved cysteine motif. Phylogenetic analysis was done using maximum parsimony, maximum likelihood and Bayesian analyses, and all resulted in similar topologies with slight differences. The genus Bothrops displayed two separate lineages. The reconciliation of gene trees and species tree indicated eight to nine duplications and 23 to 29 extinctions depending on the gene tree used. Our results together with previously published data indicate that the ancestral ß-defensin-like gene may have three exons in vertebrates and that their evolution occurred according to a birth-and-death model.