Convergent and parallel evolution in a voltage-gated sodium channel underlies TTX-resistance in the Greater Blue-ringed Octopus: Hapalochlaena lunulata

Coevolution with toxic prey produces functional trade-offs in sodium channels of predatory snakes

Seemingly unrelated traits often share the same underlying molecular mechanisms, potentially generating a pleiotropic relationship whereby selection shaping one trait can simultaneously compromise another. While such functional trade-offs are expected to influence evolutionary outcomes, their actual relevance in nature is masked by obscure links between genotype, phenotype, and fitness. Here, we describe functional trade-offs that likely govern a key adaptation and coevolutionary dynamics in a predator-prey system. Several garter snake (Thamnophis spp.) populations have evolved resistance to tetrodotoxin (TTX), a potent chemical defense in their prey, toxic newts (Taricha spp.). Snakes achieve TTX resistance through mutations occurring at toxin-binding sites in the pore of snake skeletal muscle voltage-gated sodium channels (NaV1.4). We hypothesized that these mutations impair basic NaV functions, producing molecular trade-offs that should ultimately scale up to compromised organismal performance. We investigate biophysical costs in two snake species with unique and independently evolved mutations that confer TTX resistance. We show electrophysiological evidence that skeletal muscle sodium channels encoded by toxin-resistant alleles are functionally compromised. Furthermore, skeletal muscles from snakes with resistance genotypes exhibit reduced mechanical performance. Lastly, modeling the molecular stability of these sodium channel variants partially explains the electrophysiological and muscle impairments. Ultimately, adaptive genetic changes favoring toxin resistance appear to negatively impact sodium channel function, skeletal muscle strength, and organismal performance. These functional trade-offs at the cellular and organ levels appear to underpin locomotor deficits observed in resistant snakes and may explain variation in the population-level success of toxin-resistant alleles across the landscape, ultimately shaping the trajectory of snake-newt coevolution.

G × E interactions as a basis for toxicological uncertainty

To transfer toxicological findings from model systems, e.g. animals, to humans, standardized safety factors are applied to account for intra-species and inter-species variabilities. An alternative approach would be to measure and model the actual compound-specific uncertainties. This biological concept assumes that all observed toxicities depend not only on the exposure situation (environment = E), but also on the genetic (G) background of the model (G × E). As a quantitative discipline, toxicology needs to move beyond merely qualitative G × E concepts. Research programs are required that determine the major biological variabilities affecting toxicity and categorize their relative weights and contributions. In a complementary approach, detailed case studies need to explore the role of genetic backgrounds in the adverse effects of defined chemicals. In addition, current understanding of the selection and propagation of adverse outcome pathways (AOP) in different biological environments is very limited. To improve understanding, a particular focus is required on modulatory and counter-regulatory steps. For quantitative approaches to address uncertainties, the concept of "genetic" influence needs a more precise definition. What is usually meant by this term in the context of G × E are the protein functions encoded by the genes. Besides the gene sequence, the regulation of the gene expression and function should also be accounted for. The widened concept of past and present "gene expression" influences is summarized here as Ge. Also, the concept of "environment" needs some re-consideration in situations where exposure timing (Et) is pivotal: prolonged or repeated exposure to the insult (chemical, physical, life style) affects Ge. This implies that it changes the model system. The interaction of Ge with Et might be denoted as Ge × Et. We provide here general explanations and specific examples for this concept and show how it could be applied in the context of New Approach Methodologies (NAM).

Historical Perspective of the Characterization of Conotoxins Targeting Voltage-Gated Sodium Channels

Marine toxins have potent actions on diverse sodium ion channels regulated by transmembrane voltage (voltage-gated ion channels) or by neurotransmitters (nicotinic acetylcholine receptor channels). Studies of these toxins have focused on varied aspects of venom peptides ranging from evolutionary relationships of predator and prey, biological actions on excitable tissues, potential application as pharmacological intervention in disease therapy, and as part of multiple experimental approaches towards an understanding of the atomistic characterization of ion channel structure. This review examines the historical perspective of the study of conotoxin peptides active on sodium channels gated by transmembrane voltage, which has led to recent advances in ion channel research made possible with the exploitation of the diversity of these marine toxins.

In Silico Analysis of Tetrodotoxin Binding in Voltage-Gated Sodium Ion Channels from Toxin-Resistant Animal Lineages

Multiple animal species have evolved resistance to the neurotoxin tetrodotoxin (TTX) through changes in voltage-gated sodium ion channels (VGSCs). Amino acid substitutions in TTX-resistant lineages appear to be positionally convergent with changes in homologous residues associated with reductions in TTX block. We used homology modeling coupled with docking simulations to test whether positionally convergent substitutions generate functional convergence at the level of TTX–channel interactions. We found little evidence that amino acids at convergent positions generated similar patterns among TTX-resistant animal lineages across several metrics, including number of polar contacts, polar contact position, and estimates of binding energy. Though binding energy values calculated for TTX docking were reduced for some TTX-resistant channels, not all TTX-resistant channels and not all of our analyses returned reduced binding energy values for TTX-resistant channels. Our results do not support a simple model of toxin resistance where a reduced number of bonds between TTX and the channel protein prevents blocking. Rather models that incorporate flexibility and movement of the protein overall may better describe how homologous substitutions in the channel cause changes in TTX block.

The road not taken: Evolution of tetrodotoxin resistance in the Sierra garter snake (Thamnophis couchii) by a path less traveled

The repeated evolution of tetrodotoxin (TTX) resistance provides a model for testing hypotheses about the mechanisms of convergent evolution. This poison is broadly employed as a potent antipredator defense, blocking voltage-gated sodium channels (Na_v ) in muscles and nerves, paralyzing and sometimes killing predators. Resistance in taxa bearing this neurotoxin and a few predators appears to come from convergent replacements in specific Na_v residues that interact with TTX. This stereotyped genetic response suggests molecular and phenotypic evolution may be constrained and predictable. Here, we investigate the extent of mechanistic convergence in garter snakes (Thamnophis) that prey on TTX-bearing newts (Taricha) by examining the physiological and genetic basis of TTX resistance in the Sierra garter snake (Th. couchii). We characterize variation in this predatory adaptation across populations at several biological scales: whole-animal TTX resistance; skeletal muscle resistance, functional genetic variation in three Na_v encoding loci; and levels of gene expression for one of these loci. We found Th. couchii possess extensive geographic variation in resistance at the whole-animal and skeletal muscle levels. As in other Thamnophis, resistance at both levels is highly correlated, suggesting convergence across the biological levels linking organism to organ. However, Th. couchii shows no functional variation in Na_v loci among populations or difference in candidate gene expression. Local variation in TTX resistance in Th. couchii cannot be explained by the same relationship between genotype and phenotype seen in other taxa. Thus, historical contingencies may lead different species of Thamnophis down alternative routes to local adaptation.

Convergent evolution of toxin resistance in animals

Convergence is the phenomenon whereby similar phenotypes evolve independently in different lineages. One example is resistance to toxins in animals. Toxins have evolved many times throughout the tree of life. They disrupt molecular and physiological pathways in target species, thereby incapacitating prey or deterring a predator. In response, molecular resistance has evolved in many species exposed to toxins to counteract their harmful effects. Here, we review current knowledge on the convergence of toxin resistance using examples from a wide range of toxin families. We explore the evolutionary processes and molecular adaptations driving toxin resistance. However, resistance adaptations may carry a fitness cost if they disrupt the normal physiology of the resistant animal. Therefore, there is a trade‐off between maintaining a functional molecular target and reducing toxin susceptibility. There are relatively few solutions that satisfy this trade‐off. As a result, we see a small set of molecular adaptations appearing repeatedly in diverse animal lineages, a phenomenon that is consistent with models of deterministic evolution. Convergence may also explain what has been called ‘autoresistance’. This is often thought to have evolved for self‐protection, but we argue instead that it may be a consequence of poisonous animals feeding on toxic prey. Toxin resistance provides a unique and compelling model system for studying the interplay between trophic interactions, selection pressures and the molecular mechanisms underlying evolutionary novelties.

An Updated Review of Tetrodotoxin and Its Peculiarities

Tetrodotoxin (TTX) is a crystalline, weakly basic, colorless organic substance and is one of the most potent marine toxins known. Although TTX was first isolated from pufferfish, it has been found in numerous other marine organisms and a few terrestrial species. Moreover, tetrodotoxication is still an important health problem today, as TTX has no known antidote. TTX poisonings were most commonly reported from Japan, Thailand, and China, but today the risk of TTX poisoning is spreading around the world. Recent studies have shown that TTX-containing fish are being found in other regions of the Pacific and in the Indian Ocean, as well as the Mediterranean Sea. This review aims to summarize pertinent information available to date on the structure, origin, distribution, mechanism of action of TTX and analytical methods used for the detection of TTX, as well as on TTX-containing organisms, symptoms of TTX poisoning, and incidence worldwide.

Genome Analysis of Lagocephalus sceleratus: Unraveling the Genomic Landscape of a Successful Invader

The Tetraodontidae family encompasses several species which attract scientific interest in terms of their ecology and evolution. The silver-cheeked toadfish (Lagocephalus sceleratus) is a well-known “invasive sprinter” that has invaded and spread, in less than a decade, throughout the Eastern and part of the Western Mediterranean Sea from the Red Sea through the Suez Canal. In this study, we built and analysed the first near-chromosome level genome assembly of L. sceleratus and explored its evolutionary landscape. Through a phylogenomic analysis, we positioned L. sceleratus closer to T. nigroviridis, compared to other members of the family, while gene family evolution analysis revealed that genes associated with the immune response have experienced rapid expansion, providing a genetic basis for studying how L. sceleratus is able to achieve highly successful colonisation. Moreover, we found that voltage-gated sodium channel (NaV 1.4) mutations previously connected to tetrodotoxin resistance in other pufferfishes are not found in L. sceleratus, highlighting the complex evolution of this trait. The high-quality genome assembly built here is expected to set the ground for future studies on the species biology.

Heterologous functional expression of ascidian Nav1 channels and close relationship with the evolutionary ancestor of vertebrate Nav channels

Voltage-gated sodium channels (Nav1s) are responsible for the initiation and propagation of action potentials in neurons, muscle, and endocrine cells. Many clinically used drugs such as local anesthetics and antiarrhythmics inhibit Nav1s, and a variety of inherited human disorders are caused by mutations in Nav1 genes. Nav1s consist of the main α subunit and several auxiliary β subunits. Detailed information on the structure–function relationships of Nav1 subunits has been obtained through heterologous expression experiments and analyses of protein structures. The basic properties of Nav1s, including their gating and ion permeation, were classically described in the squid giant axon and other invertebrates. However, heterologous functional expression of Nav1s from marine invertebrates has been unsuccessful. Ascidians belong to the Urochordata, a sister group of vertebrates, and the larval central nervous system of ascidians shows a similar plan to that of vertebrates. Here, we report the biophysical properties of ascidian Ciona Nav1 (CiNav1a) heterologously expressed in Xenopus oocytes. CiNav1a exhibited tetrodotoxin-insensitive sodium currents with rapid gating kinetics of activation and inactivation. Furthermore, consistent with the fact that the Ciona genome lacks orthologous genes to vertebrate β subunits, the human β1 subunit did not influence the gating properties when coexpressed with CiNav1a. Interestingly, CiNav1a contains an ankyrin-binding motif in the II–III linker, which can be targeted to the axon initial segment of mammalian cortical neurons. Our findings provide a platform to gain insight into the evolutionary and biophysical properties of Nav1s, which are important for the development of targeted therapeutics.

Adaptive venom evolution and toxicity in octopods is driven by extensive novel gene formation, expansion, and loss

BACKGROUND

Cephalopods represent a rich system for investigating the genetic basis underlying organismal novelties. This diverse group of specialized predators has evolved many adaptations including proteinaceous venom. Of particular interest is the blue-ringed octopus genus (Hapalochlaena), which are the only octopods known to store large quantities of the potent neurotoxin, tetrodotoxin, within their tissues and venom gland.

FINDINGS

To reveal genomic correlates of organismal novelties, we conducted a comparative study of 3 octopod genomes, including the Southern blue-ringed octopus (Hapalochlaena maculosa). We present the genome of this species and reveal highly dynamic evolutionary patterns at both non-coding and coding organizational levels. Gene family expansions previously reported in Octopus bimaculoides (e.g., zinc finger and cadherins, both associated with neural functions), as well as formation of novel gene families, dominate the genomic landscape in all octopods. Examination of tissue-specific genes in the posterior salivary gland revealed that expression was dominated by serine proteases in non-tetrodotoxin-bearing octopods, while this family was a minor component in H. maculosa. Moreover, voltage-gated sodium channels in H. maculosa contain a resistance mutation found in pufferfish and garter snakes, which is exclusive to the genus. Analysis of the posterior salivary gland microbiome revealed a diverse array of bacterial species, including genera that can produce tetrodotoxin, suggestive of a possible production source.

CONCLUSIONS

We present the first tetrodotoxin-bearing octopod genome H. maculosa, which displays lineage-specific adaptations to tetrodotoxin acquisition. This genome, along with other recently published cephalopod genomes, represents a valuable resource from which future work could advance our understanding of the evolution of genomic novelty in this family.

The chemistry and biology of guanidine secondary metabolites

Covering: 2017-2019Guanidine natural products isolated from microorganisms, marine invertebrates and terrestrial plants, amphibians and spiders, represented by non-ribosomal peptides, guanidine-bearing polyketides, alkaloids, terpenoids and shikimic acid derived, are the subject of this review. The topics include the discovery of new metabolites, total synthesis of natural guanidine compounds, biological activity and mechanism-of-action, biosynthesis and ecological functions.

Structures of N-Hydroxy-Type Tetrodotoxin Analogues and Bicyclic Guanidinium Compounds Found in Toxic Newts

The biosynthesis of tetrodotoxin (TTX, 1), a potent neurotoxin widely distributed in marine and terrestrial metazoans, remains unresolved. A significant issue has been identifying intermediates and shunt products associated with the biosynthetic pathway of TTX. We investigated TTX biosynthesis by screening and identifying new TTX-related compounds from Cynops ensicauda popei and Taricha granulosa. Mass spectrometry (MS)-guided screening identified two new N-hydroxy TTX analogues in newts: 1-hydroxy-8-epiTTX (2) and 1-hydroxy-8-epi-5,11-dideoxyTTX (3, previously reported as 1-hydroxy-5,11-dideoxyTTX). We prepared a new analogue, 8-epi-5,11-dideoxyTTX (4), from 3 via N-OH reduction and confirmed the presence of 4 in T. granulosa using hydrophilic interaction liquid chromatography (HILIC)-LCMS. The presence of 8-epi-type TTX analogues in both Cynops and Taricha supports a branched biosynthetic pathway of terrestrial TTX, which produces 6- and 8-epimers. In addition, new bicyclic guanidinium compounds Tgr-238 (5) and Tgr-240 (6) were identified as putative shunt products of our proposed TTX biosynthesis pathway. A structural analysis of Cep-228A (7), another bicyclic compound, was performed using NMR. Based on the structures of 5-7 and their analogues, we propose a model of the shunt and metabolic pathways of the terrestrial TTX biosynthesis.

The gene-rich genome of the scallop Pecten maximus

BACKGROUND

The king scallop, Pecten maximus, is distributed in shallow waters along the Atlantic coast of Europe. It forms the basis of a valuable commercial fishery and plays a key role in coastal ecosystems and food webs. Like other filter feeding bivalves it can accumulate potent phytotoxins, to which it has evolved some immunity. The molecular origins of this immunity are of interest to evolutionary biologists, pharmaceutical companies, and fisheries management.

FINDINGS

Here we report the genome assembly of this species, conducted as part of the Wellcome Sanger 25 Genomes Project. This genome was assembled from PacBio reads and scaffolded with 10X Chromium and Hi-C data. Its 3,983 scaffolds have an N50 of 44.8 Mb (longest scaffold 60.1 Mb), with 92% of the assembly sequence contained in 19 scaffolds, corresponding to the 19 chromosomes found in this species. The total assembly spans 918.3 Mb and is the best-scaffolded marine bivalve genome published to date, exhibiting 95.5% recovery of the metazoan BUSCO set. Gene annotation resulted in 67,741 gene models. Analysis of gene content revealed large numbers of gene duplicates, as previously seen in bivalves, with little gene loss, in comparison with the sequenced genomes of other marine bivalve species.

CONCLUSIONS

The genome assembly of P. maximus and its annotated gene set provide a high-quality platform for studies on such disparate topics as shell biomineralization, pigmentation, vision, and resistance to algal toxins. As a result of our findings we highlight the sodium channel gene Nav1, known to confer resistance to saxitoxin and tetrodotoxin, as a candidate for further studies investigating immunity to domoic acid.