The Brazilian Atlantic Bushmaster Lachesis (Linnaeus, 1766) Mitogenome With Insights On Snake Evolution And Divergence (Serpentes: Viperidae: Crotalinae).

This study presents the first complete mitogenome of the Brazilian Atlantic bushmaster Lachesis with insights into snake evolution. The total length was 17,177 bp, consisting of 13 PCGs, 22 tRNAs, two rRNAs and a duplicate control region (CRs). Almost all genes were encoded by the heavy-strand, except for the ND6 gene and eight tRNAs (tRNA-Gln, Ala, Asn, Cys, Tyr, Ser[TGA anticodon], Glu, Pro). Only ATG, ATA, and ATC were starting codons for protein-coding sequences. Stop codons mainly were TAA, AGA, AGG, and TAG; whereas ND1, ND3, and CYTB terminated with incomplete stop codons. Phylogeny retrieved Lachesis within the Crotalinae as the sister group of Agkistrodon; and the Lachesis+Agkistrodon clade as the sister group of (Sistrurus+Crotalus)+Bothrops. The tree supports Crotalinae, Viperinae, and Azemiopinae in the Viperidae family, being sister taxa of Colubridae+(Elapidae+Psammophiidae). The mean genetic distance across 15 snake families and 57 nucleotide sequences was 0.37. The overall mean value of genetic distance across the Crotalinae was 0.23, with Lachesis muta exhibiting the shortest distance of 0.2 with Agkistrodon piscivorus, Protobothrops dabieshanensis and P. flavoviridis and the greatest 0.25 with Gloydius blomhoffii, Trimeresurus albolabris, S. miliarius, and Deinagkistrodon acutus. The complete Atlantic L. muta mitogenome presented herein is only the third annotated mitogenome from more than 430 described Brazilian snake species.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium.

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

RESUMEN: El veneno de las serpientes puede variar entre y dentro de las especies. Mientras algunos grupos de viperidos del Nuevo Mundo­como las cascabeles­han sido bien estudiadas, muy poco se sabe acerca del veneno de las nauyacas de frío (Cerrophidion) que se encuentran en las zonas altas de Mesoamérica. Comparadas con las extensamente estudiadas cascabeles, que estan ampliamente distribuidas, las poblaciones de Cerrophidion, aisladas en montañas, pueden poseer trayectorias evolutivas y diferenciación en su veneno unicos. En el presente trabajo, describimos el transcriptoma de las glándulas de veneno de poblaciones de C. petlalcalensis, C. tzotzilorum, y C. godmani de México, y un individuo de C. sasai de Costa Rica. Exploramos la variación en la expresión de toxinas en Cerrophidion y la evolución en las secuencias geneticas en C. godmani específicamente. El transcriptoma de la glándula de veneno de Cerrophidion esta compuesto principalmente de Metaloproteinasas de Veneno de Serpiente, Fosfolipasas A[Formula: see text] (PLA[Formula: see text]s), y Serin Proteasas de Veneno de Serpiente. Cerrophidion petlalcalensis presenta poca variación intraespecífica; sin embargo, los transcriptomas de la glandula de veneno de C. godmani y C. tzotzilorum difieren significativamente entre poblaciones geográficamente aisladas. Curiosamente, la variación intraespecífica estuvo atribuida principalmente a la expresión de las toxinas ya que no encontramos señales de selección en las toxinas de C. godmani. Adicionalmente, encontramos miotoxinas similares a PLA[Formula: see text] en todas las especies excepto C. petlalcalensis, y PLA[Formula: see text]s similares a crotoxina en la población sureña de C. godmani. Nuestros resultados demuestran la presencia de variacion intraespecífica presente en el veneno de C. godmani y C. tzotzilorum. Las toxinas de Cerrophidion godmani muestran poca evidencia de selección direccional, y la variación en la secuencias de las toxinas es consistente con evolucion bajo un modelo de equilibrio de mutación-deriva. Algunos individuos de C. godmani de la población del sur potencialmente tienen un veneno neurotóxico dada la presencia de PLA[Formula: see text]s similares a la crotoxina, sin embargo, se necesita más evidencia para corroborar esta hipótesis.

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.

Phylogenetic relationships between three Korean pit viper Gloydius (Serpentes: Crotalinae) species using mitochondrial DNA genes.

BACKGROUND: Molecular phylogenetic studies of the Asian pit viper genus Gloydius have been widely published in Asia, but Korea population have not been conducted till date. OBJECTIVE: This study aimed to analyze the phylogenetic relationships of three Gloydius species (G. saxatilis, G. brevicaudus, and G. ussuriensis) from Korea with other Gloydius species, based on Cytochrome b and ND4. METHODS: We compared 160 samples representing the three species with those of 17 reference species and their phylogenetic status and genetic diversity were analyzed with concatenated sequences of two mitochondrial DNA. RESULTS: Korean G. brevicaudus and G. saxatilis showed high haplotype diversity and relatively low and moderate nucleotide diversity, respectively. Although G. ussuriensis showed high genetic diversity, it was low in the Baengnyeong Island population. The phylogenetic tree represented two major lineages. One major lineage comprised G. ussuriensis, G. tsushimaensis, G. blomhoffii, and G. brevicaudus. The Chinese G. ussuriensis belonged to the same clade as the Korean G. ussuriensis and was closely related to the Baengnyeong Island population. Moreover, G. tsushimaensis was closely related to G. ussuriensis from southwestern Korean and Jeju Island populations. The other major lineage comprised the remaining 12 species and G. saxatilis. Korean G. saxatilis was closely related to G. saxatilis, G. shedanoensis, and G. intermedius from China. CONCLUSION: The phylogenetic status of the Korean Gloydius species in comparison with the other Gloydius species was identified. We suggesting the conservation management unit for the Baengnyeong Island population, while the current conservation status of Korean G. saxatilis is suggested to be revised to a higher level.

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.

Identification of a Novel Cathelicidin from the Deinagkistrodon acutus Genome with Antibacterial Activity by Multiple Mechanisms.

The abuse of antibiotics and the consequent increase of drug-resistant bacteria constitute a serious threat to human health, and new antibiotics are urgently needed. Research shows that antimicrobial peptides produced by natural organisms are potential substitutes for antibiotics. Based on Deinagkistrodonacutus (known as five-pacer viper) genome bioinformatics analysis, we discovered a new cathelicidin antibacterial peptide which was called FP-CATH. Circular dichromatic analysis showed a typical helical structure. FP-CATH showed broad-spectrum antibacterial activity. It has antibacterial activity to Gram-negative bacteria and Gram-positive bacteria including methicillin-resistant Staphylococcus aureus (MRSA). The results of transmission electron microscopy (TEM) and scanning electron microscopy (SEM) showed that FP-CATH could cause the change of bacterial cell integrity, having a destructive effect on Gram-negative bacteria and inducing Gram-positive bacterial surface formation of vesicular structure. FP-CATH could bind to LPS and showed strong binding ability to bacterial DNA. In vivo, FP-CATH can improve the survival rate of nematodes in bacterial invasion experiments, and has a certain protective effect on nematodes. To sum up, FP-CATH is likely to play a role in multiple mechanisms of antibacterial action by impacting bacterial cell integrity and binding to bacterial biomolecules. It is hoped that the study of FP-CATH antibacterial mechanisms will prove useful for development of novel antibiotics.

Whole-exome sequencing and genome-wide evolutionary analyses identify novel candidate genes associated with infrared perception in pit vipers.

Pit vipers possess a unique thermal sensory system consisting of facial pits that allow them to detect minute temperature fluctuations within their environments. Biologists have long attempted to elucidate the genetic basis underlying the infrared perception of pit vipers. Early studies have shown that the TRPA1 gene is the thermal sensor associated with infrared detection in pit vipers. However, whether genes other than TRPA1 are also involved in the infrared perception of pit vipers remains unknown. Here, we sequenced the whole exomes of ten snake species and performed genome-wide evolutionary analyses to search for novel candidate genes that might be involved in the infrared perception of pit vipers. We applied both branch-length-comparison and selection-pressure-alteration analyses to identify genes that specifically underwent accelerated evolution in the ancestral lineage of pit vipers. A total of 47 genes were identified. These genes were significantly enriched in the ion transmembrane transporter, stabilization of membrane potential, and temperature gating activity functional categories. The expression levels of these candidate genes in relevant nerve tissues (trigeminal ganglion, dorsal root ganglion, midbrain, and cerebrum) were also investigated in this study. We further chose one of our candidate genes, the potassium channel gene KCNK4, as an example to discuss its possible role in the infrared perception of pit vipers. Our study provides the first genome-wide survey of infrared perception-related genes in pit vipers via comparative evolutionary analyses and reveals valuable candidate genes for future functional studies.

Discovery of the Gene Encoding a Novel Small Serum Protein (SSP) of Protobothrops flavoviridis and the Evolution of SSPs.

Small serum proteins (SSPs) are low-molecular-weight proteins in snake serum with affinities for various venom proteins. Five SSPs, PfSSP-1 through PfSSP-5, have been reported in Protobothrops flavoviridis ("habu", Pf) serum so far. Recently, we reported that the five genes encoding these PfSSPs are arranged in tandem on a single chromosome. However, the physiological functions and evolutionary origins of the five SSPs remain poorly understood. In a detailed analysis of the habu draft genome, we found a gene encoding a novel SSP, SSP-6. Structural analysis of the genes encoding SSPs and their genomic arrangement revealed the following: (1) SSP-6 forms a third SSP subgroup; (2) SSP-5 and SSP-6 were present in all snake genomes before the divergence of non-venomous and venomous snakes, while SSP-4 was acquired only by venomous snakes; (3) the composition of paralogous SSP genes in snake genomes seems to reflect snake habitat differences; and (4) the evolutionary emergence of SSP genes is probably related to the physiological functions of SSPs, with an initial snake repertoire of SSP-6 and SSP-5. SSP-4 and its derivative, SSP-3, as well as SSP-1 and SSP-2, appear to be venom-related and were acquired later.

Trait differentiation and modular toxin expression in palm-pitvipers.

BACKGROUND: Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A2 (PLA2) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris). RESULTS: We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA2 complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification. CONCLUSION: Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits.

Phylogeny and toxicological assessments of two Porthidium lansbergii lansbergii morphotypes from the caribbean region of Colombia.

After a snakebite accident, species identification is of vital importance. However, the existence of intraspecific differences in the body coloration patterns of venomous snakes can generate confusion and delay a convenient and effective treatment. This is the situation for Porthidium lansbergii lansbergii from Colombia, for which two distinctive color morphs occur, and the relationship of these morphs with venom toxicity is unknown. Therefore, venom samples from specimens of these two morphs were collected from the Colombian Caribbean region, and their protein profiles compared. Likewise, their venom functional activities were evaluated in vitro and in vivo in BALB/C mice. Additionally, using sequences of the mitochondrial cytochrome b (Cyt-b) gene, the relationship between these Colombian P. lansbergii lansbergii morphotypes was investigated, and their phylogenetic positions were determined for the first time using Bayesian inference. Despite the noticeable coloration divergence between the individuals analyzed, similar protein profiles of their venoms were observed. Additionally, neither their lethality nor biochemical activities were notably different. In general, both venoms were highly proteolytic, lacked a coagulant effect in vitro, and extended the clotting time due to the action of venom components, such as disintegrins and proteases, that induce defibrination. These results agreed with the result of our phylogenetic analysis, suggesting that the two chromatic morphs do not represent isolated populations. The phylogenetic analyses also supported the currently recognized P. lansbergii lansbergii subspecies as a monophyletic complex. In conclusion, the results of this investigation suggest similar clinical manifestations regardless of body coloration after a P. lansbergii lansbergii envenomation, and pools can therefore be used for antivenom development, medical treatments, and further research efforts.

Evolutionary history and postglacial colonization of an Asian pit viper (Gloydius halys caucasicus) into Transcaucasia revealed by phylogenetic and phylogeographic analyses.

It has been generally acknowledged that glacial climates at the time of the Pleistocene altered the patterns of species distributions, prompting latitudinal and altitudinal distribution shifts in several species, including poikilothermic species commonly known for their thermal sensitivity. However, the historical phylogeographic patterns of such species have remained largely unknown. Here, we present the historical biogeographic, phylogenetic, and phylogeographic relationships of the Caucasian pit viper, G. h. caucasicus, based on two mtDNA (cyt b and ND4) and one nDNA (c-mos) genes. This pit viper represents the westernmost member of the Crotalinae subfamily in the Palearctic and occurs in a variety of habitats, from 30 m to 3,000 m above sea level. In Iran, it is distributed on the northern and southern slopes of the Alborz Mountains, rendering it a target for phylogenetic and phylogeographic studies of a terrestrial poikilothermic animal. Our study identified four Iranian lineages of G. h. caucasicus along the northeastern to northwestern slopes of the Alborz Mountains and southern Azerbaijan (Talysh Mountains). Diversification of the Iranian lineages highlights population expansion and subsequent isolation into four plausible refugial areas during the Quaternary paleo-climatic oscillations, confirmed by our molecular dating and historical biogeographic analyses. The results of coalescence-based simulations support the incursion of the species from northeastern Iran to the western end of the Alborz, and then toward Transcaucasia via two directions: northern and southern slopes of the Alborz Mountains. Furthermore, our results clearly implied that G. h. caucasicus should be elevated to species rank and further referred to as G. caucasicus (Nikolsky, 1916).

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin.

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based ß-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

Tracing the evolution of venom phospholipases A2 in Gloydius strauchii and related pitvipers: A tale of two acidic isozymes.

Two acidic Asp49-PLA2s with Glu6 substitution and a neutral Lys49-PLA (designated Gst-K49) were cloned from G. strauchii venom glands, their full amino acid sequences were deduced. The predominant acidic PLA2 (designated Gst-E6a) contains 124 residues and the M18W30 substitutions, while the minor acidic PLA2 (designated Gst-E6b) contains 122 residues and the V18A30 substitutions. Their sequences are most similar to those of the respective orthologous PLA2s of G. intermedius venom. Gst-E6a and Gst-E6b appear to be paralogs and possibly have different predatory targets or functions. The LC-MS/MS results indicate the presence of only three PLA2 gene products in the crude venom, the relative expression levels were in the order of Gst-E6a â‰« Gst-E6b > Gst-K49, as confirmed by qPCR results. In contrast to other Gloydius, G. strauchii venom does not contain neurotoxic or basic anticoagulant Asp49-PLA2s, but Gst-K49 is the first Lys49-PLA2 identified in Gloydius venoms. However, its venom content is relatively low and its pI value 7.3 is much lower than those of other Lys49-PLA2s and. The Lys49-PLA2 genes appear to regress in the venom of most of Gloydius and related rattlesnake, and this evolutionary regression occurred before the dispersal of Asian pitvipers to the New World.

Serine protease isoforms in Gloydius intermedius venom: Full sequences, molecular phylogeny and evolutionary implications.

Nine distinct venom serine proteases (vSPs) of Gloydius intermedius were studied by transcriptomic, sub-proteomic and phylogenetic analyses. Their complete amino acid sequences were deduced after Expression Sequence Tag (EST) analyses followed by cDNA cloning and sequencing. These vSPs appear to be paralogs and contain the catalytic triads and 1-4 potential N-glycosylation sites. Their relative expression levels evaluated by qPCR were grossly consistent with their EST hit-numbers. The major vSPs were purified by HPLC and their N-terminal sequences matched well to the deduced sequences, while fragments of the minor vSPs were detected by LC-MS/MS identification. Specific amidolytic activities of the fractions from HPLC and anion exchange separation were assayed using four chromogenic substrates, respectively. Molecular phylogenetic tree based on the sequences of these vSPs and their orthologs revealed six major clusters, one of them covered four lineages of plasminogen activator like vSPs. N-glycosylation patterns and variations for the vSPs are discussed. The high sequence similarities between G. intermedius vSPs and their respective orthologs from American pitvipers suggest that most of the isoforms evolved before Asian pitvipers migrated to the New World. Our results also indicate that the neurotoxic venoms contain more kallikrein-like vSPs and hypotensive components than the hemorrhagic venoms. SIGNIFICANCE: Full sequences and expression levels of nine paralogous serine proteases (designated as GiSPs) of Gloydius intermedius venom have been studied. A kallikrein-like enzyme is most abundant and four isoforms homologous to venom plasminogen-activators are also expressed in this venom. Taken together, the present and previous data demonstrate that the neurotoxic G. intermedius venoms contain more hypotensive vSPs relative to other hemorrhagic pitviper venoms and the pitviper vSPs are highly versatile and diverse. Their structure-function relationships remain to be explored and compared. A novel, simplified phylogenetic tree based on the sequences of GiSPs and their closely related orthologs from other pitvipers reveals six major subtypes and offers a better understanding of vSP duplication and evolution in pitvipers of both the Old and New Worlds. It is well known that specific vSPs are potential therapeutic or diagnostic agents that target the plasma proteins or coagulation factors. Our results not only render deeper insights into the variation and evolution of vSPs, but may help to choose right venoms for the development of better therapeutic leads.

Venomics of Tropidolaemus wagleri, the sexually dimorphic temple pit viper: Unveiling a deeply conserved atypical toxin arsenal.

Tropidolaemus wagleri (temple pit viper) is a medically important snake in Southeast Asia. It displays distinct sexual dimorphism and prey specificity, however its venomics and inter-sex venom variation have not been thoroughly investigated. Applying reverse-phase HPLC, we demonstrated that the venom profiles were not significantly affected by sex and geographical locality (Peninsular Malaya, insular Penang, insular Sumatra) of the snakes. Essentially, venoms of both sexes share comparable intravenous median lethal dose (LD50) (0.56-0.63 µg/g) and cause neurotoxic envenomation in mice. LCMS/MS identified six waglerin forms as the predominant lethal principles, comprising 38.2% of total venom proteins. Fourteen other toxin-protein families identified include phospholipase A2, serine proteinase, snaclec and metalloproteinase. In mice, HPLC fractions containing these proteins showed insignificant contribution to the overall venom lethality. Besides, the unique elution pattern of approximately 34.5% of non-lethal, low molecular mass proteins (3-5 kDa) on HPLC could be potential biomarker for this primitive crotalid species. Together, the study unveiled the venom proteome of T. wagleri that is atypical among many pit vipers as it comprises abundant neurotoxic peptides (waglerins) but little hemotoxic proteinases. The findings also revealed that the venom is relatively well conserved intraspecifically despite the drastic morphological differences between sexes.

The complete mitochondrial genome sequence of Gloydius shedaoensis (Squamata: Viperidae).

Gloydius shedaoensis is an insular and vulnerable pitviper that is endemic to Snake Island, northeastern China. In this study, we successfully sequenced mitochondrial genomes of two individuals of G. shedaoensis. The complete mitochondrial genomes of G. shedaoensis are circular molecular with 17 222 and 17 221 bp in length respectively, which both contain 2 ribosomal RNA (rRNA) genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes, an origin of light-strand replication (OL) and two non-coding control regions. Compared with previously published mitochondrial genomes of Gloydius species, the base composition and gene arrangement are rather conservative. A Bayesian phylogenetic tree using the complete mitochondrial genomes of all viper species available showed a consistent result with previous studies.

The complete mitochondrial genome of Gloydius intermedius (Squamata: Viperidae: Crotalinae) from China.

The mitochondrial genome sequence of Gloydius intermedius is analyzed and presented for the first time. The genome was 17, 226 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 2 control region. The overall base composition was A (32.4%), C (28.8%), T (25.9%), and G (12.9%). The base compositions clearly presented the A-C skew, which was most obvious in the protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of G. intermedius. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees, indicating a close phylogenetic affinity of the 13 Crotalinae species. It appeared that no less than two major phyletic lineages were present in Crotalinae. The main clades within the Crotalinae supported are: A clade including the Protobothrops. A clade (G. brevicaudus, G. ussuriensis, G. intermedius, G. saxatilis) with the Ovophis appeared as the sister taxon to Protobothrops and was supported by bootstrap values of 88%. The four Gloydius species formed a paraphyletic group with the high bootstrap value (100 %) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of Gloydius saxatilis (Squamata: Viperidae: Crotalinae).

The mitochondrial genome sequence of Gloydius saxatilis is analyzed and presented for the public for the first time. The genome was 17,218 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 2 control regions. The overall base composition was A (32.3%), C (28.9%), T (25.8%), and G (13.0%). The base compositions presented clearly the A-C skew, which was most obviously in the protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of G. saxatilis. Mitochondrial genomes analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees, indicating a close phylogenetic affinity of the thirteen Crotalinae species. It appeared that no less than two major phyletic lineages were present in Crotalinae. The main clades within the Crotalinae supported are: A clade including the Protobothrops. A clade (G. brevicaudus, G. ussuriensis, G. intermedius, G. saxatilis) with the Ovophis as the sister taxon to Protobothrops and was supported by bootstrap values of 88%. The four Gloydius species formed a paraphyletic group with the high bootstrap value (100%) in all examinations.