Distinct regulatory networks control toxin gene expression in elapid and viperid snakes.

BACKGROUND: Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris). RESULTS: Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis. CONCLUSIONS: Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Proteomic comparisons of venoms of long-term captive and recently wild-caught Eastern brown snakes (Pseudonaja textilis) indicate venom does not change due to captivity.

UNLABELLED: Snake venom is a highly variable phenotypic character, and its variation and rapid evolution are important because of human health implications. Because much snake antivenom is produced from captive animals, understanding the effects of captivity on venom composition is important. Here, we have evaluated toxin profiles from six long-term (LT) captive and six recently wild-caught (RC) eastern brown snakes, Pseudonaja textilis, utilizing gel electrophoresis, HPLC-MS, and shotgun proteomics. We identified proteins belonging to the three-finger toxins, group C prothrombin activators, Kunitz-type serine protease inhibitors, and phospholipases A2, among others. Although crude venom HPLC analysis showed LT snakes to be higher in some small molecular weight toxins, presence/absence patterns showed no correlation with time in captivity. Shotgun proteomics indicated the presence of similar toxin families among individuals but with variation in protein species. Although no venom sample contained all the phospholipase A2 subunits that form the textilotoxin, all did contain both prothrombin activator subunits. This study indicates that captivity has limited effects on venom composition, that venom variation is high, and that venom composition may be correlated to geographic distribution. BIOLOGICAL SIGNIFICANCE: Through proteomic comparisons, we show that protein variation within LT and RC groups of snakes (Pseudonaja textilis) is high, thereby resulting in no discernible differences in venom composition between groups. We utilize complementary techniques to characterize the venom proteomes of 12 individual snakes from our study area, and indicate that individuals captured close to one another have more similar venom gel electrophoresis patterns than those captured at more distant locations. These data are important for understanding natural variation in and potential effects of captivity on venom composition.