JVATiTD - Articles

Official publication of CEVAP/UNESP

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Leandro de Mattos Pereira1,2, Elisa Alves Messias1, Bruna Pereira Sorroche1, Angela das Neves Oliveira1, Lidia Maria Rebolho Batista Arantes1, Ana Carolina de Carvalho1, Anita Mitico Tanaka-Azevedo3, Kathleen Fernandes Grego3, André Lopes Carvalho1, Matias Eliseo Melendez1,4,5 [ + show more ]

J Venom Anim Toxins incl Trop Dis, 2020, 26:e20190058
Received: 24 September 2019 | Accepted: 16 September 2020 | Published online: 21 October 2020
Collection: Snake venoms: from production to bioprospecting


Background: Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design. Methods: RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment. Results: We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin. Conclusion: It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.


Keywords: Bothrops jararaca; Venom gland; Transcriptome; Biotechnological application; Stonustoxin; Verrucotoxin.

Full Article PDF