Structures of Azemiops feae venom phospholipases and cys-rich-secretory protein and implications for taxonomy and toxinology

Venom Proteomics of Trimeresurus gracilis, a Taiwan-Endemic Pitviper, and Comparison of Its Venom Proteome and VEGF and CRISP Sequences with Those of the Most Related Species

Trimeresurus gracilis is an endemic alpine pitviper in Taiwan with controversial phylogeny, and its venom proteome remains unknown. In this study, we conducted a proteomic analysis of T. gracilis venom using high-performance liquid chromatography-tandem mass spectrometry and identified 155 toxin proteoforms that belong to 13 viperid venom toxin families. By searching the sequences of trypsin-digested peptides of the separated HPLC fractions against the NCBI database, T. gracilis venom was found to contain 40.3% metalloproteases (SVMPs), 15.3% serine proteases, 6.6% phospholipases A2, 5.0% L-amino acid oxidase, 4.6% Cys-rich secretory proteins (CRISPs), 3.2% disintegrins, 2.9% vascular endothelial growth factors (VEGFs), 1.9% C-type lectin-like proteins, and 20.2% of minor toxins, nontoxins, and unidentified peptides or compounds. Sixteen of these proteoforms matched the toxins whose full amino-acid sequences have been deduced from T. gracilis venom gland cDNA sequences. The hemorrhagic venom of T. gracilis appears to be especially rich in PI-class SVMPs and lacks basic phospholipase A2. We also cloned and sequenced the cDNAs encoding two CRISP and three VEGF variants from T. gracilis venom glands. Sequence alignments and comparison revealed that the PI-SVMP, kallikrein-like proteases, CRISPs, and VEGF-F of T. gracilis and Ovophis okinavensis are structurally most similar, consistent with their close phylogenetic relationship. However, the expression levels of some of their toxins were rather different, possibly due to their distinct ecological and prey conditions.

Differences between Two Groups of Burmese Vipers (Viperidae: Azemiops) in the Proteomic Profiles, Immunoreactivity and Biochemical Functions of Their Venoms

Two recently revised Azemiops snakes with apparent differences in their external appearances and skeletal morphologies but unclear genetic boundaries have been proposed. Some researchers have refrained from using the newly proposed taxonomy because these two “species” might be two clades corresponding to different geographical populations of Azemiops feae. To improve the understanding of the kinship of these two Burmese viper groups, more of their characteristics should be explored in depth. We performed a comparative analysis of the proteomic profiles and biochemical activities of snake venoms from these two groups (Sichuan A. feae and Zhejiang A. feae) and evaluated the immunorecognition capacity of commercial antivenoms toward them. Eight protein families were identified in venoms from these two groups, while phospholipase B was only detected in venom from Sichuan A. feae. These protein families displayed varying degrees of differences in relative abundance between venoms, and phospholipase A₂ (Sichuan A. feae: 57.15%; Zhejiang A. feae: 65.94%) was the predominated component. Gloydius brevicaudus antivenom exhibited the strongest capacity to immunologically recognize these two venoms, but this was mainly limited to components with high molecular masses, some of which differed between venoms. Additionally, Zhejiang A. feae venom was more toxic than Sichuan A. feae venom, and the venoms expressed remarkable differences in enzymatic activities, probably resulting from the variation in the relative abundance of specific protein families. Our findings unveil differences between the two Burmese viper groups in terms of proteomic profiles, immunoreactivity, and the biochemical functions of their venoms. This information will facilitate the management of snakebites caused by these snakes.

De Novo Genome Assembly Highlights the Role of Lineage-Specific Duplications in the Evolution of Venom in Fea’s Viper

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A₂ (PLA₂) proteins expressed by genes that likely arose from lineage-specific PLA₂ gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.

Structural and bioinformatic analyses of Azemiops venom serine proteases reveal close phylogeographic relationships to pitvipers from eastern China and the New World

The semi-fossil and pit-less Azemiops feae is possibly the most primitive crotalid species. Here, we have cloned and sequenced cDNAs encoding four serine proteases (vSPs) from the venom glands of Chinese A. feae. Full amino-acid sequences of the major vSP (designated as AzKNa) and three minor vSPs (designated as AzKNb, AzKNc and Az-PA) were deduced. Using Protein-BLAST search, the ten most-similar vSPs for each Azemiops vSP have been selected for multiple sequence alignment, and all the homologs are crotalid vSPs. The results suggest that the A. feae vSPs are structurally most like those of eastern-Chinese Gloydius, Viridovipera, Protobothrops and North American pitvipers, and quite different from more-specialized vSPs such as Agkistrodon venom Protein-C activators. The vSPs from Chinese A. feae and those from Vietnamese A. feae show significant sequence variations. AzKNa is acidic and contains six potential N-glycosylation sites and its surface-charge distribution differs greatly from that of AzKNb, as revealed by 3D-modeling. AzKNb and AzKNc do not contain N-glycosylation sites although most of their close homologs contain one or two. Az-PA belongs to the plasminogen-activator subtype with a conserved N²⁰-glycosylation site. The evolution of this subtype of vSPs in Azemiops and related pitvipers has been traced by phylogenetic analysis.

Novel Bradykinin-Potentiating Peptides and Three-Finger Toxins from Viper Venom: Combined NGS Venom Gland Transcriptomics and Quantitative Venom Proteomics of the Azemiops feae Viper

Feae's viper Azemipos feae belongs to the Azemiopinae subfamily of the Viperidae family. The effects of Viperidae venoms are mostly coagulopathic with limited neurotoxicity manifested by phospholipases A2. From A. feae venom, we have earlier isolated azemiopsin, a novel neurotoxin inhibiting the nicotinic acetylcholine receptor. To characterize other A. feae toxins, we applied label-free quantitative proteomics, which revealed 120 unique proteins, the most abundant being serine proteinases and phospholipases A2. In total, toxins representing 14 families were identified, among which bradykinin-potentiating peptides with unique amino acid sequences possessed biological activity in vivo. The proteomic analysis revealed also basal (commonly known as non-conventional) three-finger toxins belonging to the group of those possessing neurotoxic activity. This is the first indication of the presence of three-finger neurotoxins in viper venom. In parallel, the transcriptomic analysis of venom gland performed by Illumina next-generation sequencing further revealed 206 putative venom transcripts. Together, the study unveiled the venom proteome and venom gland transciptome of A. feae, which in general resemble those of other snakes from the Viperidae family. However, new toxins not found earlier in viper venom and including three-finger toxins and unusual bradykinin-potentiating peptides were discovered.

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

Two acidic Asp⁴⁹-PLA₂s with Glu⁶ substitution and a neutral Lys⁴⁹-PLA (designated Gst-K49) were cloned from G. strauchii venom glands, their full amino acid sequences were deduced. The predominant acidic PLA₂ (designated Gst-E6a) contains 124 residues and the M¹⁸W³⁰ substitutions, while the minor acidic PLA₂ (designated Gst-E6b) contains 122 residues and the V¹⁸A³⁰ substitutions. Their sequences are most similar to those of the respective orthologous PLA₂s 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 PLA₂ 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 Asp⁴⁹-PLA₂s, but Gst-K49 is the first Lys⁴⁹-PLA₂ identified in Gloydius venoms. However, its venom content is relatively low and its pI value 7.3 is much lower than those of other Lys⁴⁹-PLA₂s and. The Lys⁴⁹-PLA₂ 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.