Acessibilidade / Reportar erro

Gene sequence analysis of toxins from the spider Phoneutria nigriventer revealed an intronless feature

Abstract

Background:

Phoneutria nigriventer spider venom contains several cysteine-rich peptide toxins that act on different ion channels. Despite extensive studies on its venom and description of cDNA sequences of several of its toxin precursors, the gene structure of these toxins remains unknown.

Methods:

Genomic regions encoding the precursors of three previously characterized P. nigriventer toxins - PnTx1, PnTx2-5 and PnTx4(5-5) - were amplified by PCR using specific primers. PCR fragments were cloned and sequenced. Obtained sequences were compared with their corresponding cDNA sequences.

Results:

The size of PCR fragments obtained and sequences corresponding to genomic regions encoding for the toxin precursors matched their cDNA sequences.

Conclusions:

Despite a few nucleotide substitutions in the genomic regions encoding for the toxin precursors when compared with cDNA sequences, the results of the present work indicate that P. nigriventer toxins do not contain introns in their genes sequences.

Keywords:
Phoneutria; Spider toxins; Toxin genes

Background

Phoneutria nigriventer is one of the largest existing spiders from the suborder Araneomorphae (RTA clade; family Ctenidae) and one of the few in the world that can cause harm to humans [11. de Lima ME, Figueiredo SG, Matavel A, Nunes KP, da Silva CN, Almeida F de M, et al. Phoneutria nigriventer Venom and Toxins: A Review. Spider Venoms. 2015. p. 1-14. ,22. Peigneur S, de Lima ME, Tytgat J. Phoneutria nigriventer venom: A pharmacological treasure. Toxicon [Internet]. Elsevier Ltd; 2018;151:96-110. Available from: https://doi.org/10.1016/j.toxicon.2018.07.008
https://doi.org/10.1016/j.toxicon.2018.0...
]. They are wandering, solitary and aggressive spiders, relying on strength and venom toxicity for defense and prey capture rather than using silk webs. P. nigriventer venom contains several neurotoxic peptides that act on ion channels and chemical receptors of vertebrates and invertebrates [33. Gomez MV, Kalapothakis E, Guatimosim C, Prado MAM. Phoneutria nigriventer Venom : A Cocktail of Toxins That Affect Ion Channels. Cell Mol Neurobiol. 2002;22:579-88. ].

Although this venom has been studied for over 40 years and cDNA sequences of several toxin precursors have been obtained [44. Diniz MRV, Paiva ALB, Guerra-Duarte C, Nishiyama MY, Mudadu MA, De Oliveira U, et al. An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches. PLoS One. 2018;13:1-29. -1010. Diniz MR V, Paine MJI, Diniz CR, Theakston RDG, Crampton JM. Sequence of the cDNA coding for the lethal neurotoxin Tx1 from the Brazilian “armed” spider Phoneutria nigriventer predicts the synthesis and processing of a preprotoxin. J Biol Chem. 1993;268:15340-2. ], up to this moment there is no investigation on the genome or even gene structure of any of its toxins.

The earliest spiders dates back to about 300 million years ago. Due to this long existence, together with the high species diversity in this group, spiders have been studied as an interesting evolution model. As an important trait for their evolutionary success, the study of venom and toxin evolution in the molecular level may contribute in elucidating the complex history of spiders [1111. King GF. The wonderful world of spiders: Preface to the special Toxicon issue on spider venoms. Toxicon. 2004;43:471-5. ]. Morphological and behavioral data have been traditionally used to infer phylogeny, but genomic and transcriptomic molecular data are recently challenging previous assumptions on the tree of life of spiders [1212. Garrison NL, Rodriguez J, Agnarsson I, Coddington JA, Griswold CE, Hamilton CA, et al. Spider phylogenomics: Untangling the Spider Tree of Life. PeerJ. 2016;2016. ,1313. Fernández R, Kallal RJ, Dimitrov D, Ballesteros JA, Arnedo MA, Giribet G, et al. Phylogenomics, Diversification Dynamics, and Comparative Transcriptomics across the Spider Tree of Life. Curr Biol. 2018;28:1489-1497.e5. ]. The genetic architecture of their toxin genes has shown to be variable in different spider major groups - such as Araneomorphs [1414. Krapcho KJ, Kral RM, Vanwagenen BC, Eppler KG, Morgan TK. Characterization and cloning of insecticidal peptides from the primitive weaving spider Diguetia canities. Insect Biochem Mol Biol. 1995;25:991-1000. ] and Mygalomorphs [1515. Qiao P, Zuo XP, Chai ZF, Ji YH. The cDNA and genomic DNA organization of a novel toxin SHT-I from spider Ornithoctonus huwena. Acta Biochim Biophys Sin (Shanghai). 2004;36:656-60. -1818. Pineda SS, Wilson D, Mattick JS, King GF. The lethal toxin from Australian funnel-web spiders is encoded by an intronless gene. PLoS One. 2012;7. ] -, and genes encoding for venom toxins have been described both with and without introns. Therefore, the description of Phoneutria toxins gene structure, as a member of the RTA clade (the most diverse group within spiders), is a relevant contribution to the field.

Several novel components of high molecular mass have been described in P. nigriventer venom, but cysteine-rich peptide toxins are by far its most abundant component [44. Diniz MRV, Paiva ALB, Guerra-Duarte C, Nishiyama MY, Mudadu MA, De Oliveira U, et al. An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches. PLoS One. 2018;13:1-29. ]. PnTx1 (μ-ctenitoxin-Pn1a), PnTx2-5 (δ-ctenitoxin-Pn2c) and PnTx4(5-5) (ɣ-ctenitoxin-Pn1a) are among the best characterized and most abundant cysteine-rich peptide toxins in P. nigriventer venom. PnTx1 was demonstrated to inhibit sodium channel currents [1919. Martin-Moutot N, Haro L De, Santos RG Dos, Mori Y, Seagar M. Phoneutria nigriventer ω-Phonetoxin IIA: A new tool for anti-calcium channel autoantibody assays in Lambert-Eaton myasthenic syndrome. Neurobiol Dis. 2006;22:57-63. ,2020. Silva AO, Peigneur S, Diniz MRV, Tytgat J, Beirão PSL. Inhibitory effect of the recombinant Phoneutria nigriventer Tx1 toxin on voltage-gated sodium channels. Biochimie [Internet]. Elsevier Masson SAS; 2012;94:2756-63. Available from: http://dx.doi.org/10.1016/j.biochi.2012.08.016
http://dx.doi.org/10.1016/j.biochi.2012....
] and has shown great neurotoxicity, inducing tail elevation, excitation, salivation, spastic paralysis and lethality in mice [2121. Diniz MRV, Theakston RDG, Crampton JM, Nascimento Cordeiro M do, Pimenta AMC, De Lima ME, et al. Functional expression and purification of recombinant Tx1, a sodium channel blocker neurotoxin from the venom of the Brazilian “armed” spider, Phoneutria nigriventer. Protein Expr Purif. 2006;50:18-24. ]. PnTx2-5 acts on sodium channels as well, being able to inhibit sodium channel inactivation [2222. Matavel A, Fleury C, Oliveira LC, Molina F, De Lima ME, Cruz JS, et al. Structure and activity analysis of two spider toxins that alter sodium channel inactivation kinetics. Biochemistry. 2009;48:3078-88. ]. It is also one of the P. nigriventer’s most toxic venom components to mice [2323. Cordeiro M do N, Diniz CR, do Carmo Valentim A, von Eickstedt VRD, Gilroy J, Richardson M. The purification and amino acid sequences of four Tx2 neurotoxins from the venom of the Brazilian “armed” spider Phoneutria nigriventer (Keys). FEBS Lett. 1992;310:153-6. ] and can induce penile erection, hypersalivation and death by respiratory distress or pulmonary edema [2424. Yonamine CM, Troncone LRP, Camillo MAP. Blockade of neuronal nitric oxide synthase abolishes the toxic effects of Tx2-5, a lethal Phoneutria nigriventer spider toxin. Toxicon . 2004;44:169-72. ]. Unlike these two first toxins, PnTx4(5-5) is toxic only to insects [2525. De Figueiredo SG, De Lima ME, Cordeiro MN, Diniz CR , Patten D, Halliwell RF, et al. Purification and amino acid sequence of a highly insecticidal toxin from the venom of the Brazilian spider Phoneutria nigriventer which inhibits NMDA-evoked currents in rat hippocampal neurones. Toxicon . 2001;39:309-17. ], showing a remarkable effect on insect sodium channels inactivation [2626. Paiva ALB, Matavel A, Peigneur S, Cordeiro M do N, Tytgat J, Diniz MR V, et al. Differential effects of the recombinant toxin PnTx4(5-5) from the spider Phoneutria nigriventer on mammalian and insect sodium channels. Biochimie. 2016;121:326-35. ]. This toxin also seems to hold a biotechnological potential use as a neuroprotective [2727. Silva FR, Batista EML, Gomez M V., Kushmerick C, Da Silva JF, Cordeiro MN, et al. The Phoneutria nigriventer spider toxin, PnTx4-5-5, promotes neuronal survival by blocking NMDA receptors. Toxicon [Internet]. Elsevier Ltd; 2016;112:16-21. Available from: http://dx.doi.org/10.1016/j.toxicon.2016.01.056
http://dx.doi.org/10.1016/j.toxicon.2016...
] and analgesic drug lead [2828. Oliveira CFB, Alves DP, Emerich BL, Figueiredo SG de, Cordeiro M do N, Borges MH, et al. Antinociceptive effect of PnTx4(5-5), a peptide from Phoneutria nigriventer spider venom, in rat models and the involvement of glutamatergic system. J Venom Anim Toxins Incl Trop Dis. 2019;25:1-9. ].

Considering the relevance of P. nigriventer toxins and due to their biotechnological, medical and evolutionary importance, we have investigated the gene structure of these three sodium channel modulators toxins from P. nigriventer venom, searching for the presence or absence of introns.

Methods

Venom glands and genomic DNA (gDNA) were obtained from P. nigriventer adult spiders maintained at Ezequiel Dias Foundation in Belo Horizonte, Brazil (Sisgen #A26E945).

Isolation of RNA and cDNA synthesis

Total RNA was extracted from pooled venom glands of five adult specimens using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s protocol. Subsequently, total RNA was used to synthesize cDNA first strand using Super Script First-Strand Synthesis System for RT-PCR kit (Invitrogen, USA), following the manufacturer’s protocol.

Isolation of genomic DNA

Genomic DNA (gDNA) was extracted from the leg muscular tissue of an adult P. nigriventer specimen, as described by Fan and Gulley [2929. Killeen AA, Fan H, Gulley ML. DNA Extraction from Fresh or Frozen Tissues. Mol Pathol Protoc. 2003;49:5-10. ]. After spider euthanasia, two legs were isolated from the base of the cephalothorax and ground with a chilled, sterile mortar and pestle in liquid nitrogen. After pulverization, the tissue was added to 450 µL of extraction buffer (0.01 M NaCl, 20 mM Tris- HCl pH 8, 1 mM EDTA, 1% SDS and 300 µg/mL proteinase K) and incubated at 55°C for 3h. One volume of phenol:chloroform:isoamyl alcohol (25:24:1) was added and the aqueous phase was collected after centrifugation. DNA was recovered by adding 50 µL of 3 M NaOAc pH 5.2, followed by ethanol precipitation and re-suspended in 100 µL of TE buffer (Tris 10mM, EDTA 1mM).

Amplification and sequencing of toxin genes

Sequences of gDNA and cDNA encoding for toxins were amplified by PCR using specific primers (Table 1). Primers sequences for PnTx1 and PnTx2-5 precursors were designed using information from cDNA sequences deposited in GenBank database (accession numbers X73155.1 and AF014463.1, respectively). Primers for PnTx4(5-5) precursor and mature sequences were designed using a cDNA sequence obtained from a P. nigriventer venom glands transcriptome study [44. Diniz MRV, Paiva ALB, Guerra-Duarte C, Nishiyama MY, Mudadu MA, De Oliveira U, et al. An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches. PLoS One. 2018;13:1-29. ].

Table 1.
Primer sequences used for amplification of gDNA and cDNA sequences of each toxin.

PCR reactions contained 100 ng of gDNA or cDNA, 1x PCR buffer, 1 μM dNTPs, 1 μM of each primer and 1 U Platinum Taq DNA polymerase High Fidelity (Invitrogen, USA). Cycling conditions were: 5 min at 94ºC followed by 35 cycles at 94°C for 45 s, 50°C for 45 s, 72°C for 1 min, and a final cycle of 72ºC for 7 min. PCR products were subjected to electrophoresis on 2% agarose gel. gDNA fragments extracted from the agarose gel were purified using Illustra GFX PCR DNA and Gel Band Purification kit (GE Healthcare, USA) and cloned into a pGEM-T Easy vector (Promega, USA). Positive clones were sequenced in an ABI PRISMTM 3700 DNA Automatic Sequencer using the standard M13 reverse primer and Big Dye Terminator v3.1 Cycle sequencing kit (Applied Biosystems, USA).

Analysis of intron presence in cysteine-rich peptide neurotoxins in available spider genomes

Sequenced spider genomes available in NCBI Genome public database (accession number for each species are as follows: Loxosceles reclusa - JJRW01; Stegodyphus mimosarum - AZAQ01; Araneous ventricosus - BGPR01; Dysdera sylvatica - QLNU01; Trichonephila clavipes - MWRG01; Acanthoscurria geniculate - AZMS01; Pardosa pseudoannulata - SBLA01) were used to search for cysteine-rich peptide toxins sequences and evaluate the presence of introns within their genes.

Since it was not possible to find a specific cysteine-rich peptide toxin sequence published for each of these spiders species, multi FASTA archives were built with known toxin sequences from other spider species, annotated as inhibitor cysteine knot (ICK), knottin or cysteine-rich peptide toxin retrieved from ArachnoServer [3030. Pineda SS, Chaumeil P, Kunert A, Thang MWC, Li L, Nuhn M, et al. Spider toxin database ArachnoServer 3.0: an online resource for automated discovery , analysis and annotation of spider toxins. Bioinformatics. 2017;7-8. ], Knottin database [3131. Postic G, Gracy J, Périn C, Chiche L, Gelly JC. KNOTTIN: The database of inhibitor cystine knot scaffold after 10 years, toward a systematic structure modeling. Nucleic Acids Res. 2018;46:D454-8. ] and GenBank [3232. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, et al. GenBank. Nucleic Acids Res. 2013;41:36-42. ] to be used as queries to interrogate the genomes. Genomes were used individually as databases to run BLASTn and tBLASTn [3333. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403-10. ], with default parameters. Results were manually inspected, considering hits with e-value < 1e-05, searching for breaks in the obtained alignment, which would indicate the discontinuity of the sequence, pointing to the presence of introns.

Results

Using information from cDNA sequences, we designed specific primers in order to amplify genomic regions encoding for the precursors of three well characterized toxins from the spider P. nigriventer, namely PnTx1, PnTx2-5, and PnTx4(5-5). The forward primers are complementary to the beginning of the sequence region encoding for the toxins signal peptides. Reverse primers were designed to anneal to the 3’ UTR regions.

The results showed that the size of PCR fragments corresponding to amplified gDNA encoding for precursor sequences of PnTx1 and PnTx2-5 matched the amplified cDNA fragments, which were expected to be 396 bp and 271 bp, respectively (Figure 1). This indicates that these two toxins do not contain introns in their precursor gene sequences.

Figure 1.
Agarose 2% gel showing the comparison between the fragments of amplified gDNA and cDNA sequences encoding for the toxins. 1: 100 bp; 2: PnTx1 gDNA; 3: PnTx1 cDNA; 4: PnTx2-5 gDNA; 5: PnTx2-5 cDNA; 6: PnTx4 (5-5) gDNA mature sequence; 7: PnTx4(5-5) cDNA mature sequence.

Despite several attempts, we could not obtain any amplification for PnTx4(5-5) using the forward primers annealing to the sequence region encoding for the toxin signal peptide with gDNA as a template. This lack of amplification may be due to the existence of extensive nucleotide variation in DNA sequences encoding for toxins in spiders [3434. Vassilevski AA, Kozlov SA, Grishin E V. Molecular diversity of spider venom. Biochem. 2009;74:1505-34. ]. Indeed, it has been demonstrated that PnTx4(5-5) is one of the most expressed toxins in P. nigriventer venom glands, presenting several isoforms [44. Diniz MRV, Paiva ALB, Guerra-Duarte C, Nishiyama MY, Mudadu MA, De Oliveira U, et al. An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches. PLoS One. 2018;13:1-29. ]. Thus, we have designed a new forward primer to amplify only the mature sequence of PnTx4(5-5). With this new approach, amplification was observed (Fig. 1) and as for the other tested toxins, gDNA and cDNA amplified fragments presented the same size (181 bp). This indicates that at least the gDNA region encoding for PnTx4(5-5) mature sequence is also intronless.

In order to confirm if the amplified gDNA sequences of PnTx1 and PnTx2-5 matched the cDNA sequences deposited in the GenBank database (accession numbers X73155.1 and AF014463.1, respectively), the PCR products from gDNA were cloned and sequenced. Each sequence obtained for PnTx1 and PnTx2-5 presented one non-synonymous nucleotide substitution in the propeptide region, leading to one amino acid alteration when compared with the corresponding cDNA sequence deposited in GenBank (Figure 2). Furthermore, PnTx1 gDNA sequence also presented two synonymous substitutions in the toxin mature sequence (Fig. 2A) and PnTx2-5 presented one synonymous substitution in the signal peptide sequence (Fig. 2B). The gDNA sequences obtained for PnTx1 and PnTx2-5 were deposited in GenBank database under the accession numbers MN851275 and MN851276, respectively.

Figure 2.
Alignment of cDNA and gDNA sequences of (A) PnTx-1 and (B) PnTx2-5. The regions corresponding to signal peptide, propeptide, mature toxin and 3’ ÚTR are schematically shown above the sequences in white, gray, black and dark gray boxes, respectively. Primer regions are underlined in cDNA sequences. Nucleotide substitutions are indicated in red letters. The translated protein sequences are shown in black letters between the cDNA and gDNA nucleotide sequences and amino acid substitutions are indicated in orange letters. Hashtag indicates stop codons. cDNA sequences for PnTx-1 and PnTx2-5 were obtained from GenBank accession numbers X73155.1 and AF014463.1, respectively.

Our results were compared with other published data concerning the presence of introns in spider cysteine-rich peptide toxins (Figure 3). Available spider genome sequences were also used for an in-silico analysis, using a multiFASTA archive containing several cysteine-rich peptide toxin sequences. We did not find evidence of introns in the retrieved alignments with the genomic sequences of A. geniculata (Mygalomorphae), since all alignments were contiguous with the queries. Although most toxin sequences found in Araneomaphae genomes seem to be contiguous, at least three sequences in A. ventricosus, two sequences in P. pseudoannulata, two sequences in S. mimosarum and two sequences in T. clavipes showed evidence of possible introns, but a more in-depth analysis is required to confirm this.

Figure 3.
Summary of the presence of introns in spider toxin genes. According to the reviewed literature, all toxin genes from Mygalomorphae spiders studied until now are intronless. In the Araneomorphae group, toxin genes are structured both with and without introns.

Discussion

Spiders are an evolutionary successful group, with a high number of species adapted to different environments. This diversity of adaptations relies on essential molecules produced by them, which compose venom and silk. Studying the gene structure of these molecules may help to elucidate how spiders have evolved and adapted, disclosing possible mechanisms for generating molecular diversity [3535. Garb JE, Sharma PP, Ayoub NA. Recent progress and prospects for advancing arachnid genomics. Curr Opin Insect Sci [Internet]. Elsevier Inc; 2018;25:51-7. Available from: http://dx.doi.org/10.1016/j.cois.2017.11.005
http://dx.doi.org/10.1016/j.cois.2017.11...
].

In the present work we have compared the cDNA sequences (derived from mRNA expressed in P. nigriventer venom glands), and the genomic sequences (retrieved from muscular tissue of the spiders’ legs) of three of the main toxins from P. nigriventer venom, in order to elucidate this Ctenidae toxin gene structure. By comparing cDNA and gDNA sequences of toxins PnTx1 and PnTx2-5, both synonymous and non-synonymous point mutations leading to amino acid substitutions were found in the signal peptide, propeptide and in the mature toxin sequence. These nucleotide and amino acid substitutions were expected since spider toxins are known to have many isoforms, forming groups of related sequences differing by point mutations, even within a single spider specimen, so-called combinatorial libraries [3434. Vassilevski AA, Kozlov SA, Grishin E V. Molecular diversity of spider venom. Biochem. 2009;74:1505-34. ]. In P. nigriventer, the existence of extensive variation in signal peptide, propeptide and mature sequences for toxins was also already demonstrated [44. Diniz MRV, Paiva ALB, Guerra-Duarte C, Nishiyama MY, Mudadu MA, De Oliveira U, et al. An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches. PLoS One. 2018;13:1-29. ]. In addition, these are abundant toxins that are transcribed at a high frequency, which has been associated with elevated mutation rates, so called transcription-associated mutation (TAM) [3636. Aguilera A. The connection between transcription and genomic instability. EMBO J. 2002;21:195-201. ].

The most important finding of the present work is the description, for the first time, of the intronless genomic structure of some important P. nigriventer cysteine-rich peptide toxins. Regarding the presence of introns in toxins genes, 11 spider species have been analyzed to date (including this study) (Fig. 3). All spider genes encoding for cysteine-rich peptide toxins from infra-order Mygalomorphae analyzed so far are intronless [1515. Qiao P, Zuo XP, Chai ZF, Ji YH. The cDNA and genomic DNA organization of a novel toxin SHT-I from spider Ornithoctonus huwena. Acta Biochim Biophys Sin (Shanghai). 2004;36:656-60. -1818. Pineda SS, Wilson D, Mattick JS, King GF. The lethal toxin from Australian funnel-web spiders is encoded by an intronless gene. PLoS One. 2012;7. ]. Among the Araneomorphae, the gene for a cysteine-rich peptide toxin found in Diguetia canities venom, µ-diguetoxin-Dc1a, has an intron-exon structure [1414. Krapcho KJ, Kral RM, Vanwagenen BC, Eppler KG, Morgan TK. Characterization and cloning of insecticidal peptides from the primitive weaving spider Diguetia canities. Insect Biochem Mol Biol. 1995;25:991-1000. ] whereas two-domain neurotoxins from Cheiracanthium punctorium and Oxyopes lineatus do not present introns in their toxin gene sequences [3737. Sachkova MY, Slavokhotova AA, Grishin E V., Vassilevski AA. Structure of the yellow sac spider Cheiracanthium punctorium genes provides clues to evolution of insecticidal two-domain knottin toxins. Insect Mol Biol. 2014;23:527-38. ,3838. Sachkova MY, Slavokhotova AA, Grishin E V., Vassilevski AA. Genes and evolution of two-domain toxins from lynx spider venom. FEBS Lett [Internet]. Federation of European Biochemical Societies; 2014;588:740-5. Available from: http://dx.doi.org/10.1016/j.febslet.2014.01.018
http://dx.doi.org/10.1016/j.febslet.2014...
], like we have shown in the present report for P. nigriventer toxin genes (Fig. 3).

The analysis of other toxin classes in other Araneomorphae species showed a diverse scenario. Latrotoxins genes were found to be intronless in Latrodectus tridemciguttatus [3939. Danilevich VN, Grishin EV. The Genes Encoding Black Widow Spider Neurotoxins are Intronless. Russ J Bioorganic Chem. 2000;26:838-43. ], but introns were present in these type of toxins genes in Parasteatoda tepidariorum [4040. Gendreau KL, Haney RA, Schwager EE, Wierschin T, Stanke M, Richards S, et al. House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity. BMC Genomics; 2017;18:1-14. ] and Latrodectus hesperus [4141. Bhere KV, Haney RA, Ayoub NA, Garb JE. Gene structure, regulatory control, and evolution of black widow venom latrotoxins. FEBS Lett. 2014;588:3891-7. ]. Sphingomyelinases D from Loxosceles arizonica [4242. Binford GJ, Cordes MHJ, Wells MA. Sphingomyelinase D from venoms of Loxosceles spiders: Evolutionary insights from cDNA sequences and gene structure. Toxicon . 2005;45:547-60. ] also presented introns in its genes.

To the best of our knowledge, only seven spider genomes have been published to date: Acanthoscurria geniculata (Mygalomorphae) and Stegodyphus mimosarum (Araneomorphae) [4343. Sanggaard KW, Bechsgaard JS, Fang X, Duan J, Dyrlund TF, Gupta V, et al. Spider genomes provide insight into composition and evolution of venom and silk. Nat Commun. 2014;5. ]; Trichonephila clavipes (Araneomorphae) [4444. Babb PL, Lahens NF, Correa-Garhwal SM, Nicholson DN, Kim EJ, Hogenesch JB, et al. The Nephila clavipes genome highlights the diversity of spider silk genes and their complex expression. Nat Genet [Internet]. Nature Publishing Group; 2017;49:895-903. Available from: http://dx.doi.org/10.1038/ng.3852
http://dx.doi.org/10.1038/ng.3852...
]; Parasteatoda tepidariorum (Araneomorphae) [4545. Schwager EE, Sharma PP, Clarke T, Leite DJ, Wierschin T, Pechmann M, et al. The house spider genome reveals an ancient whole-genome duplication during arachnid evolution. BMC Biol; 2017;15:1-27. ]; Pardosa pseudoannulata (Araneomorpha) [4646. Yu N, Li J, Liu M, Huang L, Bao H, Yang Z, et al. Genome sequencing and neurotoxin diversity of a wandering spider &lt;em&gt;Pardosa pseudoannulata&lt;/em&gt; (pond wolf spider). bioRxiv. 2019;747147. ]; Araneus ventricosus (Araneomorphae) [4747. Kono N, Nakamura H, Ohtoshi R, Moran DAP, Shinohara A, Yoshida Y, et al. Orb-weaving spider Araneus ventricosus genome elucidates the spidroin gene catalogue. Sci Rep [Internet]. Springer US; 2019;9:1-13. Available from: http://dx.doi.org/10.1038/s41598-019-44775-2
http://dx.doi.org/10.1038/s41598-019-447...
] and Dysdera silvatica (Araneomorphae) [4848. Sánchez-Herrero JF, Frías-López C, Escuer P, Hinojosa-Alvarez S, Arnedo MA, Sánchez-Gracia A, et al. The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates. Gigascience. 2019;8:1-9. ]. In addition, the 5000 arthropod genome initiative (i5K), which is committed to sequencing and analyzing 5000 high-priority arthropods [4949. Evans JD, Brown SJ, Hackett KJJ, Robinson G, Richards S, Lawson D, et al. The i5K initiative: Advancing arthropod genomics for knowledge, human health, agriculture, and the environment. J Hered. 2013;104:595-600. ], made available a draft of Loxosceles reclusa genome. Introns were described as present in all these spider genome analysis, but whether introns are present or not specifically in cystein-rich peptide toxin genes was not informed. As previously mentioned, latrotoxins found in the genome of the common house spider P. tepidariorum did present introns in their gene structure [4040. Gendreau KL, Haney RA, Schwager EE, Wierschin T, Stanke M, Richards S, et al. House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity. BMC Genomics; 2017;18:1-14. ].

To increase the body of information regarding the presence of introns in cysteine-rich peptide toxins, we attempted to use the available spider genomic data to analyze the structure of cysteine-rich peptides in toxin gene sequences. We did not find evidence of introns in the genomic sequences of A. geniculata (Mygalomorphae), corroborating the previous findings in the literature for this group. Regarding Araneomaphae spiders, the in-silico preliminary analysis did not show evidence of intron presence in cysteine-rich peptide toxins in two species (L. reclusa and D. sylvatica). However, despite most of the toxin gene sequences being apparently intronless, in the other four analyzed species (A. ventricosus, P. pseudoannulata, S. mimosarum and T. clavipes) the discontinuity of some alignments pointed to the possible presence of introns in a few toxins gene sequences.

Here, our experimental results showed that P. nigriventer do not present introns in the cysteine-rich peptide toxin sequences analyzed, corroborating most of the other experimental and in-silico findings that indicate the absence of introns in this class of spider toxins. However, since in the genomic data analysis we found evidence of the occurrence of introns in some toxin gene sequences from other Araneomorphae spiders, we cannot discard the possibility that other Phoneutria toxins can also present introns in their sequences. For instance, the co-existence of gene copies both with and without introns has been demonstrated for the elongation factor-1α (EF-1α) in Salticidae spiders from the genus Habbronatus [5050. Hedin MC, Maddison WP. Phylogenetic utility and evidence for multiple copies of Elongation factor-1α in the spider genus Habronattus (Araneae: Salticidae). Mol Biol Evol. 2001;18:1512-21. ].

Pineda et al. [1818. Pineda SS, Wilson D, Mattick JS, King GF. The lethal toxin from Australian funnel-web spiders is encoded by an intronless gene. PLoS One. 2012;7. ], in their work describing the gene structure of toxins from Australian funnel-web spiders, hypothesized that Mygalomorphae spiders had lost their toxin gene introns through their evolutionary history, whereas Araneomorphae kept them from a common ancestor. However, considering that different families of toxins have different features concerning the presence of introns and no consensus could be found among Araneomorphae, the evolutionary history of spider peptide neurotoxins seems to be more complex and remains to be further clarified.

Venom toxins are essential for spider survival and a high level of toxin expression in the venom glands is required. Gene expression increase may occur through the polyploidization of venom gland tissue cells [5151. Rasch EM, Connelly BA. Genome Size and Endonuclear DNA Replication in Spiders. 2005;214:209-14. ]. Intron loss is another possible mechanism for increasing gene expression, since it has been reported that highly expressed genes tend to lose introns more frequently and gain introns more rarely than genes with low expression levels [5252. Vinogradov AE. Compactness of human housekeeping genes: selection for economy or genomic design? 2004;20. ]. Furthermore, the absence of introns has also been related to an increase of mutation rates [5353. Minkevich IG, Patrushev LI.Genomic Noncoding Sequences and the Size of Eukaryotic Cell Nucleus as Important Factors of Gene Protection from Chemical Mutagens. 2007;33:474-7. ], which in turn may contribute to the high variability of toxin sequences and the emergence of new spider toxins [3838. Sachkova MY, Slavokhotova AA, Grishin E V., Vassilevski AA. Genes and evolution of two-domain toxins from lynx spider venom. FEBS Lett [Internet]. Federation of European Biochemical Societies; 2014;588:740-5. Available from: http://dx.doi.org/10.1016/j.febslet.2014.01.018
http://dx.doi.org/10.1016/j.febslet.2014...
]. On the other hand, the presence of introns in toxin genes may modulate the mutation rates of each exon separated by them [5454. Olivera BM, Walker C, Cartier GE, Hooper D, Santos AD, Schoenfeld R, et al. Speciation of cone snails and interspecific hyperdivergence of their venom peptides. Potential evolutionary significance of introns. Ann. N. Y. Acad. Sci. 1999. p. 223-37.], constituting a different mechanism for toxin evolution and diversification. Besides that, as most of spider toxin DNA sequences come from venom gland transcriptome analysis , the roles of alternative splicing, gene duplication and other regulatory controls in generating venom molecular diversity still need to be further studied.

Conclusions

In this work, we have investigated the structure of the genes encoding for three known sodium channel modulator toxins from P. nigriventer venom and the results indicate that P. nigriventer toxins do not contain introns in their genes sequences. However, since we also found evidence of toxin genes with and without introns in genomes of other Araneomorphae spiders, we cannot rule out the possibility of the presence of introns in other P. nigriventer toxin genes. This can only be confirmed after sequencing its genome.

As the majority of spider toxin DNA sequences come from venom glands transcriptomes studies, the gene structure of toxins from most spider toxins are still obscure. Thus, we believe that our results can contribute to future studies on understanding the mechanisms underlying spider venom molecular diversity, as well as the evolutionary aspects of spider toxins.

Abbreviations

gDNA: genomic DNA; ICK: inhibitor cysteine knot; PnTx1: μ-ctenitoxin-Pn1a; PnTx2-5: δ-ctenitoxin-Pn2c; PnTx4(5-5): ɣ-ctenitoxin-Pn1a; TAM: transcription-associated mutation.

References

  • 1. de Lima ME, Figueiredo SG, Matavel A, Nunes KP, da Silva CN, Almeida F de M, et al. Phoneutria nigriventer Venom and Toxins: A Review. Spider Venoms. 2015. p. 1-14.
  • 2. Peigneur S, de Lima ME, Tytgat J. Phoneutria nigriventer venom: A pharmacological treasure. Toxicon [Internet]. Elsevier Ltd; 2018;151:96-110. Available from: https://doi.org/10.1016/j.toxicon.2018.07.008
    » https://doi.org/10.1016/j.toxicon.2018.07.008
  • 3. Gomez MV, Kalapothakis E, Guatimosim C, Prado MAM. Phoneutria nigriventer Venom : A Cocktail of Toxins That Affect Ion Channels. Cell Mol Neurobiol. 2002;22:579-88.
  • 4. Diniz MRV, Paiva ALB, Guerra-Duarte C, Nishiyama MY, Mudadu MA, De Oliveira U, et al. An overview of Phoneutria nigriventer spider venom using combined transcriptomic and proteomic approaches. PLoS One. 2018;13:1-29.
  • 5. Cardoso FC, Pacífico LG, Carvalho DC, Victória JMN, Neves ALG, Chávez-Olórtegui C, et al. Molecular cloning and characterization of Phoneutria nigriventer toxins active on calcium channels. Toxicon. 2003;41:755-63.
  • 6. Penaforte CL, Prado VF, Prado MAM, Romano-Silva MA, Guimarães PEM, De Marco L, et al. Molecular cloning of cDNAs encoding insecticidal neurotoxic peptides from the spider Phoneutria nigriventer. Toxicon. 2000;38:1443-9.
  • 7. Matavel A, Cruz JS, Penaforte CL, Araújo DAM, Kalapothakis E, Prado VF, et al. Electrophysiological characterization and molecular identification of the Phoneutria nigriventer peptide toxin PnTx2-61. FEBS Lett. 2002;523:219-23.
  • 8. Kushmerick C, Kalapothakis E, Beirão PSL, Penaforte CL, Prado VF, Cruz JS, et al. Phoneutria nigriventer toxin Tx3-1 blocks A-type K+ currents controlling Ca2+ oscillation frequency in GH3 cells. J Neurochem. 1999;72:1472-81.
  • 9. Kalapothakis E, Penaforte CL, Leão RM, Cruz JS, Prado VF, Cordeiro MN, et al. Cloning, cDNA sequence analysis and patch clamp studies of a toxin from the venom of the armed spider (Phoneutria nigriventer). Toxicon. 1998;36:1971-80.
  • 10. Diniz MR V, Paine MJI, Diniz CR, Theakston RDG, Crampton JM. Sequence of the cDNA coding for the lethal neurotoxin Tx1 from the Brazilian “armed” spider Phoneutria nigriventer predicts the synthesis and processing of a preprotoxin. J Biol Chem. 1993;268:15340-2.
  • 11. King GF. The wonderful world of spiders: Preface to the special Toxicon issue on spider venoms. Toxicon. 2004;43:471-5.
  • 12. Garrison NL, Rodriguez J, Agnarsson I, Coddington JA, Griswold CE, Hamilton CA, et al. Spider phylogenomics: Untangling the Spider Tree of Life. PeerJ. 2016;2016.
  • 13. Fernández R, Kallal RJ, Dimitrov D, Ballesteros JA, Arnedo MA, Giribet G, et al. Phylogenomics, Diversification Dynamics, and Comparative Transcriptomics across the Spider Tree of Life. Curr Biol. 2018;28:1489-1497.e5.
  • 14. Krapcho KJ, Kral RM, Vanwagenen BC, Eppler KG, Morgan TK. Characterization and cloning of insecticidal peptides from the primitive weaving spider Diguetia canities. Insect Biochem Mol Biol. 1995;25:991-1000.
  • 15. Qiao P, Zuo XP, Chai ZF, Ji YH. The cDNA and genomic DNA organization of a novel toxin SHT-I from spider Ornithoctonus huwena. Acta Biochim Biophys Sin (Shanghai). 2004;36:656-60.
  • 16. Jiang L, Chen J, Peng L, Zhang Y, Xiong X, Liang S. Genomic organization and cloning of novel genes encoding toxin-like peptides of three superfamilies from the spider Orinithoctonus huwena. Peptides. 2008;29:1679-84.
  • 17. Tang X, Zhang Y, Hu W, Xu D, Tao H, Yang X, et al. Molecular diversification of peptide toxins from the tarantula haplopelma hainanum (Ornithoctonus hainana) venom based on transcriptomic, peptidomic, and genomic analyses. J Proteome Res. 2010;9:2550-64.
  • 18. Pineda SS, Wilson D, Mattick JS, King GF. The lethal toxin from Australian funnel-web spiders is encoded by an intronless gene. PLoS One. 2012;7.
  • 19. Martin-Moutot N, Haro L De, Santos RG Dos, Mori Y, Seagar M. Phoneutria nigriventer ω-Phonetoxin IIA: A new tool for anti-calcium channel autoantibody assays in Lambert-Eaton myasthenic syndrome. Neurobiol Dis. 2006;22:57-63.
  • 20. Silva AO, Peigneur S, Diniz MRV, Tytgat J, Beirão PSL. Inhibitory effect of the recombinant Phoneutria nigriventer Tx1 toxin on voltage-gated sodium channels. Biochimie [Internet]. Elsevier Masson SAS; 2012;94:2756-63. Available from: http://dx.doi.org/10.1016/j.biochi.2012.08.016
    » http://dx.doi.org/10.1016/j.biochi.2012.08.016
  • 21. Diniz MRV, Theakston RDG, Crampton JM, Nascimento Cordeiro M do, Pimenta AMC, De Lima ME, et al. Functional expression and purification of recombinant Tx1, a sodium channel blocker neurotoxin from the venom of the Brazilian “armed” spider, Phoneutria nigriventer. Protein Expr Purif. 2006;50:18-24.
  • 22. Matavel A, Fleury C, Oliveira LC, Molina F, De Lima ME, Cruz JS, et al. Structure and activity analysis of two spider toxins that alter sodium channel inactivation kinetics. Biochemistry. 2009;48:3078-88.
  • 23. Cordeiro M do N, Diniz CR, do Carmo Valentim A, von Eickstedt VRD, Gilroy J, Richardson M. The purification and amino acid sequences of four Tx2 neurotoxins from the venom of the Brazilian “armed” spider Phoneutria nigriventer (Keys). FEBS Lett. 1992;310:153-6.
  • 24. Yonamine CM, Troncone LRP, Camillo MAP. Blockade of neuronal nitric oxide synthase abolishes the toxic effects of Tx2-5, a lethal Phoneutria nigriventer spider toxin. Toxicon . 2004;44:169-72.
  • 25. De Figueiredo SG, De Lima ME, Cordeiro MN, Diniz CR , Patten D, Halliwell RF, et al. Purification and amino acid sequence of a highly insecticidal toxin from the venom of the Brazilian spider Phoneutria nigriventer which inhibits NMDA-evoked currents in rat hippocampal neurones. Toxicon . 2001;39:309-17.
  • 26. Paiva ALB, Matavel A, Peigneur S, Cordeiro M do N, Tytgat J, Diniz MR V, et al. Differential effects of the recombinant toxin PnTx4(5-5) from the spider Phoneutria nigriventer on mammalian and insect sodium channels. Biochimie. 2016;121:326-35.
  • 27. Silva FR, Batista EML, Gomez M V., Kushmerick C, Da Silva JF, Cordeiro MN, et al. The Phoneutria nigriventer spider toxin, PnTx4-5-5, promotes neuronal survival by blocking NMDA receptors. Toxicon [Internet]. Elsevier Ltd; 2016;112:16-21. Available from: http://dx.doi.org/10.1016/j.toxicon.2016.01.056
    » http://dx.doi.org/10.1016/j.toxicon.2016.01.056
  • 28. Oliveira CFB, Alves DP, Emerich BL, Figueiredo SG de, Cordeiro M do N, Borges MH, et al. Antinociceptive effect of PnTx4(5-5), a peptide from Phoneutria nigriventer spider venom, in rat models and the involvement of glutamatergic system. J Venom Anim Toxins Incl Trop Dis. 2019;25:1-9.
  • 29. Killeen AA, Fan H, Gulley ML. DNA Extraction from Fresh or Frozen Tissues. Mol Pathol Protoc. 2003;49:5-10.
  • 30. Pineda SS, Chaumeil P, Kunert A, Thang MWC, Li L, Nuhn M, et al. Spider toxin database ArachnoServer 3.0: an online resource for automated discovery , analysis and annotation of spider toxins. Bioinformatics. 2017;7-8.
  • 31. Postic G, Gracy J, Périn C, Chiche L, Gelly JC. KNOTTIN: The database of inhibitor cystine knot scaffold after 10 years, toward a systematic structure modeling. Nucleic Acids Res. 2018;46:D454-8.
  • 32. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, et al. GenBank. Nucleic Acids Res. 2013;41:36-42.
  • 33. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403-10.
  • 34. Vassilevski AA, Kozlov SA, Grishin E V. Molecular diversity of spider venom. Biochem. 2009;74:1505-34.
  • 35. Garb JE, Sharma PP, Ayoub NA. Recent progress and prospects for advancing arachnid genomics. Curr Opin Insect Sci [Internet]. Elsevier Inc; 2018;25:51-7. Available from: http://dx.doi.org/10.1016/j.cois.2017.11.005
    » http://dx.doi.org/10.1016/j.cois.2017.11.005
  • 36. Aguilera A. The connection between transcription and genomic instability. EMBO J. 2002;21:195-201.
  • 37. Sachkova MY, Slavokhotova AA, Grishin E V., Vassilevski AA. Structure of the yellow sac spider Cheiracanthium punctorium genes provides clues to evolution of insecticidal two-domain knottin toxins. Insect Mol Biol. 2014;23:527-38.
  • 38. Sachkova MY, Slavokhotova AA, Grishin E V., Vassilevski AA. Genes and evolution of two-domain toxins from lynx spider venom. FEBS Lett [Internet]. Federation of European Biochemical Societies; 2014;588:740-5. Available from: http://dx.doi.org/10.1016/j.febslet.2014.01.018
    » http://dx.doi.org/10.1016/j.febslet.2014.01.018
  • 39. Danilevich VN, Grishin EV. The Genes Encoding Black Widow Spider Neurotoxins are Intronless. Russ J Bioorganic Chem. 2000;26:838-43.
  • 40. Gendreau KL, Haney RA, Schwager EE, Wierschin T, Stanke M, Richards S, et al. House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity. BMC Genomics; 2017;18:1-14.
  • 41. Bhere KV, Haney RA, Ayoub NA, Garb JE. Gene structure, regulatory control, and evolution of black widow venom latrotoxins. FEBS Lett. 2014;588:3891-7.
  • 42. Binford GJ, Cordes MHJ, Wells MA. Sphingomyelinase D from venoms of Loxosceles spiders: Evolutionary insights from cDNA sequences and gene structure. Toxicon . 2005;45:547-60.
  • 43. Sanggaard KW, Bechsgaard JS, Fang X, Duan J, Dyrlund TF, Gupta V, et al. Spider genomes provide insight into composition and evolution of venom and silk. Nat Commun. 2014;5.
  • 44. Babb PL, Lahens NF, Correa-Garhwal SM, Nicholson DN, Kim EJ, Hogenesch JB, et al. The Nephila clavipes genome highlights the diversity of spider silk genes and their complex expression. Nat Genet [Internet]. Nature Publishing Group; 2017;49:895-903. Available from: http://dx.doi.org/10.1038/ng.3852
    » http://dx.doi.org/10.1038/ng.3852
  • 45. Schwager EE, Sharma PP, Clarke T, Leite DJ, Wierschin T, Pechmann M, et al. The house spider genome reveals an ancient whole-genome duplication during arachnid evolution. BMC Biol; 2017;15:1-27.
  • 46. Yu N, Li J, Liu M, Huang L, Bao H, Yang Z, et al. Genome sequencing and neurotoxin diversity of a wandering spider &lt;em&gt;Pardosa pseudoannulata&lt;/em&gt; (pond wolf spider). bioRxiv. 2019;747147.
  • 47. Kono N, Nakamura H, Ohtoshi R, Moran DAP, Shinohara A, Yoshida Y, et al. Orb-weaving spider Araneus ventricosus genome elucidates the spidroin gene catalogue. Sci Rep [Internet]. Springer US; 2019;9:1-13. Available from: http://dx.doi.org/10.1038/s41598-019-44775-2
    » http://dx.doi.org/10.1038/s41598-019-44775-2
  • 48. Sánchez-Herrero JF, Frías-López C, Escuer P, Hinojosa-Alvarez S, Arnedo MA, Sánchez-Gracia A, et al. The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates. Gigascience. 2019;8:1-9.
  • 49. Evans JD, Brown SJ, Hackett KJJ, Robinson G, Richards S, Lawson D, et al. The i5K initiative: Advancing arthropod genomics for knowledge, human health, agriculture, and the environment. J Hered. 2013;104:595-600.
  • 50. Hedin MC, Maddison WP. Phylogenetic utility and evidence for multiple copies of Elongation factor-1α in the spider genus Habronattus (Araneae: Salticidae). Mol Biol Evol. 2001;18:1512-21.
  • 51. Rasch EM, Connelly BA. Genome Size and Endonuclear DNA Replication in Spiders. 2005;214:209-14.
  • 52. Vinogradov AE. Compactness of human housekeeping genes: selection for economy or genomic design? 2004;20.
  • 53. Minkevich IG, Patrushev LI.Genomic Noncoding Sequences and the Size of Eukaryotic Cell Nucleus as Important Factors of Gene Protection from Chemical Mutagens. 2007;33:474-7.
  • 54. Olivera BM, Walker C, Cartier GE, Hooper D, Santos AD, Schoenfeld R, et al. Speciation of cone snails and interspecific hyperdivergence of their venom peptides. Potential evolutionary significance of introns. Ann. N. Y. Acad. Sci. 1999. p. 223-37.
  • Availability of data and materials

    All data generated or analyzed during this study are included in this article.
  • Funding

    Not applicable.
  • Ethics approval

    The authors declare that the spiders used in this work were collected with approval of the Brazilian Institute of the Environment and Renewable Natural Resources (IBAMA). The project was also registered in National System for Management of Genetic Heritage and Associated Traditional Knowledge (SisGen): #A26E945.
  • Consent for publication

    Not applicable.

Publication Dates

  • Publication in this collection
    30 Apr 2020
  • Date of issue
    2020

History

  • Received
    21 Oct 2019
  • Accepted
    18 Mar 2020
Centro de Estudos de Venenos e Animais Peçonhentos (CEVAP/UNESP) Av. Universitária, 3780, Fazenda Lageado, Botucatu, SP, CEP 18610-034, Brasil, Tel.: +55 14 3880-7693 - Botucatu - SP - Brazil
E-mail: editorial.jvatitd@unesp.br