A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function

A plethora of rodents: Rattlesnake predators generate unanticipated patterns of venom resistance in a grassland ecosystem

Predation has the potential to impart strong selective pressures on organisms within their environments, resulting in adaptive changes in prey that minimize risk of predation. Pressures from venomous snakes present an exceptional challenge to prey, as venom represents a unique chemical arsenal evolutionarily tailored to incapacitate prey. In response, venom resistance has been detected in various snake prey species, and to varying degrees. This study analyzes venom resistance in an eastern Colorado grassland habitat, where the Prairie Rattlesnake (Crotalus viridis) and Desert Massasauga Rattlesnake (Sistrurus tergeminus edwardsii) co-occur with a suite of grassland rodents. We test for venom resistance across rodent and snake pairings using two geographically distant field sites to determine the role of 1) predation pressure and trophic ecology, and 2) sympatric and allopatric patterns of venom resistance. Resistance was measured using serum-based metalloproteinase inhibition assays to determine potential inhibition of proteolytic activity, augmented by median lethal dose (LD50) assays on rodent species to assess toxicity of crude venoms. Resistance is present in several rodent species, with strong resistance present in populations of Eastern Woodrat (Neotoma floridana), Ord's Kangaroo Rat (Dipodomys ordii), and Northern Grasshopper Mouse (Onychomys leucogaster). Resistance is less developed in other species, including the House Mouse (Mus musculus) and Plains Pocket Mouse (Perognathus flavescens). An unexpected differential is present, where Lincoln County Kangaroo Rats are highly resistant to venom of co-occurring Prairie Rattlesnakes yet are sensitive to an allopatric population of Prairie Rattlesnakes in Weld County. Lincoln Co. Northern Grasshopper Mice also demonstrate extremely elevated resistance to Weld Co. Prairie Rattlesnake venoms, and they may possess resistance mechanisms for myotoxin a, an abundant component of Weld Co. C. v viridis venoms. This study illustrates the complexity of venom resistance in biological communities that can exist when incorporating multiple species interactions. Future studies aimed at characterizing resistance mechanisms at the molecular level will provide a more detailed physiological context for understanding mechanisms by which resistance to venoms occurs.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

Exploring the Evolvability of Plant Specialized Metabolism: Uniqueness Out Of Uniformity and Uniqueness Behind Uniformity

The huge structural diversity exhibited by plant specialized metabolites has primarily been considered to result from the catalytic specificity of their biosynthetic enzymes. Accordingly, enzyme gene multiplication and functional differentiation through spontaneous mutations have been established as the molecular mechanisms that drive metabolic evolution. Nevertheless, how plants have assembled and maintained such metabolic enzyme genes and the typical clusters that are observed in plant genomes, as well as why identical specialized metabolites often exist in phylogenetically remote lineages, is currently only poorly explained by a concept known as convergent evolution. Here, we compile recent knowledge on the co-presence of metabolic modules that are common in the plant kingdom but have evolved under specific historical and contextual constraints defined by the physicochemical properties of each plant specialized metabolite and the genetic presets of the biosynthetic genes. Furthermore, we discuss a common manner to generate uncommon metabolites (uniqueness out of uniformity) and an uncommon manner to generate common metabolites (uniqueness behind uniformity). This review describes the emerging aspects of the evolvability of plant specialized metabolism that underlie the vast structural diversity of plant specialized metabolites in nature.

Venom-gland transcriptomics of the Malayan pit viper (Calloselasma rhodostoma) for identification, classification, and characterization of venom proteins

The Malayan pit viper (Calloselasma rhodostoma) is a hemotoxic snake widely found in Southeast Asia and is responsible for the majority of poisoning cases in this region, including Thailand. However, a comprehensive knowledge of the venom protein profile and classification, as well as novel venom proteins, of this viper is still limited. Recently, the detailed composition of several snake venoms has been discovered through the use of transcriptome analysis. Therefore, the aim of this study was to employ a next-generation sequencing platform and bioinformatics analysis to undertake venom-gland de novo transcriptomics of Malayan pit vipers. Furthermore, 21,272 functional coding genes were identified from 36,577 transcripts, of which 314 transcripts were identified as toxin proteins, accounting for 61.41% of total FPKM, which can be categorized into 22 toxin gene families. The most abundant are snake venom metalloproteinase kistomin (P0CB14) and zinc metalloproteinase/disintegrin (P30403), which account for 60.47% of total toxin FPKM and belong to the SVMP toxin family, followed by snake venom serine protease 1 (O13059) and Snaclec rhodocetin subunit beta (P81398), which account for 6.84% and 5.50% of total toxin FPKM and belong to the snake venom serine protease (SVSP) and Snaclec toxin family, respectively. Amino acid sequences of the aforementioned toxins were compared with those identified in other important medical hemotoxic snakes from Southeast Asia, including the Siamese Russell's viper (Daboia siamensis) and green pit viper (Trimeresurus albolabris), in order to analyze their protein homology. The results demonstrated that ranges of 58%-62%, 31%-60%, and 48%-59% identity was observed among the SVMP, Snaclec, and SVSP toxin families, respectively. Understanding the venom protein profile and classification is essential in interpreting clinical symptoms during human envenomation and developing potential therapeutic applications. Moreover, the variability of toxin families and amino acid sequences among related hemotoxic snakes found in this study suggests the use and development of universal antivenom for the treatment of envenomating patients is still challenging.