Unique structural characteristics and evolution of a cluster of venom phospholipase A2 isozyme genes of Protobothrops flavoviridis snake

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

Molecular diversity and evolutionary trends of cysteine-rich peptides from the venom glands of Chinese spider Heteropoda venatoria

Heteropoda venatoria in the family Sparassidae is highly valued in pantropical countries because the species feed on domestic insect pests. Unlike most other species of Araneomorphae, H. venatoria uses the great speed and strong chelicerae (mouthparts) with toxin glands to capture the insects instead of its web. Therefore, H. venatoria provides unique opportunities for venom evolution research. The venom of H. venatoria was explored by matrix-assisted laser desorption/ionization tandem time-of-flight and analyzing expressed sequence tags. The 154 sequences coding cysteine-rich peptides (CRPs) revealed 24 families based on the phylogenetic analyses of precursors and cysteine frameworks in the putative mature regions. Intriguingly, four kinds of motifs are first described in spider venom. Furthermore, combining the diverse CRPs of H. venatoria with previous spider venom peptidomics data, the structures of precursors and the patterns of cysteine frameworks were analyzed. This work revealed the dynamic evolutionary trends of venom CRPs in H. venatoria: the precursor has evolved an extended mature peptide with more cysteines, and a diminished or even vanished propeptides between the signal and mature peptides; and the CRPs evolved by multiple duplications of an ancestral ICK gene as well as recruitments of non-toxin genes.

Alternative mRNA Splicing in Three Venom Families Underlying a Possible Production of Divergent Venom Proteins of the Habu Snake, Protobothrops flavoviridis

Snake venoms are complex mixtures of toxic proteins encoded by various gene families that function synergistically to incapacitate prey. A huge repertoire of snake venom genes and proteins have been reported, and alternative splicing is suggested to be involved in the production of divergent gene transcripts. However, a genome-wide survey of the transcript repertoire and the extent of alternative splicing still remains to be determined. In this study, the comprehensive analysis of transcriptomes in the venom gland was achieved by using PacBio sequencing. Extensive alternative splicing was observed in three venom protein gene families, metalloproteinase (MP), serine protease (SP), and vascular endothelial growth factors (VEGF). Eleven MP and SP genes and a VEGF gene are expressed as a total of 81, 61, and 8 transcript variants, respectively. In the MP gene family, individual genes are transcribed into different classes of MPs by alternative splicing. We also observed trans-splicing among the clustered SP genes. No other venom genes as well as non-venom counterpart genes exhibited alternative splicing. Our results thus indicate a potential contribution of mRNA alternative and trans-splicing in the production of highly variable transcripts of venom genes in the habu snake.

Unique structure (construction and configuration) and evolution of the array of small serum protein genes of Protobothrops flavoviridis snake

The nucleotide sequence of Protobothrops flavoviridis (Pf) 30534 bp genome segment which contains genes encoding small serum proteins (SSPs) was deciphered. The genome segment contained five SSP genes (PfSSPs), PfSSP-4, PfSSP-5, PfSSP-1, PfSSP-2, and PfSSP-3 in this order and had characteristic configuration and constructions of the particular nucleotide sequences inserted. Comparison between the configurations of the inserted chicken repeat-1 (CR1) fragments of P. flavoviridis and Ophiophagus hannah (Oh) showed that the nucleotide segment encompassing from PfSSP-1 to PfSSP-2 was inverted. The inactive form of PfSSP-1, named PfSSP-1δ(Ψ), found in the intergenic region (I-Reg) between PfSSP-5 and PfSSP-1 had also been destroyed by insertions of the plural long interspersed nuclear elements (LINEs) and DNA transposons. The L2 LINE inserted into the third intron or the particular repetitive sequences inserted into the second intron structurally divided five PfSSPs into two subgroups, the Long SSP subgroup of PfSSP-1, PfSSP-2 and PfSSP-5 or the Short SSP subgroup of PfSSP-3 and PfSSP-4. The mathematical analysis also showed that PfSSPs of the Long SSP subgroup evolved alternately in an accelerated and neutral manner, whereas those of the Short SSP subgroup evolved in an accelerated manner. Moreover, the ortholog analysis of SSPs of various snakes showed that the evolutionary emerging order of SSPs was as follows: SSP-5, SSP-4, SSP-2, SSP-1, and SSP-3. The unique interpretation about accelerated evolution and the novel idea that the transposable elements such as LINEs and DNA transposons are involved in maintaining the host genome besides its own transposition natures were proposed.

The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes

Here we use a chromosome-level genome assembly of a prairie rattlesnake (Crotalus viridis), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies.

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

Squamate reptiles challenge paradigms of genomic repeat element evolution set by birds and mammals

Broad paradigms of vertebrate genomic repeat element evolution have been largely shaped by analyses of mammalian and avian genomes. Here, based on analyses of genomes sequenced from over 60 squamate reptiles (lizards and snakes), we show that patterns of genomic repeat landscape evolution in squamates challenge such paradigms. Despite low variance in genome size, squamate genomes exhibit surprisingly high variation among species in abundance (ca. 25–73% of the genome) and composition of identifiable repeat elements. We also demonstrate that snake genomes have experienced microsatellite seeding by transposable elements at a scale unparalleled among eukaryotes, leading to some snake genomes containing the highest microsatellite content of any known eukaryote. Our analyses of transposable element evolution across squamates also suggest that lineage-specific variation in mechanisms of transposable element activity and silencing, rather than variation in species-specific demography, may play a dominant role in driving variation in repeat element landscapes across squamate phylogeny.

Large-scale patterns of genomic repeat element evolution have been studied mainly in birds and mammals. Here, the authors analyze the genomes of over 60 squamate reptiles and show high variation in repeat elements compared to mammals and birds, and particularly high microsatellite seeding in snakes.

Population Genomic Analysis of a Pitviper Reveals Microevolutionary Forces Underlying Venom Chemistry

Venoms are among the most biologically active secretions known, and are commonly believed to evolve under extreme positive selection. Many venom gene families, however, have undergone duplication, and are often deployed in doses vastly exceeding the LD₅₀ for most prey species, which should reduce the strength of positive selection. Here, we contrast these selective regimes using snake venoms, which consist of rapidly evolving protein formulations. Though decades of extensive studies have found that snake venom proteins are subject to strong positive selection, the greater action of drift has been hypothesized, but never tested. Using a combination of de novo genome sequencing, population genomics, transcriptomics, and proteomics, we compare the two modes of evolution in the pitviper, Protobothrops mucrosquamatus. By partitioning selective constraints and adaptive evolution in a McDonald–Kreitman-type framework, we find support for both hypotheses: venom proteins indeed experience both stronger positive selection, and lower selective constraint than other genes in the genome. Furthermore, the strength of selection may be modulated by expression level, with more abundant proteins experiencing weaker selective constraint, leading to the accumulation of more deleterious mutations. These findings show that snake venoms evolve by a combination of adaptive and neutral mechanisms, both of which explain their extraordinarily high rates of molecular evolution. In addition to positive selection, which optimizes efficacy of the venom in the short term, relaxed selective constraints for deleterious mutations can lead to more rapid turnover of individual proteins, and potentially to exploration of a larger venom phenotypic space.

The Deep Origin and Recent Loss of Venom Toxin Genes in Rattlesnakes

The genetic origin of novel traits is a central but challenging puzzle in evolutionary biology. Among snakes, phospholipase A2 (PLA2)-related toxins have evolved in different lineages to function as potent neurotoxins, myotoxins, or hemotoxins. Here, we traced the genomic origin and evolution of PLA2 toxins by examining PLA2 gene number, organization, and expression in both neurotoxic and non-neurotoxic rattlesnakes. We found that even though most North American rattlesnakes do not produce neurotoxins, the genes of a specialized heterodimeric neurotoxin predate the origin of rattlesnakes and were present in their last common ancestor (∼22 mya). The neurotoxin genes were then deleted independently in the lineages leading to the Western Diamondback (Crotalus atrox) and Eastern Diamondback (C. adamanteus) rattlesnakes (∼6 mya), while a PLA2 myotoxin gene retained in C. atrox was deleted from the neurotoxic Mojave rattlesnake (C. scutulatus; ∼4 mya). The rapid evolution of PLA2 gene number appears to be due to transposon invasion that provided a template for non-allelic homologous recombination.

Insights into the Evolution of a Snake Venom Multi-Gene Family from the Genomic Organization of Echis ocellatus SVMP Genes

The molecular events underlying the evolution of the Snake Venom Metalloproteinase (SVMP) family from an A Disintegrin And Metalloproteinase (ADAM) ancestor remain poorly understood. Comparative genomics may provide decisive information to reconstruct the evolutionary history of this multi-locus toxin family. Here, we report the genomic organization of Echis ocellatus genes encoding SVMPs from the PII and PI classes. Comparisons between them and between these genes and the genomic structures of Anolis carolinensis ADAM28 and E. ocellatus PIII-SVMP EOC00089 suggest that insertions and deletions of intronic regions played key roles along the evolutionary pathway that shaped the current diversity within the multi-locus SVMP gene family. In particular, our data suggest that emergence of EOC00028-like PI-SVMP from an ancestral PII(e/d)-type SVMP involved splicing site mutations that abolished both the 3′ splice AG acceptor site of intron 12* and the 5′ splice GT donor site of intron 13*, and resulted in the intronization of exon 13* and the consequent destruction of the structural integrity of the PII-SVMP characteristic disintegrin domain.

Interisland variegation of venom [Lys(49)]phospholipase A2 isozyme genes in Protobothrops genus snakes in the southwestern islands of Japan

Protobothrops tokarensis (Pt), a Crotalinae snake, inhabits only Takarajima and Kodakarajima islands of the Tokara Islands located in the immediate north of Amami-Oshima island of Japan. Kodakarajima P. tokarensis venom gland cDNA library gave four types of phospholipase A2 (PLA2) cDNAs encoding neutral [Asp(49)]PLA2, basic [Asp(49)]PLA2, highly basic [Asp(49)]PLA2, and [Lys(49)]PLA2. As the amino acid sequences encoded by their open reading frames (ORFs) were identical to those of PLA2, PLA-B, PLA-N, and BPI (a [Lys(49)]PLA2), respectively, from Amami-Oshima P. flavoviridis (Pf) venom, they were named PtPLA2, PtPLA-B, PtPLA-N, and PtBPI. Chromatography of P. tokarensis venom gave three PLA2 isozymes, PtPLA2, PtPLA-B, and PtBPI. However, BPII and BPIII ([Lys(49)]PLA2s) expressed in Amami-Oshima P. flavoviridis venom were not found in P. tokarensis venom. Genomic polymerase chain reaction (PCR) for P. tokarensis liver DNAs with the unique primers gave PtBPI gene. Notably it was found that LINE (long interspersed nuclear element)-1 fragment is inserted into second intron of PtBPI gene. The LINE-1 fragment may prevent duplication of PtBPI gene and thus formation of plural [Lys(49)]PLA2 genes in P. tokarensis genome. The interisland variegation of venom [Lys(49)]PLA2 isozyme genes in Protobothrops genus snakes in the southwestern islands of Japan is discussed.

Snake venoms are integrated systems, but abundant venom proteins evolve more rapidly

Background

While many studies have shown that extracellular proteins evolve rapidly, how selection acts on them remains poorly understood. We used snake venoms to understand the interaction between ecology, expression level, and evolutionary rate in secreted protein systems. Venomous snakes employ well-integrated systems of proteins and organic constituents to immobilize prey. Venoms are generally optimized to subdue preferred prey more effectively than non-prey, and many venom protein families manifest positive selection and rapid gene family diversification. Although previous studies have illuminated how individual venom protein families evolve, how selection acts on venoms as integrated systems, is unknown.

Results

Using next-generation transcriptome sequencing and mass spectrometry, we examined microevolution in two pitvipers, allopatrically separated for at least 1.6 million years, and their hybrids. Transcriptomes of parental species had generally similar compositions in regard to protein families, but for a given protein family, the homologs present and concentrations thereof sometimes differed dramatically. For instance, a phospholipase A₂ transcript comprising 73.4 % of the Protobothrops elegans transcriptome, was barely present in the P. flavoviridis transcriptome (<0.05 %). Hybrids produced most proteins found in both parental venoms. Protein evolutionary rates were positively correlated with transcriptomic and proteomic abundances, and the most abundant proteins showed positive selection. This pattern holds with the addition of four other published crotaline transcriptomes, from two more genera, and also for the recently published king cobra genome, suggesting that rapid evolution of abundant proteins may be generally true for snake venoms. Looking more broadly at Protobothrops, we show that rapid evolution of the most abundant components is due to positive selection, suggesting an interplay between abundance and adaptation.

Conclusions

Given log-scale differences in toxin abundance, which are likely correlated with biosynthetic costs, we hypothesize that as a result of natural selection, snakes optimize return on energetic investment by producing more of venom proteins that increase their fitness. Natural selection then acts on the additive genetic variance of these components, in proportion to their contributions to overall fitness. Adaptive evolution of venoms may occur most rapidly through changes in expression levels that alter fitness contributions, and thus the strength of selection acting on specific secretome components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1832-6) contains supplementary material, which is available to authorized users.

The finding of a group IIE phospholipase A2 gene in a specified segment of Protobothrops flavoviridis genome and its possible evolutionary relationship to group IIA phospholipase A2 genes

The genes encoding group IIE phospholipase A₂, abbreviated as IIE PLA₂, and its 5' and 3' flanking regions of Crotalinae snakes such as Protobothrops flavoviridis, P. tokarensis, P. elegans, and Ovophis okinavensis, were found and sequenced. The genes consisted of four exons and three introns and coded for 22 or 24 amino acid residues of the signal peptides and 134 amino acid residues of the mature proteins. These IIE PLA₂s show high similarity to those from mammals and Colubridae snakes. The high expression level of IIE PLA₂s in Crotalinae venom glands suggests that they should work as venomous proteins. The blast analysis indicated that the gene encoding OTUD3, which is ovarian tumor domain-containing protein 3, is located in the 3' downstream of IIE PLA₂ gene. Moreover, a group IIA PLA₂ gene was found in the 5' upstream of IIE PLA₂ gene linked to the OTUD3 gene (OTUD3) in the P. flavoviridis genome. It became evident that the specified arrangement of IIA PLA₂ gene, IIE PLA₂ gene, and OTUD3 in this order is common in the genomes of humans to snakes. The present finding that the genes encoding various secretory PLA₂s form a cluster in the genomes of humans to birds is closely related to the previous finding that six venom PLA₂ isozyme genes are densely clustered in the so-called NIS-1 fragment of the P. flavoviridis genome. It is also suggested that venom IIA PLA₂ genes may be evolutionarily derived from the IIE PLA₂ gene.

Restriction and recruitment-gene duplication and the origin and evolution of snake venom toxins

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive "just-so story" in evolutionary biology.

The rise of operon-like gene clusters in plants

Gene clusters are common features of prokaryotic genomes also present in eukaryotes. Most clustered genes known are involved in the biosynthesis of secondary metabolites. Although horizontal gene transfer is a primary source of prokaryotic gene cluster (operon) formation and has been reported to occur in eukaryotes, the predominant source of cluster formation in eukaryotes appears to arise de novo or through gene duplication followed by neo- and sub-functionalization or translocation. Here we aim to provide an overview of the current knowledge and open questions related to plant gene cluster functioning, assembly, and regulation. We also present potential research approaches and point out the benefits of a better understanding of gene clusters in plants for both fundamental and applied plant science.

Tracing monotreme venom evolution in the genomics era

The monotremes (platypuses and echidnas) represent one of only four extant venomous mammalian lineages. Until recently, monotreme venom was poorly understood. However, the availability of the platypus genome and increasingly sophisticated genomic tools has allowed us to characterize platypus toxins, and provides a means of reconstructing the evolutionary history of monotreme venom. Here we review the physiology of platypus and echidna crural (venom) systems as well as pharmacological and genomic studies of monotreme toxins. Further, we synthesize current ideas about the evolution of the venom system, which in the platypus is likely to have been retained from a venomous ancestor, whilst being lost in the echidnas. We also outline several research directions and outstanding questions that would be productive to address in future research. An improved characterization of mammalian venoms will not only yield new toxins with potential therapeutic uses, but will also aid in our understanding of the way that this unusual trait evolves.

Integrated "omics" profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus

BACKGROUND

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.

RESULTS

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 adulthood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus.

CONCLUSIONS

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.

Structural characteristics and evolution of the Protobothrops elegans pancreatic phospholipase A2 gene in contrast with those of Protobothrops genus venom phospholipase A2 genes

The nucleotide sequence of the gene encoding Protobothrops elegans (Crotalinae) pancreatic phospholipase A(2) (PLA(2)), abbreviated PePancPLA(2), was determined by means of inverted PCR techniques. Since its deduced amino acid sequence contains a pancreatic loop and shows high similarity to that of Laticauda semifasciata (Elapinae) group IB pancreatic PLA(2), PePancPLA(2) is classified into group IB PLA(2). The nucleotide sequences of the PePancPLA(2) gene, the L. semifasciata group IB pancreatic PLA(2) gene, and the L. semifasciata group IA venom PLA(2) gene are similar to one another but greatly dissimilar to those of Protobothrops genus (Crotalinae) group II venom PLA(2) genes, suggesting that the Elapinae group IB PLA(2) gene and the group IA PLA(2) gene appeared after Elapinae was established, and that the Crotalinae group II venom PLA(2) genes came into existence before Elapinae and Crotalinae diverged. A phylogenetic analysis of their amino acid sequences confirms this.

First draft of the genomic organization of a PIII-SVMP gene

The evolutionary pathway of highly toxic proteins expressed in snake venom glands from proteins without toxic function and expressed in non-parotid tissues remains poorly understood. Here we examine gene structure of a representative of a venom protein with an ADAMs metalloproteinase evolutionary origin. The structure of the 15,652 bp Echis ocellatus pre-pro EOC00089-like PIII-SVMP gene was assembled from PCR-amplified sequences of overlapping genomic fragments. The gene comprises 12 exons interrupted by 11 introns. In a homology model of the EOC00089-like protein, the insertion of introns interrupting coding regions lie just after or between secondary structure elements. Long interspersed nuclear retroelements (LINE) L2/CR1 and RTE/Bov-B, short interspersed nuclear retroelements SINE/Sauria, and a hobo-activator DNA (Charlie, hAT) transposon were identified within introns 1, 3, 7 and 8. Pairwise amino acid sequence comparisons between EOC00089-like PIII-SVMP and its closest orthologs, ADAM28, from a mammal, Homo sapiens, and the lizard, Anolis carolinensis, showed that the ORFs of these three proteins share 42%/59%, 49%/69%, and 48%/65% (identity/similarity), respectively. The protein-coding positions interrupted by each of the 11 introns of the Echis PIII-SVMP gene are entirely conserved in the A. carolinensis and human ADAM28 genes. However, the lizard and the human ADAM28 genes contain 5 introns not present in the E. ocellatus gene. Furthermore, Echis and Anolis introns exhibit quantitatively and qualitatively distinctions in their inserted retroelements. These findings identify introns as possible key elements in the recruitment and amplification process of SVMPs into the venom gland of extant snakes. Ongoing reptile genome sequencing projects may shed light on this intriguing aspect of the emergence and evolution of venom toxin genes. Furthermore, the organization of the PIII-SVMP reported here provides a genomic explanation for the emergence of dimeric disintegrin subunits encoded by short messengers.

Sequencing the genome of the Burmese python (Python molurus bivittatus) as a model for studying extreme adaptations in snakes

The Consortium for Snake Genomics is in the process of sequencing the genome and creating transcriptomic resources for the Burmese python. Here, we describe how this will be done, what analyses this work will include, and provide a timeline.

A proposal to sequence the genome of a garter snake (Thamnophis sirtalis)

Here we develop an argument in support of sequencing a garter snake (Thamnophis sirtalis) genome, and outline a plan to accomplish this. This snake is a common, widespread, nonvenomous North American species that has served as a model for diverse studies in evolutionary biology, physiology, genomics, behavior and coevolution. The anole lizard is currently the only genome sequence available for a non-avian reptile. Thus, the garter snake at this time would be the first available snake genome sequence and as such would provide much needed comparative representation of non-avian reptilian genomes, and would also allow critical new insights for vertebrate comparative genomic studies. We outline the major areas of discovery that the availability of the garter snake genome would enable, and describe a plan for whole-genome sequencing.

Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing

BACKGROUND

A long term research goal of venomics, of applied importance for improving current antivenom therapy, but also for drug discovery, is to understand the pharmacological potential of venoms. Individually or combined, proteomic and transcriptomic studies have demonstrated their feasibility to explore in depth the molecular diversity of venoms. In the absence of genome sequence, transcriptomes represent also valuable searchable databases for proteomic projects.

RESULTS

The venom gland transcriptomes of 8 Costa Rican taxa from 5 genera (Crotalus, Bothrops, Atropoides, Cerrophidion, and Bothriechis) of pitvipers were investigated using high-throughput 454 pyrosequencing. 100,394 out of 330,010 masked reads produced significant hits in the available databases. 5.165,220 nucleotides (8.27%) were masked by RepeatMasker, the vast majority of which corresponding to class I (retroelements) and class II (DNA transposons) mobile elements. BLAST hits included 79,991 matches to entries of the taxonomic suborder Serpentes, of which 62,433 displayed similarity to documented venom proteins. Strong discrepancies between the transcriptome-computed and the proteome-gathered toxin compositions were obvious at first sight. Although the reasons underlaying this discrepancy are elusive, since no clear trend within or between species is apparent, the data indicate that individual mRNA species may be translationally controlled in a species-dependent manner. The minimum number of genes from each toxin family transcribed into the venom gland transcriptome of each species was calculated from multiple alignments of reads matched to a full-length reference sequence of each toxin family. Reads encoding ORF regions of Kazal-type inhibitor-like proteins were uniquely found in Bothriechis schlegelii and B. lateralis transcriptomes, suggesting a genus-specific recruitment event during the early-Middle Miocene. A transcriptome-based cladogram supports the large divergence between A. mexicanus and A. picadoi, and a closer kinship between A. mexicanus and C. godmani.

CONCLUSIONS

Our comparative next-generation sequencing (NGS) analysis reveals taxon-specific trends governing the formulation of the venom arsenal. Knowledge of the venom proteome provides hints on the translation efficiency of toxin-coding transcripts, contributing thereby to a more accurate interpretation of the transcriptome. The application of NGS to the analysis of snake venom transcriptomes, may represent the tool for opening the door to systems venomics.

Discovery of highly divergent repeat landscapes in snake genomes using high-throughput sequencing

We conducted a comprehensive assessment of genomic repeat content in two snake genomes, the venomous copperhead (Agkistrodon contortrix) and the Burmese python (Python molurus bivittatus). These two genomes are both relatively small (∼1.4 Gb) but have surprisingly extensive differences in the abundance and expansion histories of their repeat elements. In the python, the readily identifiable repeat element content is low (21%), similar to bird genomes, whereas that of the copperhead is higher (45%), similar to mammalian genomes. The copperhead's greater repeat content arises from the recent expansion of many different microsatellites and transposable element (TE) families, and the copperhead had 23-fold greater levels of TE-related transcripts than the python. This suggests the possibility that greater TE activity in the copperhead is ongoing. Expansion of CR1 LINEs in the copperhead genome has resulted in TE-mediated microsatellite expansion ("microsatellite seeding") at a scale several orders of magnitude greater than previously observed in vertebrates. Snakes also appear to be prone to horizontal transfer of TEs, particularly in the copperhead lineage. The reason that the copperhead has such a small genome in the face of so much recent expansion of repeat elements remains an open question, although selective pressure related to extreme metabolic performance is an obvious candidate. TE activity can affect gene regulation as well as rates of recombination and gene duplication, and it is therefore possible that TE activity played a role in the evolution of major adaptations in snakes; some evidence suggests this may include the evolution of venom repertoires.

Estimating genetic variability in non-model taxa: a general procedure for discriminating sequence errors from actual variation

Genetic variation is the driving force of evolution and as such is of central interest for biologists. However, inadequate discrimination of errors from true genetic variation could lead to incorrect estimates of gene copy number, population genetic parameters, phylogenetic relationships and the deposition of gene and protein sequences in databases that are not actually present in any organism. Misincorporation errors in multi-template PCR cloning methods, still commonly used for obtaining novel gene sequences in non-model species, are difficult to detect, as no previous information may be available about the number of expected copies of genes belonging to multi-gene families. However, studies employing these techniques rarely describe in any great detail how errors arising in the amplification process were detected and accounted for. Here, we estimated the rate of base misincorporation of a widely-used PCR-cloning method, using a single copy mitochondrial gene from a single individual to minimise variation in the template DNA, as 1.62×10(-3) errors per site, or 9.26×10(-5) per site per duplication. The distribution of errors among sequences closely matched that predicted by a binomial distribution function. The empirically estimated error rate was applied to data, obtained using the same methods, from the Phospholipase A(2) toxin family from the pitviper Ovophis monticola. The distribution of differences detected closely matched the expected distribution of errors and we conclude that, when undertaking gene discovery or assessment of genetic diversity using this error-prone method, it will be informative to empirically determine the rate of base misincorporation.