Evolution Of Toxins Research Articles

Reconstructing the Origins of the Somatostatin and Allatostatin-C Signaling Systems Using the Accelerated Evolution of Biodiverse Cone Snail Toxins.

Somatostatin and its related peptides (SSRPs) form an important family of hormones with diverse physiological roles. The ubiquitous presence of SSRPs in vertebrates and several invertebrate deuterostomes suggests an ancient origin of the SSRP signaling system. However, the existence of SSRP genes outside of deuterostomes has not been established, and the evolutionary history of this signaling system remains poorly understood. Our recent discovery of SSRP-like toxins (consomatins) in venomous marine cone snails (Conus) suggested the presence of a related signaling system in mollusks and potentially other protostomes. Here, we identify the molluscan SSRP-like signaling gene that gave rise to the consomatin family. Following recruitment into venom, consomatin genes experienced strong positive selection and repeated gene duplications resulting in the formation of a hyperdiverse family of venom peptides. Intriguingly, the largest number of consomatins was found in worm-hunting species (>400 sequences), indicating a homologous system in annelids, another large protostome phylum. Consistent with this, comprehensive sequence mining enabled the identification of SSRP-like sequences (and their corresponding orphan receptor) in annelids and several other protostome phyla. These results established the existence of SSRP-like peptides in many major branches of bilaterians and challenge the prevailing hypothesis that deuterostome SSRPs and protostome allatostatin-C are orthologous peptide families. Finally, having a large set of predator–prey SSRP sequences available, we show that although the cone snail’s signaling SSRP-like genes are under purifying selection, the venom consomatin genes experience rapid directional selection to target receptors in a changing mix of prey.

Toxic Relationships and Arms-Race Coevolution Revisited.

Toxin evolution in animals is one of the most fascinating and complex subjects of scientific inquiry today. Gaining an understanding of toxins poses a multifaceted challenge given the diverse modes of acquisition, evolutionary adaptations, and abiotic components that affect toxin phenotypes. Here, we highlight some of the main genetic and ecological factors that influence toxin evolution and discuss the role of antagonistic interactions and coevolutionary dynamics in shaping the direction and extent of toxicity and resistance in animals. We focus on toxic Pacific newts (family Salamandridae, genus Taricha) as a system to investigate and better evaluate the widely distributed toxin they possess, tetrodotoxin (TTX), and the hypothesized model of arms-race coevolution with snake predators that is used to explain phenotypic patterns of newt toxicity. Finally, we propose an alternative coevolutionary model that incorporates TTX-producing bacteria and draws from an elicitor-receptor concept to explain TTX evolution and ecology.

Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins

BackgroundThe explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes.ResultsTranscriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae – Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae – Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms.ConclusionsWe found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator–prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Auxiliary interfaces support the evolution of specific toxin-antitoxin pairing.

Toxin-antitoxin (TA) systems are a large family of genes implicated in the regulation of bacterial growth and its arrest in response to attacks. These systems encode nonsecreted toxins and antitoxins that specifically pair, even when present in several paralogous copies per genome. Salmonella enterica serovar Typhimurium contains three paralogous TacAT systems that block bacterial translation. We determined the crystal structures of the three TacAT complexes to understand the structural basis of specific TA neutralization and the evolution of such specific pairing. In the present study, we show that alteration of a discrete structural add-on element on the toxin drives specific recognition by their cognate antitoxin underpinning insulation of the three pairs. Similar to other TA families, the region supporting TA-specific pairing is key to neutralization. Our work reveals that additional TA interfaces beside the main neutralization interface increase the safe space for evolution of pairing specificity.

Insights into the genome architecture and evolution of Shiga toxin encoding bacteriophages of Escherichia coli

BackgroundA total of 179 Shiga toxin-producing Escherichia coli (STEC) complete genomes were analyzed in terms of serotypes, prophage coding regions, and stx gene variants and their distribution. We further examined the genetic diversity of Stx-converting phage genomes (Stx phages), focusing on the lysis-lysogeny decision and lytic cassettes.ResultsWe show that most STEC isolates belong to non-O157 serotypes (73 %), regardless the sources and geographical regions. While the majority of STEC genomes contain a single stx gene (61 %), strains containing two (35 %), three (3 %) and four (1 %) stx genes were also found, being stx2 the most prevalent gene variant. Their location is exclusively found in intact prophage regions, indicating that they are phage-borne. We further demonstrate that Stx phages can be grouped into four clusters (A, B, C and D), three subclusters (A1, A2 and A3) and one singleton, based on their shared gene content. This cluster distribution is in good agreement with their predicted virion morphologies. Stx phage genomes are highly diverse with a vast number of 1,838 gene phamilies (phams) of related sequences (of which 677 are orphams i.e. unique genes) and, although having high mosaicism, they are generally organized into three major transcripts. While the mechanisms that guide lysis–lysogeny decision are complex, there is a strong selective pressure to maintain the stx genes location close to the lytic cassette composed of predicted SAR-endolysin and pin-holin lytic proteins. The evolution of STEC Stx phages seems to be strongly related to acquiring genetic material, probably from horizontal gene transfer events.ConclusionsThis work provides novel insights on the genetic structure of Stx phages, showing a high genetic diversity throughout the genomes, where the various lysis-lysogeny regulatory systems are in contrast with an uncommon, but conserved, lytic system always adjacent to stx genes.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Identification of Novel Toxin Genes from the Stinging Nettle Caterpillar Parasa lepida (Cramer, 1799): Insights into the Evolution of Lepidoptera Toxins.

Simple SummaryMany caterpillar species can produce toxins that cause harmful reactions to humans, varying from mild irritation to death. Currently, there is very limited knowledge about caterpillar toxin diversity, because only a few species have been investigated. We used the transcriptome technique to identify candidate toxin genes from the nettle caterpillar Parasa lepida (Cramer, 1799). It is a common pest of oil palm, coconut, and mango in South and South-East Asia, which can cause severe pain and allergic responses to those in contact with them. We reported 168 candidate toxin genes. Most of them are members of the toxin genes families commonly recruited in animal venoms such as serine protease and serine protease inhibitors. However, we identified 21 novel genes encoding knottin-like peptides expressed at a high level in the transcriptome. Their predicted 3D structures are similar to neurotoxins in scorpion and tarantula. Our study suggests that P. lepida venom contains diverse toxin proteins that potentially cause allergic reactions and pain. This study sheds light on the hidden diversity of toxin proteins in caterpillar lineage, which could be future fruitful new drug sources.Many animal species can produce venom for defense, predation, and competition. The venom usually contains diverse peptide and protein toxins, including neurotoxins, proteolytic enzymes, protease inhibitors, and allergens. Some drugs for cancer, neurological disorders, and analgesics were developed based on animal toxin structures and functions. Several caterpillar species possess venoms that cause varying effects on humans both locally and systemically. However, toxins from only a few species have been investigated, limiting the full understanding of the Lepidoptera toxin diversity and evolution. We used the RNA-seq technique to identify toxin genes from the stinging nettle caterpillar, Parasa lepida (Cramer, 1799). We constructed a transcriptome from caterpillar urticating hairs and reported 34,968 unique transcripts. Using our toxin gene annotation pipeline, we identified 168 candidate toxin genes, including protease inhibitors, proteolytic enzymes, and allergens. The 21 P. lepida novel Knottin-like peptides, which do not show sequence similarity to any known peptide, have predicted 3D structures similar to tarantula, scorpion, and cone snail neurotoxins. We highlighted the importance of convergent evolution in the Lepidoptera toxin evolution and the possible mechanisms. This study opens a new path to understanding the hidden diversity of Lepidoptera toxins, which could be a fruitful source for developing new drugs.

ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics.

Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide. Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.

Rapid molecular evolution of Spiroplasma symbionts of Drosophila.

Spiroplasma is a genus of Mollicutes whose members include plant pathogens, insect pathogens and endosymbionts of animals. Spiroplasma phenotypes have been repeatedly observed to be spontaneously lost in Drosophila cultures, and several studies have documented a high genomic turnover in Spiroplasma symbionts and plant pathogens. These observations suggest that Spiroplasma evolves quickly in comparison to other insect symbionts. Here, we systematically assess evolutionary rates and patterns of Spiroplasma poulsonii , a natural symbiont of Drosophila. We analysed genomic evolution of sHy within flies, and sMel within in vitro culture over several years. We observed that S. poulsonii substitution rates are among the highest reported for any bacteria, and around two orders of magnitude higher compared with other inherited arthropod endosymbionts. The absence of mismatch repair loci mutS and mutL is conserved across Spiroplasma , and likely contributes to elevated substitution rates. Further, the closely related strains sMel and sHy (>99.5 % sequence identity in shared loci) show extensive structural genomic differences, which potentially indicates a higher degree of host adaptation in sHy, a protective symbiont of Drosophila hydei. Finally, comparison across diverse Spiroplasma lineages confirms previous reports of dynamic evolution of toxins, and identifies loci similar to the male-killing toxin Spaid in several Spiroplasma lineages and other endosymbionts. Overall, our results highlight the peculiar nature of Spiroplasma genome evolution, which may explain unusual features of its evolutionary ecology.

Variation in venom composition in the Australian funnel-web spiders Hadronyche valida.

Mygalomorph venom properties and active components, which have importance in medicine, agronomy, venomics, ecology and evolution, have been widely studied, but only a small fraction have been characterised. Several studies have shown inter-individual variation in the composition of venom peptides based on ontogeny, sexual dimorphism, season and diet. However, intra-individual variation in venom composition, which could play a key role in the evolution, diversification and function of toxins, is poorly understood. In this study, we demonstrate significant intra- and inter-individual variation in venom composition in the Australian funnel-web spider Hadronyche valida, highlighting that individuals show different venom profiles over time. Fourteen (four juvenile and ten adult females) funnel-web spiders, maintained under the same environmental conditions and diet, were milked a total of four times, one month apart. We then used reversed-phase high performance liquid chromatography/electrospray ionisation mass spectrometry to generate venom fingerprints containing the retention time and molecular weights of the different toxin components in the venom. Across all individuals, we documented a combined total of 83 individual venom components. Only 20% of these components were shared between individuals. Individuals showed variation in the composition of venom peptides, with some components consistently present over time, while others were only present at specific times. When individuals were grouped using the Jaccard clustering index and Kernel Principal Component Analysis, spiders formed two distinct clusters, most likely due to their origin or time of collection. This study contributes to the understanding of variation in venom composition at different levels (intra-individual, and intra- and inter-specific) and considers some of the mechanisms of selection that may contribute to venom diversification within arachnids. In addition, inter-specific variation in venom composition can be highly useful as a chemotaxonomic marker to identify funnel-web species.

Toxin variation among salamander populations: discussing potential causes and future directions.

Amphibians produce defensive chemicals which provide protection against both predators and infections. Within species, populations can differ considerably in the composition and amount of these chemical defenses. Studying intraspecific variation in toxins and linking it to environmental variables may help us to identify the selective drivers of toxin evolution, such as predation pressure and infection risk. Recently, there has been a renewed interest in the unique toxins produced by salamanders from the genus Salamandra: the samandarines. Despite this attention, intraspecific variation has largely been ignored within Salamandra-species. The aim of this study was to investigate whether geographic variation in profiles of samandarines exists, by sampling 4 populations of Salamandra atra over its range in the Dinaric Alps. In addition, we preliminary explored whether potential variation could be explained by predation (counting the number of snake species) and infection risk (cultivation and genomic analyses of collected soil samples). Salamanders from the 4 populations differed in toxin composition and in the size of their poison glands, although not in overall toxin quantity. Nor predation nor infection risk could explain this variation, as populations barely differed in these variables. Sampling over a much broader geographic range, using better estimators for predation and infection risk, will contribute to an improved understanding of how environment may shape variation in chemical defenses. Nevertheless, as the 4 populations of S. atra did differ in their toxin profiles, we propose that this species provides an interesting opportunity for further ecological and evolutionary studies on amphibian toxins.

How a Scorpion Toxin Selectively Captures a Prey Sodium Channel: The Molecular and Evolutionary Basis Uncovered.

The growing resistance of insects to chemical pesticides is reducing the effectiveness of conventional methods for pest control and thus, the development of novel insecticidal agents is imperative. Scorpion toxins specific for insect voltage-gated sodium channels (Navs) have been considered as one of the most promising insecticide alternatives due to their host specificity, rapidly evoked toxicity, biodegradability, and the lack of resistance. However, they have not been developed for uses in agriculture and public health, mainly because of a limited understanding of their molecular and evolutionary basis controlling their phylogenetic selectivity. Here, we show that the traditionally defined insect-selective scorpion toxin LqhIT2 specifically captures a prey Nav through a conserved trapping apparatus comprising a three-residue-formed cavity and a structurally adjacent leucine. The former serves as a detector to recognize and bind a highly exposed channel residue conserved in insects and spiders, two major prey items for scorpions; and the latter subsequently seizes the "moving" voltage sensor via hydrophobic interactions to reduce activation energy for channel opening, demonstrating its action in an enzyme-like manner. Based on the established toxin-channel interaction model in combination with toxicity assay, we enlarged the toxic spectrum of LqhIT2 to spiders and certain other arthropods. Furthermore, we found that genetic background-dependent cavity shapes determine the species selectivity of LqhIT2-related toxins. We expect that the discovery of the trapping apparatus will improve our understanding of the evolution and design principle of Nav-targeted toxins from a diversity of arthropod predators and accelerate their uses in pest control.

The evolutionary dynamics of venom toxins made by insects and other animals

Animal venoms are recognised as unique biological systems in which to study molecular evolution. Venom use has evolved numerous times among the insects, and insects today use venom to capture prey, defend themselves from predators, or to subdue and modulate host responses during parasitism. However, little is known about most insect venom toxins or the mode and tempo by which they evolve. Here, I review the evolutionary dynamics of insect venom toxins, and argue that insects offer many opportunities to examine novel aspects of toxin evolution. The key questions addressed are: How do venomous animals evolve from non-venomous animals, and how does this path effect the composition and pharmacology of the venom? What genetic processes (gene duplication, co-option, neofunctionalisation) are most important in toxin evolution? What kinds of selection pressures are acting on toxin-encoding genes and their cognate targets in envenomated animals? The emerging evidence highlights that venom composition and pharmacology adapts quickly in response to changing selection pressures resulting from new ecological interactions, and that such evolution occurs through a stunning variety of genetic mechanisms. Insects offer many opportunities to investigate the evolutionary dynamics of venom toxins due to their evolutionary history rich in venom-related adaptations, and their quick generation time and suitability for culture in the laboratory.

Subtyping analysis reveals new variants and accelerated evolution of Clostridioides difficile toxin B

Clostridioides difficile toxins (TcdA and TcdB) are major exotoxins responsible for C. difficile infection (CDI) associated diseases. The previously reported TcdB variants showed distinct biological features, immunoactivities, and potential pathogenicity in disease progression. Here, we performed global comparisons of amino acid sequences of both TcdA and TcdB from 3,269 C. difficile genomes and clustered them according to the evolutionary relatedness. We found that TcdB was much diverse and could be divided into eight subtypes, of which four were first described. Further analysis indicates that the tcdB gene undergoes accelerated evolution to maximize diversity. By tracing TcdB subtypes back to their original isolates, we found that the distribution of TcdB subtypes was not completely aligned with the phylogeny of C. difficile. These findings suggest that the tcdB genes not only frequently mutate, but also continuously transfer and exchange among C. difficile strains.

Venom without glands: Novel methods to investigate toxin diversity, function and evolution in ribbon worms (Nemertea)

Venom without glands: Novel methods to investigate toxin diversity, function and evolution in ribbon worms (Nemertea)

Detection of Abrin-Like and Prepropulchellin-Like Toxin Genes and Transcripts Using Whole Genome Sequencing and Full-Length Transcript Sequencing of Abrus precatorius.

The sequenced genome and the leaf transcriptome of a near relative of Abrus pulchellus and Abrus precatorius was analyzed to characterize the genetic basis of toxin gene expression. From the high-quality genome assembly, a total of 26 potential coding regions were identified that contain genes with abrin-like, pulchellin-like, and agglutinin-like homology, with full-length transcripts detected in leaf tissue for 9 of the 26 coding regions. All of the toxin-like genes were identified within only five isolated regions of the genome, with each region containing 1 to 16 gene variants within each genomic region (<1 Mbp). The Abrus precatorius cultivar sequenced here contains genes which encode for proteins that are homologous to certain abrin and prepropulchellin genes previously identified, and we observed substantial diversity of genes and predicted gene products in Abrus precatorius and previously characterized toxins. This suggests diverse toxin repertoires within Abrus, potentially the results of rapid toxin evolution.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents.

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Exploring Toxin Evolution: Venom Protein Transcript Sequencing and Transcriptome-Guided High-Throughput Proteomics.

Studying animal toxin evolution requires sequences of these proteins and peptides, and transcript sequences allow for the construction of cladograms and evaluation of selection pressures from nonsynonymous and synonymous nucleotide mutation ratios. In addition, these translated sequences can be useful as custom databases for peptide identifications within venoms and for better proteomic quantification. Obtaining these transcripts is achieved by sequencing cDNA originating from venom gland tissue or venom. This chapter provides the methodology for (1) targeted sequencing of transcripts from a single venom protein family (RNA isolation and 3'RACE [rapid amplification of cDNA ends]), (2) generation of a venom gland transcriptome with next-generation sequencing (NGS) technology (de novo transcriptome assembly, toxin transcript identification, quantification, and positive selection analysis), and (3) combined high-throughput proteomics to identify secreted venom components. Transcriptomics has become fundamental for studying toxin evolution, but it creates many challenges for scientists who are unfamiliar with working with RNA, managing large NGS datasets and executing the required programs, particularly considering that there is an overabundance of available software in this field and not all perform optimally for venom gland transcriptome assembly. This chapter provides one pipeline for the integration of both low- and high-throughput transcriptomics with proteomics to characterize venoms.

The draft genome of Actinia tenebrosa reveals insights into toxin evolution.

Sea anemones have a wide array of toxic compounds (peptide toxins found in their venom) which have potential uses as therapeutics. To date, the majority of studies characterizing toxins in sea anemones have been restricted to species from the superfamily, Actinioidea. No highly complete draft genomes are currently available for this superfamily, however, highlighting our limited understanding of the genes encoding toxins in this important group. Here we have sequenced, assembled, and annotated a draft genome for Actinia tenebrosa. The genome is estimated to be approximately 255 megabases, with 31,556 protein‐coding genes. Quality metrics revealed that this draft genome matches the quality and completeness of other model cnidarian genomes, including Nematostella, Hydra, and Acropora. Phylogenomic analyses revealed strong conservation of the Cnidaria and Hexacorallia core‐gene set. However, we found that lineage‐specific gene families have undergone significant expansion events compared with shared gene families. Enrichment analysis performed for both gene ontologies, and protein domains revealed that genes encoding toxins contribute to a significant proportion of the lineage‐specific genes and gene families. The results make clear that the draft genome of A. tenebrosa will provide insight into the evolution of toxins and lineage‐specific genes, and provide an important resource for the discovery of novel biological compounds.

Toxins from scratch? Diverse, multimodal gene origins in the predatory robber fly Dasypogon diadema indicate a dynamic venom evolution in dipteran insects.

BackgroundVenoms and the toxins they contain represent molecular adaptations that have evolved on numerous occasions throughout the animal kingdom. However, the processes that shape venom protein evolution are poorly understood because of the scarcity of whole-genome data available for comparative analyses of venomous species.ResultsWe performed a broad comparative toxicogenomic analysis to gain insight into the genomic mechanisms of venom evolution in robber flies (Asilidae). We first sequenced a high-quality draft genome of the hymenopteran hunting robber fly Dasypogon diadema, analysed its venom by a combined proteotranscriptomic approach, and compared our results with recently described robber fly venoms to assess the general composition and major components of asilid venom. We then applied a comparative genomics approach, based on 1 additional asilid genome, 10 high-quality dipteran genomes, and 2 lepidopteran outgroup genomes, to reveal the evolutionary mechanisms and origins of identified venom proteins in robber flies.ConclusionsWhile homologues were identified for 15 of 30 predominant venom protein in the non-asilid genomes, the remaining 15 highly expressed venom proteins appear to be unique to robber flies. Our results reveal that the venom of D. diadema likely evolves in a multimodal fashion comprising (i) neofunctionalization after gene duplication, (ii) expression-dependent co-option of proteins, and (iii) asilid lineage-specific orphan genes with enigmatic origin. The role of such orphan genes is currently being disputed in evolutionary genomics but has not been discussed in the context of toxin evolution. Our results display an unexpected dynamic venom evolution in asilid insects, which contrasts the findings of the only other insect toxicogenomic evolutionary analysis, in parasitoid wasps (Hymenoptera), where toxin evolution is dominated by single gene co-option. These findings underpin the significance of further genomic studies to cover more neglected lineages of venomous taxa and to understand the importance of orphan genes as possible drivers for venom evolution.