Molecular genetics of ecological diversification: duplication and rapid evolution of toxin genes of the venomous gastropod Conus

Potential Ancestral Conoidean Toxins in the Venom Cocktail of the Carnivorous Snail Raphitoma purpurea (Montagu, 1803) (Neogastropoda: Raphitomidae)

Venomous marine gastropods of the superfamily Conoidea possess a rich arsenal of toxins, including neuroactive toxins. Venom adaptations might have played a fundamental role in the radiation of conoideans; nevertheless, there is still no knowledge about the venom of the most diversified family of the group: Raphitomidae Bellardi, 1875. In this study, transcriptomes were produced from the carcase, salivary glands, and proximal and distal venom ducts of the northeastern Atlantic species Raphitoma purpurea (Montagu, 1803). Using a gut barcoding approach, we were also able to report, for the first time, molecular evidence of a vermivorous diet for the genus. Transcriptomic analyses revealed over a hundred putative venom components (PVC), including 69 neurotoxins. Twenty novel toxin families, including some with high levels of expansion, were discovered. No significant difference was observed between the distal and proximal venom duct secretions. Peptides related to cone snail toxins (Cerm06, Pgam02, and turritoxin) and other venom-related proteins (disulfide isomerase and elevenin) were retrieved from the salivary glands. These salivary venom components may constitute ancestral adaptations for venom production in conoideans. Although often neglected, salivary gland secretions are of extreme importance for understanding the evolutionary history of conoidean venom.

Prey Shifts Drive Venom Evolution in Cone Snails

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of antivenom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared with venom of nonpiscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Proteomic analysis of the venom of Conus flavidus from Red Sea reveals potential pharmacological applications

BACKGROUND

Venomous marine cone snails produce unique neurotoxins called conopeptides or conotoxins, which are valuable for research and drug discovery. Characterizing Conus venom is important, especially for poorly studied species, as these tiny and steady molecules have considerable potential as research tools for detecting new pharmacological applications. In this study, a worm-hunting cone snail, Conus flavidus inhabiting the Red Sea coast were collected, dissected and the venom gland extraction was subjected to proteomic analysis to define the venom composition, and confirm the functional structure of conopeptides.

RESULTS

Analysis of C. flavidus venom identified 117 peptide fragments and assorted them to conotoxin precursors and non-conotoxin proteins. In this procedure, 65 conotoxin precursors were classified and identified to 16 conotoxin precursors and hormone superfamilies. In the venom of C. flavidus, the four conotoxin superfamilies T, A, O2, and M were the most abundant peptides, accounting for 75.8% of the total conotoxin diversity. Additionally, 19 non-conotoxin proteins were specified in the venom, as well as several potentially biologically active peptides with putative applications.

CONCLUSION

Our research displayed that the structure of the C. flavidus-derived proteome is similar to other Conus species and includes toxins, ionic channel inhibitors, insulin-like peptides, and hyaluronidase. This study provides a foundation for discovering new conopeptides from C. flavidus venom for pharmaceutical use.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Mammalian neurotoxins, Blarina paralytic peptides, cause hyperpolarization of human T-type Ca channel hCav3.2 activation

Among the rare venomous mammals, the short-tailed shrew Blarina brevicauda has been suggested to produce potent neurotoxins in its saliva to effectively capture prey. Several kallikrein-like lethal proteases have been identified, but the active substances of B. brevicauda remained unclear. Here, we report Blarina paralytic peptides (BPPs) 1 and 2 isolated from its submaxillary glands. Synthetic BPP2 showed mealworm paralysis and a hyperpolarization shift (-11 mV) of a human T-type Ca2+ channel (hCav3.2) activation. The amino acid sequences of BPPs were similar to those of synenkephalins, which are precursors of brain opioid peptide hormones that are highly conserved among mammals. However, BPPs rather resembled centipede neurotoxic peptides SLPTXs in terms of disulfide bond connectivity and stereostructure. Our results suggested that the neurotoxin BPPs were the result of convergent evolution as homologs of nontoxic endogenous peptides that are widely conserved in mammals. This finding is of great interest from the viewpoint of the chemical evolution of vertebrate venoms.

Genomic basis of multiphase evolution driving divergent selection of zinc-finger homeodomain genes

Gene families divergently evolve and become adapted as different genes with specific structures and functions in living organisms. We performed comprehensive structural and functional analyses of Zinc-finger homeodomain genes (ZF-HDs), including Mini zinc-finger genes (MIFs) and Zinc-finger with homeodomain genes (ZHDs), displaying competitive functions each other. Intensive annotation updates for 90 plant genomes verified that most MIFs (MIF-Is) exhibited distinct motif compositions from ZHDs, although some MIFs (MIF-Zs) contained ZHD-specific motifs. Phylogenetic analyses suggested that MIF-Zs and ZHDs originated from the same ancestral gene, whereas MIF-Is emerged from a distinct progenitor. We used a gene-editing system to identify a novel function of MIF-Is in rice: regulating the surface material patterns in anthers and pollen through transcriptional regulation by interacting ZHDs. Kingdom-wide investigations determined that (i) ancestral MIFs diverged into MIF-Is and MIF-Zs in the last universal common ancestor, (ii) integration of HD into the C-terminal of MIF-Zs created ZHDs after emergence of green plants and (iii) MIF-Is and ZHDs subsequently expanded independently into specific plant lineages, with additional formation of MIF-Zs from ZHDs. Our comprehensive analysis provides genomic evidence for multiphase evolution driving divergent selection of ZF-HDs.

Whole Genome Duplication and Gene Evolution in the Hyperdiverse Venomous Gastropods

The diversity of venomous organisms and the toxins they produce have been increasingly investigated, but taxonomic bias remains important. Neogastropods, a group of marine predators representing almost 22% of the known gastropod diversity, evolved a wide range of feeding strategies, including the production of toxins to subdue their preys. However, whether the diversity of these compounds is at the origin of the hyperdiversification of the group and how genome evolution may correlate with both the compounds and species diversities remain understudied. Among the available gastropods genomes, only eight, with uneven quality assemblies, belong to neogastropods. Here, we generated chromosome-level assemblies of two species belonging to the Tonnoidea and Muricoidea superfamilies (Monoplex corrugatus and Stramonita haemastoma). The two obtained high-quality genomes had 3 and 2.2 Gb, respectively, and 92-89% of the total assembly conformed 35 pseudochromosomes in each species. Through the analysis of syntenic blocks, Hox gene cluster duplication, and synonymous substitutions distribution pattern, we inferred the occurrence of a whole genome duplication event in both genomes. As these species are known to release venom, toxins were annotated in both genomes, but few of them were found in homologous chromosomes. A comparison of the expression of ohnolog genes (using transcriptomes from osphradium and salivary glands in S. haemastoma), where both copies were differentially expressed, showed that most of them had similar expression profiles. The high quality of these genomes makes them valuable reference in their respective taxa, facilitating the identification of genome-level processes at the origin of their evolutionary success.

Coordinated adaptations define the ontogenetic shift from worm- to fish-hunting in a venomous cone snail

Marine cone snails have attracted researchers from all disciplines but early life stages have received limited attention due to difficulties accessing or rearing juvenile specimens. Here, we document the culture of Conus magus from eggs through metamorphosis to reveal dramatic shifts in predatory feeding behaviour between post-metamorphic juveniles and adult specimens. Adult C. magus capture fish using a set of paralytic venom peptides combined with a hooked radular tooth used to tether envenomed fish. In contrast, early juveniles feed exclusively on polychaete worms using a unique "sting-and-stalk" foraging behaviour facilitated by short, unbarbed radular teeth and a distinct venom repertoire that induces hypoactivity in prey. Our results demonstrate how coordinated morphological, behavioural and molecular changes facilitate the shift from worm- to fish-hunting in C. magus, and showcase juvenile cone snails as a rich and unexplored source of novel venom peptides for ecological, evolutionary and biodiscovery studies.

The deep-rooted origin of disulfide-rich spider venom toxins

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name 'Adi Shakti', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.

Comprehensive Analysis of Copy Number Variations on Glycoside Hydrolase 45 Genes among Different Bursaphelenchus xylophilus Strains

Bursaphelenchus xylophilus is considered the most dangerous quarantine pest in China. It causes enormous economic and ecological losses in many countries from Asia and Europe. The glycoside hydrolase 45 gene family has been demonstrated in early studies to contribute to the cell wall degradation ability of B. xylophilus during its infection. However, the copy number variation (CNV) of the GH45 gene and its association with B. xylophilus pathogenicity were not fully elucidated. In this study, we found that the GH45 gene with two copies is the most predominant type among 259 B. xylophilus strains collected from China and Japan. Additionally, 18 strains are identified as GH45 genes with a single copy, and only two strains are verified to have three copies. Subsequent expression analysis and inoculation test suggest that the copy numbers of the GH45 gene are correlated with gene expression as well as the B. xylophilus pathogenicity. B. xylophilus strains with more copies of the GH45 gene usually exhibit more abundant expression and cause more severe wilt symptoms on pine trees. The aforementioned results indicated the potential regulatory effects of CNV in B. xylophilus and provided novel information to better understand the molecular pathogenesis of this devastating pest.

Venomics survey of six myrmicine ants provides insights into the molecular and structural diversity of their peptide toxins

Among ants, Myrmicinae represents the most speciose subfamily. The venom composition previously described for these social insects is extremely variable, with alkaloids predominant in some genera while, conversely, proteomics studies have revealed that some myrmicine ant venoms are peptide-rich. Using integrated transcriptomic and proteomic approaches, we characterized the venom peptidomes of six ants belonging to the different tribes of Myrmicinae. We identified a total of 79 myrmicitoxins precursors which can be classified into 38 peptide families according to their mature sequences. Myrmicine ant venom peptidomes showed heterogeneous compositions, with linear and disulfide-bonded monomers as well as dimeric toxins. Several peptide families were exclusive to a single venom whereas some were retrieved in multiple species. A hierarchical clustering analysis of precursor signal sequences led us to divide the myrmicitoxins precursors into eight families, including some that have already been described in other aculeate hymenoptera such as secapin-like peptides and voltage-gated sodium channel (Na_V) toxins. Evolutionary and structural analyses of two representatives of these families highlighted variation and conserved patterns that might be crucial to explain myrmicine venom peptide functional adaptations to biological targets.

Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA₂s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

Proteomic analysis of Red Sea Conus taeniatus venom reveals potential biological applications

Background:

Diverse and unique bioactive neurotoxins known as conopeptides or conotoxins are produced by venomous marine cone snails. Currently, these small and stable molecules are of great importance as research tools and platforms for discovering new drugs and therapeutics. Therefore, the characterization of Conus venom is of great significance, especially for poorly studied species.

Methods:

In this study, we used bioanalytical techniques to determine the venom profile and emphasize the functional composition of conopeptides in Conus taeniatus, a neglected worm-hunting cone snail.

Results:

The proteomic analysis revealed that 84.0% of the venom proteins were between 500 and 4,000 Da, and 16.0% were > 4,000 Da. In C. taeniatus venom, 234 peptide fragments were identified and classified as conotoxin precursors or non-conotoxin proteins. In this process, 153 conotoxin precursors were identified and matched to 23 conotoxin precursors and hormone superfamilies. Notably, the four conotoxin superfamilies T (22.87%), O1 (17.65%), M (13.1%) and O2 (9.8%) were the most abundant peptides in C. taeniatus venom, accounting for 63.40% of the total conotoxin diversity. On the other hand, 48 non-conotoxin proteins were identified in the venom of C. taeniatus. Moreover, several possibly biologically active peptide matches were identified, and putative applications of the peptides were assigned.

Conclusion:

Our study showed that the composition of the C. taeniatus-derived proteome is comparable to that of other Conus species and contains an effective mix of toxins, ionic channel inhibitors and antimicrobials. Additionally, it provides a guidepost for identifying novel conopeptides from the venom of C. taeniatus and discovering conopeptides of potential pharmaceutical importance.

Venom duct origins of prey capture and defensive conotoxins in piscivorous Conus striatus

The venom duct origins of predatory and defensive venoms has not been studied for hook-and-line fish hunting cone snails despite the pharmacological importance of their venoms. To better understand the biochemistry and evolution of injected predatory and defensive venoms, we compared distal, central and proximal venom duct sections across three specimens of C. striatus (Pionoconus) using proteomic and transcriptomic approaches. A total of 370 conotoxin precursors were identified from the whole venom duct transcriptome. Milked defensive venom was enriched with a potent cocktail of proximally expressed inhibitory α-, ω- and μ-conotoxins compared to milked predatory venom. In contrast, excitatory κA-conotoxins dominated both the predatory and defensive venoms despite their distal expression, suggesting this class of conotoxin can be selectively expressed from the same duct segment in response to either a predatory or defensive stimuli. Given the high abundance of κA-conotoxins in the Pionoconus clade, we hypothesise that the κA-conotoxins have evolved through adaptive evolution following their repurposing from ancestral inhibitory A superfamily conotoxins to facilitate the dietary shift to fish hunting and species radiation in this clade.

Copy Number Variations of Glycoside Hydrolase 45 Genes in Bursaphelenchus xylophilus and Their Impact on the Pathogenesis of Pine Wilt Disease

The pine wood nematode Bursaphelenchus xylophilus parasitizes millions of pine trees worldwide each year, causing severe wilt and the death of host trees. Glycoside hydrolase 45 genes of B. xylophilus are reported to have been acquired by horizontal gene transfer from fungi and are responsible for cell wall degradation during nematode infection. Previous studies ignored the possibility of copy number variations of such genes. In this study, we determined that two of the glycoside hydrolase 45 genes evolved to maintain multiple copies with distinct expression levels, enabling the nematode to infect a variety of pine hosts. Additionally, tandem repeat variations within coding regions were also detected between different copies of glycoside hydrolase 45 genes that could result in changes in protein sequences and serve as an effective biological marker to detect copy number variations among different B. xylophilus populations. Consequently, we were able to further identify the copy number variations of glycoside hydrolase 45 genes among B. xylophilus strains with different virulence. Our results provide new insights into the pathogenicity of B. xylophilus, provide a practical marker to genotype copy number variations and may aid in population classification.

How the Toxin got its Toxicity

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules—“toxins”—share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as “exochemicals.” There has been debate within the field of toxinology concerning the evolutionary pathways leading to the “recruitment” of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel “environment,” i.e., a range of novel partners for molecular interaction. After stressing the distinction between “activity” and “function,” we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of “neofunctionalization”—simply the origins of a novel molecular function by any process—and “recruitment”—the “weaponization” of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of “birth-and-death.” Thus, recruitment leads to a form of “downwards causation,” in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Curses or Cures: A Review of the Numerous Benefits Versus the Biosecurity Concerns of Conotoxin Research

Conotoxins form a diverse group of peptide toxins found in the venom of predatory marine cone snails. Decades of conotoxin research have provided numerous measurable scientific and societal benefits. These include their use as a drug, diagnostic agent, drug leads, and research tools in neuroscience, pharmacology, biochemistry, structural biology, and molecular evolution. Human envenomations by cone snails are rare but can be fatal. Death by envenomation is likely caused by a small set of toxins that induce muscle paralysis of the diaphragm, resulting in respiratory arrest. The potency of these toxins led to concerns regarding the potential development and use of conotoxins as biological weapons. To address this, various regulatory measures have been introduced that limit the use and access of conotoxins within the research community. Some of these regulations apply to all of the ≈200,000 conotoxins predicted to exist in nature of which less than 0.05% are estimated to have any significant toxicity in humans. In this review we provide an overview of the many benefits of conotoxin research, and contrast these to the perceived biosecurity concerns of conotoxins and research thereof.

Presence-absence polymorphisms of single-copy genes in the stony coral Acropora digitifera

Background

Despite the importance of characterizing genetic variation among coral individuals for understanding phenotypic variation, the correlation between coral genomic diversity and phenotypic expression is still poorly understood.

Results

In this study, we detected a high frequency of genes showing presence–absence polymorphisms (PAPs) for single-copy genes in Acropora digitifera. Among 10,455 single-copy genes, 516 (5%) exhibited PAPs, including 32 transposable element (TE)-related genes. Five hundred sixteen genes exhibited a homozygous absence in one (102) or more than one (414) individuals (n = 33), indicating that most of the absent alleles were not rare variants. Among genes showing PAPs (PAP genes), roughly half were expressed in adults and/or larvae, and the PAP status was associated with differential expression among individuals. Although 85% of PAP genes were uncharacterized or had ambiguous annotations, 70% of these genes were specifically distributed in cnidarian lineages in eumetazoa, suggesting that these genes have functional roles related to traits related to cnidarians or the family Acroporidae or the genus Acropora. Indeed, four of these genes encoded toxins that are usually components of venom in cnidarian-specific cnidocytes. At least 17% of A. digitifera PAP genes were also PAPs in A. tenuis, the basal lineage in the genus Acropora, indicating that PAPs were shared among species in Acropora.

Conclusions

Expression differences caused by a high frequency of PAP genes may be a novel genomic feature in the genus Acropora; these findings will contribute to improve our understanding of correlation between genetic and phenotypic variation in corals.

Venomic Interrogation Reveals the Complexity of Conus striolatus Venom

Given the complexity of cone snail venoms, high throughput venomics approaches are required to fully investigate venom composition, envenomation strategies, and evolutionary trajectories. This study describes 158 conotoxins in the venom transcriptome of the little studied C. striolatus from the fish hunting clade Pionoconus. Despite similar gene superfamily distributions along the venom duct, only 18 common transcripts were identified between distal, central, and proximal venom duct transcriptomes. Proteomic analysis of the injected predatory venom collected from the same individual revealed an ~18-fold enhanced complexity at the proteomic level, consistent with complex post-translational modifications and variable venom peptide processing occurring in the venom duct. Overall, C. striolatus venom was dominated by M, O1, O2, and A gene superfamily conotoxins and conkunitzins, which are potential modulators of sodium, calcium, and potassium channels. Conkunitzins and gene superfamily A peptides dominated the proximal over the distal duct, the M and O1 gene superfamily peptides were distributed along the full length of the duct, while the O2 gene superfamily peptides dominated the distal duct. Interestingly, the predatory injected venom of C. striolatus was dominated by peptides from gene superfamilies M, O1, O2, A, and conkunitzins, suggesting the predatory venom of C. striolatus may arise at multiple sites along the venom duct.

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

Conotoxin Diversity in the Venom Gland Transcriptome of the Magician's Cone, Pionoconus magus

The transcriptomes of the venom glands of two individuals of the magician’s cone, Pionoconus magus, from Okinawa (Japan) were sequenced, assembled, and annotated. In addition, RNA-seq raw reads available at the SRA database from one additional specimen of P. magus from the Philippines were also assembled and annotated. The total numbers of identified conotoxin precursors and hormones per specimen were 118, 112, and 93. The three individuals shared only five identical sequences whereas the two specimens from Okinawa had 30 sequences in common. The total number of distinct conotoxin precursors and hormones for P. magus was 275, and were assigned to 53 conotoxin precursor and hormone superfamilies, two of which were new based on their divergent signal region. The superfamilies that had the highest number of precursors were M (42), O1 (34), T (27), A (18), O2 (17), and F (13), accounting for 55% of the total diversity. The D superfamily, previously thought to be exclusive of vermivorous cones was found in P. magus and contained a highly divergent mature region. Similarly, the A superfamily alpha 4/3 was found in P. magus despite the fact that it was previously postulated to be almost exclusive of the genus Rhombiconus. Differential expression analyses of P. magus compared to Chelyconus ermineus, the only fish-hunting cone from the Atlantic Ocean revealed that M and A2 superfamilies appeared to be more expressed in the former whereas the O2 superfamily was more expressed in the latter.

Lack of Signal for the Impact of Conotoxin Gene Diversity on Speciation Rates in Cone Snails

Understanding why some groups of organisms are more diverse than others is a central goal in macroevolution. Evolvability, or the intrinsic capacity of lineages for evolutionary change, is thought to influence disparities in species diversity across taxa. Over macroevolutionary time scales, clades that exhibit high evolvability are expected to have higher speciation rates. Cone snails (family: Conidae, $>$900 spp.) provide a unique opportunity to test this prediction because their toxin genes can be used to characterize differences in evolvability between clades. Cone snails are carnivorous, use prey-specific venom (conotoxins) to capture prey, and the genes that encode venom are known and diversify through gene duplication. Theory predicts that higher gene diversity confers a greater potential to generate novel phenotypes for specialization and adaptation. Therefore, if conotoxin gene diversity gives rise to varying levels of evolvability, conotoxin gene diversity should be coupled with macroevolutionary speciation rates. We applied exon capture techniques to recover phylogenetic markers and conotoxin loci across 314 species, the largest venom discovery effort in a single study. We paired a reconstructed timetree using 12 fossil calibrations with species-specific estimates of conotoxin gene diversity and used trait-dependent diversification methods to test the impact of evolvability on diversification patterns. Surprisingly, we did not detect any signal for the relationship between conotoxin gene diversity and speciation rates, suggesting that venom evolution may not be the rate-limiting factor controlling diversification dynamics in Conidae. Comparative analyses showed some signal for the impact of diet and larval dispersal strategy on diversification patterns, though detection of a signal depended on the dataset and the method. If our results remain true with increased taxonomic sampling in future studies, they suggest that the rapid evolution of conid venom may cause other factors to become more critical to diversification, such as ecological opportunity or traits that promote isolation among lineages.

TOXIFY: a deep learning approach to classify animal venom proteins

In the era of Next-Generation Sequencing and shotgun proteomics, the sequences of animal toxigenic proteins are being generated at rates exceeding the pace of traditional means for empirical toxicity verification. To facilitate the automation of toxin identification from protein sequences, we trained Recurrent Neural Networks with Gated Recurrent Units on publicly available datasets. The resulting models are available via the novel software package TOXIFY, allowing users to infer the probability of a given protein sequence being a venom protein. TOXIFY is more than 20X faster and uses over an order of magnitude less memory than previously published methods. Additionally, TOXIFY is more accurate, precise, and sensitive at classifying venom proteins.

Effects of Predator-Prey Interactions on Predator Traits: Differentiation of Diets and Venoms of a Marine Snail

Species interactions are fundamental ecological forces that can have significant impacts on the evolutionary trajectories of species. Nonetheless, the contribution of predator-prey interactions to genetic and phenotypic divergence remains largely unknown. Predatory marine snails of the family Conidae exhibit specializations for different prey items and intraspecific variation in prey utilization patterns at geographic scales. Because cone snails utilize venom to capture prey and venom peptides are direct gene products, it is feasible to examine the evolution of genes associated with changes in resource utilization. Here, we compared feeding ecologies and venom duct transcriptomes of individuals from three populations of Conus miliaris, a species that exhibits geographic variation in prey utilization and dietary breadth, in order to determine the extent to which dietary differences are correlated with differences in venom composition, and if expanded niche breadth is associated with increased variation in venom composition. While populations showed little to no overlap in resource utilization, taxonomic richness of prey was greatest at Easter Island. Changes in dietary breadth were associated with differences in expression patterns and increased genetic differentiation of toxin-related genes. The Easter Island population also exhibited greater diversity of toxin-related transcripts, but did not show increased variance in expression of these transcripts. These results imply that differences in dietary breadth contribute more to the structural and regulatory differentiation of venoms than differences in diet.

Thinking twice about the evolution of photosynthesis

Sam Granick opened his seminal 1957 paper titled 'Speculations on the origins and evolution of photosynthesis' with the assertion that there is a constant urge in human beings to seek beginnings (I concur). This urge has led to an incessant stream of speculative ideas and debates on the evolution of photosynthesis that started in the first half of the twentieth century and shows no signs of abating. Some of these speculative ideas have become commonplace, are taken as fact, but find little support. Here, I review and scrutinize three widely accepted ideas that underpin the current study of the evolution of photosynthesis: first, that the photochemical reaction centres used in anoxygenic photosynthesis are more primitive than those in oxygenic photosynthesis; second, that the probability of acquiring photosynthesis via horizontal gene transfer is greater than the probability of losing photosynthesis; and third, and most important, that the origin of anoxygenic photosynthesis pre-dates the origin of oxygenic photosynthesis. I shall attempt to demonstrate that these three ideas are often grounded in incorrect assumptions built on more assumptions with no experimental or observational support. I hope that this brief review will not only serve as a cautionary tale but also that it will open new avenues of research aimed at disentangling the complex evolution of photosynthesis and its impact on the early history of life and the planet.

Early Archean origin of Photosystem II

Photosystem II is a photochemical reaction center that catalyzes the light-driven oxidation of water to molecular oxygen. Water oxidation is the distinctive photochemical reaction that permitted the evolution of oxygenic photosynthesis and the eventual rise of eukaryotes. At what point during the history of life an ancestral photosystem evolved the capacity to oxidize water still remains unknown. Here, we study the evolution of the core reaction center proteins of Photosystem II using sequence and structural comparisons in combination with Bayesian relaxed molecular clocks. Our results indicate that a homodimeric photosystem with sufficient oxidizing power to split water had already appeared in the early Archean about a billion years before the most recent common ancestor of all described Cyanobacteria capable of oxygenic photosynthesis, and well before the diversification of some of the known groups of anoxygenic photosynthetic bacteria. Based on a structural and functional rationale, we hypothesize that this early Archean photosystem was capable of water oxidation to oxygen and had already evolved protection mechanisms against the formation of reactive oxygen species. This would place primordial forms of oxygenic photosynthesis at a very early stage in the evolutionary history of life.

The Evolutionary History and Functional Divergence of Trehalase (treh) Genes in Insects

Trehalases (treh) have been found in different organisms, such as bacteria, fungi, yeast, nematodes, insects, vertebrates, and plants. Their biochemical properties are extremely variable and not yet fully understood. Gene expression patterns have shown differences among insect species suggesting a potential functional diversification of trehalase enzymes during their evolution. A second gene family encoding for enzymes with hypothetical trehalase activity has been repeatedly annotated in insect genome as acid trehalases/acid trehalase-like (ath), but its functional role is still not clear. The currently available large amount of genomic data from many insect species may enable a better understanding of the evolutionary history, phylogenetic relationships and possible roles of trehalase encoding genes in this taxon. The aim of the present study is to infer the evolutionary history of trehalases and acid trehalase genes in insects and analyze the trehalase functional divergence during their evolution, combining phylogenetic and genomic synteny/colinearity analyses.

Phylogenetic Comparative Methods can Provide Important Insights into the Evolution of Toxic Weaponry

The literature on chemical weaponry of organisms is vast and provides a rich understanding of the composition and mechanisms of the toxins and other components involved. However, an ecological or evolutionary perspective has often been lacking and is largely limited to (1) molecular evolutionary studies of particular toxins (lacking an ecological view); (2) comparisons across different species that ignore phylogenetic relatedness (lacking an evolutionary view); or (3) descriptive studies of venom composition and toxicology that contain post hoc and untested ecological or evolutionary interpretations (a common event but essentially uninformative speculation). Conveniently, comparative biologists have prolifically been developing and using a wide range of phylogenetic comparative methods that allow us to explicitly address many ecological and evolutionary questions relating to venoms and poisons. Nevertheless, these analytical tools and approaches are rarely used and poorly known by biological toxinologists and toxicologists. In this review I aim to (1) introduce phylogenetic comparative methods to the latter audience; (2) highlight the range of questions that can be addressed using them; and (3) encourage biological toxinologists and toxicologists to either seek out adequate training in comparative biology or seek collaboration with comparative biologists to reap the fruits of a powerful interdisciplinary approach to the field.

Whole-Genome Sequencing of Chinese Yellow Catfish Provides a Valuable Genetic Resource for High-Throughput Identification of Toxin Genes

Naturally derived toxins from animals are good raw materials for drug development. As a representative venomous teleost, Chinese yellow catfish (Pelteobagrus fulvidraco) can provide valuable resources for studies on toxin genes. Its venom glands are located in the pectoral and dorsal fins. Although with such interesting biologic traits and great value in economy, Chinese yellow catfish is still lacking a sequenced genome. Here, we report a high-quality genome assembly of Chinese yellow catfish using a combination of next-generation Illumina and third-generation PacBio sequencing platforms. The final assembly reached 714 Mb, with a contig N50 of 970 kb and a scaffold N50 of 3.65 Mb, respectively. We also annotated 21,562 protein-coding genes, in which 97.59% were assigned at least one functional annotation. Based on the genome sequence, we analyzed toxin genes in Chinese yellow catfish. Finally, we identified 207 toxin genes and classified them into three major groups. Interestingly, we also expanded a previously reported sex-related region (to ≈6 Mb) in the achieved genome assembly, and localized two important toxin genes within this region. In summary, we assembled a high-quality genome of Chinese yellow catfish and performed high-throughput identification of toxin genes from a genomic view. Therefore, the limited number of toxin sequences in public databases will be remarkably improved once we integrate multi-omics data from more and more sequenced species.

Conotoxin Diversity in Chelyconus ermineus (Born, 1778) and the Convergent Origin of Piscivory in the Atlantic and Indo-Pacific Cones

The transcriptome of the venom duct of the Atlantic piscivorous cone species Chelyconus ermineus (Born, 1778) was determined. The venom repertoire of this species includes at least 378 conotoxin precursors, which could be ascribed to 33 known and 22 new (unassigned) protein superfamilies, respectively. Most abundant superfamilies were T, W, O1, M, O2, and Z, accounting for 57% of all detected diversity. A total of three individuals were sequenced showing considerable intraspecific variation: each individual had many exclusive conotoxin precursors, and only 20% of all inferred mature peptides were common to all individuals. Three different regions (distal, medium, and proximal with respect to the venom bulb) of the venom duct were analyzed independently. Diversity (in terms of number of distinct members) of conotoxin precursor superfamilies increased toward the distal region whereas transcripts detected toward the proximal region showed higher expression levels. Only the superfamilies A and I3 showed statistically significant differential expression across regions of the venom duct. Sequences belonging to the alpha (motor cabal) and kappa (lightning-strike cabal) subfamilies of the superfamily A were mainly detected in the proximal region of the venom duct. The mature peptides of the alpha subfamily had the α4/4 cysteine spacing pattern, which has been shown to selectively target muscle nicotinic-acetylcholine receptors, ultimately producing paralysis. This function is performed by mature peptides having a α3/5 cysteine spacing pattern in piscivorous cone species from the Indo-Pacific region, thereby supporting a convergent evolution of piscivory in cones.

Transcriptional fates of human-specific segmental duplications in brain

Despite the importance of duplicate genes for evolutionary adaptation, accurate gene annotation is often incomplete, incorrect, or lacking in regions of segmental duplication. We developed an approach combining long-read sequencing and hybridization capture to yield full-length transcript information and confidently distinguish between nearly identical genes/paralogs. We used biotinylated probes to enrich for full-length cDNA from duplicated regions, which were then amplified, size-fractionated, and sequenced using single-molecule, long-read sequencing technology, permitting us to distinguish between highly identical genes by virtue of multiple paralogous sequence variants. We examined 19 gene families as expressed in developing and adult human brain, selected for their high sequence identity (average >99%) and overlap with human-specific segmental duplications (SDs). We characterized the transcriptional differences between related paralogs to better understand the birth-death process of duplicate genes and particularly how the process leads to gene innovation. In 48% of the cases, we find that the expressed duplicates have changed substantially from their ancestral models due to novel sites of transcription initiation, splicing, and polyadenylation, as well as fusion transcripts that connect duplication-derived exons with neighboring genes. We detect unannotated open reading frames in genes currently annotated as pseudogenes, while relegating other duplicates to nonfunctional status. Our method significantly improves gene annotation, specifically defining full-length transcripts, isoforms, and open reading frames for new genes in highly identical SDs. The approach will be more broadly applicable to genes in structurally complex regions of other genomes where the duplication process creates novel genes important for adaptive traits.

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.

Genetic Diversity of Seven Cattle Breeds Inferred Using Copy Number Variations

Copy number variations (CNVs) comprise deletions, duplications, and insertions found within the genome larger than 50 bp in size. CNVs are thought to be primary role-players in breed formation and adaptation. South Africa boasts a diverse ecology with harsh environmental conditions and a broad spectrum of parasites and diseases that pose challenges to livestock production. This has led to the development of composite cattle breeds which combine the hardiness of Sanga breeds and the production potential of the Taurine breeds. The prevalence of CNVs within these respective breeds of cattle and the prevalence of CNV regions (CNVRs) in their diversity, adaptation and production is however not understood. This study therefore aimed to ascertain the prevalence, diversity, and correlations of CNVRs within cattle breeds used in South Africa. Illumina Bovine SNP50 data and PennCNV were utilized to identify CNVRs within the genome of 287 animals from seven cattle breeds representing Sanga, Taurine, Composite, and cross breeds. Three hundred and fifty six CNVRs of between 36 kb to 4.1 Mb in size were identified. The null hypothesis that one CNVR loci is independent of another was tested using the GENEPOP software. One hunded and two and seven of the CNVRs in the Taurine and Sanga/Composite cattle breeds demonstrated a significant (p ≤ 0.05) association. PANTHER overrepresentation analyses of correlated CNVRs demonstrated significant enrichment of a number of biological processes, molecular functions, cellular components, and protein classes. CNVR genetic variation between and within breed group was measured using phiPT which allows intra-individual variation to be suppressed and hence proved suitable for measuring binary CNVR presence/absence data. Estimate PhiPT within and between breed variance was 2.722 and 0.518 respectively. Pairwise population PhiPT values corresponded with breed type, with Taurine Holstein and Angus breeds demonstrating no between breed CNVR variation. Phylogenetic trees were drawn. CNVRs primarily clustered animals of the same breed type together. This study successfully identified, characterized, and analyzed 356 CNVRs within seven cattle breeds. CNVR correlations were evident, with many more correlations being present among the exotic Taurine breeds. CNVR genetic diversity of Sanga, Taurine and Composite breeds was ascertained with breed types exposed to similar selection pressures demonstrating analogous incidences of CNVRs.

Targeted Sequencing of Venom Genes from Cone Snail Genomes Improves Understanding of Conotoxin Molecular Evolution

To expand our capacity to discover venom sequences from the genomes of venomous organisms, we applied targeted sequencing techniques to selectively recover venom gene superfamilies and nontoxin loci from the genomes of 32 cone snail species (family, Conidae), a diverse group of marine gastropods that capture their prey using a cocktail of neurotoxic peptides (conotoxins). We were able to successfully recover conotoxin gene superfamilies across all species with high confidence (> 100× coverage) and used these data to provide new insights into conotoxin evolution. First, we found that conotoxin gene superfamilies are composed of one to six exons and are typically short in length (mean = ∼85 bp). Second, we expanded our understanding of the following genetic features of conotoxin evolution: 1) positive selection, where exons coding the mature toxin region were often three times more divergent than their adjacent noncoding regions, 2) expression regulation, with comparisons to transcriptome data showing that cone snails only express a fraction of the genes available in their genome (24-63%), and 3) extensive gene turnover, where Conidae species varied from 120 to 859 conotoxin gene copies. Finally, using comparative phylogenetic methods, we found that while diet specificity did not predict patterns of conotoxin evolution, dietary breadth was positively correlated with total conotoxin gene diversity. Overall, the targeted sequencing technique demonstrated here has the potential to radically increase the pace at which venom gene families are sequenced and studied, reshaping our ability to understand the impact of genetic changes on ecologically relevant phenotypes and subsequent diversification.

Animal toxins for channelopathy treatment

Ion channels are transmembrane proteins that allow passive flow of ions inside and/or outside of cells or cell organelles. Except mutations lead to nonfunctional protein production or abolished receptor entrance on the membrane surface an altered channel may have two principal conditions that can be corrected. The channel may conduct fewer ions through (loss-of-function mutations) or too many ions (gain-of-function mutations) compared to a normal channel. Toxins from animal venoms are specialised molecules that are generally oriented toward interactions with ion channels. This is a result of long coevolution between predators and their prey. On the molecular level, toxins activate or inhibit ion channels, so they are ideal molecules for restoring conductance in mutated channels. Another aspect of this long coevolution is that a broad variety of toxins have been fine tuned to recognize the channels of different species, keeping many amino acids substitution among sequences. Many peptide ligands with high selectivity to specific receptor subtypes have been isolated from animal venoms, some of which are absolutely non-toxic to humans and mammalians. It is expected that molecules that are selective to each known receptor can be found in animal venoms, but the pool of toxins currently does not override all receptors described as being involved in channelopathies. Modern investigating methods have enhanced the search process for selective ligands. One prominent method is a site-directed mutagenesis of existing toxins to change the selectivity or/and affinity to the selected receptor, which has shown positive results. This article is part of the Special Issue entitled 'Channelopathies.'

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

Linking neuroethology to the chemical biology of natural products: interactions between cone snails and their fish prey, a case study

From a biological perspective, a natural product can be defined as a compound evolved by an organism for chemical interactions with another organism including prey, predator, competitor, pathogen, symbiont or host. Natural products hold tremendous potential as drug leads and have been extensively studied by chemists and biochemists in the pharmaceutical industry. However, the biological purpose for which a natural product evolved is rarely addressed. By focusing on a well-studied group of natural products-venom components from predatory marine cone snails-this review provides a rationale for why a better understanding of the evolution, biology and biochemistry of natural products will facilitate both neuroscience and the potential for drug leads. The larger goal is to establish a new sub-discipline in the broader field of neuroethology that we refer to as "Chemical Neuroethology", linking the substantial work carried out by chemists on natural products with accelerating advances in neuroethology.

Selection To Increase Expression, Not Sequence Diversity, Precedes Gene Family Origin and Expansion in Rattlesnake Venom

Gene duplication is the primary mechanism leading to new genes and phenotypic novelty, but the proximate evolutionary processes underlying gene family origin, maintenance, and expansion are poorly understood. Although sub- and neofunctionalization provide clear long-term advantages, selection does not act with foresight, and unless a redundant gene copy provides an immediate fitness advantage, the copy will most likely be lost. Many models for the evolution of genes immediately following duplication have been proposed, but the robustness and applicability of these models is unclear because of the lack of data at the population level. We used qPCR, protein expression data, genome sequencing, and hybrid enrichment to test three competing models that differ in whether selection favoring the spread of duplicates acts primarily on expression level or sequence diversity for specific toxin-encoding loci in the eastern diamondback rattlesnake (Crotalus adamanteus). We sampled 178 individuals and identified significant inter- and intrapopulation variation in copy number, demonstrated that copy number was significantly and positively correlated with protein expression, and found little to no sequence variation across paralogs in all populations. Collectively, these results demonstrate that selection for increased expression, not sequence diversity, was the proximate evolutionary process underlying gene family origin and expansion, providing data needed to resolve the debate over which evolutionary processes govern the fates of gene copies immediately following duplication.

Small Packages, Big Returns: Uncovering the Venom Diversity of Small Invertebrate Conoidean Snails

Venomous organisms used in research were historically chosen based on size and availability. This opportunity-driven strategy created a species bias in which snakes, scorpions, and spiders became the primary subjects of venom research. Increasing technological advancements have enabled interdisciplinary studies using genomics, transcriptomics, and proteomics to expand venom investigation to animals that produce small amounts of venom or lack traditional venom producing organs. One group of non-traditional venomous organisms that have benefitted from the rise of -omic technologies is the Conoideans. The Conoidean superfamily of venomous marine snails includes, the Terebridae, Turridae (s.l), and Conidae. Conoidea venom is used for both predation and defense, and therefore under strong selection pressures. The need for conoidean venom peptides to be potent and specific to their molecular targets has made them important tools for investigating cellular physiology and bioactive compounds that are beneficial to improving human health. A convincing case for the potential of Conoidean venom is made with the first commercially available conoidean venom peptide drug Ziconotide (Prialt®), an analgesic derived from Conus magus venom that is used to treat chronic pain in HIV and cancer patients. Investigation of conoidean venom using -omics technology provides significant insights into predator-driven diversification in biodiversity and identifies novel compounds for manipulating cellular communication, especially as it pertains to disease and disorders.

Venom Insulins of Cone Snails Diversify Rapidly and Track Prey Taxa

A specialized insulin was recently found in the venom of a fish-hunting cone snail, Conus geographus Here we show that many worm-hunting and snail-hunting cones also express venom insulins, and that this novel gene family has diversified explosively. Cone snails express a highly conserved insulin in their nerve ring; presumably this conventional signaling insulin is finely tuned to the Conus insulin receptor, which also evolves very slowly. By contrast, the venom insulins diverge rapidly, apparently in response to biotic interactions with prey and also possibly the cones' own predators and competitors. Thus, the inwardly directed signaling insulins appear to experience predominantly purifying sele\ction to target an internal receptor that seldom changes, while the outwardly directed venom insulins frequently experience directional selection to target heterospecific insulin receptors in a changing mix of prey, predators and competitors. Prey insulin receptors may often be constrained in ways that prevent their evolutionary escape from targeted venom insulins, if amino-acid substitutions that result in escape also degrade the receptor's signaling functions.

Intragenome Diversity of Gene Families Encoding Toxin-like Proteins in Venomous Animals

The evolution of venoms is the story of how toxins arise and of the processes that generate and maintain their diversity. For animal venoms these processes include recruitment for expression in the venom gland, neofunctionalization, paralogous expansions, and functional divergence. The systematic study of these processes requires the reliable identification of the venom components involved in antagonistic interactions. High-throughput sequencing has the potential of uncovering the entire set of toxins in a given organism, yet the existence of non-venom toxin paralogs and the misleading effects of partial census of the molecular diversity of toxins make necessary to collect complementary evidence to distinguish true toxins from their non-venom paralogs. Here, we analyzed the whole genomes of two scorpions, one spider and one snake, aiming at the identification of the full repertoires of genes encoding toxin-like proteins. We classified the entire set of protein-coding genes into paralogous groups and monotypic genes, identified genes encoding toxin-like proteins based on known toxin families, and quantified their expression in both venom-glands and pooled tissues. Our results confirm that genes encoding toxin-like proteins are part of multigene families, and that these families arise by recruitment events from non-toxin genes followed by limited expansions of the toxin-like protein coding genes. We also show that failing to account for sequence similarity with non-toxin proteins has a considerable misleading effect that can be greatly reduced by comparative transcriptomics. Our study overall contributes to the understanding of the evolutionary dynamics of proteins involved in antagonistic interactions.

Is Mutation Random or Targeted?: No Evidence for Hypermutability in Snail Toxin Genes

Ever since Luria and Delbruck, the notion that mutation is random with respect to fitness has been foundational to modern biology. However, various studies have claimed striking exceptions to this rule. One influential case involves toxin-encoding genes in snails of the genus Conus, termed conotoxins, a large gene family that undergoes rapid diversification of their protein-coding sequences by positive selection. Previous reconstructions of the sequence evolution of conotoxin genes claimed striking patterns: (1) elevated synonymous change, interpreted as being due to targeted "hypermutation" in this region; (2) elevated transversion-to-transition ratios, interpreted as reflective of the particular mechanism of hypermutation; and (3) much lower rates of synonymous change in the codons encoding several highly conserved cysteine residues, interpreted as strong position-specific codon bias. This work has spawned a variety of studies on the potential mechanisms of hypermutation and on causes for cysteine codon bias, and has inspired hypermutation hypotheses for various other fast-evolving genes. Here, I show that all three findings are likely to be artifacts of statistical reconstruction. First, by simulating nonsynonymous change I show that high rates of dN can lead to overestimation of dS. Second, I show that there is no evidence for any of these three patterns in comparisons of closely related conotoxin sequences, suggesting that the reported findings are due to breakdown of statistical methods at high levels of sequence divergence. The current findings suggest that mutation and codon bias in conotoxin genes may not be atypical, and that random mutation and selection can explain the evolution of even these exceptional loci.

Dietary breadth is positively correlated with venom complexity in cone snails

BACKGROUND

Although diet is believed to be a major factor underlying the evolution of venom, few comparative studies examine both venom composition and diet across a radiation of venomous species. Cone snails within the family, Conidae, comprise more than 700 species of carnivorous marine snails that capture their prey by using a cocktail of venomous neurotoxins (conotoxins or conopeptides). Venom composition across species has been previously hypothesized to be shaped by (a) prey taxonomic class (i.e., worms, molluscs, or fish) and (b) dietary breadth. We tested these hypotheses under a comparative phylogenetic framework using ecological data from past studies in conjunction with venom duct transcriptomes sequenced from 12 phylogenetically disparate cone snail species, including 10 vermivores (worm-eating), one molluscivore, and one generalist.

RESULTS

We discovered 2223 unique conotoxin precursor peptides that encoded 1864 unique mature toxins across all species, >90 % of which are new to this study. In addition, we identified two novel gene superfamilies and 16 novel cysteine frameworks. Each species exhibited unique venom profiles, with venom composition and expression patterns among species dominated by a restricted set of gene superfamilies and mature toxins. In contrast with the dominant paradigm for interpreting Conidae venom evolution, prey taxonomic class did not predict venom composition patterns among species. We also found a significant positive relationship between dietary breadth and measures of conotoxin complexity.

CONCLUSIONS

The poor performance of prey taxonomic class in predicting venom components suggests that cone snails have either evolved species-specific expression patterns likely as a consequence of the rapid evolution of conotoxin genes, or that traditional means of categorizing prey type (i.e., worms, mollusc, or fish) and conotoxins (i.e., by gene superfamily) do not accurately encapsulate evolutionary dynamics between diet and venom composition. We also show that species with more generalized diets tend to have more complex venoms and utilize a greater number of venom genes for prey capture. Whether this increased gene diversity confers an increased capacity for evolutionary change remains to be tested. Overall, our results corroborate the key role of diet in influencing patterns of venom evolution in cone snails and other venomous radiations.

From Mollusks to Medicine: A Venomics Approach for the Discovery and Characterization of Therapeutics from Terebridae Peptide Toxins

Animal venoms comprise a diversity of peptide toxins that manipulate molecular targets such as ion channels and receptors, making venom peptides attractive candidates for the development of therapeutics to benefit human health. However, identifying bioactive venom peptides remains a significant challenge. In this review we describe our particular venomics strategy for the discovery, characterization, and optimization of Terebridae venom peptides, teretoxins. Our strategy reflects the scientific path from mollusks to medicine in an integrative sequential approach with the following steps: (1) delimitation of venomous Terebridae lineages through taxonomic and phylogenetic analyses; (2) identification and classification of putative teretoxins through omics methodologies, including genomics, transcriptomics, and proteomics; (3) chemical and recombinant synthesis of promising peptide toxins; (4) structural characterization through experimental and computational methods; (5) determination of teretoxin bioactivity and molecular function through biological assays and computational modeling; (6) optimization of peptide toxin affinity and selectivity to molecular target; and (7) development of strategies for effective delivery of venom peptide therapeutics. While our research focuses on terebrids, the venomics approach outlined here can be applied to the discovery and characterization of peptide toxins from any venomous taxa.

Target-Driven Positive Selection at Hot Spots of Scorpion Toxins Uncovers Their Potential in Design of Insecticides

Positive selection sites (PSSs), a class of amino acid sites with an excess of nonsynonymous to synonymous substitutions, are indicators of adaptive molecular evolution and have been detected in many protein families involved in a diversity of biological processes by statistical approaches. However, few studies are conducted to evaluate their functional significance and the driving force behind the evolution (i.e., agent of selection). Scorpion α-toxins are a class of multigene family of peptide neurotoxins affecting voltage-gated Na(+ )(Nav) channels, whose members exhibit differential potency and preference for insect and mammalian Nav channels. In this study, we undertook a systematical molecular dissection of nearly all the PSSs newly characterized in the Mesobuthus α-toxin family and a two-residue insertion ((19)AlaPhe(20)) located within a positively selected loop via mutational analysis of α-like MeuNaTxα-5, one member affecting both insect and mammalian Nav channels. This allows to identify hot-spot residues on its functional face involved in interaction with the receptor site of Nav channels, which comprises two PSSs (Ile(40) and Leu(41)) and the small insertion, both located on two spatially separated functional loops. Mutations at these hot-spots resulted in a remarkably decreased anti-mammalian activity in MeuNaTxα-5 with partially impaired or enhanced insecticide activity, suggesting the potential of PSSs in designing promising candidate insecticides from scorpion α-like toxins. Based on an experiment-guided toxin-channel complex model and high evolutionary variability in the receptor site of predators and prey of scorpions, we provide new evidence for target-driven adaptive evolution of scorpion toxins to deal with their targets' diversity.

Age-related association of venom gene expression and diet of predatory gastropods

Background

Venomous organisms serve as wonderful systems to study the evolution and expression of genes that are directly associated with prey capture. To evaluate the relationship between venom gene expression and prey utilization, we examined these features among individuals of different ages of the venomous, worm-eating marine snail Conus ebraeus. We determined expression levels of six genes that encode venom components, used a DNA-based approach to evaluate the identity of prey items, and compared patterns of venom gene expression and dietary specialization.

Results

C. ebraeus exhibits two major shifts in diet with age—an initial transition from a relatively broad dietary breadth to a narrower one and then a return to a broader diet. Venom gene expression patterns also change with growth. All six venom genes are up-regulated in small individuals, down-regulated in medium-sized individuals, and then either up-regulated or continued to be down-regulated in members of the largest size class. Venom gene expression is not significantly different among individuals consuming different types of prey, but instead is coupled and slightly delayed with shifts in prey diversity.

Conclusion

These results imply that changes in gene expression contribute to intraspecific variation of venom composition and that gene expression patterns respond to changes in the diversity of food resources during different growth stages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0592-5) contains supplementary material, which is available to authorized users.