BIOLOGICAL SIGNIFICANCE: Advents in proteomics and bioinformatics have vigorously propelled the scientific discoveries of toxins from various lineages of venomous snakes. The Malayan pit viper, Calloselasma rhodostoma, is a medically important species in Southeast Asia as its bite can cause envenomation, while the venom is also a source of bioactive compounds for drug discovery. Detailed profiling of the venom, however, is inadequate possibly due to the complex nature of the venom and technical limitation in separating the constituents into details. Integrating a multi-step fractionation method, this study successfully revealed a comprehensive and quantitative profile of the composition of the venom of this medically important venomous snake. The relative abundance of the various venom proteins is determined in a global profile, providing useful information for understanding the pathogenic roles of the different toxins in C. rhodostoma envenomation. Notably, the principal hemotoxins were identified in great details, including the variety of toxin subunits and isoforms. The findings indicate that these toxins are the principal targets for effective antivenom neutralization, and should be addressed in the production of a pan-regional polyspecific antivenom. In addition, minor toxin components not reported previously in the venom were also detected in this study, enriching the current toxin database for the venomous snakes.
Methods: We have employed high-throughput RNA-Seq technology to uncover the transcriptome changes of P. monodon hepatopancreas when challenged with VpAHPND. The shrimps were challenged with VpAHPND through immersion method with dissected hepatopancreas samples for the control group (APm-CTL) and treatment group at 3 (APm-T3), 6 (APm-T6), and 24 (APm-T24) hours post-AHPND infection sent for RNA-Seq. The transcriptome de novo assembly and Unigene expression determination were conducted using Trinity, Tgicl, Bowtie2, and RSEM software. The differentially expressed transcripts were functionally annotated mainly through COG, GO, and KEGG databases.
Results: The sequencing reads generated were filtered to obtain 312.77 Mb clean reads and assembled into 48662 Unigenes. Based on the DEGs pattern identified, it is inferred that the PAMPs carried by VpAHPND or associated toxins are capable of activating PRRs, which leads to subsequent pathway activation, transcriptional modification, and antibacterial responses (Phagocytosis, AMPs, proPO system). DAMPs are released in response to cell stress or damage to further activate the sequential immune responses. The comprehensive interactions between VpAHPND, chitin, GbpA, mucin, chitinase, and chitin deacetylase were postulated to be involved in bacterial colonization or antibacterial response.
Conclusions: The outcomes of this research correlate the different stages of P. monodon immune response to different time points of AHPND infection. This finding supports the development of biomarkers for the detection of early stages of VpAHPND colonization in P. monodon through host immune expression changes. The potential genes to be utilized as biomarkers include but not limited to C-type lectin, HMGB1, IMD, ALF, serine proteinase, and DSCAM.
METHODS: The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.
RESULTS AND DISCUSSION: The toxin transcripts showed high redundancy (41-82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin's fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.