Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom
© Leão et al; licensee BioMed Central Ltd. 2009
Received: 13 October 2008
Accepted: 16 March 2009
Published: 16 March 2009
Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization.
A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A2 (PLA2s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA2) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens.
Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.
The coral snake (genus Micrurus) is the most abundant, diverse and representative member of the family Elapidae in the New World. It has a wide geographic distribution which covers the southwest United States, Central America, and southern Argentina . Compared to the family Viperidae, the numbers of accidents caused by the coral snake are not great. Coral snakes are not aggressive and only attack when threatened. Still, when an accident occurs, the symptoms are usually severe, leading to death by asphyxia after only 5 or 6 hours due to strong neurotoxic effects .
In Brazil, cases of envenoming by coral snakes are caused mainly by Micrurus corallinus and Micrurus frontalis, species inhabiting in highly populated areas in the Central, South and Southeast regions. Many of their characteristics, such as ophiophagous diet, fossorial habit and living in tropical latitudes, make it difficult to obtain and keep them in captivity. This limitation in maintenance, the small size of their venom glands and, consequently, low production of venom have been the major factors hindering the production of Brazilian anti-elapidic serum. Furthermore, the Butantan Institute, in Sao Paulo, utilizes almost all of the venom obtained to produce the anti-elapidic serum, limiting biochemical studies . In fact, the quantity of venom available for the generation of serum is not enough to supply the national needs. While the incidence of accidents is small when compared to that for other genera, the wide geographic dispersion of Micrurus and the lethality of its bite require the serum to be distributed all over the country, raising its demand.
Nowadays, the transcriptomic analysis of venom glands to obtain a general profile of the toxins composing the venom is a common experimental approach, applied especially in snakes of the families Viperidae [4–8] and Colubridae . Nevertheless, there are no systematic transcriptome reports of this kind for the family Elapidae, neither from American coral snakes nor from African-Asian species. Moreover, considering the difficulties in obtaining the venom even for antiserum production, a complete set of the most abundant cDNAs from the venom glands of coral snake species could open possibilities for alternative ways of immunization. One of them is the production of recombinant proteins for immunization, which is an obvious choice, since horse hyper-immunization demands large amounts of proteins and requires robust expression systems, such as those using E. coli. However, prokaryotic systems fail to express complex disulfide-bonded proteins, as in the case of snake toxins. Several of our attempts to produce M. corallinus recombinant toxins involved complicated refolding and eventually resulted in time-consuming procedures [10, 11]. A more elegant approach is the DNA immunization with plasmids expressing exogenous DNA in animal cells, thus stimulating immunological responses . For snake toxins, Harrison and colleagues  produced high titers of antibodies after DNA immunization using snake toxin cDNAs, showing a potential systemic treatment of injuries caused by Viperidae snake bites. They then immunized mice with a synthetic chimeric DNA containing seven different epitopes from snake venom metalloproteases expressed in the venom glands of the viper Echis ocellatus .
The dataset of Micrurus corallinus transcriptome presented here aims not only to describe the putative toxins expressed by this sinuous Elapidae snake, but also to be a source of potential antigens for an alternative production of antiserum. Therefore, some cDNAs meeting eligible criteria, such as abundance, were selected and injected into mice for a preliminary evaluation of their capacity to produce an immune response as DNA immunogens.
Results and Discussion
The EST databank of M. corallinus contains 1438 sequences grouped into 611 clusters using CAP3 software with 98% identity. Since the three-finger proteins (3FTx) present in M. corallinus venom are very diverse and a few substitutions of residues could change the clade in which they are classified, some of the clusters were inspected and manually regrouped to avoid joining different EST products in the same cluster.
The 611 clusters of contiguous sequences include 483 formed by individual clones (singlets) and 128 clusters with two or more sequences (contigs). After the application of a filter against sequences of ribosomal RNA, mitochondrial DNA and other contaminants, 5% of the clones were removed for the subsequent statistical analysis.
The analysis using Blast revealed that 485 clusters (corresponding to 1250 clones, or 92% of the total) produced significant similarities with the database (e-value < 10-05), and that 103 clusters (107 clones or 8% of total) could not be identified. The information obtained was separated into two general functions of venom gland transcripts: toxins (toxic functions) and other transcripts involved in all cellular functions (non-toxins). Toxin transcripts totaled 46% of all clones in 93 clusters. Non-toxins also totaled 46% of all clones but in 392 clusters. Therefore, the non-toxins are less redundant (1.6 clone/cluster) than toxins (6.7 clones/cluster), which means that each toxin is, on average, four times more transcribed than a cellular transcript.
Analysis of the putative toxins
Three-Finger Toxins (3FTx)
The 3FTx are the most abundant group of possible toxins, comprising the majority of expressed sequences (~52% of the total toxin clones). These proteins are small, with less than 80 amino acids and generally have four or five disulfide bonds. They have very different pharmacological effects, but are generally associated with post-synaptic blockade.
The analysis of the 3FTx revealed that 326 clones (~24% of total) could be grouped into 55 clusters. Many clusters had to be recreated manually after the initial grouping made by the CAP3 program, because it was noted that some initial contigs were actually grouping different variants. Therefore, each initial contig was inspected to identify the SNPs (single nucleotide polymorphisms) common to more than one sequence in the contig. When these common SNPs were found, the quality of their base pairs at each point was evaluated, and the respective sequences were removed and rearranged in a new cluster. It was not possible to determine if a single animal shows all these isoforms or if they are expressed by distinct animals, since the cDNA library was constructed using 10 coral snakes.
Two groups are much more abundant than the rest, representing the majority of the 3FTx (Figure 2). Both contain sequences that are very similar to α-neurotoxin isoforms previously described in Micrurus corallinus . The first group, named here Atg1 (selected as antigen 1, as discussed later), was similar to a cDNA called neurotoxin homolog 8 (nxh8), previously isolated by Ho et al. (1995), which is a peculiar form that differs from most 3FTx. After a 21-amino acid signal peptide, the mature form consists of 65 amino acids, along with the possible presence of an extra disulfide bond in the first loop. Most short 3FTx have only four very conserved disulfide bridges but few unrelated toxins show this extra pair of cysteines (Figure 2). The second group, named Atg2 (selected as antigen 2, as discussed later), refers to more typical 3FTx, the neurotoxin homologs 7/3/1 (nxh7, nxh3, nxh1 cDNAs), as previously described . Some of the clusters showed differences in the signal peptide region and some amino acids in the mature protein region. The mature protein consists of 57 amino acids. Comparing these two groups, Atg1 displays just 42% similarity with the Atg2 group. The rest of the 3FTx, some of them also highly transcribed, represent new proteins with an identity value of no more than 50% to the sequences of 3FTx in the databanks (Figure 2).
Phospholipase A2 (PLA2)
With 120 amino acids and 6 or 7 disulfide bonds, PLA2 are ubiquitous enzymes in venoms, composed of alpha-helix and some beta-sheet structures . Some PLA2s are pre-synaptic neurotoxins or β-neurotoxins with specific targets in neuromuscular junctions , promoting an inhibition in the release of neuromediators. Some others act like a myotoxin that depolarizes muscular cell membranes, leading to necrosis. PLA2 may exert neurotoxic actions by binding to nerve cell membranes and catalyzing phospholipid hydrolysis, with the production of lysophospholipids and free fatty acids .
The most abundant cluster obtained was MCOR0036C, with 196 clones. This cluster shows 100% identity with the PLA2 sequence of Micrurus corallinus from the database (AY157830). This PLA2 was expressed in E. coli and immunologically characterized . It has, in position 49, one aspartate (Asp49) intimately involved in Ca2+ binding, having an essential function in the catalytic activity of this enzyme. This residue is highly conserved in catalytically functional PLA2, but in many Viperidae myotoxic PLA2s, it is substituted by a lysine residue (Lys49) . No such latter form was detected in this transcriptome.
C-type lectins are non-enzymatic proteins present in the venom of snakes of all families. Like plant lectins, some of the venom lectins are able to bind to carbohydrates. In the presence of Ca2+, the C-type lectins begin several biological processes, such as adhesion, agglutination, endocytosis and neutralization of the pathogen. The lectins could be divided into seven different groups according to their structural characteristics . The C-type lectins that have lost their capability of binding to carbohydrates are called lectin-like, because they still maintain the structural characteristics in common with true lectins. Venom lectins may act as agonists or sometimes antagonists in platelet aggregation and affect thrombosis and homeostasis by activating and inhibiting specific receptors in platelet membranes .
The C-type lectins found here were also abundant toxins. Sixty-five clones show similarity with these proteins and were grouped into 13 clusters, representing approximately 5% of the total transcripts expressed in the tissue. Many of the clusters correspond to full-lengths ORFs that could be aligned with other C-type lectins from venom [see Additional file 2].
C-type lectins are ubiquitous components of the venom, but are very diversified in structure-function terms, and thus, it is interesting to understand how the C-type lectins from M. corallinus fit in the phylogeny of lectins. A Bayesian phylogenetic analysis revealed that the proteins of snake venom set in a unique group, distinct from the animal physiologic lectins [see Additional file 2]. Within the group of venom lectins, it is possible to see a clear distinction between the true C-type lectins (carbohydrate-binding) and their structural homologs, the lectin-like (probabilities 0.81 and 1.0, respectively). All the M. corallinus lectins are arranged inside the true C-type lectin group, together with the other elapids lectins. The only exception is the cluster MCOR0067C (clearly distinct from the others in alignment), which is more similar, in terms of identity, to the short physiologic proteins related to lithostatin than to the venom lectins (true or not). Recently, the same kind of sequence was observed in the B. insularis transcriptome (BINS0004C). Thus, the finding of another lectin with those characteristics supports the possibility that this lectin could be a third group of venom lectins or that it has a specialized function in the venom gland.
For the organization inside the branch of true lectins, the majority of clusters are grouped together with the lectin of Lapemis hardwick, an elapid of the Hydrophiinae subfamily (sea snakes). Only one cluster, MCOR0090C, is grouped with the real lectins of other elapids and viperids. Either way, these components seem to be well diversified in M. corallinus.
Natriuretic Peptide (NP)
NPs have been described in snake venoms, where their precursors are organized in different ways. In the Elapidae, a 110- to 140-amino acid ORF contains a signal peptide, a short propeptide and the NP with a C-terminal extension resembling that of vertebrate atrial and brain NPs (A/BNPs) . In contrast, in the Viperidae (Crotalinae subfamily), the precursor is a more complex molecule with a 180- to 270-amino acid ORF containing a signal peptide, a region with a variable number of the pharmacologically relevant bradykinin-potentiating peptides (BPPs) (from one to seven in different snake species) followed by a long intervening linker sequence and the C-type NP (CNP), without a C-terminal extension . In the venom of the dipsadid snake Philodryas olfersii, it is similar to the Viperidae BPP/CNP precursor, including the signal peptide, the linker and the CNP (without C-terminal extension), but lacking the entire BPPs region . The BPPs are important inhibitors of the angiotensin-converting enzyme that served as basis for the design of the anti-hypertensive drug captopril.
Here we found 5 clusters, totaling 14 clones. The cDNA representing the cluster MCOR0001C (6 clones) is very similar to the one previously characterized , except for the absence of the extra cysteine observed in that cDNA. In fact, resequencing the original clone from Ho and colleagues, it was clear that the sequence described earlier was incorrect (data not shown), and that the correct one corresponds to that of the MCOR0001C cluster. All the natriuretic peptides described represent a structural ring with 17 amino acids formed by an expanding disulfide bond with a few amino acids in the N-terminal. In the C-terminal region, the type B and C of mammalian NPs have 4 to 7 amino acids. The natriuretic peptide of M. corallinus has some peculiarities: a) the peptide has an unusually longer C-terminal extension; and b) this possible natriuretic peptide is flanked by 2 homologous regions with unknown function .
It is difficult to determine the importance of such molecule that is abundant, with part of the pro-domain and the C-terminal end, but not showing the NP. One clue could be the presence, even in the deleted form, of one of the repetitive regions, which has been suggested to have typical elements of active biological peptides . This precursor may be chosen to produce a new peptide based on this segment. It is interesting to note that the BPPs/CNPs precursors, whose structures have been known for a while, are now revealing new active peptides in their pro-peptide regions. Graham and colleagues, in 2005 , found bradykinin inhibitory peptides in the same precursor of the BPPs. Recently, Wagstaff and colleagues described that the conserved spacing region from some precursors generates peptides (EKW) capable of inhibiting the metalloprotease in the same venom, in a supposed mechanism of protection against proteolysis in the gland .
L-Amino Acid Oxidase
These are larger enzymes of 58 kDa and act basically on L-amino acids, converting them to keto acids and generating H2O2, which seems to contribute to the inhibition of platelet aggregation and to cause cellular injury . It was not possible to obtain a full-length clone from our dataset, probably due to the very long cDNA (2.8 kb) . Two clusters, MCOR0200S and MCOR0263S, each one containing just one clone (singlet), were found in the M. corallinus databank. L-amino acid oxidase has already been reported in elapids (ex.: AAY89682), but not in Micrurus.
In snake venoms, the scaffold of Kunitz-type serine protease inhibitors occurs as two functionally distinct molecules: the non-neurotoxic serine protease inhibitor, and the neurotoxic ones, also called dendrotoxins, lacking major inhibitory action on proteases. The former ones are capable of inhibiting trypsin and chymotrypsin, being common in both the Viperidae and Elapidae . The dendrotoxins are present in the venom of some elapids , such as the black mamba (Dendroaspis) and the kraits (Bungarus). They exert their toxic action by increasing the release of neurotransmitters, such as acetylcholine, from pre-synaptic membranes of cholinergic synapses. The specific targets of dendrotoxins are high-voltage-activated Ca2+ channels (L-, N-, and P-type) in neuronal cell membranes, which are blocked .
Two low expressed clusters, MCOR0110C and MCOR0611S, were found in this database. Both seem to be more similar to the inhibitors than to the neurotoxins, matching several Viperidae inhibitors, although the degree of conservation between these two types is low, and further analysis is needed in order to define their exact function.
Venom metalloproteases are normally grouped according to the presence of a specific domain. All of them have a proteolytic domain characterized by the presence of a catalytic site conserving the motif HEXXHXXGXXH, responsible for the binding of a metal ion, generally Zn2+. After the catalytic domain, there may be a disintegrin domain, capable of binding to platelet integrins or to endothelial cells. This region could contain the RGD motif (RGD disintegrins) or have its sequence changed to another, such as ECD (non-RGD disintegrins), which would change its specificity. The C-terminal region of non-RGD disintegrins could also have a cysteine-rich domain with unknown function. These enzymes are very abundant in the family Viperidae but are not so common and diversified in the Elapidae.
The metalloprotease found here is represented by one partial cluster, MCOR0063C, containing just 4 clones (1123 bp). It is similar to cobrin found in Naja naja, a species of the family Elapidae from the Old World. It is an example of a P-III type of metalloprotease, i.e., the sequence shows a non-RGD disintegrin type followed by a cysteine-rich domain .
Snake venom serine proteases act on elements of the coagulation cascade and are believed to occur as trypsin-like enzymes, widely described in the Viperidae, or as the more complex factor X-like proteins found in some terrestrial Elapidae . Recently, some trypsin-like enzymes were described in the Elapidae, e.g., ABN72544, showing that the factor X-like is not the only serine protease fold in this family.
In the transcriptome of M. corallinus, we did not find either type, but only a partial singlet cluster (MCOR0160S), showing more similarity to non-venom serine proteases from mammals (~52% identity), than to snake ones (~39%). The main hits are with prostasins from basal mammals such as platypus (Monotremata) and opossum (Metatheria). This is a channel-activating protease involved in the growth of epithelial cells. The portion matching MCOR0160S corresponds to a part of the peptidase S1 domain, indicating that proteolytic activity should be present. Since the sequence is partial, it is not possible to define the presence of extra domains, such as the transmembrane or the propeptide C-terminal end. Due to the ubiquitous presence of serine protease in snake venoms, we include this one among the possible toxins from M. corallinus, although further confirmation is still needed.
Like their homologs in mammals, responsible for regulating neuronal differentiation, NGFs (neurotrophic factor family) of venom act as low potential agonists of TrkA receptors (tyrosine-kinase), competing with endogenous NGF in binding to the receptor and influencing the development of cholinergic neurons. Besides inducing a typical growth of fibers in cell culture, snake venom and human NGFs also exhibit non-neuronal effects, such as the induction of plasma overflow or the release of histamines from blood cells . The effect of these factors in envenoming probably renders the bite site more susceptible to other components of the venom and distribute the substances that are difficult to infiltrate the target tissue .
This component was also found in our library, although less abundant than in other transcriptomes. Probably due to the length of the sequence, the only cluster obtained (MCOR0149S) is partial.
WAP proteins (whey acidic proteins), or waprin or nawaprin (described first in Naja nigrocollis), play a possible role as a protease inhibitor. It is similar to elafin (elastase inhibitor of human leukocytes), clatrin-like (calcium transport inhibitor) and other extra cellular protease inhibitors . The WAP domain generally consists of 50 amino acid residues, with 8 conserved cysteine residues forming 4 disulfide bridges. Although the cysteines residues are conserved, the inter-cysteine segments are totally different from those of other WAP family members. Among the ESTs of M. corallinus, we obtained only one cluster, MCOR0526S, containing just one clone coding for a WAP protein. It is a partial sequence without signal peptide and part of the mature protein.
General considerations on M. corallinus transcriptome
Whereas other snake families have been well investigated through the transcriptomic approach, Elapidae cDNA sequences were obtained only from timely cloning efforts. The toxins in the M. corallinus transcriptome presented here represent approximately 46% (625 clones grouped in 93 clusters) of all the information. This very high expression of a functionally related group of protein is expected since this tissue is specialized in toxin synthesis, and is similar to the results observed from other snake families such as the Viperidae [5–7] and Colubridae .
The general profile of toxins revealed a complex diversity of toxins (10 different classes of toxins). Protein profiles of Micrurus venoms (e.g., [10, 41]) usually show a large number of proteins of around 6–24 kDa (mainly 3FTx with various activities and PLA2) with minor bands at higher molecular mass. Although PLA2 and 3FTx make up the vast majority of ESTs, 8 other classes were found, including proteins with a predicted higher molecular mass. All these classes have been previously observed in at least one species of the Elapidae, but as far as we know, not in the same species or genus. Olamendi-Portugal and colleagues  investigated the proteome set of a Central America coral snake (Micrurus surinamensis) and found that the main group of toxins corresponded to 3FTx. However, they also identified PLA2, other low-molecular mass proteins and L-amino acid oxidase. Together, these data suggest that the complexity of venom components in Micrurus is readily comparable with that of the family Viperidae. Of course, other datasets are needed to extrapolate these Micrurus profiles.
Nevertheless, the neurotoxic action of pre- and post-synaptic toxins still predominate. Considering this action, M. corallinus venom is frequently compared and contrasted with that of M. frontalis, because both are the major species in terms of occurrence and accidents in Brazil, but the latter shows an almost exclusively post-synaptic venom, whereas the former has both pharmacological actions . The presence of elevated levels of PLA2 messengers in M. corallinus, almost equal to the 3FTx levels, is in agreement with this, suggesting that the responsibility for this action may be associated with the level of this component. Unfortunately, we still do not have the M. frontallis transcriptome data in order to compare them.
Besides the existence of diverse classes of toxins, it is worth noting the large number of different cDNAs in each class. 3FTx, for example, are represented by 52 differently assembled cDNAs. It is true that many forms may result from the assembly of non-overlapping segments or correspond to cDNAs with few nucleotide polymorphisms (with many silent substitutions), which may be associated with individual variation resulting from the pool of the venom glands used. Nevertheless, exclusively considering the amino acid sequences from the full-length 3FTx, 10 of the 28 full-length 3FTx proteins observed are less than 80% conserved relative to any other, and 7 of them are less than 50% conserved relative to any other. If we exclude the signal peptide from the analysis, percentages are obtained that reflect even more diversity. This indicates the great variability of these sequences in a single species. It is surprising that a species with a very specialized diet (feeding mainly on other small snakes) shows this diversity. However, in the proteome of M. surinamensi various subtypes of 3FTx are present. Five of the most abundant, including short- and long-chain 3FTx, were functionally characterized and shown to be acetylcholine receptor with different affinities and peculiar specificities. M. surinamensi is also a diet specialist, feeding on fish, but many of their toxins seem to be active also in mammals . In a family of short toxins subjected to high evolutionary pressures on structure-function relationships, this may represent very different kinds of activities, such as several receptor specificities.
Identification of putative antigenic candidates
Considering the problems in the production of M. corallinus antiserum for humans due to the difficulties in obtaining its venom, alternative ways to obtain their antigens are important. The direct use of M. corallinus recombinant proteins as immunogens have been tested before [10, 11] and although they were shown to produce immune responses indicating the recognition of the native proteins, very complicated steps were required for protein refolding and to obtain active toxins. Therefore, we decided to evaluate the feasibility of using cDNAs derived from this transcriptome initiative as immunogens for producing an antiserum. By analyzing the dataset, we confirmed the predominant toxins in this venom and suggested some proteins that could represent good antigenic candidates.
The proposed candidates should code for proteins that are abundant in the venom. Since 3FTx and PLA2 account for 85% of toxin ESTs, they were chosen for this first evaluation of immune response through DNA immunization. The 3FTx candidates were selected based on the abundance of each transcript, with aim of choosing the toxins that would be abundant in the venom. Among them, we considered only the most divergent forms, indicated by a neighbor-joining tree, avoiding the selection of closely related isoforms. Thus, we aimed to select diverse abundant (potentially important) cDNAs of 3FTx, which are the major toxins, and one representative of PLA2, rather than to choose the most immunogenic ones, as used elsewhere .
The 35 clusters of 3FTx with full length sequences were aligned and self-compared, giving rise to groups of clusters showing more similarity (Figure 2) [see Additional file 1]. One group containing 5 clusters and 45 clones was represented by MCOR0604C, here named Atg1, and is similar to the nxh8 (AJ344067). The second group (Atg2) is the most abundant 3FTx set, containing 117 clones, represented by the MCOR0599C cluster, which is similar to the neurotoxin homologs 7/3/1 (nxh1, nxh3, nxh7). Among the remaining 3FTx, two subgroups were defined: Atg3 included 29 clones (represented by the MCOR0039C cluster) and Atg4 had 27 clones (represented by the MCOR0100C cluster).
The obvious fifth candidate (Atg5) was the MCOR0036C cluster, representing the most abundant cluster obtained (196 clones), corresponding to the putative M. corallinus PLA2.
DNA immunization and immune response
The five selected candidates were subcloned in the mammalian vector, pSecTag2A and the DNA injected in mice. In parallel, some constructions of the same cDNAs cloned in prokaryotic expression systems were used to produce recombinant antigens (data not shown) for the evaluation of the antisera.
To evaluate the antiserum raised against the two other 3FTx antigens, ELISA plates were coated with the 3FTx recombinant proteins available (Atg1 and Atg2) and exposed to the sera obtained from Atg3 and Atg4 DNA immunization. Therefore, a positive response would indicate that there is cross reaction between the selected antigenic candidates. In this test, the sera of mice immunized with Atg3 or Atg4 showed better recognition of the recombinant proteins corresponding to Atg2 than those corresponding to the Atg1 [see Additional file 3]. This indicates that Atg1 has less conserved epitopes with the other antigens analyzed, which is in accordance with its sequence properties that showed MCOR0604C (nxh8) as a truly distinctive type of 3FTx. Besides, the low cross-reaction titers even with the most conserved Atg2 suggests that the distinct 3FTx isoforms selected, alone, may not elicit a response sufficient to neutralize the venom. For this reason, a set of these toxins should be used together to ensure complete protection.
The immunological response against Atg5 (PLA2) only showed low reaction titers (1:160), but was also capable of recognizing the recombinant proteins in Western blots (Figure 5).
Therefore, our results indicate that immunization with DNA is able to generate an immunological response against recombinant proteins. Meanwhile, the venom was poorly recognized by these sera in ELISA. Perhaps, optimizing the immunization strategy would improve the recognition. For practical reasons, in this first evaluation we utilized here a simple intra-muscular protocol, by injecting the DNA in mice in their anterior tibia muscle or quadriceps. Epidermal injection should improve the immune response, as previously demonstrated , particularly if associated with biobalistic protocols. This procedure has already been shown in several studies to be the most efficient to induce an immunological response of the Th2 type (IL-4, IL-5, IL10), efficient to generate more antibodies, which is desirable for hyperimmune serum production. Nevertheless, the reaction titers of anti-toxin sera are not directly related to its efficiency; if some important epitopes were neutralized, the general toxicity of the venom could be blocked even with low-titer serum, as clearly shown before .
Years of DNA immunization studies have shown that, unlike with inactivated antigen vaccines or subunit vaccines, genetic vaccines result in an antigenic presentation via MHC class I and MHC class II molecules, which mimics the resulting process of natural infection to activate T lymphocytes CD4+, CD8+ and the production of antibodies. The different types of induced immunological response for DNA administration clearly justify their application in the infectious disease field . The question is how to derive this specific answer for the purpose of generating a hyperimmune serum, and less cellular protection. Some attempts to achieve this objective are in progress. Besides the few DNA immunization works referring to snake toxins, some efforts are being made to investigate the generation of hyperimmune serum against bacterial or viral diseases. Fischer and colleagues in 2003 , for example, using a cationic lipid composition together with DNA, generated reasonable titers of antibodies against rabies virus in horses. Herrmann and colleagues in 2006  were capable to produce an antiserum in rabbits against anthrax that was capable of neutralizing the deleterious effects of the Bacillus anthracis toxin. Thus, quoting Harrison et al., 2000, : "The application of DNA-based methodologies to the development of therapeutic antivenom represents the first major conceptual change in antivenom production in over a century and has the potential to provide a more cost-effective, less hazardous and more immunologically specific therapy than those used currently to treat envenoming by snakes"
The transcriptome analysis of M. corallinus provides a large profile of Elapidae toxin cDNAs. Ten classes of possible toxins were found, representing a great diversity of toxins for a venom believed to be almost exclusively neurotoxic. Nevertheless, the possible neurotoxins (3FTx and PLA2) are in fact the majority, totaling 85% of toxin transcripts. The possible post-synaptic components (3FTx) are very diverse in terms of sequences, possibly aiming to achieve different kinds of receptors. The pre-synaptic component (PLA2), in contrast, is more conserved, with the main transcript being represented by 196 of the 1438 ESTs analyzed here. Nevertheless, the high expression of both types of possible neurotoxins is in agreement with the known presence of pre- and post-synaptic activities in the venom of this species.
From a biotechnological point of view, this transcriptome set represents a library of naturally selected templates, especially for molecules acting on nicotinic receptors or ion channels which may be useful for pharmacological purposes. Regarding envenomation treatment, the utilization of genetic immunization based on the survey of transcriptome data carried out here was shown to be feasible for generating immune responses, although more optimization is still needed. Although far from therapeutic application, these findings represent the first steps in the production of an alternative anti-M. corallinus venom serum.
cDNA library preparation
Two cDNA libraries were constructed in 1992, with RNA poli(A+) extracted from 10 specimens of M. corallinus . The cDNAs were divided in two fractions (400 to 600 bp and >600 bp) and the transcripts were linked to λgt11D (Pharmacia, USA) pre-digested phage with EcoRI/NotI restriction enzymes. The phage library was amplified after infection in E. coli Y1088.
To isolate bacteriophage plaques, first a colony of E. coli Y1088 was inoculated in 5 mL of NZCYM culture media containing 0.2% of maltose . After the addiction of 0.01 M MgSO4, the absorbance was accomplished until reaching 2.0 and then kept at 4°C. Dilutions of phage libraries were mixed with 100 μL of this E. coli Y1088 preparation and incubated for 20 minutes at 37°C. NZCYM culture media containing top agarose (0.7% agarose) was then added to the mixture and it was distributed over a bottom agar (1.5% agar) plate prepared with NZCYM culture media. After 12–16 hours at 37°C, the phage plaques were ready to be isolated.
Amplification of phage cDNA insert by PCR and DNA sequencing
Each lytic colony was collected using a Pasteur pipette tip and inserted into a microcentrifuge tube containing 200 μl of SM media . After 1 hour of incubation with agitation at room temperature, 10 μl of this solution was used to amplify the cDNA inserted in the genetic material of each phage with Taq polymerase enzyme (Platinum Taq DNA Polymerase – Invitrogen) and λgt11F and R primers by PCR, in a 96 well plate. Eight microliters of the amplified product was then chemically cleaned by the addition of 5 units of Exonuclease I (GE Healthcare), 1 unit of Shrimp Alkaline Phosphatase (SAP) (GE Healthcare) for 1 hour at 37°C and then inactivated for 15 minutes at 80°C. From the purified PCR product, 300 ng was sequenced with BigDye2 dideoxy-terminators (Applied Biosystems) and λgt11F primer on an ABI 3100 sequencer.
Cluster assembly and identification
After sequencing the DNA, the electropherogram files were analyzed in a semi-automatic way, as described . The Phred program was used http://www.phrap.com to remove bad quality sequences (window length of 75 bases with 75% of standard quality < 25). Then, adapter and phage sequences were removed by the CrossMatch Program. An examination was carried out manually and sequences below 150 bp were discarded. ESTs were then assembled in clusters of contiguous sequences using the CAP3 program , set for 98% or more of base identity in a high-quality region. The cluster sequences were searched in the GeneBank NCBI database with BLASTX and BLASTN algorithms. Non-identified sequences and those with unpredicted function were checked for the presence of signal peptide by SignalP 3.0 program http://www.cbs.dtu.dk/services/SignalP/. A final annotation table was generated in Microsoft Excel format containing all the relevant information about clusters. ESTs sequences were deposited in Genbank dbEST under accession numbers FL589790 to FL591230.
Antigenic candidates selection, cloning and immunization
Five cDNAs were selected after an analysis of the Micrurus corallinus transcriptomic databank and named Atg1 to Atg5 (Antigens 1 to 5), as described in the Results. Specific oligos were drawn for each antigenic candidate. The mature protein region was amplified by PCR and cloned, first, in pGEM-T vector, then digested and subcloned in the eukaryotic expression vector, pSecTag2A (Invitrogen). This vector allows the secretion of the protein to the medium due to the secretion signal from the kappa chain of immunoglobulins. The expressed protein also contains some tags like a poly-histidine (His) and the myc epitopes (c-myc). To obtain the large amount of DNA needed for the immunization process, maxipreps were made with the Qiagen Plasmid Purification kit (Qiagen). They were then quantified and its purity was determined (OD 260/280) by spectrophotometry.
The DNA immunization of the candidates was made in five groups of BALB/c mice, females of homogenous weight (18–22 g) and age (6–8 weeks old). The experiments in animals was approved by the Ethics Committee of Butantan Institute, under protocol 258/06, according to international guidelines. Each group of five animals was immunized with an antigenic candidate. The animals were anesthetized i.p. with approximately 200 μL of anesthetic (ketamine 500 μL/10 mL and xilazine 1000 μL/10 mL). The injection of the DNA vaccine was intramuscular (IM). Each animal received 100 μg of DNA per dose and one dose every 15 days for a total of three doses (300 μg/animal). The same procedure was used to immunize a group of five animals with a pool of three immunogens (33 μg of Atg1, 33 μg of Atg2 and 33 μg of Atg5 per dose).
Immune response analysis
Mice were bled individually from the retro-orbital plexus to detect antibodies in the blood through ELISA (Enzyme-linked Immunosorbent Assay). An initial bleeding was done before the immunizations, in order to obtain the pre-immune serum and the final bleeding was made after the third dose. Recombinant proteins from 3 of the 5 candidates (Atg1, Atg2, Atg5) were produced in E. coli, as His tagged fusion proteins, and purified by Ni2+ affinity chromatography, as described previously [10, 11]. These proteins (1 μg) were used for the ELISA coating in NaHCO3 and Na2CO3 0,05 M pH 9,6 buffer. After blocking it with 10% milk in PBS, dilutions of the sera were prepared in ranges from 1:20 to 1:10240 and incubated during 1.5 hour with a 1% BSA + 0.1% PBS-T solution, in 37°C. The plates were then incubated with a secondary antibody, IgG anti-mouse conjugated with peroxidase for 1 hour in 37°C and revealed with 100 μl of 0.4 mg/ml OPD (ortofenilenodiamine), 0.05% H2O2 in citrate-phosphate buffer (0.1 M citric acid + 0.2 M sodium phosphate). The reaction was blocked with the addition of 50 μl of 8N H2SO4 after 15 minutes. The experiments were realized in triplicates.
This work was supported by grants from the Brazilian agency Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP). We thank Giuseppe Puorto and Otávio A. V. Marques for obtaining the specimens in 1992, and Marcelo Bento Soares for helping in the generation of the cDNA library used in this work. We are grateful to Henrique R. Ramos for his help in the immunization experiments. Dr. A. Leyva provided English editing of the manuscript.
- Dowling HG, Duellman WE: Systematic Herpetology: A synopsis of families and higher categories. 1978, New York: Hiss PublicationsGoogle Scholar
- Vital Brazil O: [Neurotoxic snake venoms]. AMB Rev Assoc Med Bras. 1980, 26 (6): 212-218.Google Scholar
- Ho PL, Soares MB, Yamane T, Raw I: Reverse Biology applied to Micrurus corallinus, a South American Coral Snake. Toxicol Toxin Reviews. 1995, 14 (3): 10-Google Scholar
- Junqueira-de-Azevedo ILM, Ho PL: A survey of gene expression and diversity in the venom glands of the pitviper snake Bothrops insularis through the generation of expressed sequence tags (ESTs). Gene. 2002, 299 (1–2): 279-291. 10.1016/S0378-1119(02)01080-6.View ArticleGoogle Scholar
- Junqueira-de-Azevedo IL, Ching AT, Carvalho E, Faria F, Nishiyama MY, Ho PL, Diniz MR: Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution. Genetics. 2006, 173 (2): 877-889. 10.1534/genetics.106.056515.PubMed CentralView ArticlePubMedGoogle Scholar
- Qinghua L, Xiaowei Z, Wei Y, Chenji L, Yijun H, Pengxin Q, Xingwen S, Songnian H, Guangmei Y: A catalog for transcripts in the venom gland of the Agkistrodon acutus: identification of the toxins potentially involved in coagulopathy. Biochem Biophys Res Commun. 2006, 341 (2): 522-531. 10.1016/j.bbrc.2006.01.006.View ArticlePubMedGoogle Scholar
- Wagstaff SC, Harrison RA: Venom gland EST analysis of the saw-scaled viper, Echis ocellatus, reveals novel alpha9beta1 integrin-binding motifs in venom metalloproteinases and a new group of putative toxins, renin-like aspartic proteases. Gene. 2006, 377: 21-32. 10.1016/j.gene.2006.03.008.View ArticlePubMedGoogle Scholar
- Jia Y, Cantu BA, Sanchez EE, Perez JC: Complementary DNA sequencing and identification of mRNAs from the venomous gland of Agkistrodon piscivorus leucostoma. Toxicon. 2008, 51 (8): 1457-1466. 10.1016/j.toxicon.2008.03.028.PubMed CentralView ArticlePubMedGoogle Scholar
- Ching AT, Rocha MM, Paes Leme AF, Pimenta DC, de Fatima DFM, Serrano SM, Ho PL, Junqueira-de-Azevedo IL: Some aspects of the venom proteome of the Colubridae snake Philodryas olfersii revealed from a Duvernoy's (venom) gland transcriptome. FEBS Lett. 2006, 580 (18): 4417-4422. 10.1016/j.febslet.2006.07.010.View ArticlePubMedGoogle Scholar
- Prieto da Silva AR, Yamagushi IK, Morais JF, Higashi HG, Raw I, Ho PL, Oliveira JS: Cross reactivity of different specific Micrurus antivenom sera with homologous and heterologous snake venoms. Toxicon. 2001, 39 (7): 949-953. 10.1016/S0041-0101(00)00233-6.View ArticlePubMedGoogle Scholar
- de Oliveira UC, Assui A, da Silva AR, de Oliveira JS, Ho PL: Cloning and characterization of a basic phospholipase A2 homologue from Micrurus corallinus (coral snake) venom gland. Toxicon. 2003, 42 (3): 249-255. 10.1016/S0041-0101(03)00139-9.View ArticlePubMedGoogle Scholar
- Tang DC, DeVit M, Johnston SA: Genetic immunization is a simple method for eliciting an immune response. Nature. 1992, 356 (6365): 152-154. 10.1038/356152a0.View ArticlePubMedGoogle Scholar
- Harrison RA: Development of venom toxin-specific antibodies by DNA immunisation: rationale and strategies to improve therapy of viper envenoming. Vaccine. 2004, 22 (13–14): 1648-1655. 10.1016/j.vaccine.2003.09.046.View ArticlePubMedGoogle Scholar
- Wagstaff SC, Laing GD, Theakston RD, Papaspyridis C, Harrison RA: Bioinformatics and multiepitope DNA immunization to design rational snake antivenom. PLoS Med. 2006, 3 (6): e184-10.1371/journal.pmed.0030184.PubMed CentralView ArticlePubMedGoogle Scholar
- Fry BG, Vidal N, Norman JA, Vonk FJ, Scheib H, Ramjan SF, Kuruppu S, Fung K, Hedges SB, Richardson MK, et al: Early evolution of the venom system in lizards and snakes. Nature. 2006, 439 (7076): 584-588. 10.1038/nature04328.View ArticlePubMedGoogle Scholar
- Nirthanan S, Gwee MC: Three-finger alpha-neurotoxins and the nicotinic acetylcholine receptor, forty years on. J Pharmacol Sci. 2004, 94 (1): 1-17. 10.1254/jphs.94.1.View ArticlePubMedGoogle Scholar
- Silveira de Oliveira J, Rossan de Brandao Prieto da Silva A, Soares MB, Stephano MA, de Oliveira Dias W, Raw I, Ho PL: Cloning and characterization of an alpha-neurotoxin-type protein specific for the coral snake Micrurus corallinus. Biochem Biophys Res Commun. 2000, 267 (3): 887-891. 10.1006/bbrc.1999.2033.View ArticlePubMedGoogle Scholar
- Dufton MJ, Hider RC: Classification of phospholipases A2 according to sequence. Evolutionary and pharmacological implications. Eur J Biochem. 1983, 137 (3): 545-551. 10.1111/j.1432-1033.1983.tb07860.x.View ArticlePubMedGoogle Scholar
- Harris JB: Phospholipases in snake venoms and their effects on nerve and muscle. Snake Toxins, International Encyclopedia of Pharmacology and Therapeutics. Edited by: Harvey AL. 1991, New York: Pergamon Press, 91-129.Google Scholar
- Mebs D, Ownby CL: Myotoxic components of snake venoms: their biochemical and biological activities. Pharmacol Ther. 1990, 48 (2): 223-236. 10.1016/0163-7258(90)90081-C.View ArticlePubMedGoogle Scholar
- Rossetto O, Morbiato L, Caccin P, Rigoni M, Montecucco C: Presynaptic enzymatic neurotoxins. J Neurochem. 2006, 97 (6): 1534-1545. 10.1111/j.1471-4159.2006.03965.x.View ArticlePubMedGoogle Scholar
- Maraganore JM, Heinrikson RL: The lysine-49 phospholipase A2 from the venom of Agkistrodon piscivorus piscivorus. Relation of structure and function to other phospholipases A2. J Biol Chem. 1986, 261 (11): 4797-4804.PubMedGoogle Scholar
- Fujimi TJ, Kariya Y, Tsuchiya T, Tamiya T: Nucleotide sequence of phospholipase A(2) gene expressed in snake pancreas reveals the molecular evolution of toxic phospholipase A(2) genes. Gene. 2002, 292 (1–2): 225-231. 10.1016/S0378-1119(02)00682-0.View ArticlePubMedGoogle Scholar
- Doley R, Tram NN, Reza MA, Kini RM: Unusual accelerated rate of deletions and insertions in toxin genes in the venom glands of the pygmy copperhead (Austrelaps labialis) from Kangaroo island. BMC Evol Biol. 2008, 8: 70-10.1186/1471-2148-8-70.PubMed CentralView ArticlePubMedGoogle Scholar
- Weis WI, Taylor ME, Drickamer K: The C-type lectin superfamily in the immune system. Immunol Rev. 1998, 163: 19-34. 10.1111/j.1600-065X.1998.tb01185.x.View ArticlePubMedGoogle Scholar
- Ogawa T, Chijiwa T, Oda-Ueda N, Ohno M: Molecular diversity and accelerated evolution of C-type lectin-like proteins from snake venom. Toxicon. 2005, 45 (1): 1-14. 10.1016/j.toxicon.2004.07.028.View ArticlePubMedGoogle Scholar
- Ho PL, Soares MB, Maack T, Gimenez I, Puorto G, Furtado MF, Raw I: Cloning of an unusual natriuretic peptide from the South American coral snake Micrurus corallinus. Eur J Biochem. 1997, 250 (1): 144-149. 10.1111/j.1432-1033.1997.00144.x.View ArticlePubMedGoogle Scholar
- Murayama N, Hayashi MA, Ohi H, Ferreira LA, Hermann VV, Saito H, Fujita Y, Higuchi S, Fernandes BL, Yamane T, et al: Cloning and sequence analysis of a Bothrops jararaca cDNA encoding a precursor of seven bradykinin-potentiating peptides and a C-type natriuretic peptide. Proc Natl Acad Sci USA. 1997, 94 (4): 1189-1193. 10.1073/pnas.94.4.1189.PubMed CentralView ArticlePubMedGoogle Scholar
- Graham RL, Graham C, McClean S, Chen T, O'Rourke M, Hirst D, Theakston D, Shaw C: Identification and functional analysis of a novel bradykinin inhibitory peptide in the venoms of New World Crotalinae pit vipers. Biochem Biophys Res Commun. 2005, 338 (3): 1587-1592. 10.1016/j.bbrc.2005.10.130.View ArticlePubMedGoogle Scholar
- Wagstaff SC, Favreau P, Cheneval O, Laing GD, Wilkinson MC, Miller RL, Stocklin R, Harrison RA: Molecular characterisation of endogenous snake venom metalloproteinase inhibitors. Biochem Biophys Res Commun. 2008, 365 (4): 650-656. 10.1016/j.bbrc.2007.11.027.View ArticlePubMedGoogle Scholar
- Nathan I, Dvilansky A, Yirmiyahu T, Aharon M, Livne A: Impairment of platelet aggregation by Echis colorata venom mediated by L-amino acid oxidase or H2O2. Thromb Haemost. 1982, 48 (3): 277-282.PubMedGoogle Scholar
- Raibekas AA, Massey V: Primary structure of the snake venom L-amino acid oxidase shows high homology with the mouse B cell interleukin 4-induced Fig 1 protein. Biochem Biophys Res Commun. 1998, 248 (3): 476-478. 10.1006/bbrc.1998.9024.View ArticlePubMedGoogle Scholar
- Zupunski V, Kordis D, Gubensek F: Adaptive evolution in the snake venom Kunitz/BPTI protein family. FEBS Lett. 2003, 547 (1–3): 131-136. 10.1016/S0014-5793(03)00693-8.View ArticlePubMedGoogle Scholar
- Harvey AL, Karlsson E: Protease inhibitor homologues from mamba venoms: facilitation of acetylcholine release and interactions with prejunctional blocking toxins. Br J Pharmacol. 1982, 77 (1): 153-161.PubMed CentralView ArticlePubMedGoogle Scholar
- Schweitz H, Heurteaux C, Bois P, Moinier D, Romey G, Lazdunski M: Calcicludine, a venom peptide of the Kunitz-type protease inhibitor family, is a potent blocker of high-threshold Ca2+ channels with a high affinity for L-type channels in cerebellar granule neurons. Proc Natl Acad Sci USA. 1994, 91 (3): 878-882. 10.1073/pnas.91.3.878.PubMed CentralView ArticlePubMedGoogle Scholar
- Fox JW, Serrano SM: Insights into and speculations about snake venom metalloproteinase (SVMP) synthesis, folding and disulfide bond formation and their contribution to venom complexity. FEBS J. 2008, 275 (12): 3016-3030. 10.1111/j.1742-4658.2008.06466.x.View ArticlePubMedGoogle Scholar
- Kini RM: Serine proteases affecting blood coagulation and fibrinolysis from snake venoms. Pathophysiol Haemost Thromb. 2005, 34 (4–5): 200-204. 10.1159/000092424.View ArticlePubMedGoogle Scholar
- Kashima S, Soares AM, Roberto PG, Pereira JO, Astolfi-Filho S, Cintra AO, Fontes MR, Giglio JR, de Castro Franca S: cDNA sequence and molecular modeling of a nerve growth factor from Bothrops jararacussu venomous gland. Biochimie. 2002, 84 (7): 675-680. 10.1016/S0300-9084(02)01429-3.View ArticlePubMedGoogle Scholar
- Kostiza T, Meier J: Nerve growth factors from snake venoms: chemical properties, mode of action and biological significance. Toxicon. 1996, 34 (7): 787-806. 10.1016/0041-0101(96)00023-2.View ArticlePubMedGoogle Scholar
- Torres AM, Wong HY, Desai M, Moochhala S, Kuchel PW, Kini RM: Identification of a novel family of proteins in snake venoms. Purification and structural characterization of nawaprin from Naja nigricollis snake venom. J Biol Chem. 2003, 278 (41): 40097-40104. 10.1074/jbc.M305322200.View ArticlePubMedGoogle Scholar
- Olamendi-Portugal T, Batista CV, Restano-Cassulini R, Pando V, Villa-Hernandez O, Zavaleta-Martinez-Vargas A, Salas-Arruz MC, Rodriguez de la Vega RC, Becerril B, Possani LD: Proteomic analysis of the venom from the fish eating coral snake Micrurus surinamensis: novel toxins, their function and phylogeny. Proteomics. 2008, 8 (9): 1919-1932. 10.1002/pmic.200700668.View ArticlePubMedGoogle Scholar
- Harrison RA, Moura-Da-Silva AM, Laing GD, Wu Y, Richards A, Broadhead A, Bianco AE, Theakston RD: Antibody from mice immunized with DNA encoding the carboxyl-disintegrin and cysteine-rich domain (JD9) of the haemorrhagic metalloprotease, Jararhagin, inhibits the main lethal component of viper venom. Clin Exp Immunol. 2000, 121 (2): 358-363. 10.1046/j.1365-2249.2000.01287.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Fischer L, Minke J, Dufay N, Baudu P, Audonnet JC: Rabies DNA vaccine in the horse: strategies to improve serological responses. Vaccine. 2003, 21 (31): 4593-4596. 10.1016/S0264-410X(03)00504-8.View ArticlePubMedGoogle Scholar
- Herrmann JE, Wang S, Zhang C, Panchal RG, Bavari S, Lyons CR, Lovchik JA, Golding B, Shiloach J, Lu S: Passive immunotherapy of Bacillus anthracis pulmonary infection in mice with antisera produced by DNA immunization. Vaccine. 2006, 24 (31–32): 5872-5880. 10.1016/j.vaccine.2006.04.065.View ArticlePubMedGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: A laboratory manual. New York. 1989Google Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.