Integrated “omics” profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus
© Durban et al.; licensee BioMed Central Ltd. 2013
Received: 12 November 2012
Accepted: 14 March 2013
Published: 10 April 2013
Understanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus using a combined proteomics and transcriptomics approach.
Proteomic analysis showed that the ontogenetic shift in the venom composition of C. s. simus is essentially characterized by a gradual reduction in the expression of serine proteinases and PLA2 molecules, particularly crotoxin, a β-neurotoxic heterodimeric PLA2, concominantly with an increment of PI and PIII metalloproteinases at age 9–18 months. Comparison of the transcriptional activity of the venom glands of neonate and adult C. s. simus specimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulhood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus.
Existing snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112–125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snake's venom gland, are further challenges for future research.
The presence of a venom system is a shared derived character of the advanced snakes [1–3]. Venom represents an adaptive trophic trait in snake evolution [4–7] that has played a central role in the origin of the advanced snakes during the Cenozoic era [8, 9]. The diversity of modern snakes appeared during the Paleocene period of the Cenozoic Era, approximately 54–64 Mya, following the split of the Pareatidae from the remaining Caenophidians . In Occidental culture, Francesco Redi (1626–1697), court physician to Ferdinando II de’ Medici, Grand Duke of Tuscany and his successor, Cosimo III, is credited for discovering how vipers produce venom and inject it into their prey . One century latter, Abbé Gasparo Ferdinando Felice Fontana (1730–1805), first director of the Museum of Physics and Natural History in Florence, performed experiments on the venom of the European viper. His classic text published in 1781  is regarded the Opera Prima of modern toxinology .
Research on venoms has been continuously advanced by technological developments . Particularly for the past decade, and fueled by the application of "omic" technologies [14–18], the field of molecular toxinology has experienced a sustained exponential growth. In-depth venom proteomics and transcriptomic analyses have recently become available for a number of venomous snake lineages, revealing a vast unexplored source of evolutionary innovation [2, 19, 20]. Snake venoms are well documented as having different venom compositions and toxicity based on taxonomic or geographical locations . Inter- and intraspecific individual, gender-specific, regional, and seasonal variations in venom toxin composition may reflect local adaptations conferring fitness advantages to individuals specializing on different prey, or phylogenetic carry on. Understanding the processes that drive the evolution of snake venom variability is a topic of intense research interest in molecular and evolutionary toxinology. Recent studies of the molecular basis of adaptations have sought to understand the relative importance of gene regulation effects as determinants of venom phenotype [6, 22–24]. Most significantly, a number of snakes show age-related (ontogenetic) changes in venom composition [25–32]. Surprisingly, despite previous suggestions that ontogenetic changes in venom are prey-related , juvenile Dusky Pigmy rattlesnakes, Sistrurus miliarius barbouri, raised on different diets showed similar albeit highly-variable venom compositions by the end of a 26-months study, suggesting little effect of diet on the overall make-up of venom in snakes this age or younger . Over the same period shifts in venom composition occurred in females raised on different prey in all diet treatments with the magnitude of those changes strongly related to diet. This work provided evidence that venom composition is somewhat plastic in both juvenile and adult S. m. barbouri and that, at least in adults, prey consumed may influence the relative abundance of possibly functionally-distinct classes of venom proteins . However, the molecular mechanisms that underly age-related changes in venom remain elusive.
In this work we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus. Biogeographical data indicate a basal cladogenesis in the Central American C. simus clade, dating back to the late Miocene/early Pliocene (6.4-6.7 Mya) . Neonate and juvenile C. s. simus venoms contain large relative amounts of crotoxin, a heterodimeric PLA2 molecule exhibiting presynaptic β-neurotoxicity, along with low concentration of hemorrhagic snake venom metalloproteinases (SVMPs) . By contrast, adult C. s. simus venom presents a high content of SVMP and is largely devoid of neurotoxic activity [29, 35, 36]. Juvenile and adult C. s. simus venom phenotypes broadly correspond, respectively, to type II (low metalloproteinase activity and high toxicity, LD50 <1 μg/g mouse body weight) and type I (high levels of SVMPs and low toxicity, LD50 >1 μg/g mouse body weight) venoms defined by Mackessy [37, 38]. Here we investigate the transition from type II to type I venom phenotype in C. s. simus through an integrated "omics" approach involving proteomic survey of time-course venom composition variation, from birth through adulthood, and Next Generation sequencing of the venom gland transcriptomes and microRNAs of neonate and adult specimens. With this experimental design we intended to determine whether venom gene expression is transcriptionally controlled by developmental stage-dependent factors, or whether the ontogenetic changes in venom composition involve regulation at the post-transcriptional level.
Results and discussion
Time-resolved proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus
A previous proteomic survey of the venom of the Central American rattlesnake, C. s. simus, laid the foundation for understanding the changes in toxin composition and overall pharmacological features between adult (predominantly hemorrhagic) and newborn (mainly neurotoxic) snakes . This study revealed prominent stage-dependent protein expression changes between the pooled venoms from newborn and adult specimens from Costa Rica, characterized by a shift from a type II to a type I venom phenotype . Such a conspicuous venom composition change has been hypothesized to be related to variations in the size of prey consumed by snakes of different ages and the variable requirements for immobilizing and digesting them [37, 38]. It has also been suggested that the significance of venom adaptation to specific diets represents a trade-off between the metabolic cost of venom production and increasing foraging efficiency [7, 39]. In the present work we have sought to establish the molecular mechanism responsible for the reported ontogenetic shift in venom protein composition.
Analysis of the venom gland transcriptomes of neonate and adult C. s. simus. Retroelements
454 sequencing statistics
Mean length read (nt)
Max. length read (nt)
Assembled reads (%)
Number of contigs
Mean length contig (nt)
Max. length contig (nt)
Contig sequences were inspected for repetitive elements using Repeat Masker. 97,204 bases (1.93% of adult C. s. simus venom gland transcriptome) and 65,742 nucleotides (1.94% of the newborn venom gland transcriptome) were masked with N characters, a large part of them comprising Short and Long INterspersed repetitive Elements (SINEs and LINEs) retroelements (Additional file 1: Tables S1 and S2). Transposable elements (TE) that propagate within the host genome via RNA intermediates occupy a large fraction of eukaryotic genomes. Their mobility and amplification represent a major source of genomic variation [50, 51]. Retrotransposable elements have been reported in the transcriptomes of Bothrops insularis (4.1% of ESTs) , Lachesis muta (0.3%) , and Philodryas olfersii (4.1%) , in PLA2 genes from the venom gland of Vipera ammodytes[55, 56] and Protobothrops flavoviridis[57, 58], and in an E. ocellatus PIII-SVMP gene . Although their functional relevance in the venom gland remains unknown, transposable elements appear to be overrepresented in UTRs of mRNAs of rapidly evolving genes , suggesting that they have played a role in the diversification and expansion of these gene families [61, 62].
Sauria SINE have been characterized in all major lineages of squamate reptiles [62, 63], and phylogenetic analysis of E. ocellatus Sauria SINEs  indicated that their origin correlates with the time frame of the evolution of the snake venom system. Sauria SINEs may have arisen by a fusion of a tRNA-related sequence with the 3' end of a LINE  more than 200 million years ago and are uniquely widespread in lepidosaurian genomes . SINEs have no protein coding capacity, and their retrotransposition depends on a "target-primed reverse transcription" by autonomous partner LINEs, that encode an endonuclease for cleaving the genomic integration site and a reverse transcriptase to copy the RNA to DNA [65, 66]. Since Sauria SINEs share an identical 3' sequence with Bov-B LINEs, it has been proposed  that they utilize in trans the enzymatic machinery of Bov-B LINEs for their mobility and dispersal throughout the genome. In squamates and turtles, CR1 and L2 LINES are also partners of diverse SINEs [63, 67].
An overview of the landscape of retrotransposable elements reported in the sauropsida taxon, a sister group of mammals that includes all extant reptiles and birds, has recently emerged from analysis of the draft genomes of the red jungle fowl, Gallus gallus, and the green anole, Anolis carolinensis. Whereas a single type of LINE, CR1, comprises over 80% of all interspersed repeats in the chicken genome (200,000 copies; 9% of the chicken genome), approximately 30% of the A. carolinensis genome is composed of a wide diversity of LINE and SINE mobile element families. Since SINEs are among the largest multigene families in reptilian genomes, they may act as sites for homologous recombination events between dispersed SINEs, resulting in a variety of genetic rearrangements, including duplication, deletion and translocation, that likely represent mechanisms that generates genetic diversity in reptilian genomes . The prevalence of transposable elements in untranslated regions of mRNAs of recently expanded gene classes suggested that TEs could affect gene expression through the donation of transcriptional regulatory signals . The indistinguishable distribution of TEs in the venom gland transcriptomes of neonate and adult C. s. simus (Additional file 1: Tables S1 and S2) argues against this type of transcriptional regulation to explain the ontogenetic shift in venom composition.
Comparison of toxin-coding transcript distribution in newborn and adult C. s. simustranscriptomes provides clues for streamlining their divergent venom phenotypes
RPKM (Reads per Kilobase per Million mapped reads)-normalized contigs and singletons from the venom gland transcriptomes of newborn and adult C. s. simus aligned to a reference snake venom toxin ORF
Estimation of minimum number of toxin genes transcribed into the venom gland transcriptomes of newborn (N) and adult (A) C. s. simus
N + A
Distinct venom gland miRNA sets in newborn and adult snakes: modulators of the Central American rattlesnake's ontogenetic venom composition variation?
A. Summary of candidate miRNA sequencing statistics. B. Clustering and miRBase annotation statistics
Total raw reads
Prinseq quality filtered reads
FastX too-short reads
FastX adapter-only reads
Mean length read (nt)
BLAST miRBase hits
% Anolis hits
Concluding remarks and perspectives
The most important concept that emerges from our results is the possibility that miRNAs play a role in both the ontogenetic shift observed in certain venoms and in the plasticity of venoms from adult snakes underlying adaptations to changing environments. An important caveat of our current understanding of miRNA target recognition and post-transcriptional venom gland transcriptome regulation is the absence of mapeable snake genomes. Nonetheless, almost half of the most abundant miRNAs sequenced from the venom gland of both, newborn and adult, C. s. simus were orthologous to sequences in A. carolinensis. More significant is the fact that a C. s. simus newborn venom gland-exclusive miRNA and a venom gland miRNA uniquely found in adult snakes target, respectively, mRNAs encoding a PIII-SVMP and the B-subunit of crotoxin/Mojave toxin. Relevant to this point, Mojave toxin is a neurotoxic heterodimeric PLA2, whose translation into the venom proteome is ontogenetically regulated  (Figure 1). The noncovalent association of two dissimilar subunits, the small acidic, nonenzymatic, and nontoxic A-subunit with the basic and weakly toxic PLA2 B-subunit increases the lethal potency of the uncomplexed crotoxin B-subunit by enabling the toxin to reach its specific site of action at the neuromuscular junction [40, 41, 89–91]. Hence, miRNAs targeting a crotoxin subunit messenger would serve the goal of eliminating the generation of the neurotoxic heterodimer. Age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulthood can potentially explain what is observed in the age-dependent proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus illustrated in Figure 1.
Large-scale proteomic analysis, performed for one miRNA (miR-223) in only one cell type (murine neutrophils), revealed that although some proteins were repressed by 50-80%, miR-223 typically had more modest effects , suggesting that perhaps other miRNAs in their endogenous context have targets for which protein output is more dramatically repressed. Clearly, a challenge ahead in molecular toxinology is to design laboratory experiments to uncover the impact, molecular details, and evolution of the regulatory miRNA-target interactions that shape venom phenotype during snake development. In this regard, it has been proposed that a single origin of venom in Squamata, the order of reptiles including lizards and snakes, dates back roughly 200 Mya to the Late Triassic/Early Jurassic period [1, 2, 93]. Existing snake venom toxins are the result of recruitment events by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112–125 Mya of snake evolution. Given the central role that diet has played in the adaptive radiation of snakes , venom thus represents a key adaptation that has played an important role in the diversification of snakes. Our findings here reported support the view that understanding the basis for the phenotypic diversity of snake venoms requires a deep understanding of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results, though restricted to individual specimens from a single species, suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation and impact of these miRNAs on patterns of expression in the snake's venom gland, are further challenges for future research.
Proteomics assessment of the ontogenetic shift of venom composition in C. simus
Venoms from a large number (>25) of adult Crotalus s. simus specimens and from 11 captive-born siblings (from an age of 8-weeks to 21-months) kept at the serpentarium of the Instituto Clodomiro Picado, University of Costa Rica in San José were collected by snake biting on a parafilm-wrapped jar (adults and juveniles) or by aspiration from the fangs of neonates using an Eppendorf pipette. Crude venoms were centrifuged at low speed to remove cells and debris, lyophilized, weighed on a microbalance, and stored at −20°C until used. Venom proteins were separated by reverse-phase HPLC using a Teknokroma Europa C18 (0.4 cm × 25 cm, 5 mm particle size, 300 Å pore size) column and an Agilent LC 1100 High Pressure Gradient System equipped with DAD detector and micro-Auto-sampler. The flow-rate was set to 1 ml/min and the column was developed with a linear gradient of 0.1% TFA in water (solution A) and acetonitrile (solution B), isocratically (5% B) for 10 min, followed by 5–25 % B for 20 min, 25-45% B for 120 min, and 45-70% for 40 min. Isolated proteins were characterized as described .
Snake venom gland cDNA synthesis and 454 sequencing
Venom glands of an 8-week-old and an adult C. s. simus specimens were removed 3 days after venom milking, when transcription is maximal , from anesthetized snakes using fine forceps and immediately placed in RNAlater™ solution (Qiagen). About 50 mg of tissue were disrupted and homogenized by a rotor-stator homogenizer, and total RNA was isolated using RNeasy Mini kit (Qiagen), quantified in a NanoDrop ND_spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) through the A260/A280 ratio, and quality-checked on an agarose gel discerning the 28S and 18S bands of ribosomal RNA. First strand cDNA was synthesized using RevertAid™ H Minus First Strand cDNA Synthesis Kit (Fermentas), which selectively transcribes full-length polyadenylated mRNA. The manufacturer's recommendations were followed except where specified. Approximately 5 μg of total RNA was used as starting material. In order to avoid polymerase slippage, a modified 3' 54-mer adaptor (5' GAGCTAGTTCTGGAG(T)16VN), which includes a type IIs enzyme (GsuI) site (underlined), was used for first-strand synthesis. This modified oligonucleotide effectively converts the long run of adenosine residues at the polyA tail into a sequence that causes fewer problems for dideoxy sequencing chemistry, and thus the resulting cDNA libraries were enriched in 3'-end-transcripts. To avoid internal cuts, the cDNA was hemimethylated by adding 5-methyl-dCTP to the dNTPs mix. The first strand cDNA was used as template for second strand synthesis by E. coli DNA Polymerase I and RNase H. Double strand (ds) cDNA was precipitated with ethanol and the pellet was resuspended in 70 μL of nuclease-free water and subjected to enzymatic digestion with GsuI for 4 hours at 30°C. The enzyme was then inactivated at 65°C for 20 minutes and the digested cDNA was stored at −20°C. For 454 pyrosequencing, the GS FLX General DNA Library Preparation Method Manual workflow (Roche Diagnostics) was followed. To this end, 3 μg of final non-normalized cDNA library were sheared by nebulization into small fragments. The fragment ends were polished and short A/B adaptors were ligated onto both ends, providing priming regions to support both emulsion amplification and the pyrosequencing process. A biotin tag on the B adaptor allowed immobilization of the dscDNA library fragments onto streptavidin-conjugated magnetic beads and the subsequent isolation of the library of single strand cDNA sequencing templates. Each of the eight cDNA libraries was tagged with a unique 10-base sequence (MID, Multiplex IDentifier) that is recognized by the sequencing analysis software, allowing for automated sorting of MID-containing reads. Barcoded libraries were simultaneously sequenced in a Genome sequencing FLX Titanium System (Roche Applied Science) at Life Sequencing S.L. (Parc Científic Universitat de Valencia, Paterna, Valencia, Spain; http://www.lifesequencing.com) using the method developed by Margulies et al. .
Bioinformatic analysis of the 454 reads
454 data were analyzed using the workflow developed in  to identify sequences of toxin molecules by similarity search against nucleotide databases, which includes available NGS algorithms and in-house scripts. An initial quality test step of both the raw reads provided by the Genome Sequencer FLX System prior to the assembly process and the longer contig sequences obtained after running the 454 Newbler assembler program (version 2.6) (Titanium chemistry) was run using the Prinseq program (standalone version 0.17.4) . Interspersed repeats and low complexity DNA sequences were masked from the transcript reads using RepeatMasker (version 3.2.9) . RepeatMasker is available from the Institute for Systems Biology (http://www.systemsbiology.org) addressing http://repeatmasker.org. The program screens DNA sequences for interspersed repeats and low complexity DNA sequences included in the Repbase database (http://www.girinst.org). Repbase, a comprehensive database of repetitive element consensus sequences (update of 20 September 2011), operates as a service of the Genetic Information Research Institute (http://www.girinst.org). Data and computational resources for the Pre-Masked Genomes page is provided courtesy of the UCSC Genome Bioinformatics group (http://genome.ucsc.edu). Masked contig sequences were translated into the 6 possible reading frames and blasted against the non-redundant NCBI database (http://blast.ncbi.nlm.nih.gov, release of February 2012) and the UniProtKB/Swiss-Prot Toxin Annotation Program database (http://us.expasy.org/sprot/tox-prot), using BlastX and BlastN  algorithms (version 2.2.24), specifying a cut-off value of e-03 and BLOSUM62 as scoring matrix. Snake venom gland-specific transcripts were selected from best BLAST-hit descriptions identifying GenBank entries belonging to the taxonomic suborder Serpentes. This taxonomic group is represented by 44,141 nucleotide records comprising entries from 2,396 different species. A second filtering round was carried out using a list of keywords (including the acronyms of all known toxin protein families described so far to distinguish putative snake venom toxins from non-toxin (ribosomal, mitochondrial, nuclear, etc.) ordinary proteins . The phylogenetically nearest top-hit full-length sequence was designated as the reference sequence onto which all toxin family-specific reads were aligned to create a multiple alignment using COBALT . The multiple alignment was then parsed to create an assembled (consensus) toxin sequence in which each amino acid position is supported by at least four reads. We considered two contigs as different if their nucleotide sequences depart in more number of positions than expected from a sequencing error rate of 1.5%, and the same mutated residues were found in at least two other reads. Positions where two or more amino-acids fulfilled this criterium were annotated as variable residues suggesting the occurrence of different alleles (isoforms) of the protein.
Possible differential expression of venom proteins between adult and the newborn individuals was assessed with the non-parametric NOISeq-sim algorithm  using the following parameters recommended for counts without replicates in the NOISeq-sim manual were used: k (counts equal to zero) = 0.5; nss (number of samples to be simulated) ≥ 5; pnr (percentage of the total sequencing depth) = 0.2; v (variability in the total sequencing depth of the simulated sample) = 0.02. Reads were normalized by the length of the phylogenetically nearest sequence and a threshold of 0.9 was allowed.
Comparison of newborn and adult C. s. simusvenom transcriptomes
To assess the degree of similarity between transcripts synthesized by newborn an adult venom glands, a Perl script was written that aligned singletons and Newbler-assembled contigs onto the open reading frame (ORF) of the phylogenetically nearest protein sequence used as reference. The aligned nucleotide sequences were re-assembled with MIRA (http://www.chevreux.org/projects_mira.html) to infer the minimum number of different assemblies. These resulting sequences from newborn and adult individuals were compiled into a single FASTA file, translated into protein sequences, and manually inspected to discard possible mispaired BLAST annotations of local regions due to incorrect frame translations. Protein sequences were then clustered with CD-HIT (standalone version 4.5.7 built on 27th February 2012) [102, 103] to identify protein sequences shared between newborn and adult transcriptomes.
Snake venom gland microRNA profiling
Total RNA from neonate and adult venom glands was used to isolate small RNA libraries using a RNeasy Mini Kit following the manufacture's (Qiagen) instructions. Samples were quantified using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) through the 260/280 absorbance relation. 5 microliter of neonate (258.6 ng/microliter) and adult (490.7 ng/microliter) size (<200 nt)-enriched microRNA (miRNA) library were sequenced on an Ion Torrent Personal Genome Machine (PGM™) Sequencing platform at Life Sequencing S.L.
Bioinformatic processing of the Ion-Torrent miRNA reads
Using Prinseq , the Ion-Torrent miRNA reads were filtered out according to nucleotide length (min_length 15 nt; max_length 40 nt) and the presence of either the 30-mer IonA or the 23-mer P1 adapter sequences. Adapter sequences were removed using the fasxt_clipper tool from the Fastx-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit). The resulting reads were clustered using CD-HIT-454  setting an identity threshold of 0.98 and the accurate mode option g1. HIT-454 was designed to identify artificial duplicates from raw 454 sequencing reads, including exact duplicates and near identical duplicates. Script cdhit-cluster-consensus (v.13) of the CD-HIT suite of programs was run to derive consensus sequences for each of the miRNA clusters. This program is currently maintained by Dr. Li's group (http://weizhong-lab.ucsd.edu). The relative abundances of miRNAs with frequency count ≥ 100 were normalized by scaling the number of reads clustered by CD-HIT-454 to the total number of processed reads (220,548 for newborn + 389,064 for adult). To identify unique and shared sequences between the newborn and adult datasets, these newborn and adult miRNA clusters were compared between themselves using CD-HIT. Differential expression between miRNAs in these two datasets was assessed with NOISeq-sim  allowing a threshold of 0.9. miRNAs comprising ≥ 100 reads were subjected to BLAST analysis against miRBase (release 18, November 2011, which included 21643 mature miRNA products from 168 species) (http://www.mirbase.org) .
The search for miRNA targets
The precise rules and energetics for pairing between a miRNA and its mRNA target sites are not completely understood [104, 105]. A variety of computational algorithms aimed at predicting miRNA target genes incorporate rules for miRNA-mRNA interactions such as base pairing complementarity and favourable miRNA-mRNA duplex thermodynamics. Current prediction methods are diverse, both in approach and performance . We have employed the freely available target prediction, position-weighted local alignment miRanda algorithm (standalone version 3.3a) [80, 81] and the MapMi webserver (version 1.5.0-build 01, release March 2012) (http://www.ebi.ac.uk/enright-srv/MapMi)  to identify candidate miRNA target sites in 3'-UTR regions of 454 reads by base complementarity, and putative miRNA loci, respectively. MapMi is a tool designed to locate miRNA precursor sequences in existing genomic sequences, using potential mature miRNA sequences as input. miRanda is an algorithm for finding genomic targets for microRNAs that incorporates current biological knowledge on target rules and computes optimal sequence complementarity between a set of mature microRNAs and a given mRNA using a weighted dynamic programming algorithm. In addition, miRanda uses an estimate of the free energy of formation of the miRNA:mRNA duplex as a secondary filter . The following parameters were set to reduce the estimated false positives to ≤ 9% [80, 81]: total Score >100; ΔG of the intermediate duplex < −19 Kcal/mol; and output of 2 or more 454-contig targets.
The raw 454 GS FLX Titanium reads of C. s. simus adult and neonate venom gland transcriptomes, and the Ion Torrent PGM reads (miRNA sequences) of adult and neonate C. s. simus have been archived as Standard Flowgram Format (sff) files with the NCBI Sequence Read Archive (SRA) (http://www.ncbi.nlm.nih.gov/sra?term=SRA051956) under accession codes SRX143982 and SRX143985, SRX143983 and SRX143984, respectively.
Procedures for snake maintenance, sacrifice, and gland extraction in this study followed the Quality Management Protocols from the Instituto Clodomiro Picado, University of Costa Rica, and comply the Animal Welfare Law #7451, chapters II and III, and the Biodiversity Law, Decree 7788, Republic of Costa Rica. Research permission was allowed under Resolution ACG- SINAC- PI-012-2010 (Ministry of Environment of Costa Rica).
This work has been financed by grant BFU2010-17373 from the Ministerio de Economía y Competitividad, Madrid, Spain; PROMETEO/2010/005 from the Generalitat Valenciana; and projects 741-B2-093 from the Vicerrectoría de Investigación, Universidad de Costa Rica, CRUSA-CSIC (2009CR0021), and CYTED (206AC0281). We also acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI). JD is grateful to Eva Alloza (CIPF) and José Afonso Guerra-Assunçâo (EBI) for useful comments and advice regarding miRNA sequence analysis.
- Fry BG, Vidal N, Norman JA, Vonk FJ, Scheib H, Ramjan SF, Kuruppu S, Fung K, Hedges SB, Richardson MK, Hodgson WC, Ignjatovic V, Summerhayes R, Kochva E: Early evolution of the venom system in lizards and snakes. Nature. 2006, 439: 584-588. 10.1038/nature04328.PubMedGoogle Scholar
- Fry BG, Casewell NR, Wüster W, Vidal N, Young B, Jackson TN: The structural and functional diversification of the Toxicofera reptile venom system. Toxicon. 2012, 60: 434-448. 10.1016/j.toxicon.2012.02.013.PubMedGoogle Scholar
- Vonk FJ, Admiraal JF, Jackson K, Reshef R, de Bakker MA, Vanderschoot K, van den Berge I, van Atten M, Burgerhout E, Beck A, Mirtschin PJ, Kochva E, Witte F, Fry BG, Woods AE, Richardson MK: Evolutionary origin and development of snake fangs. Nature. 2008, 454: 630-633. 10.1038/nature07178.PubMedGoogle Scholar
- Greene HW: Dietary correlates of the origin and radiation of snakes. Am Zool. 1983, 23: 431-441.Google Scholar
- Daltry JC, Wüster W, Thorpe RS: Diet and snake venom evolution. Nature. 1996, 379: 537-540. 10.1038/379537a0.PubMedGoogle Scholar
- Gibbs HL, Mackessy SP: Functional basis of a molecular adaptation: prey-specific toxic effects of venom from Sistrurus rattlesnakes. Toxicon. 2008, 53: 672-679.Google Scholar
- Barlow A, Pook CE, Harrison RA, Wüster W: Coevolution of diet and prey-specific venom activity supports the role of selection in snake venom evolution. Proc Biol Sci. 2009, 276: 2443-2449. 10.1098/rspb.2009.0048.PubMed CentralPubMedGoogle Scholar
- Vidal N, Hedges SB: The molecular evolutionary tree of lizards, snakes, and amphisbaenians. C Rendus Biol. 2009, 332: 129-139. 10.1016/j.crvi.2008.07.010.Google Scholar
- Vidal N, Rage J-C, Couloux A, Hedges SB: Snakes (Serpentes). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford Univ. Press, 390-397.Google Scholar
- Redi F: Osservazioni Intorno alle Vipere. Edited by: All’Insegna Della Stella F. 1664, a digitalized version is freely accessible in http://archive.org/details/osservazioniint00redigoogGoogle Scholar
- Fontana F: Traité sur le vénin de la vipere, sur les poisons amaricains, sur le laurier-cerise et sur quelques autres poisons végetaux. 1781, Florence: GibelinGoogle Scholar
- Hawgood BJ: Abbé Felice Fontana (1730–1805): founder of modern toxinology. Toxicon. 1995, 33: 591-601. 10.1016/0041-0101(95)00006-8.PubMedGoogle Scholar
- Calvete JJ: Venomics, what else?. Toxicon. 2012, 60: 427-433. 10.1016/j.toxicon.2012.05.012.PubMedGoogle Scholar
- de Azevedo ILM J, Diniz MRV, Ho PL: Venom gland transcriptomic analysis. Animal Toxins: State of the Art. Perspectives in Health and Biotechnology. Edited by: De Lima ME, Pimenta AMC, Martin-Euclaire MF, Zingali RB, Rochat H. 2009, Belo Horizonte: Editora UFMG, 693-713.Google Scholar
- Durban J, Juárez P, Angulo Y, Lomonte B, Flores-Diaz M, Alape-Girón A, Sasa M, Sanz L, Gutiérrez JM, Dopazo J, Conesa A, Calvete JJ: Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing. BMC Genomics. 2011, 12: 259-10.1186/1471-2164-12-259.PubMed CentralPubMedGoogle Scholar
- Rokyta DR, Wray KP, Lemmon AR, Lemmon EM, Caudle SB: A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes. Toxicon. 2011, 57: 657-671. 10.1016/j.toxicon.2011.01.008.PubMedGoogle Scholar
- Calvete JJ: Proteomic tools against the neglected pathology of snake bite envenoming. Exp Rev Proteomics. 2011, 8: 739-758. 10.1586/epr.11.61.Google Scholar
- Calvete JJ: Snake venomics: from the inventory of toxins to biology. Toxicon. 2013, in pressGoogle Scholar
- Fry BG, Roelants K, Champagne DE, Scheib H, Tyndall JD, King GF, Nevalainen TJ, Norman JA, Lewis RJ, Norton RS, Renjifo C, de la Vega RC: The toxicogenomic multiverse: convergent recruitment of proteins into animal venoms. Annu Rev Genomics Hum Genet. 2009, 10: 483-511. 10.1146/annurev.genom.9.081307.164356.PubMedGoogle Scholar
- Fry BG, Scheib H, De LM, Junqueira De Azevedo I, Silva DA, Casewell NR: Novel transcripts in the maxillary venom glands of advanced snakes. Toxicon. 2012, 59: 696-708. 10.1016/j.toxicon.2012.03.005.PubMedGoogle Scholar
- Chippaux JP, Williams V, White J: Snake venom variability: methods of study, results, and interpretation. Toxicon. 1991, 29: 1279-1303. 10.1016/0041-0101(91)90116-9.PubMedGoogle Scholar
- Gibbs HL, Sanz L, Calvete JJ: Snake population venomics: proteomics-based analyses of individual variation reveals significant gene regulation effects on venom protein expression in Sistrurus rattlesnakes. J Mol Evol. 2009, 68: 113-125. 10.1007/s00239-008-9186-1.PubMedGoogle Scholar
- Gibbs HL, Sanz L, Chiucchi JE, Farrell TM, Calvete JJ: Proteomic analysis of ontogenetic and diet-related changes in venom composition of juvenile and adult Dusky Pigmy rattlesnakes (Sistrurus miliarius barbouri). J Proteomics. 2011, 74: 2169-2179. 10.1016/j.jprot.2011.06.013.PubMedGoogle Scholar
- Chiucchi JE, Gibbs HL: Similarity of contemporary and historical gene flow among highly fragmented populations of an endangered rattlesnake. Mol Ecol. 2010, 19: 5345-5358. 10.1111/j.1365-294X.2010.04860.x.PubMedGoogle Scholar
- Mackessy SP: Venom ontogeny in the Pacific rattlesnakes, Crotalus viridis helleri and C. viridis oreganus. Copeia. 1988, 1988: 92-101. 10.2307/1445927.Google Scholar
- Guércio RAP, Shevchenko A, Shevchenko A, López-Lozano JL, Paba J, Sousa MV, Ricart CAO: Ontogenetic variations in the venom proteome of the Amazonian snake Bothrops atrox. Proteome Sci. 2006, 4: 11-10.1186/1477-5956-4-11.PubMed CentralPubMedGoogle Scholar
- Mackessy SP, Sixberry NM, Heyborne WH, Fritts T: Venom of the Brown Treesnake, Boiga irregularis: ontogenetic shifts and taxa-specific toxicity. Toxicon. 2006, 47: 537-548. 10.1016/j.toxicon.2006.01.007.PubMedGoogle Scholar
- Alape-Girón A, Sanz L, Escolano J, Flores-Díaz M, Madrigal M, Sasa M, Calvete JJ: Snake venomics of the lancehead pitviper Bothrops asper: geographic, individual, and ontogenetic variations. J Proteome Res. 2008, 7: 3556-3571. 10.1021/pr800332p.PubMedGoogle Scholar
- Calvete JJ, Sanz L, Cid P, De La Torre P, Flores-Diaz M, Dos Santos MC, Borges A, Bremo A, Angulo Y, Lomonte B, Alape-Girón A, Gutiérrez JM: Snake venomics of the Central American rattlesnake Crotalus simus and the South American Crotalus durissus complex points to neurotoxicity as an adaptive paedomorphic trend along Crotalus dispersal in South America. J Proteome Res. 2010, 9: 528-544. 10.1021/pr9008749.PubMedGoogle Scholar
- Zelanis A, Travaglia-Cardoso SR, Furtado MFD: Ontogenetic changes in Bothrops insularis (Serpentes: Viperidae) snake venom and its biological implications. South Am J Herpetol. 2008, 3: 45-50.Google Scholar
- Zelanis A, Tashima AK, Pinto AF, Leme AF, Stuginski DR, Furtado MF, Sherman NE, Ho PL, Fox JW, Serrano SM: Bothrops jararaca venom proteome rearrangement upon neonate to adult transition. Proteomics. 2011, 11: 4218-4228. 10.1002/pmic.201100287.PubMedGoogle Scholar
- Madrigal M, Sanz L, Flores-Diaz M, Sasa M, Núñez V, Alape-Girón A, Calvete JJ: Snake venomics across genus Lachesis. Ontogenetic changes in the venom composition of L. stenophrys and comparative proteomics of the venoms of adult L. melanocephala and L. acrochorda. J Proteomics. in pressGoogle Scholar
- Andrade DV, Abe AS: Relationship of venom ontogeny and diet in Bothrops. Herpetologica. 1999, 55: 200-204.Google Scholar
- Wüster W, Ferguson JE, Quijada-Mascareñas JA, Pook CE, Salomão MG, Thorpe RS: Tracing an invasion: landbridges, refugia and the phylogeography of the Neotropical rattlesnake (Serpentes: Viperidae: Crotalus durissus). Mol Ecol. 2005, 14: 1095-1108. 10.1111/j.1365-294X.2005.02471.x.PubMedGoogle Scholar
- Gutiérrez JM, Dos Santos MC, Furtado MF, Rojas G: Biochemical and pharmacological similarities between the venoms of newborn Crotalus durissus durissus and adult Crotalus durissus terrificus rattlesnakes. Toxicon. 1991, 29: 1273-1277. 10.1016/0041-0101(91)90201-2.PubMedGoogle Scholar
- Saravia P, Rojas E, Arce V, Guevara C, López JC, Chaves E, Velásquez R, Rojas G, Gutiérrez JM: Geographic and ontogenic variability in the venom of the neotropical rattlesnake Crotalus durissus: Pathophysiological and therapeutic implications. Rev Biol Trop. 2002, 50: 337-346.PubMedGoogle Scholar
- Mackessy SP: Venom composition in rattlesnakes: trends and biological significance. The Biology of Rattlesnakes. Edited by: Hayes WK, Beaman KR, Cardwell MD, Bush SP. 2008, Loma Linda, California: Loma Linda University Press, 495-510.Google Scholar
- Mackessy SP: Evolutionary trends in venom composition in the Western Rattlesnakes (Crotalus viridis sensu lato): Toxicity vs. tenderizers. Toxicon. 2010, 55: 1463-1474. 10.1016/j.toxicon.2010.02.028.PubMedGoogle Scholar
- McCue MD: Cost of producing venom in three North American pitviper species. Copeia. 2006, 2006: 818-825. 10.1643/0045-8511(2006)6[818:COPVIT]2.0.CO;2.Google Scholar
- Bon C: Multicomponent neurotoxic phospholipases A2. Venom phospholipase A2 enzymes: structure, function and mechanism. Edited by: Kini RM. 1997, Chichester: Wiley, 269-286.Google Scholar
- Faure G, Xu H, Saul FA: Crystal structure of crotoxin reveals key residues involved in the stability and toxicity of this potent heterodimeric β-neurotoxin. J Mol Biol. 2011, 412: 176-191. 10.1016/j.jmb.2011.07.027.PubMedGoogle Scholar
- Faure G, Saul F: Crystallographic characterization of functional sites of crotoxin and ammodytoxin, potent β-neurotoxins from Viperidae venom. Toxicon. 2012, 60: 531-538. 10.1016/j.toxicon.2012.05.009.PubMedGoogle Scholar
- Gutiérrez JM: Clinical toxicology of snakebite in Central America. Handbook of clinical toxicology of animal venoms and poisons. Edited by: Meier J, White J. 1995, Boca Raton, Florida: CRC, 645-665.Google Scholar
- Massey DJ, Calvete JJ, Sánchez EE, Sanz L, Richards K, Curtis R, Boesen K: Venom variability and envenoming severity outcomes of the Crotalus scutulatus scutulatus (Mojave rattlesnake) from Southern Arizona. J Proteomics. 2012, 75: 2576-2587. 10.1016/j.jprot.2012.02.035.PubMedGoogle Scholar
- Glenn JL, Sraight RC, Wolfe MC, Hardy DL: Geographical variation in Crotalus scutulatus scutulatus (Mojave rattlesnake) venom properties. Toxicon. 1983, 21: 119-130.PubMedGoogle Scholar
- Glenn JL, Sraight RC: Intergradation of two different venom populations of the Mojave rattlesnake (Crotalus scutulatus scutulatus) in Arizona. Toxicon. 1989, 27: 411-418. 10.1016/0041-0101(89)90203-1.PubMedGoogle Scholar
- Gutiérrez JM, Lomonte B, León G, Alape-Girón A, Flores-Díaz M, Sanz L, Angulo Y, Calvete JJ: Snake venomics and antivenomics: Proteomic tools in the design and control of antivenoms for the treatment of snakebite envenoming. J Proteomics. 2009, 72: 165-182. 10.1016/j.jprot.2009.01.008.PubMedGoogle Scholar
- Williams DJ, Gutiérrez JM, Calvete JJ, Wüster W, Ratanabanangkoon K, Paiva O, Brown NI, Casewell NR, Harrison RA, Rowley PD, O'Shea M, Jensen SD, Winkel KD, Warrell DA: Ending the drought: new strategies for improving the flow of affordable, effective antivenoms in Asia and Africa. J Proteomics. 2011, 74: 1735-1767. 10.1016/j.jprot.2011.05.027.PubMedGoogle Scholar
- Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM: Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 2007, 8: R143-10.1186/gb-2007-8-7-r143.PubMed CentralPubMedGoogle Scholar
- Biémont C, Vieira C: Junk DNA as an evolutionary force. Nature. 2006, 443: 521-524. 10.1038/443521a.PubMedGoogle Scholar
- Castoe TA, Hall KT, Guibotsy Mboulas ML, Gu W, de Koning AP, Fox SE, Poole AW, Vemulapalli V, Daza JM, Mockler T, Smith EN, Feschotte C, Pollock DD: Discovery of highly divergent repeat landscapes in snake genomes using high-throughput sequencing. Genome Biol Evol. 2011, 3: 641-653. 10.1093/gbe/evr043.PubMed CentralPubMedGoogle Scholar
- Junqueira-de-Azevedo ILM, Ho PL: A survey of gene expression and diversity in the venom glands of the pit viper snake Bothrops insularis through the generation of expressed sequence tags (ESTs). Gene. 2002, 299: 279-291. 10.1016/S0378-1119(02)01080-6.Google Scholar
- Junqueira-de-Azevedo ILM, Ching ATC, Carvalho E, Faria F, Nishiyama MY, Ho PL, Diniz MRV: Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution. Genetics. 2006, 172: 877-889.Google Scholar
- Ching AT, Rocha MM, Leme AFP, Pimenta DC, Furtado MFD, Serrano SM, Ho PL, de Azevedo IL J: Some aspects of the venom proteome of the Colubridae snake Philodryas olfersii revealed from a Duvernoy’s (venom) gland transcriptome. FEBS Lett. 2006, 580: 4417-4422. 10.1016/j.febslet.2006.07.010.PubMedGoogle Scholar
- Kordis D, Gubenšek F: Bov-B long interspersed repeated DNA (LINE) sequences are present in Vipera ammodytes phospholipase A2 genes and in genomes of Viperidae snakes. Eur J Biochem. 1997, 246: 772-779. 10.1111/j.1432-1033.1997.00772.x.PubMedGoogle Scholar
- Kordis D, Gubenšek F: The Bov-B lines found in Vipera ammodytes toxic PLA2 genes are widespread in snake genomes. Toxicon. 1998, 36: 1585-1590. 10.1016/S0041-0101(98)00150-0.PubMedGoogle Scholar
- Ikeda N, Chijiwa T, Matsubara K, Oda-Ueda N, Hattori S, Matsuda Y, Ohno M: Unique structural characteristics and evolution of a cluster of venom phospholipase A2 isozyme genes of Protobothrops flavoviridis snake. Gene. 2010, 461: 15-25. 10.1016/j.gene.2010.04.001.PubMedGoogle Scholar
- Nobuhisa I, Ogawa T, Deshimaru M, Chijiwa T, Nakashima K-I, Chuman Y, Shimohigashi Y, Fukumaki Y, Hattori S, Ohno M: Retrotransposable CR1-like elements in Crotalinae snake genomes. Toxicon. 1998, 36: 915-920. 10.1016/S0041-0101(97)00104-9.PubMedGoogle Scholar
- Sanz L, Harrison RA, Calvete JJ: First draft of the genomic organization of a PIII-SVMP gene. Toxicon. 2012, 60: 455-469. 10.1016/j.toxicon.2012.04.331.PubMedGoogle Scholar
- van de Lagemaat LN, Landry JR, Mager DL, Medstrand P: Transposable elements in mammals promote regulatory variation and diversification of genes with specialized functions. Trends Genet. 2003, 19: 530-536. 10.1016/j.tig.2003.08.004.PubMedGoogle Scholar
- Medstrand P, Van de Lagemaat LN, Dunn CA, Landry JR, Svenback D, Mager DL: Impact of transposable elements on the evolution of mammalian gene regulation. Cytogenet Genome Res. 2005, 110: 342-345. 10.1159/000084966.PubMedGoogle Scholar
- Piskurek O, Austin CC, Okada N: Sauria SINEs: novel short interspersed transposable elements that are widespread in reptile genomes. J Mol Evol. 2006, 62: 630-644. 10.1007/s00239-005-0201-5.PubMedGoogle Scholar
- Piskurek O, Okada N: Poxviruses as possible vectors for horizontal transfer of retroposons from reptiles to mammals. Proc Natl Acad Sci USA. 2007, 104: 12046-12051. 10.1073/pnas.0700531104.PubMed CentralPubMedGoogle Scholar
- Oshima K, Okada N: SINEs and LINEs: symbionts of eukaryotic genomes with a common tail. Cytogenet Genome Res. 2005, 110: 475-490. 10.1159/000084981.Google Scholar
- Luan DD, Korman MH, Jakubczak JL, Eickbush TH: Reverse transcription of R2Bm RNA is primed by a nick at the chromosomal target site: a mechanism for non-LTR retrotransposition. Cell. 1993, 72: 595-605. 10.1016/0092-8674(93)90078-5.PubMedGoogle Scholar
- Deininger PL, Batzer MA: Mammalian retroelements. Genome Res. 2002, 12: 1455-1465. 10.1101/gr.282402.PubMedGoogle Scholar
- Piskurek O, Nishihara H, Okada N: The evolution of two partner LINE/SINE families and a full-length chromodomain-containing Ty3/Gypsy LTR element in the first reptilian genome of Anolis carolinensis. Gene. 2009, 441: 111-118. 10.1016/j.gene.2008.11.030.PubMedGoogle Scholar
- International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.Google Scholar
- Alföldi J, Di Palma F, Grabherr M, Williams C, Kong L, Mauceli E, Russell P, Lowe CB, Glor RE, Jaffe JD, Ray DA, Boissinot S, Shedlock AM, Botka C, Castoe TA, Colbourne JK, Fujita MK, Moreno RG, ten Hallers BF, Haussler D, Heger A, Heiman D, Janes DE, Johnson J, de Jong PJ, Koriabine MY, Lara M, Novick PA, Organ CL, Peach SE: The genome of the green anole lizard and a comparative analysis with birds and mammals. Nature. 2011, 477: 587-591. 10.1038/nature10390.PubMed CentralPubMedGoogle Scholar
- Kazazian HHJ: Mobile elements: drivers of genome evolution. Science. 2004, 303: 1626-1632. 10.1126/science.1089670.PubMedGoogle Scholar
- Cline M, Smooth M, Cerami E, Kushinski A, Landys N, Workman C, Christmas R, Avila-Campilo I, Creech M, Gross B, Hanspers K, Isserlin R, Kelley R, Killcoyne S, Lotia S, Maere S, Morris J, Ono K, Pavlovic V, Pico AR, Vailaya A, Wang PL, Adler A, Conklin BR, Hood L, Kuiper M, Sander C, Schmulevich I, Schwikowski B, Warner GJ, Ideker T, Bader GD: Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007, 2: 2366-2382. 10.1038/nprot.2007.324.PubMed CentralPubMedGoogle Scholar
- Sharp PA: The centrality of RNA. Cell. 2009, 136: 577-580. 10.1016/j.cell.2009.02.007.PubMedGoogle Scholar
- Licatalosi DD, Darnell RB: RNA processing and its regulation: global insights into biological networks. Nat Rev Genet. 2010, 11: 75-87.PubMed CentralPubMedGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355. 10.1038/nature02871.PubMedGoogle Scholar
- He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.PubMedGoogle Scholar
- Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136: 215-233. 10.1016/j.cell.2009.01.002.PubMed CentralPubMedGoogle Scholar
- Höck H, Meister G: The Argonaute protein family. Genome Biol. 2008, 9: 210-10.1186/gb-2008-9-2-210.PubMed CentralPubMedGoogle Scholar
- Wheeler BM, Heimberg AM, Moy VN, Sperling EA, Holstein TW, Heber S, Peterson KJ: The deep evolution of metazoan microRNAs. Evol Dev. 2009, 11: 145-150.Google Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acid Res. 2011, 39: D152-D157. 10.1093/nar/gkq1027.PubMed CentralPubMedGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5: R1-10.1186/gb-2003-5-1-r1.PubMed CentralPubMedGoogle Scholar
- John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2: e363-10.1371/journal.pbio.0020363.PubMed CentralPubMedGoogle Scholar
- Wuchty S, Fontana W, Hofacker IL, Schuster P: Complete suboptimal folding of RNA and the stability of secondary structures. Biopolymers. 1999, 49: 145-165. 10.1002/(SICI)1097-0282(199902)49:2<145::AID-BIP4>3.0.CO;2-G.PubMedGoogle Scholar
- Guerra-Assunção JA, Enright AJ: MapMi: automated mapping of microRNA loci. BMC Bioinformatics. 2012, 11: 133-Google Scholar
- Chen K, Rajewsky N: Deep conservation of microRNA-target relationships and 3' UTR motifs in vertebrates, flies, and nematodes. Cold Spring Harb Symp Quant Biol. 2006, 71: 149-156. 10.1101/sqb.2006.71.039.PubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.PubMedGoogle Scholar
- Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27: 91-105. 10.1016/j.molcel.2007.06.017.PubMed CentralPubMedGoogle Scholar
- Ma JB, Yuan YR, Meister G, Pei Y, Tuschi T, Patel J: Structural basis for 5'-end-specific recognition of guide RNA by the A. fulgidus Piwi protein. Nature. 2005, 434: 666-670. 10.1038/nature03514.PubMed CentralPubMedGoogle Scholar
- Parker JS, Roe SM, Barford D: Structural insight into mRNA recognition from Piwi domain-siRNA guide complex. Nature. 2005, 434: 663-666. 10.1038/nature03462.PubMed CentralPubMedGoogle Scholar
- Bon C, Changeux JP, Jeng TW, Fraenkel-Conrat H: Postsynaptic effects of crotoxin and its isolated subunits. Eur J Biochem. 1979, 99: 471-481. 10.1111/j.1432-1033.1979.tb13278.x.PubMedGoogle Scholar
- Faure G, Copic A, Le Porrier S, Gubensek F, Bon C, Krizaj I: Crotoxin acceptor protein isolated from Torpedo electric organ: binding properties to crotoxin by surface plasmon resonance. Toxicon. 2003, 41: 509-517. 10.1016/S0041-0101(02)00394-X.PubMedGoogle Scholar
- Pereañez JA, Gómez ID, Patiño AC: Relationship between the structure and the enzymatic activity of crotoxin complex and its phospholipase A2 subunit: An in silico approach. J Mol Graph Model. 2012, 35: 36-42.PubMedGoogle Scholar
- Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455: 64-71. 10.1038/nature07242.PubMed CentralPubMedGoogle Scholar
- Vidal N, Hedges SB: The phylogeny of squamate reptiles (lizards, snakes, and amphisbaenians) inferred from nine nuclear protein-coding genes. Comptes Rendus Biologies. 2005, 328: 1000-1008. 10.1016/j.crvi.2005.10.001.PubMedGoogle Scholar
- Paine MJ, Desmond HP, Theakston RDG, Crampton JM: Gene expression in Echis carinatus (carpet viper) venom glands following milking. Toxicon. 1992, 30: 379-386. 10.1016/0041-0101(92)90534-C.PubMedGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen Y, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Ho CH, Irzyk GP, Jando SC, Alenquer MLI, Jarvie TP, Jirage KB, Kim J, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.PubMed CentralPubMedGoogle Scholar
- Schmieder R, Edwards R: Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011, 27: 863-864. 10.1093/bioinformatics/btr026.PubMed CentralPubMedGoogle Scholar
- Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. http://repeatmasker.org,
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedGoogle Scholar
- Papadopoulos JS, Agarwala R: COBALT: constraint-based alignment tool for multiple protein sequences. Bioinformatics. 2007, 23: 1073-1079. 10.1093/bioinformatics/btm076.PubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.PubMedGoogle Scholar
- Tarazona S, García-Alcalde F, Dopazo J, Ferrer A, Conesa A: Differential expression in RNA-seq: a matter of depth. Genome Res. 2011, 21: 2213-2223. 10.1101/gr.124321.111.PubMed CentralPubMedGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.PubMedGoogle Scholar
- Niu BF, Fu LM, Sun SL, Li WZ: Artificial and natural duplicates in pyrosequencing reads of metagenomic data. BMC Bioinformatics. 2010, 11: 187-10.1186/1471-2105-11-187.PubMed CentralPubMedGoogle Scholar
- Stefani G, Slack FJ: A ‘pivotal’ new rule for microRNA-mRNA interactions. Nature Struct & Mol Biol. 2012, 19: 265-266. 10.1038/nsmb.2256.Google Scholar
- Chi SW, Hannon GJ, Darnell RB: An alternative mode of microRNA target recognition. Nature Struct & Mol Biol. 2012, 19: 321-328. 10.1038/nsmb.2230.Google Scholar
- Seelbach M, Schwanhausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455: 58-63. 10.1038/nature07228.Google 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.