Next Article in Journal
Aflatoxin B1 Degradation by Ery4 Laccase: From In Vitro to Contaminated Corn
Next Article in Special Issue
De Novo Venom Gland Transcriptome Assembly and Characterization for Calloselasma rhodostoma (Kuhl, 1824), the Malayan Pit Viper from Malaysia: Unravelling Toxin Gene Diversity in a Medically Important Basal Crotaline
Previous Article in Journal
Biocontrol Capabilities of Bacillus subtilis E11 against Aspergillus flavus In Vitro and for Dried Red Chili (Capsicum annuum L.)
Previous Article in Special Issue
Anti-Metalloproteases: Production and Characterization of Polyclonal IgG Anti-F2 Fraction Antibodies Purified from the Venom of the Snake Bitis arietans
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 15
Issue 5
10.3390/toxins15050309
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring Toxin Genes of Myanmar Russell’s Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics

by
Khin Than Yee
1,
Jason Macrander
2,
Olga Vasieva
3,4 and
Ponlapat Rojnuckarin
5,*
1
Department of Medical Research, Ministry of Health, Yangon 11191, Myanmar
2
Department of Biology, Florida Southern College, Lakeland, FL 33801, USA
3
Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, UK
4
BioSynthetic Machines, Inc., Chicago, IL 60062, USA
5
Excellence Center in Translational Hematology, Division of Hematology, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
*
Author to whom correspondence should be addressed.
Toxins 2023, 15(5), 309; https://doi.org/10.3390/toxins15050309
Submission received: 24 January 2023 / Revised: 3 April 2023 / Accepted: 20 April 2023 / Published: 26 April 2023
(This article belongs to the Special Issue Recent Updates in Venomics and Applications)
Browse Figures
Review Reports Versions Notes

Abstract

:
The Russell’s viper (Daboia siamensis) is a medically important venomous snake in Myanmar. Next-generation sequencing (NGS) shows potential to investigate the venom complexity, giving deeper insights into snakebite pathogenesis and possible drug discoveries. mRNA from venom gland tissue was extracted and sequenced on the Illumina HiSeq platform and de novo assembled by Trinity. The candidate toxin genes were identified via the Venomix pipeline. Protein sequences of identified toxin candidates were compared with the previously described venom proteins using Clustal Omega to assess the positional homology among candidates. Candidate venom transcripts were classified into 23 toxin gene families including 53 unique full-length transcripts. C-type lectins (CTLs) were the most highly expressed, followed by Kunitz-type serine protease inhibitors, disintegrins and Bradykinin potentiating peptide/C-type natriuretic peptide (BPP-CNP) precursors. Phospholipase A2, snake venom serine proteases, metalloproteinases, vascular endothelial growth factors, L-amino acid oxidases and cysteine-rich secretory proteins were under-represented within the transcriptomes. Several isoforms of transcripts which had not been previously reported in this species were discovered and described. Myanmar Russell’s viper venom glands displayed unique sex-specific transcriptome profiles which were correlated with clinical manifestation of envenoming. Our results show that NGS is a useful tool to comprehensively examine understudied venomous snakes.
Keywords:
venom gland transcriptomes; Russell’s viper; Daboia siamensis; next-generation sequencing; toxin genes; Venomix pipeline
Key Contribution: Myanmar Russell’s viper venom glands showed unique sex-specific transcriptomes. Sequence data are correlated with clinical envenoming and being a reference for future proteomic and toxicogenomic research.

Graphical Abstract

1. Introduction

Russell’s viper is a widespread species found in South and Southeast Asia, as well as China. Based on previous molecular [1] and morphological [2] studies, the Daboia genus has been divided into two well-defined taxa: D. russelii (west of the Bengal Bay) and D. siamensis (east of the Bengal Bay). The species designation and regional distribution align with variations observed in the clinical features of Russell’s viper envenomation [2,3]. The venom proteomes of Russell’s viper vary across their geographical distribution and have been attributed to these differences [4,5,6,7,8,9].
Snake venoms are mixtures of proteins working synergistically to produce various toxic functions. These compounds evolved to immobilize, kill and initially digest their preys [10,11]. The myriad toxin compositions and functions are a byproduct of their interrelation between natural selection and biochemical processes shaping the protein–protein or protein–cellular interactions. Venom variation is found not only across lineages, but also within individuals of the same species [12]. Additional factors influencing venom variation include age, sex, diet, season, frequency of extraction and a post-extraction process [13,14,15]. The ontogenetic change and venom gland location can also influence the venom composition in an individual snake [16].
Transcriptomic and proteomic approaches have been used to better understand venom diversity [17,18,19]. Transcriptomic studies using next-generation sequencing (NGS) techniques deliver a much broader view of the venom protein coding genes and a deeper understanding of toxin transcript abundance. Transcriptomic approaches permit the identification of unique isoforms that are candidates for drug discovery [20]. Using these techniques, multiple toxin gene copies have been identified within comprehensive datasets, aiding our understanding of how venom protein gene families have evolved. Each toxin family shows a unique evolutionary rate, which can be linked to a change in ecological opportunities, such as predation [21].
Proteomic approaches are often used to provide an expanded overview of venom composition and expression. Previous proteomic studies of Russell’s viper from different countries have shown that toxin gene families exhibit varying venom compositions throughout their distribution [4,6,7,8,22,23].
Despite the differences across their population observed in proteomic data, the Russell’s viper protein-coding toxin genes have not been previously investigated through a comprehensive transcriptomic analysis, nor have any Russell’s viper studies focused on sex-specific venom expression despite anecdotal evidence indicating size- or sex-specific envenomation symptoms in humans. The present study aims to generate a sex-specific de novo venom-gland transcriptomes of Myanmar Russell’s viper (D. siamensis) and provide a deeper insight into the correlation between the toxin composition and pathophysiology of Russell’s viper envenomation. The results of the current study can be used to validate diverse predicted toxins as well as identify novel toxin families in this species that are yet to be identified.

2. Results

2.1. De Novo Transcriptome Assembly

Sequencing of the cDNA libraries yielded a total of 78.64 and 74.44 million clean reads for male and female Myanmar Russell’s vipers, respectively (Table 1). De novo assembly using the Trinity [24] program created 88,325 contigs (N50 = 1870) comprising 77,668 Unigenes (N50 = 1397) for the male snakes and 50,858 contigs (N50 = 1379) comprising 46,884 Unigenes (N50 = 1183) for the female snakes.

2.2. Categorization of Transcripts and Gene Expression

The expression levels of toxin transcripts from both samples were compared with those of non-toxin transcripts from tissues of the same species (Bioproject PRJNA269317, Biosamples SAMN03081270) (Supplementary Tables S1 and S2).
The toxin transcripts were confirmed to be highly expressed in venom gland tissues. They were categorized according to TPM values. The individual toxin transcripts with a TPM value of more than 100 were categorized as highly expressed, a TPM value of 10 to 100 as moderately expressed and a TPM value of less than 10 as lowly expressed toxin groups. Only highly expressed groups were considered for further detailed evaluation and description in this study. Non-toxin full-length transcripts, such as glutaminyl-peptide cyclotransferase (highly expressed), tomoregulin-2 (moderately expressed) and endothelin-converting enzyme (lowly expressed), were also identified in the current transcriptome. Highly expressed transcripts were no less than 150-fold higher in the venom gland when compared to samples comprising other tissues (brain, kidney, heart, lung, spleen and liver).
The highly expressed venom transcripts were categorized into 23 toxin gene families, including 53 unique full-length sequences. The most prevalent transcripts were C-type lectins (CTLs) (26%), Kunitz-type protease inhibitors (KSPIs) (19%), disintegrins (18%) and bradykinin potentiating peptide/C-type natriuretic precursors (BPP-CNPs) (16%). The remaining 21% consisted of phospholipase A2 (PLA2), serine proteases, metalloproteinases, vascular endothelial growth factor (VEGF), L-amino acid oxidase (LAAO) and cysteine-rich secretory proteins (CRISPs) (Figure 1).
When comparing males and females as shown in Figure 2, the male sample had CTLs, KSPIs, BPP-CNPs and disintegrins as the most highly expressed toxin groups, while metalloproteinases, CTLs, BPP-CNPs and KSPIs were the top highly expressed toxin groups in female snakes. The expression levels of transcripts for most male toxin groups were higher than those in the female-mapped reads, with 1.39- to 11.68-fold differences. However, the number of transcripts for serine proteases, glutaminyl-peptide cyclotransferase-1 and metalloproteinases, were 1.16, 1.22 and 2.33 times higher in female than male venom glands, respectively. Some toxin transcripts, such as cathelicidin-related protein (Vipericidin) and techylectin-like protein (Veficolin), were only identified in male samples.

2.3. Snake Venom C-Type Lectins (CTLs) and Related Proteins

The present study revealed 23 CTL transcripts (4 full-length and 19 partial-length transcripts), including 7 different subtypes of CTLs. They are (I) newly identified Bitis gabonica C-type lectin 1 homolog, (II) D. siamensis p31 β subunit, (III) RVV-X LC1, (IV) D. siamensis p68α subunit and Microvipera lebetina Snaclec A15 homolog, (V) Dabocetin α and Microvipera lebetina Snaclec A12 homolog, (VI) RVV-X LC2 and (VII) D. siamensis p68 β subunits and D. siamensis Snaclec-7 (Supplementary Table S3). Figure 3 shows the gene tree of these sequences. RVV-X CL1 and RVV-X CL2 linked to the metalloproteinase (RVV-X) to act on factor X of the coagulation cascade.
Our sequence alignment (Supplementary Figure S1) identified that these sequences contained six conserved cysteine residues (typically eight in total for each sequence), along with a WIGL conserved motif for CTL and the WND Gla-domain binding site.

2.4. Kunitz-Type Serine Protease Inhibitors

A total of 15 Kunitz-type serine protease inhibitor (KSPI) transcripts were identified, including seven full-length and eight partial-length transcripts. Thirteen unique isoforms of KSPIs were observed (Table 2).

2.5. RTS-Disintegrins

A total of three full-length disintegrin transcripts were identified, including two isoforms of RTS-containing disintegrins which were also described as Dis1a and Dis1b in our previous paper [25]. In the gene tree analysis, Dis1a shared a high sequence similarity with jerdostatin from Protobothrops jerdonii (Q7ZZM2), while Dis1b was in the same branch as acostatin-α from D. siamensis (AUF416558) (Supplementary Figure S2). In the sequence alignment, the N-terminal sequence of Dis1b was more similar to that of ascostatin-α than jerdostatin and it possessed an RTS motif instead of an RGD motif (Figure 4).

2.6. Bradykinin-Potentiating Peptides and C-Type Natriuteric Peptide (BPP-CNP) Precursors

A total of four full-length BPP-CNP transcripts were identified in this study. In Figure 5, the repeated BPP units were present in a BPP-CNP precursor. These BPP units appeared in tripeptide metalloproteinase inhibitors. CNPs were found to be variable between these sequences. The sequence alignments of the natriuretic peptide domains of Myanmar Russell’s viper BPP-CNP precursors with those of human and other snakes are shown in Figure 6.

2.7. Phospholipase A2 (PLA2)

In the focal transcriptomes, there were five full-length and two partial-length PLA2 candidates. Five full-length transcripts contain three group II PLA2 and two group III PLA2 transcripts. Three of the group II PLA2s (acidic and basic PLA2) were highly expressed (TPM >10,000–14,000), followed by two group III PLA2 transcripts (TPM 14–35), whereas two partial-length transcripts, group IIF PLA2 homologs of Crotralus atrox and a crotoxin homolog from Protobothrops guttatus, were the lowest expressed transcripts (TPM 1–8). The highly expressed group II PLA2 transcripts were aligned with PLA2 from Eastern Russell’s viper PLA2s, DsM b1and daboiatoxin A from Myanmar, PLA2-I from China and RV-7 from Taiwan (Figure 7).

2.8. Snake Venom Serine Proteases (SVSPs)

There were 14 transcripts (2 full-length and 12 partial-length) annotated as serine proteinases that were highly expressed (TPM 448.14–5495.65). Four transcripts gave an annotation as RVV-Vγ of D. siamensis (Supplementary Figure S3), three transcripts as α-fibrinogenases and one transcript as β-fibrinogenase of D. siamensis. Four transcripts were homologous to VLSP-1 (Figure 8) and one transcript to VLSP-3 from Macrovipera lebetina. The remaining one was partially similar to VSPBF-DSABSI from D. siamensis and partially to VLSP-3 from Macrovipera lebetina. We could not find out the exact similarity due to the partial-length nature of the remaining transcripts (Table 3).

2.9. Snake Venom Metalloproteinases (SVMPs)

In the current transcriptomes, two full-length transcripts for RVV-X heavy chain (Russell’s viper factor X activator, A0A2H4Z2W1) and two partial-length DSAIP (Daboia siamensis VLAIP-like SVMP, A0A2H4Z2Y9) are the most highly expressed SVMP transcripts (>1.0–1.8 × 104 TPM). Of the remaining three full-length transcripts with TPM values of 15–47, c15329_g2_i1_M_FL and c10072_g2_i1_F_FL were homologous to the A disintegrin and metalloproteinase domain (ADAM)-containing protein 10 from Vipera anatolica senliki (A0A6G5ZVR9) (Figure 9), and c19271_g2_i1_M_FL was homologous to SVMP from Agkistrodon contortrix contortrix (A0A1W7RJN7) (Figure 10). All conceptually translated proteins contained the signal peptide, a metalloproteinase domain containing a Zn+-binding catalytic site and a disintegrin domain. There was an additional ADAM17 MPD domain containing a CxxC motif in the c19271_g2_i1_M_FL transcript.

2.10. Snake Venom Vascular Endothelial Growth Factors (svVEGFs)

In the current study, seven full-length VEGF transcripts were identified including three isoforms which comprised 136, 148 and 192 amino acids. Two transcripts were homologous to VR-1 (D. russelii) (P67861), as shown in Supplementary Figure S4. Two transcripts were related to VEGF1 (Vipera anatolica senliki) (QHR82857) and VEGFA-VIPAA (Vipera ammodytes ammodytes) (C0K3N5), as shown in Figure 11A,B.

2.11. L-Amino Acid Oxidases (LAAOs)

In this study, two full-length and one partial-length LAAO transcripts with 99.93% identity (84% query coverage, E value 0.0) to DrLAO (G8XQX1) were identified. The sequence alignment is shown in Supplementary Figure S5.

2.12. Cysteine-Rich Venom Proteins (CRISPs)

There were two full-length highly expressed cysteine-rich secretory protein (CRISP) transcripts identified in the current transcriptomes which gave the annotation of Dr-CRPK from D. russelii (ACE73567), currently termed D. siamensis (Supplementary Figure S6). These deduced amino acid sequences were also identical to those of CRISPs (MRV CRISP) (JZ980958, JZ980965 and JZ980970), which were expressed in previous cDNA libraries [26].

3. Discussion

The de novo venom gland transcriptomic analysis provides a unique profile of the toxin composition of Myanmar D. siamensis venom and its suggestive toxicogenomic diversity. Our previous study on the venom gland cDNA library revealed unique gene expression patterns [26]; however, the use of NGS technology allowed us to provide a broader examination into the variation, distribution and diversity of toxin transcripts. Using Trinity as our de novo transcriptome assembly program, we were able to assemble a robust transcriptome using high-quality (>98% Q20 for both male and female) raw reads and subsequent mapping for gene expression analysis. Our initial exploration into these transcriptomes relied on RSEM and gene expression information to identify focal toxin candidate transcripts; however, future studies would benefit from exploring additional assembly programs providing a multi-faceted approach to analyzing these transcriptomes and account for chimeric transcripts [27]. The venom proteome of D. siamensis comprised six toxin families, including serine proteinases, metalloproteinases, PLA2, LAAO, VEGF and CTLs [4]. The BPP-CNP transcripts containing SVMPIs were identified only in the transcriptomes due to the co-expression of SVMPI tripeptides and atrial natriuretic peptides in the same transcripts [28]. Only the processed products could be detected in the venom. Low-molecular-weight disintegrins (<10 kDa) were detected only in a transcriptome because a proteomic study could identify proteins larger than 10 kDa. The reason that CRISPs (27 kDa) were not found in the proteome was unclear.
We found sex-related differences in venom transcript expression in Russell’s viper, also detected in Bothrops jararaca [12], Thamnophis elegans [17] and Mesobuthus martensii [29]. These differences reflected the sex-based variations of the venom proteome and were closely associated with the biochemical and biological properties of different genders in the same species [14,30,31,32]. In B. jararaca and B. moojeni, the variability in the venom composition was associated with animal sizes, as females are larger than males [32,33]. This coincided with the sexual dimorphism of the venom donors used in the study, as females (average weight in ounces × total length in inches x girth in inches: 49.9 × 47.8 × 6.5) were larger than males (29.8 × 47 × 5.3) and more venom was collected from females. There is also a correlation between venom yield and the length of Myanmar Russell’s vipers, which is related to the severity of envenoming, time for coagulopathy development and the severity of local manifestations [34].
The transcripts corresponding to the toxin groups described below were found to be most abundant in both sexes and aligned with Myanmar Russell’s viper envenomation symptoms. Snake venom C-type lectins (snaclecs) are heterodimers of two closely related α and β subunits [35]. They target factor IX, factor X, prothrombin or α-thrombin (as pro- or anticoagulant) and platelet receptors, such as glycoprotein (Gp) Ibα, Gp Ia/IIa, or Gp VI (as pro- or antithrombotic) [36,37]. Snaclecs activate platelets and thrombus formation causing thrombo-inflammation during viper envenoming [38]. Dabocetin α together with its β subunit interacts with platelet GpIb receptor, inhibiting platelet aggregation [39]. D. siamensis Snaclec-7 is a β subunit of a GpIb-binding protein [40], and together with dabocetin α is involved in platelet inhibition. The M. lebetina Snaclec A12 homolog, originally found in D. siamensis, has a stronger sequence similarity to α than β subunits compared to other heterodimeric snaclecs [41].
Snake venom protease inhibitors primarily disrupt blood coagulation resulting in blood loss and the death of preys [42]. The largest and most widely distributed family is Kunitz-type protein inhibitors, which strongly inhibit serine proteases. They contain 58–60 amino acid residues, including three disulfide bridges [43]. Several Kunitz-type protease inhibitors from Daboia species have previously been identified.
The snake venom disintegrins are a group of small cysteine-rich proteins found in a variety of snake families. They are divided into seven groups according to disulfide patterns and sequence signatures [44]. Jerdostatin homologs (Dis1a and Dis1b) containing eight cysteines belong to group 7A. Most of these small disintegrins are synthesized from short-coding genes, instead of a proteolytic process from PII SVMPs like the majority of disintegrins [45]. The RTS-disintegrins are expressed from the RPTLN gene and have a high conserved signal peptide sequence with SVMP genes, suggesting they may share tissue-specific regulatory elements [46,47]. The disintegrin motif is responsible for binding specificity to integrin receptors [48]. Modulating disintegrin motifs is interesting for the development of antiangiogenesis, antimetastasis, antiproliferative and proapoptotic agents for cancer treatments. Lowly expressed VGD- and MLD-containing dimeric disintegrins observed in our previous study [25] were not evaluated in the current study as the candidate transcripts had TPM values lower than 100 and were therefore not analyzed.
The BPP-CNP transcripts usually carry a tripeptide repeat together with bradykinin-potentiating peptides (BPPs) and natriuretic peptide (NP) precursors [28,49]. The two tripeptides (pERW and pEKW) from Myanmar Russell’s viper completely inhibited the gelatinolytic effects of RVV-X and the fibrinogenolytic activities of daborhagin [28]. These viper-restricted inhibitory peptides protected the venom gland from autodigestion and prevented hydrolysis of other proteins in venom glands [50]. Their inhibitory activity towards the SVMP Zn2+ binding site [51] has made them popular targets for drug design in areas of neurodegenerative and cardiovascular diseases. Natriuretic peptides (NPs) in vertebrates promote salt secretion, regulating cardiovascular and renal systems. Three mammalian NPs are atrial natriuretic peptide (ANP), B-type natriuretic peptide (BNP) and C-type natriuretic peptide (CNP). ANP and BNP are released by cardiomyocytes in response to hypertension and hypervolemia, and CNP is produced by endothelial cells. The sequence alignment of various snake venom and human natriuretic peptides is shown in Figure 6, where five residues (F, D, R, I/V and L/F) are crucial for receptor binding [50]. Among MRV CNP transcripts, one had the same conserved amino acid residues, while remaining three possessed V instead of I and F instead of L residues. These variations should be explored in relation to their roles in terms of toxicity and their potential therapeutic properties. NPs have their specific niche in medical research in the treatment of heart failure and hypertension [51,52].
Phospholipases A2 (PLA2) catalyze a hydrolysis of glycerophospholipids, fatty acids and lysophospholipids and release arachidonic acids to trigger inflammation. Snake venom PLA2s (svPLA2) are small, calcium-dependent, secreted enzymes ranging from 12 to 14 kDa. The class I PLA2s are found in Elapidae and Hydrophiidae, while class II are in Crotalidae and class III are in lizard and bee venoms [53]. There are two major types in Russell’s vipers. PLA2s containing an N-terminal Asn (N) are found in species from China, Thailand, Myanmar and Pakistan, while the venoms of southern Indian and Sri Lankan species contain PLA2 with an N-terminal serine (S). They confer hypotensive with neurotoxic and myonecrotic properties, respectively [26]. Furthermore, Russell’s viper PLA2 showed nephrotoxicity in mice [54]. All of the three PLA2 transcripts from the current transcriptomes contained N-terminal Asparagine (N) and conserved Aspartic acid (D) at the 49 position. In all transcripts, 14 cysteine residues for seven conserved disulfide bonds indicated that they were GIIA PLA2s. They all possessed Thr41 which was found to be conserved in all the neurotoxic or myotoxic PLA2s. From Myanmar D. siamensis venom, DsM b1 had neurotoxic and lethal effects [55] and dabotoxin showed neurotoxic, myotoxic, cytotoxic, edema-inducing and indirect hemolytic activities [56].
Snake venom serine proteases (SVSPs) have two six-stranded β-barrels and a catalytic triad of His57, Asp102 and Ser159. Viper venoms contain a large variety of SVSPs, mostly as monomeric glycoproteins of 26 to 67 kDa. They act on almost all components of hemostatic and fibrinolytic systems and are responsible for the hemorrhage and shock after envenomation. Based on their biological roles, they are classified as activators of the fibrinolytic, procoagulant, anticoagulant and platelet-aggregating enzymes [57,58,59]. SVSPs from D. russelii are procoagulants that (i) activate factor V (RVV-V), (ii) have thrombin-like action causing defibrinogenation (Russelobin) and (iii) are αβ fibrinogenases with factor V-activation (RV-FVPα-δ). RV-FVPs showed in vivo anticoagulant effects in addition to in vitro procoagulant effects [60].
Snake venom metalloproteinases (SVMPs) are zinc endopeptidases varying from 20 to 100 kDa. They are the major components of viper venom, varying from 11% to 65% of the whole venom cocktail [61,62]. SVMPs have sequence and domain similarities to ADAM (A disintegrin and metalloproteinase) and can be grouped into three major classes (PI-PIII) according to their organization [53,63]. SVMPs are responsible for local and systemic hemorrhage, myonecrosis, the disruption of hemostasis mediated by procoagulant or anticoagulant effects or platelet aggregation, tissue necrosis and pro-inflammation [61]. These activities are mediated by metalloproteinase, disintegrin and/or disintegrin-like domains [62]. Snake venom metalloproteinases (SVMPs) isolated from Russell’s viper were RVV-X (85 kDa), Daborhagin (65 kDa), VRR-73 (73 kDa), VRH-1 (23 kDa) and RVBCMP (25 kDa). RVV-X is a well-characterized class P-IIId SVMP. It activates factor X via a single cleavage at the Arg194-Ile195 bond in factor X, which is responsible for disseminated intravascular coagulation (DIC) in snakebite patients [64]. RVV-X also binds to receptors on platelets and endothelial cells, resulting in clinical bleeding [60,65] and contributing to acute kidney injury in animal models [66]. The structure and biological activity of DSAIP still needs to be elucidated. RVV-X and daborhagin are found in both Myanmar and Indian RV species, while VRR-73, VRH-1 and RVBCMP are from only Indian RV [61]. Daborhagin-M is a P-IIIc hemorrhagic SVMP with high proteolytic activity toward fibrinogen, fibronectin and collagen [67]. These activities may be related to the severe bleeding in the patients.
A disintegrin and metalloproteinase (ADAM) family proteins in focal transcriptomes were expressed at moderate levels. Transcripts homologous to VIPAN (ADAM 10) from V. a. senliki and AGKCO (ADAM 17) from A. c. contortrix were found in Myanmar Russell’s viper. ADAM 10 played a major role in the signaling pathway and ADAM 17 in the processing of tumor necrosis factor-α and other cell-surface molecules [68]. Although ADAM 17 is considered to be the primary sheddase for TNF-α, ADAM 10 also worked as a TNF-α sheddase in certain cell types [69]. The biological functions of ADAM 10 and 17 from snakes require further investigation.
Snake venom vascular endothelial growth factors (svVEGFs) play important roles in angiogenesis (endothelial growth and proliferation), hypotension and vascular permeability by binding to VEGF receptors. svVEGFs are 25-kDa, homodimeric heparin-binding proteins with 50% identity to VEGF-A165. They have a VEGF homology domain (VHD), with eight conserved cysteines forming a cysteine knot. The C-terminal variable region corresponds to heparin and NP-1 (coreceptor) binding sites [70,71]. VR-1, a svVEGF from D. russelii, increased endothelial cell proliferation, vascular permeability and hypotension through the specific activation of endothelial KDR (kinase insert domain-containing receptor) [72]. The VR-1 isolated from D. siamensis also increased endothelial cell proliferation [73]. The markedly enhanced vascular permeability caused by svVEGFs can lead to hypotension and shock.

4. Conclusions

The de novo venom gland transcriptomic analysis provided a unique profile of venom gene expression in D. siamensis from Myanmar. The expression patterns of toxin-encoding genes showed high structural redundancy in comparison to non-toxin genes, suggesting distinct functions of the variants in each group. The expression patterns of venom genes were also sex-specific. The distribution and expression levels of the principal toxin components yielded deeper insights into the toxic syndromes of Russell’s viper envenoming, which included bleeding, hypotension and renal failure. Moreover, the full-length toxin sequences from this study would be the references for future proteomic analysis. Newly identified venom components will be a source for future structure–functional relationship studies and the development of new therapeutic modalities. Taking into account that the genomes of these species were not yet completed, this study provides rich data for a fast advancement in toxicogenomic research to study evolutionary origins and diversity in snakes.

5. Materials and Methods

5.1. Sample Preparation

Four (two male and two female) adult Russell’s viper, D. siamensis, were captured from Kyun Gyan Gon Township, Yangon, Myanmar. Venom milking was carried out prior to tissue harvesting to stimulate the toxin transcription. The snakes were allowed to rest for four days to maximize RNA transcription. The venom glands were then removed, cut into pieces smaller than 5 × 5 mm and stored in an RNAlater solution (Ambion, Austin, TX, USA) at 4 °C overnight before transferring to −80 °C.

5.2. RNA Extraction, cDNA Library Construction and Sequencing

Total RNA from venom gland tissues was extracted using Trizol reagent (Life Technologies, Carlsbad, CA, USA) or a Total RNA Purification Kit (Jena Bioscience GmbH, Jena, Germany). The sample were stored at −80 °C. Subsequently, mRNA from samples was isolated using PolyAT Tract mRNA Isolation System (Promega, Madison, WI, USA) or a FastTrack MAG mRNA Isolation Kit (Invitrogen, Carlsbad, CA, USA). Both systems used MagneSphere® technology. Subsequently, mRNA was precipitated, concentrated using 3M sodium acetate and isopropanol and stored at −80 °C. Isolated mRNA was suspended in RNase-free water and sent to Macrogen Inc. Geumchun-gu, Seoul, Republic of Korea. The mRNA sequencing libraries were prepared using the TruSeq RNA sample preparation kit (Illumina Inc., San Diego, CA, USA) with selected insert sizes of 200–400 bp. Libraries were then sequenced on the Illumina HiSeq2000 platform using (2 × 100 bp) paired-end reads.

5.3. De Novo Transcriptome Assembly and Expression Quantification

The quality of raw data or raw reads from the Illumina platform was assessed with the quality assessment software FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc; accessed on 18 November 2014) and cleaned up using Trimmomatric (0.32) (http://www.usadellab.org/cms/?page=trimmomatic; accessed on 20 November 2014) to remove adaptor sequences, nucleotide errors and low-quality sequences. Two assembly tools, Trinity [24] and CLC Genomic Workbench 7.5 (CLC bio, a QIAGEN Company, Aarhus, Denmark), were used to assemble and compare alternative pathways for de novo transcriptome assemblies. The quality assessment tool QUAST was used to evaluate the quality of the assembled transcriptomes [74], indicating that the Trinity assembly had a higher quality score than the CLC genomic assembly (Supplementary Table S4). Thus, assembled data via Trinity were used for downstream analysis.
The transcriptomes were assembled in the Trinity software (r20140717) [24]. For each transcriptome, expression values were determined as transcripts per million (TPM) and fragments per kilobase million (FKPM) were calculated in the program RSEM [75,76] as part of the Trinity package. The raw data of the assembly were uploaded to NCBI database under BioProject PRJNA545823, Biosample SAMN11938797 for males and SAMN11939170 for females. The expression levels in terms of both TPM and FKPM of toxin transcripts were mapped to transcripts from non-venom gland tissues, such as brain, kidney, heart, lung, liver and spleen (Bioproject PRJNA269317, Biosamples SAMN03081270), of D. siamensis.

5.4. Venom Gland Transcript Classification

Toxin transcripts from the two transcriptomic datasets were identified and characterized using Venomix [77]. The Venomix pipeline identified candidate toxin genes based on BLAST hits (E-value 1 × 10−20) against the ToxProt dataset [78]. A transcript was considered a “candidate” if the transcript had a lower E-value associated with a toxin than with a non-toxin protein in Uniprot. The candidate transcripts were translated into their protein sequences using Transdecoder [79] and further evaluated in ToxClassifier [80]. A predicted protein sequence with a score of >1 was considered a toxin candidate. The annotation of generated transcript sequences was aligned by standalone BLASTP (v. 2.10.0) to the protein database UniProtKB/Swiss-Prot [81] (Swissprot, Swiss Institute of Bioinformatics, Lausanne, Switzerland) of Serpentes (taxid: 8570) organism. Sequences with no matches were then searched against non-redundant protein sequences (nr) of the same organism. Proteins with the highest ranks in the BLASTP results were referenced to determine the functional annotations.
The transcripts from this study were named as “c” (component/contig), “g”’ (subgraph/trinity gene) and “i” (path sequence/trinity isoform) followed by numbers. The letter “F” stood for female, “M” for male, “PL” for partial-length and “FL” for full-length transcripts.

5.5. Multiple Sequences Alignment and Gene Tree Reconstruction

The most highly toxin transcripts (TPM values ≥ 100) were conceptually translated into protein sequences and subjected to alignment using Clustal Omega [82]. Alignments of the first-described and previously reported sequences in Russell’s viper are displayed in the manuscript and Supplementary Material, respectively. The conserved domains and functional motifs are identified and referenced in the figures. For some toxin gene families, the alignments were used to reconstruct in the program RAxML [83] using the PROTGAMMABLOSUM62 model with 1000 bootstrap replicates.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/toxins15050309/s1, Figure S1: Clustal Omega sequence alignment of the polypeptides encoded by 23 transcript candidate C-type lectins) from Myanmar D. siamensis; Figure S2: The gene tree analysis of the disintegrins from Myanmar D. siamensis, Dis1a (c13890_g2_i1_M) and Dis1b (c11797_g1_i1_F), with disintegrins from other species; Figure S3: Sequence alignment of RVV-Vγ from Daboia siamensis and deduced amino acid sequences of transcripts from the transcriptomes; Figure S4: Sequence alignment of deduced amino acid sequences of VEGF transcripts with VR-1 from D. russelii (P67861); Figure S5: Sequence alignment of DrLAO (Daboia russelii, G8XQX1.1) with three deduced LAAO transcripts obtained from Myanmar Russell’s viper transcriptome; Figure S6: Sequence alignment of deduced amino acid sequences of CRISP transcripts from the transcriptomes with cysteine-rich secretory protein Dr-CRPK from Taiwan Russell’s viper (ACE73567); Table S1: Comparison of expression levels of toxin transcripts from male sample with non-toxin transcripts from tissues of D. siamensis; Table S2: Comparison of expression levels of toxin transcripts from female sample with non-toxin transcripts from tissues of D. siamensis; Table S3: Transcripts with snake venom C-type lectins (CTLs) annotation; Table S4; QUAST reports of assembled transcriptomes by two different assemblers; Table S5: List of full-length sequences with toxin gene families of Myanmar Russell’s viper transcriptome. Reference [84] are cited in the supplementary materials.

Author Contributions

Conceptualization, P.R. and K.T.Y.; methodology, P.R., K.T.Y. and J.M.; software, J.M. and O.V.; validation, P.R., J.M. and O.V.; formal analysis, K.T.Y., J.M. and O.V.; investigation, K.T.Y.; resources, P.R.; data curation, J.M. and O.V.; writing—original draft preparation, K.T.Y.; writing—review and editing, K.T.Y., P.R., J.M. and O.V.; visualization, P.R., J.M. and O.V.; supervision, P.R. and O.V.; project administration, P.R.; funding acquisition, P.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Thailand Science Research and Innovation Fund Chulalongkorn University, CU_FRB65_hea(17)_024_30_05.

Institutional Review Board Statement

The study was approved by the Animal Care and Use Committee, Chulalongkorn University (CU-ACUC) (Protocol code 17/2552) on 28 July 2015.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data of the transcriptome assembly were uploaded to NCBI database under BioProject PRJNA545823, Biosample SAMN11938797 for males and SAMN11939170 for females.

Acknowledgments

Administrative persons and snake handlers from the Snake Farm, Myanmar Pharmaceutical Factory, Yangon, Myanmar and Lawan Chanhome, Head of Snake Farm, QSMI, Bangkok, Thailand for help in the collection of snakes and removal of venom glands.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Thorpe, R.S.; Pook, C.E.; Malhotra, A. Symptoms of envenoming. Herpetol. J. 2007, 17, 209–218. [Google Scholar]
  2. Wüster, W. The genus Daboia (Serpentes: Viperidae): Russell’s viper. Hamadryad 1998, 23, 33–40. [Google Scholar]
  3. Warrell, D.A. Snake venoms in science and clinical medicine 1. Russell’s viper: Biology, venom and treatment of bites. Trans. R. Soc. Trop. Med. Hyg. 1989, 83, 732–740. [Google Scholar] [CrossRef]
  4. Risch, M.; Georgieva, D.; von Bergen, M.; Jehmlich, N.; Genov, N.; Arni, R.K.; Betzel, C. Snake venomics of the Siamese Russell’s viper (Daboia russelli siamensis)–relation to pharmacological activities. J. Proteom. 2009, 72, 256–269. [Google Scholar] [CrossRef] [PubMed]
  5. Thomas, R.; Alfred, S.; Zaw, A.; Gentili, S.; Chataway, T.; Peh, C.A.; Hurtado, P.; White, J.; Wiese, M. A demonstration of variation in venom composition of Daboia russelii (Russell’s viper), a significantly important snake of Myanmar, by tandem mass spectrometry. Toxicon 2019, 158, S44–S45. [Google Scholar] [CrossRef]
  6. Lingam, T.M.C.; Tan, K.Y.; Tan, C.H. Proteomics and antivenom immunoprofiling of Russell’s viper (Daboia siamensis) venoms from Thailand and Indonesia. J. Venom. Anim. Toxins Incl. Trop. Dis. 2020, 26, 1–17. [Google Scholar] [CrossRef]
  7. Tan, K.Y.; Tan, N.H.; Tan, C.H. Venom proteomics and antivenom neutralization for the Chinese eastern Russell’s viper, Daboia siamensis from Guangxi and Taiwan. Sci. Rep. 2018, 8, 8545. [Google Scholar] [CrossRef]
  8. Pla, D.; Sanz, L.; Quesada-Bernat, S.; Villalta, M.; Baal, J.; Chowdhury, M.A.W.; León, G.; Gutiérrez, J.M.; Kuch, U.; Calvete, J.J. Phylovenomics of Daboia russelii across the Indian subcontinent. Bioactivities and comparative in vivo neutralization and in vitro third-generation antivenomics of antivenoms against venoms from India, Bangladesh and Sri Lanka. J. Proteom. 2019, 207, 103443. [Google Scholar] [CrossRef] [PubMed]
  9. Mukherjee, A.K.; Kalita, B.; Mackessy, S.P. A proteomic analysis of Pakistan Daboia russelii russelii venom and assessment of potency of Indian polyvalent and monovalent antivenom. J. Proteom. 2016, 144, 73–86. [Google Scholar] [CrossRef]
  10. Casewell, N.R.; Huttley, G.A.; Wüster, W. Dynamic evolution of venom proteins in squamate reptiles. Nat. Commun. 2012, 3, 1066. [Google Scholar] [CrossRef] [PubMed]
  11. Calvete, J.J.; Sanz, L.; Angulo, Y.; Lomonte, B.; Gutiérrez, J.M. Venoms, venomics, antivenomics. FEBS Lett. 2009, 583, 1736–1743. [Google Scholar] [CrossRef]
  12. Zelanis, A.; Andrade-Silva, D.; Rocha, M.M.; Furtado, M.F.; Serrano, S.M.T.; Junqueira-de-Azevedo, I.L.M.; Ho, P.L. A transcriptomic view of the proteome variability of newborn and adult Bothrops jararaca snake venoms. PLoS Negl. Trop. Dis. 2012, 6, e1554. [Google Scholar] [CrossRef] [PubMed]
  13. Barlow, A.; Pook, C.E.; Harrison, R.A.; Wüster, W. Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution. Proc. R. Soc. B Biol. Sci. 2009, 276, 2443–2449. [Google Scholar] [CrossRef] [PubMed]
  14. Menezes, M.C.; Furtado, M.F.; Travaglia-Cardoso, S.R.; Camargo, A.C.M.; Serrano, S.M.T. Sex-based individual variation of snake venom proteome among eighteen Bothrops jararaca siblings. Toxicon 2006, 47, 304–312. [Google Scholar] [CrossRef] [PubMed]
  15. Mirtschin, P.J.; Shine, R.; Nias, T.J.; Dunstan, N.L.; Hough, B.J.; Mirtschin, M. Influences on venom yield in Australian tigersnakes (Notechis scutatus) and brownsnakes (Pseudonaja textilis: Elapidae, Serpentes). Toxicon 2002, 40, 1581–1592. [Google Scholar] [CrossRef]
  16. Brahma, R.K.; McCleary, R.J.R.; Kini, R.M.; Doley, R. Venom gland transcriptomics for identifying, cataloging, and characterizing venom proteins in snakes. Toxicon 2015, 93, 1–10. [Google Scholar] [CrossRef]
  17. Schwartz, T.S.; Tae, H.; Yang, Y.; Mockaitis, K.; Van Hemert, J.L.; Proulx, S.R.; Choi, J.H.; Bronikowski, A.M. A garter snake transcriptome: Pyrosequencing, de novo assembly, and sex-specific differences. BMC Genom. 2010, 11, 694. [Google Scholar] [CrossRef]
  18. Cardoso, K.C.; Da Silva, M.J.; Costa, G.G.L.; Torres, T.T.; Del Bem, L.E.V.; Vidal, R.O.; Menossi, M.; Hyslop, S. A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu). BMC Genom. 2010, 11, 605. [Google Scholar] [CrossRef]
  19. Tan, C.H.; Tan, K.Y.; Fung, S.Y.; Tan, N.H. Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah). BMC Genom. 2015, 16, 687. [Google Scholar] [CrossRef]
  20. Ducancel, F.; Durban, J.; Verdenaud, M. Transcriptomics and venomics: Implications for medicinal chemistry. Future Med. Chem. 2014, 6, 1629–1643. [Google Scholar] [CrossRef]
  21. Barua, A.; Mikheyev, A.S. Toxin expression in snake venom evolves rapidly with constant shifts in evolutionary rates. Proc. R. Soc. B Biol. Sci. 2020, 287, 20200613. [Google Scholar] [CrossRef]
  22. Faisal, T.; Tan, K.Y.; Sim, S.M.; Quraishi, N.; Tan, N.H.; Tan, C.H. Proteomics, functional characterizatioin and antivenom neutralization of the venom of Pakistani Russell’s viper (Daboia russelii) from the wild. J. Proteom. 2018, 183, 1–3. [Google Scholar] [CrossRef]
  23. Saethang, T.; Somparn, P.; Payungporn, S.; Sriswasdi, S.; Yee, K.T.; Hodge, K.; Knepper, M.A.; Chanhome, L.; Khow, O.; Chaiyabutr, N.; et al. Identification of Daboia siamensis venom using intergrated multi-omics data. Sci. Rep. 2022, 12, 13140. [Google Scholar] [CrossRef]
  24. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef]
  25. Yee, K.T.; Tongsima, S.; Vasieva, O.; Ngamphiw, C.; Wilantho, A.; Wilkinson, M.C.; Somparn, P.; Pisitkun, T.; Rojnuckarin, P. Analysis of snake venom metalloproteinases from Myanmar Russell’s viper transcriptome. Toxicon 2018, 146, 31–41. [Google Scholar] [CrossRef] [PubMed]
  26. Yee, K.T.; Rojnuckarin, P. Complementary DNA library of Myanmar Russell’s viper (Daboia russelii siamensis) venom gland. Comp. Biochem. Physiol. Part-C Toxicol. Pharmacol. 2020, 227, 108634. [Google Scholar] [CrossRef]
  27. Holding, M.L.; Margres, M.J.; Mason, A.J.; Parkinson, C.L.; Rokyta, D.R. Evaluating the performance of de novo assembly methods for venom-gland transcriptomics. Toxins 2018, 10, 249. [Google Scholar] [CrossRef] [PubMed]
  28. Yee, K.T.; Pitts, M.; Tongyoo, P.; Rojnuckarin, P.; Wilkinson, M.C. Snake venom metalloproteinases and their peptide inhibitors from Myanmar Russell’s viper venom. Toxins 2017, 9, 15. [Google Scholar] [CrossRef]
  29. Yang, Y.; Zhu, Y.; Luo, Y.; Liu, Q.; Hua, X.; Li, J.; Gao, F.; Hofer, J.; Gao, X.; Xiao, L.; et al. Transcriptome analysis of Mesobuthus martensii revealed the differences of their toxins between females and males. Eur. Zool. J. 2022, 89, 1294–1311. [Google Scholar] [CrossRef]
  30. Lourenço, A.; Zorzella Creste, C.F.; Curtolo de Barros, L.; Delazari dos Santos, L.; Pimenta, D.C.; Barraviera, B.; Ferreira, R.S. Individual venom profiling of Crotalus durissus terrificus specimens from a geographically limited region: Crotamine assessment and captivity evaluation on the biological activities. Toxicon 2013, 69, 75–81. [Google Scholar] [CrossRef]
  31. Tan, C.H.; Tan, K.Y.; Yap, M.K.K.; Tan, N.H. Venomics of Tropidolaemus wagleri, the sexually dimorphic temple pit viper: Unveiling a deeply conserved atypical toxin arsenal. Sci. Rep. 2017, 7, 43237. [Google Scholar] [CrossRef]
  32. Amorim, F.G.; Costa, T.R.; Baiwir, D.; De Pauw, E.; Quinton, L.; Sampaio, S.V. Proteopeptidomic, functional and immunoreactivity characterization of Bothrops moojeni snake venom: Influence of snake gender on venom composition. Toxins 2018, 10, 177. [Google Scholar] [CrossRef] [PubMed]
  33. Freitas-de-Sousa, L.A.; Nachtigall, P.G.; Portes-Junior, J.A.; Holding, M.L.; Nystrom, G.S.; Ellsworth, S.A.; Guimarães, N.C.; Tioyama, E.; Ortiz, F.; Silva, B.R.; et al. Size Matters: An evaluation of the molecular basis of ontogenetic modifications in the composition of Bothrops jararacussu snake venom. Toxins 2020, 12, 791. [Google Scholar] [CrossRef]
  34. Pe, T.; Aye, B.; Myint, A.-A.; Swe, T.-N.; Warrell, D. Effect of the snake length’s and recent feeding on venom antigenaemia and severity of envenoming in Russell’s viper (Daboia russelii siamensis) bites. In Management of Snakebite and Research; WHO, Regional Office for South-East-Asia: Geneva, Switzerland, 2002; pp. 41–45. [Google Scholar]
  35. Clemetson, K.J.; Morita, T.; Manjunatha Kini, R. Scientific and standardization committee communications: Classification and nomenclature of snake venom C-type lectins and related proteins. J. Thromb. Haemost. 2009, 7, 360. [Google Scholar] [CrossRef] [PubMed]
  36. Clemetson, K.J. Snaclecs (snake C-type lectins) that inhibit or activate platelets by binding to receptors. Toxicon 2010, 56, 1236–1246. [Google Scholar] [CrossRef] [PubMed]
  37. Morita, T. Structures and functions of snake venom CLPs (C-type lectin-like proteins) with anticoagulant-, procoagulant-, and platelet-modulating activities. Toxicon 2005, 45, 1099–1114. [Google Scholar] [CrossRef]
  38. Teixeira, C.; Fernandes, C.M.; Leiguez, E.; Chudzinski-Tavassi, A.M. Inflammation Induced by platelet-activating Viperid snake venoms: Perspectives on thromboinflammation. Front. Immunol. 2019, 10, 2082. [Google Scholar] [CrossRef]
  39. Zhong, S.R.; Jin, Y.; Wu, J.B.; Chen, R.Q.; Jia, Y.H.; Wang, W.Y.; Xiong, Y.L.; Zhang, Y. Characterization and molecular cloning of dabocetin, a potent antiplatelet C-type lectin-like protein from Daboia russellii siamensis venom. Toxicon 2006, 47, 104–112. [Google Scholar] [CrossRef] [PubMed]
  40. Zhong, S.; Jin, Y.; Li, D.; Wang, W.; Xiong, Y. Molecular cloning and sequence analyses of cDNAs encoding seven C-type lectin-like protein subunits from Daboia russellii siamensis. Zool. Res. 2005, 26, 337–343. [Google Scholar]
  41. Jebali, J.; Bazaa, A.; Sarray, S.; Benhaj, K.; Karboul, A.; El Ayeb, M.; Marrakchi, N.; Gargouri, A. C-type lectin protein isoforms of Macrovipera lebetina: cDNA cloning and genetic diversity. Toxicon 2009, 53, 228–237. [Google Scholar] [CrossRef]
  42. Inagaki, H. Snake venom protease inhibitors: Enhanced identification, expanding biological function, and promising future. In Snake Venoms; Inagaki, H., Vogel, C.-W., Mukherjee, A.K., Rahmy, T.R., Gopalakrishnakone, P., Eds.; Springer: Dordrecht, The Netherlands, 2015; pp. 1–26. [Google Scholar] [CrossRef]
  43. Thakur, R.; Mukherjee, A.K. Pathophysiological significance and therapeutic applications of snake venom protease inhibitors. Toxicon 2017, 131, 37–47. [Google Scholar] [CrossRef] [PubMed]
  44. Vasconcelos, A.A.; Estrada, J.C.; David, V.; Wermelinger, L.S.; Almeida, F.C.L.; Zingali, R.B. Structure-function relationship of the disintegrin family: Sequence signature and integrin interaction. Front. Mol. Biosci. 2021, 8, 783301. [Google Scholar] [CrossRef] [PubMed]
  45. Juárez, P.; Wagstaff, S.C.; Sanz, L.; Harrison, R.A.; Calvete, J.J. Molecular cloning of Echis ocellatus disintegrins reveals non-venom-secreted proteins and a pathway for the evolution of ocellatusin. J. Mol. Evol. 2005, 63, 183–193. [Google Scholar] [CrossRef] [PubMed]
  46. Sanz-Soler, R.; Sanz, L.; Calvette, J.J. Distribution of RPTNL genes across Reptilia: Hypothesized role for RPTLN in the evolution of SVMPs. Integr. Comp. Biol. 2016, 56, 989–1003. [Google Scholar] [CrossRef]
  47. Sanz-Soler, R.; Lorente, C.; Company, B.; Sanz, L.; Juarez, P.; Perez, A.; Zhang, Y.; Jin, Y.; Chen, R.; Eble, J.A.; et al. Recombinant expression of mutants of the Frankenstein disintegrin, RTS-ocellatusin. Evidence for the independent origin of RGD and KTS/RTS disintegins. Toxicon 2012, 60, 665–675. [Google Scholar] [CrossRef]
  48. Calvete, J.J. Structure-function correlations of snake venom disintegrins. Curr. Pharm. Des. 2005, 11, 829–835. [Google Scholar] [CrossRef] [PubMed]
  49. Wagstaff, S.C.; Favreau, P.; Cheneval, O.; Laing, G.D.; Wilkinson, M.C.; Miller, R.L.; Stöcklin, R.; Harrison, R.A. Molecular characterisation of endogenous snake venom metalloproteinase inhibitors. Biochem. Biophys. Res. Commun. 2008, 365, 650–656. [Google Scholar] [CrossRef] [PubMed]
  50. Munawar, A.; Ali, S.A.; Akrem, A.; Betzel, C. Snake venom peptides: Tools of biodiscovery. Toxins 2018, 10, 474. [Google Scholar] [CrossRef] [PubMed]
  51. Averin, A.S.; Utkin, Y.N. Cardiovascular Effects of Snake Toxins: Cardiotoxicity and Cardioprotection. Acta Nat. 2021, 13, 4–14. [Google Scholar] [CrossRef]
  52. Frangieh, J.; Rima, M.; Fajloun, Z.; Henrion, D.; Sabatier, J.M.; Legros, C.; Mattei, C. Snake venom components: Tools and cures to target cardiovascular diseases. Molecules 2021, 26, 2223. [Google Scholar] [CrossRef]
  53. Coronado, M.A.; de Moraes, F.R.; Ullah, A.; Masood, R.; Mariutti, R.; Brognaro, H.; Georgieva, D.; Murakami, M.T.; Betzel, C.; Arni, R.K. Three-dimensional structures and mechanism of snake venom serine proteinases, metalloproteinases, and phospholipase A2s. In Venom Genomics and Proteomics; Gopalakrishnakone, P., Calvete, J., Eds.; Springer: Dordrecht, The Netherlands, 2016; pp. 239–267. [Google Scholar] [CrossRef]
  54. Chaisakul, J.; Khow, O.; Wiwatwarayos, K.; Rusmili, M.R.A.; Prasert, W.; Othman, I.; Abidin, S.A.Z.; Charoenpitakchai, M.; Hodgson, W.C.; Chanhome, L.; et al. A biochemical and pharmacological characterization of phospholipase A2 and metalloproteinase fractions from Eastern Russell’s viper (Daboia siamensis) venom: Two major components associated with acute kidney injury. Toxins 2021, 13, 521. [Google Scholar] [CrossRef]
  55. Tsai, I.-H.; Tasi, H.-Y.; Wang, Y.-M.; Pe, T.; Warrel, D.A. Venom phospholipases of Russell’s vipers from Myanmar and eastern India-Cloning, characterizarion and phylogeographic analysis. Biochim. Biophys. Acta. 2007, 1774, 1020–1028. [Google Scholar] [CrossRef]
  56. Thwin, M.-M.; Gopalakrishnakone, P.; Yuen, R.; Tan, C.H. A major lethal factor of the venom of Burmese Russell’s viper (Daboia russellil siamensis): Isolation, n-terminal sequencing and biological activities of daboiatoxin. Toxicon 1995, 33, 63–76. [Google Scholar] [CrossRef]
  57. Zhang, Y.; Li, W.; Gao, R.; Lu, Q.; Wang, W.; Xiong, Y. Structure-function features and effects on blood coagulation of snake venom serine proteases. Zool. Res. 1998, 19, 97–106. [Google Scholar]
  58. Kang, T.S.; Georgieva, D.; Genov, N.; Murakami, M.T.; Sinha, M.; Kumar, R.P.; Kaur, P.; Kumar, S.; Dey, S.; Sharma, S.; et al. Enzymatic toxins from snake venom: Structural characterization and mechanism of catalysis. FEBS J. 2011, 278, 4544–4576. [Google Scholar] [CrossRef]
  59. Serrano, S.M.T.; Maroun, R.C. Snake venom serine proteinases: Sequence homology vs. substrate specificity, a paradox to be solved. Toxicon 2005, 45, 1115–1132. [Google Scholar] [CrossRef]
  60. Thakur, R.; Mukherjee, A.K. A brief appraisal on Russell’s viper venom (Daboia russelii russelii) proteinases. In Snake Venoms; Inagaki, H., Vogel, C.W., Mukherjee, A., Rahmy, T., Eds.; Springer: Dordrecht, The Netherlands, 2017; pp. 123–144. [Google Scholar] [CrossRef]
  61. Gâz Florea, Ş.A.; Gâz Florea, A.; Kelemen, H.; Muntean, D.-L. Snake Venom Metalloproteinases. Acta Medica Marisiensis 2016, 62, 106–111. [Google Scholar] [CrossRef]
  62. Takeda, S.; Takeya, H.; Iwanaga, S. Snake venom metalloproteinases: Structure, function and relevance to the mammalian ADAM/ADAMTS family proteins. Biochim. Biophys. Acta Proteins Proteom. 2012, 1824, 164–176. [Google Scholar] [CrossRef] [PubMed]
  63. Chellapandi, P. Structural, functional and therapeutic aspects of snake venom metalloproteinases. Mini. Rev. Org. Chem. 2014, 11, 28–44. [Google Scholar] [CrossRef]
  64. Takeda, S.; Igarashi, T.; Mori, H. Crystal structure of RVV-X: An example of evolutionary gain of specificity by ADAM proteinases. FEBS Lett. 2007, 581, 5859–5864. [Google Scholar] [CrossRef] [PubMed]
  65. Chen, H.S.; Chen, J.M.; Lin, C.W.; Khoo, K.H.; Tsai, I.H. New insights into the functions and N-glycan structures of factor X activator from Russell’s viper venom. FEBS J. 2008, 275, 3944–3958. [Google Scholar] [CrossRef] [PubMed]
  66. Suntravat, M.; Yusuksawad, M.; Sereemaspun, A.; Pérez, J.C.; Nuchprayoon, I. Effect of purified Russell’s viper venom-factor X activator (RVV-X) on renal hemodynamics, renal functions, and coagulopathy in rats. Toxicon 2011, 58, 230–238. [Google Scholar] [CrossRef]
  67. Chen, H.S.; Tsai, H.Y.; Wang, Y.M.; Tsai, I.H. P-III hemorrhagic metalloproteinases from Russell’s viper venom: Cloning, characterization, phylogenetic and functional site analyses. Biochimie 2008, 90, 1486–1498. [Google Scholar] [CrossRef] [PubMed]
  68. Takeda, S. ADAM and ADAMTS family proteins and snake venom metalloproteinases: A structural overview. Toxins 2016, 8, 155. [Google Scholar] [CrossRef]
  69. Hikita, A.; Tanaka, N.; Yamane, S.; Ikeda, Y.; Furukawa, H.; Tohma, S.; Suzuki, R.; Tanaka, S.; Mitomi, H.; Fukui, N. Involvement of a disintegrin and metalloproteinase 10 and 17 in shedding of tumor necrosis factor-α. Biochem. Cell Biol. 2009, 87, 581–593. [Google Scholar] [CrossRef] [PubMed]
  70. Yamazaki, Y.; Matsunaga, Y.; Tokunaga, Y.; Obayashi, S.; Saito, M.; Morita, T. Snake venom vascular endothelial growth factors (VEGF-Fs) exclusively vary their structures and functions among species. J. Biol. Chem. 2009, 284, 9885–9891. [Google Scholar] [CrossRef]
  71. Ferreira, I.G.; Pucca, M.B.; de Oliveira, I.S.; Cerni, F.A.; de Cássia da Silva Jacob, B.; Arantes, E.C. Snake venom vascular endothelial growth factors (svVEGFs): Unravelling their molecular structure, functions, and research potential. Cytokine Growth Factor Rev. 2021, 60, 133–143. [Google Scholar] [CrossRef] [PubMed]
  72. Yamazaki, Y.; Takani, K.; Atoda, H.; Morita, T. Snake venom vascular endothelial growth factors (VEGFs) exhibit potent activity through their specific recognition of KDR (VEGF receptor 2). J. Biol. Chem. 2003, 278, 51985–51988. [Google Scholar] [CrossRef]
  73. Tokunaga, Y.; Yamazaki, Y.; Morita, T. Specific distribution of VEGF-F in Viperinae snake venoms: Isolation and characterization of a VEGF-F from the venom of Daboia russelli siamensis. Arch. Biochem. Biophys. 2005, 439, 241–247. [Google Scholar] [CrossRef]
  74. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef] [PubMed]
  75. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef]
  76. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed]
  77. Macrander, J.; Panda, J.; Janies, D.; Daly, M.; Reitzel, A.M. Venomix: A simple bioinformatic pipeline for identifying and characterizing toxin gene candidates from transcriptomic data. Peer J. 2018, 6, e5361. [Google Scholar] [CrossRef]
  78. Jungo, F.; Bairoch, A.; Prot, T. the toxin protein annotation program of the Swiss-Prot protein knowledgebase. Toxicon 2005, 45, 293–301. [Google Scholar] [CrossRef] [PubMed]
  79. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-Seq: Reference generation and analysis with Trinity. Nat Protoc. 2013, 8, 8. [Google Scholar] [CrossRef]
  80. Gacesa, R.; Barlow, D.J.; Long, P.F. Machine learning can differentiate venom toxins from other proteins having non-toxic physiological functions. PeerJ Comput. Sci. 2016, 2, e90. [Google Scholar] [CrossRef]
  81. UniProt Consortium. UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res. 2021, 49, D480–D489. [Google Scholar] [CrossRef] [PubMed]
  82. Madeira, F.; Pearce, M.; Tivey, A.R.N.; Basutkar, P.; Lee, J.; Edbali, O.; Madhusoodanan, N.; Kolesnikov, A.; Lopez, R. Search and sequence analysis tools services from EMBL-EBI in 2022. Nucleic Acids Res. 2022, 50, W276–W279. [Google Scholar] [CrossRef]
  83. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  84. Bembom, O.; Ivanek, R. seqLogo: Sequence logos for DNA sequence alignments. R pakage version 1.64.0. Available online: https://bioconductor.org/packages/release/bioc/html/seqLogo.html (accessed on 26 March 2022).
Toxins 15 00309 g001 550
Figure 1. Overall toxin transcript (the top ten highly expressed toxin families) distribution in Myanmar Russell’s viper venom gland.
Figure 1. Overall toxin transcript (the top ten highly expressed toxin families) distribution in Myanmar Russell’s viper venom gland.
Toxins 15 00309 g001
Toxins 15 00309 g002 550
Figure 2. Comparison of toxin transcriptome compositions (the top ten highly expressed toxin families) of male and female snakes.
Figure 2. Comparison of toxin transcriptome compositions (the top ten highly expressed toxin families) of male and female snakes.
Toxins 15 00309 g002
Toxins 15 00309 g003 550
Figure 3. The maximum likelihood C-type lectins (CTLs) gene tree including D. siamensis with CTLs from other species. Bootstrap values greater than 50 are shown. Subtypes based on sequence alignments are indicated in Roman numerals. Refer to Supplementary Table S3 for highly expressed transcript TPM values.
Figure 3. The maximum likelihood C-type lectins (CTLs) gene tree including D. siamensis with CTLs from other species. Bootstrap values greater than 50 are shown. Subtypes based on sequence alignments are indicated in Roman numerals. Refer to Supplementary Table S3 for highly expressed transcript TPM values.
Toxins 15 00309 g003
Toxins 15 00309 g004 550
Figure 4. The peptide sequence alignment of Myanmar Russell’s viper disintegrin with acostatin-α from D. siamensis and jerdostatin from Protobothrops jerdonii. Signal peptides are underlined, while cysteine, RGD and RTS binding motifs are highlighted in blue, yellow and green correspondingly. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Figure 4. The peptide sequence alignment of Myanmar Russell’s viper disintegrin with acostatin-α from D. siamensis and jerdostatin from Protobothrops jerdonii. Signal peptides are underlined, while cysteine, RGD and RTS binding motifs are highlighted in blue, yellow and green correspondingly. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Toxins 15 00309 g004
Toxins 15 00309 g005a 550Toxins 15 00309 g005b 550
Figure 5. Sequence comparison of BPP-CNP precursors from Myanmar Russell’s viper with those of snake venom metalloprotease inhibitor (A8YPR6) from Echis ocellatus (SVMI_ECHOC) and MPi-3 (A0A6B7FNM6) (VIPAA_MPi-3) and endogenous tripeptide metalloproteinase inhibitor (A0A1I9KNP8) from Vipera ammodytes ammodytes aligned by Clustal Omega. Signal peptides are underlined. BPPs are shown by wavy underlines. C-type natriuretic peptides are highlighted in gray. Metalloproteinase tripeptide inhibitors, QKW and QRW are highlighted in yellow and green, respectively. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Figure 5. Sequence comparison of BPP-CNP precursors from Myanmar Russell’s viper with those of snake venom metalloprotease inhibitor (A8YPR6) from Echis ocellatus (SVMI_ECHOC) and MPi-3 (A0A6B7FNM6) (VIPAA_MPi-3) and endogenous tripeptide metalloproteinase inhibitor (A0A1I9KNP8) from Vipera ammodytes ammodytes aligned by Clustal Omega. Signal peptides are underlined. BPPs are shown by wavy underlines. C-type natriuretic peptides are highlighted in gray. Metalloproteinase tripeptide inhibitors, QKW and QRW are highlighted in yellow and green, respectively. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Toxins 15 00309 g005aToxins 15 00309 g005b
Toxins 15 00309 g006 550
Figure 6. Sequence comparison of human natriuretic peptides (hANP (P01160), hBNP(P16860) and hCNP (P23582)) with those of different snake venoms (DENAN_DNP (P28374) from Dentroaspis angusticeps, LACMU_CNP (Q27J49) from Lachesis muta muta, CROAT_CNP (P0CV87) from Crotalus atrox and BOTJA_CNP(Q6LEM5) from Bothrops jararaca) aligned by Clustal Omega. Conserved cysteine residues are underlined, while relatively conserved amino acid residues known for NPR binding are highlighted. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues. Important residues for NP receptor binding are highlighted in yellow color with variable residues are in green and pink color.
Figure 6. Sequence comparison of human natriuretic peptides (hANP (P01160), hBNP(P16860) and hCNP (P23582)) with those of different snake venoms (DENAN_DNP (P28374) from Dentroaspis angusticeps, LACMU_CNP (Q27J49) from Lachesis muta muta, CROAT_CNP (P0CV87) from Crotalus atrox and BOTJA_CNP(Q6LEM5) from Bothrops jararaca) aligned by Clustal Omega. Conserved cysteine residues are underlined, while relatively conserved amino acid residues known for NPR binding are highlighted. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues. Important residues for NP receptor binding are highlighted in yellow color with variable residues are in green and pink color.
Toxins 15 00309 g006
Toxins 15 00309 g007 550
Figure 7. Alignment of deduced amino acid sequences of highly expressed full-length PLA2 transcripts from transcriptome with basic PLA2 DsM-b1 (A8CG82), PLA2-I (Q7ZZQ1), acidic PLA2 RV-7 (P31100) and daboiatoxin A (Q7T2R1) from D. siamensis. Signal peptide is single-underlined. Catalytic His48, Tyr52 and Asp90 are highlighted in yellow. Ca2+ binding motif Y25-GCYCGX-GG33 is dotted-underlined. α-helices are double-underlined. The region of β-strands of the β-wing are bold-underlined. The conserved Aspartate residues at 49 position are shown in red. Cysteine residues for disulfide bonds are in green. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Figure 7. Alignment of deduced amino acid sequences of highly expressed full-length PLA2 transcripts from transcriptome with basic PLA2 DsM-b1 (A8CG82), PLA2-I (Q7ZZQ1), acidic PLA2 RV-7 (P31100) and daboiatoxin A (Q7T2R1) from D. siamensis. Signal peptide is single-underlined. Catalytic His48, Tyr52 and Asp90 are highlighted in yellow. Ca2+ binding motif Y25-GCYCGX-GG33 is dotted-underlined. α-helices are double-underlined. The region of β-strands of the β-wing are bold-underlined. The conserved Aspartate residues at 49 position are shown in red. Cysteine residues for disulfide bonds are in green. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Toxins 15 00309 g007
Toxins 15 00309 g008 550
Figure 8. Sequence alignment of VLSP-1 from Macrovipera lebetinus and deduced amino acid sequences of transcripts from the transcriptomes. Signal peptide is single-underlined. Catalytic triad His57, Asp102 and Ser195 are highlighted in yellow. Cysteine residues in six disulfide bridges are highlighted in green. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Figure 8. Sequence alignment of VLSP-1 from Macrovipera lebetinus and deduced amino acid sequences of transcripts from the transcriptomes. Signal peptide is single-underlined. Catalytic triad His57, Asp102 and Ser195 are highlighted in yellow. Cysteine residues in six disulfide bridges are highlighted in green. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Toxins 15 00309 g008
Toxins 15 00309 g009a 550Toxins 15 00309 g009b 550Toxins 15 00309 g009c 550
Figure 9. Sequence alignment of VIPAN, A disintegrin and metalloproteinase domain (ADAM)-containing protein 10 from Vipera anatolica senliki (A0A6G5ZVR9) and deduced amino acid sequences of c15329_g2_i1_M_FL and c10072_g2_i1_F_FL. Signal peptide is single-underlined. Zinc dependent metalloproteinase domain is double-underlined. Disintegrin domain is dotted-underlined. Cysteine residues are highlighted in green. Zinc binding site and RCD sequence are highlighted in yellow. (*) fully conserved residues; (:) strongly similar residues.
Figure 9. Sequence alignment of VIPAN, A disintegrin and metalloproteinase domain (ADAM)-containing protein 10 from Vipera anatolica senliki (A0A6G5ZVR9) and deduced amino acid sequences of c15329_g2_i1_M_FL and c10072_g2_i1_F_FL. Signal peptide is single-underlined. Zinc dependent metalloproteinase domain is double-underlined. Disintegrin domain is dotted-underlined. Cysteine residues are highlighted in green. Zinc binding site and RCD sequence are highlighted in yellow. (*) fully conserved residues; (:) strongly similar residues.
Toxins 15 00309 g009aToxins 15 00309 g009bToxins 15 00309 g009c
Toxins 15 00309 g010a 550Toxins 15 00309 g010b 550
Figure 10. Sequence alignment of AGKCO (snake venom metalloproteinase from Agkistrodon contortrix contortrix, A0A1W7RJN7) and deduced amino acid sequences of c19271_g2_i1_M_FL. Signal peptide is single-underlined. Zinc dependent metalloproteinase domain is double-underlined. The disintegrin domain is dash-underlined. ADAM17-MPD domain is wave-underlined. CxxC motif (membrane-proximal domain) is highlighted in red. PxxP motif is highlighted in pink. Zinc binding site and EECD sequence are highlighted in yellow. Cysteine residues are highlighted in green. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Figure 10. Sequence alignment of AGKCO (snake venom metalloproteinase from Agkistrodon contortrix contortrix, A0A1W7RJN7) and deduced amino acid sequences of c19271_g2_i1_M_FL. Signal peptide is single-underlined. Zinc dependent metalloproteinase domain is double-underlined. The disintegrin domain is dash-underlined. ADAM17-MPD domain is wave-underlined. CxxC motif (membrane-proximal domain) is highlighted in red. PxxP motif is highlighted in pink. Zinc binding site and EECD sequence are highlighted in yellow. Cysteine residues are highlighted in green. (*) fully conserved residues; (:) strongly similar residues; (.) weakly similar residues.
Toxins 15 00309 g010aToxins 15 00309 g010b
Toxins 15 00309 g011 550
Figure 11. Sequence alignment of deduced amino acid sequences of VEGF transcripts with vascular endothelial growth factor 1 (VEGF1) from Vipera anatolica senliki (QHR82857) (A) and vascular endothelial growth factor A precursor (VEGFA-VIPAA) from Vipera ammodytes ammodytes (C0K3N5) (B). Signal peptide is single-underlined. Cysteine residues are in yellow color. (*) fully conserved residues; (:) strongly similar residues.
Figure 11. Sequence alignment of deduced amino acid sequences of VEGF transcripts with vascular endothelial growth factor 1 (VEGF1) from Vipera anatolica senliki (QHR82857) (A) and vascular endothelial growth factor A precursor (VEGFA-VIPAA) from Vipera ammodytes ammodytes (C0K3N5) (B). Signal peptide is single-underlined. Cysteine residues are in yellow color. (*) fully conserved residues; (:) strongly similar residues.
Toxins 15 00309 g011
Table
Table 1. The sequencing and the assembly quality of the venom gland transcriptomes of Myanmar Russell’s vipers (Daboia siamensis).
Table 1. The sequencing and the assembly quality of the venom gland transcriptomes of Myanmar Russell’s vipers (Daboia siamensis).
Trinity Assembly ParametersMaleFemale
Total raw reads82,251,23279,578,046
Total clean reads78,647,67074,444,200
Total clean nucleotides (nt)7,752,565,8677,235,655,822
Q20 percentage98.93%98.55%
GC percentage46.26%53.43%
Contigs created88,32550,858
Total assembled bases75,163,97639,200,694
Mean length (nt)850.99770.79
N5018701379
Unigenes/transcripts assembled (only longest isoform per gene)77,66846,884
Table
Table 2. Kunitz-type protease inhibitor (KSPI) transcripts from Myanmar Russell’s viper transcriptomes.
Table 2. Kunitz-type protease inhibitor (KSPI) transcripts from Myanmar Russell’s viper transcriptomes.
No.Contig IDAnnotationAccession No.SpeciesAmino Acid TPM Value
1c17840_g2_i1_M_FLKunitz-type serine protease inhibitor 1Q2ES50 Daboia russelii84109,755.22
c11654_g1_i1_F_FL16,745.60
2c6152_g1_i1_M_FLSerine protease inhibitor 1A0A6G5ZVV7Vipera anatolica senliki1481.59
c6152_g1_i2_M_FL0.76
3c3688_g1_i1_M_FLSerine protease inhibitor 4A0A6G5ZUI5Vipera anatolica senliki51426.46
c5062_g2_i1_F_FL9.09
4c12297_g1_i1_M_FLTissue factor pathway inhibitorU3TBN9Protobothrops flavoviridis3032.45
5c1956_g1_i1_M_PLSerine protease inhibitor 3A0A6G5ZW53Vipera anatolica senliki1031.48
6c12495_g1_i1_M_PLSerine protease inhibitor 6A0A6G5ZUI2Vipera anatolica senliki27256.61
c12495_g1_i2_M_PL50.11
c12495_g1_i3_M_PL74.49
c2247_g1_i2_F_PL51.7
c2247_g1_i1_F_PL33.08
7c15981_g1_i1_M_PLkunitz/BPTI-like toxinA0A6P9AX37Protobothrops guttatus13947.25
8c11558_g3_i1_F_PLkunitz/BPTI-like toxinXP_015682430.1Protobothrops mucrosquamatus13925.45
Table
Table 3. Transcripts with snake venom serine protease (SVSP) annotation.
Table 3. Transcripts with snake venom serine protease (SVSP) annotation.
No.Contig IDAnnotationAccession No.SpeciesAmino Acid TPM Value
1c20548_g3_i4_M_FLFactor V activator RVV-VγP18965Daboia siamensis2595495.65
2c11820_g1_i3_F_PLThrombin-like_enzyme_TLBm (Factor V activator RVV-Vγ)P18965D. siamensis2512205.37
3c11820_g1_i1_F_PLThrombin-like_enzyme_TLBm (Factor V activator RVV-Vγ)P18965D. siamensis2411794.30
4c11820_g1_i4_F_PLThrombin-like_enzyme_TLBm (Factor V activator RVV-Vγ)P18965D. siamensis2511355.44
5c38983_g1_i1_F_PLα-fibrinogenase-likeE5L0E3D. siamensis1053768.67
6c20548_g3_i1_M_PLα-fibrinogenase-likeE5L0E3D. siamensis1313196.39
7c20548_g3_i3_M_PLα-fibrinogenase-likeE5L0E3D. siamensis155952.41
8c11820_g1_i9_F_PLThrombin-like_enzyme_catroxobin-1 (β-fibrinogenase-like)E5L0E4D. siamensis126981.63
9c20548_g3_i5_M_FLSerine protease VLSP-1E0Y418Microvipera lebetina2653027.55
10c11820_g1_i8_F_PLThrombin-like_enzyme_TLBm (Serine protease VLSP-1)E0Y418M. lebetina2421400.97
11c11820_g1_i6_F_PLThrombin-like_enzyme_TLBm (Serine protease VLSP-1)E0Y418M. lebetina189816.73
12c11820_g1_i2_F_PLThrombin-like_enzyme_TLBm (Serine protease VLSP-1)E0Y418M. lebetina126448.44
13c11820_g1_i5_F_PLThrombin-like_enzyme_TLBm (Serine protease VLSP-3)E0Y420M. lebetina1341318.30
14c11820_g1_i7_F_PLThrombin-like_enzyme_catroxobin-1 (β-fibrinogenase-like)E5L0E4D. siamensis104919.53
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yee, K.T.; Macrander, J.; Vasieva, O.; Rojnuckarin, P. Exploring Toxin Genes of Myanmar Russell’s Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics. Toxins 2023, 15, 309. https://doi.org/10.3390/toxins15050309

AMA Style

Yee KT, Macrander J, Vasieva O, Rojnuckarin P. Exploring Toxin Genes of Myanmar Russell’s Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics. Toxins. 2023; 15(5):309. https://doi.org/10.3390/toxins15050309

Chicago/Turabian Style

Yee, Khin Than, Jason Macrander, Olga Vasieva, and Ponlapat Rojnuckarin. 2023. "Exploring Toxin Genes of Myanmar Russell’s Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics" Toxins 15, no. 5: 309. https://doi.org/10.3390/toxins15050309

APA Style

Yee, K. T., Macrander, J., Vasieva, O., & Rojnuckarin, P. (2023). Exploring Toxin Genes of Myanmar Russell’s Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics. Toxins, 15(5), 309. https://doi.org/10.3390/toxins15050309

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop