Venom evolution widespread in fishes: a phylogenetic road map for the bioprospecting of piscine venoms

Knowledge of evolutionary relationships or phylogeny allows for effective predictions about the unstudied characteristics of species. These include the presence and biological activity of an organism's venoms. To date, most venom bioprospecting has focused on snakes, resulting in six stroke and cancer treatment drugs that are nearing U.S. Food and Drug Administration review. Fishes, however, with thousands of venoms, represent an untapped resource of natural products. The first step involved in the efficient bioprospecting of these compounds is a phylogeny of venomous fishes. Here, we show the results of such an analysis and provide the first explicit suborder-level phylogeny for spiny-rayed fishes. The results, based on approximately 1.1 million aligned base pairs, suggest that, in contrast to previous estimates of 200 venomous fishes, >1,200 fishes in 12 clades should be presumed venomous. This assertion was corroborated by a detailed anatomical study examining potentially venomous structures in >100 species. The results of these studies not only alter our view of the diversity of venomous fishes, now representing >50% of venomous vertebrates, but also provide the predictive phylogeny or "road map" for the efficient search for potential pharmacological agents or physiological tools from the unexplored fish venoms.     

Evolution of venom delivery structures in madtom catfishes (Siluriformes: Ictaluridae)

The use of venom to subdue prey or deter predators has evolved multiple times in numerous animal lineages. Catfishes represent one of the most easily recognized, but least studied groups of venomous fishes. Venom glands surround spines on the dorsal and pectoral fins that serve as venom delivery structures. Species of madtom catfishes in the genus Noturus were found to each have one of four venom delivery morphologies: (1) smooth spine with no venom gland; (2) smooth spine with venom gland associated with shaft of spine; (3) serrated spine with venom gland associated with shaft of spine; and (4) serrated spine with venom gland associated with shaft of spine and posterior serrations. Analyses accounting for the phylogenetic history of Noturus species suggest that a serrated pectoral spine with a venom gland is the ancestral condition for the genus. The presence of serrations and a venom gland have been largely conserved among Noturus species, but sting morphology has changed at least five times within the genus. Four of these changes have resulted in a loss of morphological complexity, including the loss of posterior serrations, loss of venom glands associated with the posterior serrations, and one complete loss of the venom gland. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 115–129.     

SdPI, the first functionally characterized Kunitz-type trypsin inhibitor from scorpion venom

Background

Kunitz-type venom peptides have been isolated from a wide variety of venomous animals. They usually have protease inhibitory activity or potassium channel blocking activity, which by virtue of the effects on predator animals are essential for the survival of venomous animals. However, no Kunitz-type peptides from scorpion venom have been functionally characterized.

Principal Findings

A new Kunitz-type venom peptide gene precursor, SdPI, was cloned and characterized from a venom gland cDNA library of the scorpion Lychas mucronatus. It codes for a signal peptide of 21 residues and a mature peptide of 59 residues. The mature SdPI peptide possesses a unique cysteine framework reticulated by three disulfide bridges, different from all reported Kunitz-type proteins. The recombinant SdPI peptide was functionally expressed. It showed trypsin inhibitory activity with high potency (K_i = 1.6×10⁻⁷ M) and thermostability.

Conclusions

The results illustrated that SdPI is a potent and stable serine protease inhibitor. Further mutagenesis and molecular dynamics simulation revealed that SdPI possesses a serine protease inhibitory active site similar to other Kunitz-type venom peptides. To our knowledge, SdPI is the first functionally characterized Kunitz-type trypsin inhibitor derived from scorpion venom, and it represents a new class of Kunitz-type venom peptides.     

Expression pattern of three-finger toxin and phospholipase A2 genes in the venom glands of two sea snakes, Lapemis curtus and Acalyptophis peronii : comparison of evolution of these toxins in land snakes, sea kraits and sea snakes

Background  

Snake venom composition varies widely both among closely related species and within the same species, based on ecological variables. In terrestrial snakes, such variation has been proposed to be due to snakes' diet. Land snakes target various prey species including insects (arthropods), lizards (reptiles), frogs and toads (amphibians), birds (aves), and rodents (mammals), whereas sea snakes target a single vertebrate class (fishes) and often specialize on specific types of fish. It is therefore interesting to examine the evolution of toxins in sea snake venoms compared to that of land snakes.     

Histology,histochemistry, and emptying mechanism of the venom glands of some elapid snakes

The venom glands of several species of elapid snakes are described. The main venom gland consists of many tubules which usually contain large amounts of secretion product. The accessory gland surrounds the entire venom duct and is usually composed of uniform mucous epithelium. The epithelium lining the tubules of the accessory gland of Naja naja is composed of two distinct types of cells. Histochemical tests indicate that the main venom gland reacts with mercury bromphenol blue and PAS but not with alcian blue. The accessory gland reacts with PAS and alcian blue, and not with mercury bromphenol blue. Treatment of sections with sialidase demonstrates the presence of a sialomucin in the accessory gland. Stimulation of the muscles associated with the venom gland offers an indication of the venom expulsion mechanism of Bungarus caeruleus. A comparison of the venom apparatus of elapid and viperid snakes emphasizes marked differences in the internal anatomy of the venom glands, muscles associated with the gland, and arrangement of glandular components. The morphological differences and dissimilar venom expulsion mechanisms support the recent view of the polyphyletic origin of venomous snakes.     

Venomous auger snail Hastula (Impages) hectica (Linnaeus, 1758): molecular phylogeny, foregut anatomy and comparative toxinology

The >10,000 living venomous marine snail species [superfamily Conoidea (Fleming, 1822)] include cone snails (Conus), the overwhelming focus of research. Hastula hectica (Linnaeus, 1758), a venomous snail in the family Terebridae (M?rch, 1852) was comprehensively investigated. The Terebridae comprise a major monophyletic group within Conoidea. H. hectica has a striking radular tooth to inject venom that looks like a perforated spear; in Conus, the tooth looks like a hypodermic needle. H. hectica venom contains a large complement of small disulfide-rich peptides, but with no apparent overlap with Conus in gene superfamilies expressed. Although Conus peptide toxins are densely post-translationally modified, no post-translationally modified amino acids were found in any Hastula venom peptide. The results suggest that different major lineages of venomous molluscs have strikingly divergent toxinological and venom-delivery strategies.     

Purification and Functional Characterisation of Rhiminopeptidase A,a Novel Aminopeptidase from the Venom of Bitis gabonica rhinoceros

Background

Snake bite is a major neglected public health issue within poor communities living in the rural areas of several countries throughout the world. An estimated 2.5 million people are bitten by snakes each year and the cost and lack of efficacy of current anti-venom therapy, together with the lack of detailed knowledge about toxic components of venom and their modes of action, and the unavailability of treatments in rural areas mean that annually there are around 125,000 deaths worldwide. In order to develop cheaper and more effective therapeutics, the toxic components of snake venom and their modes of action need to be clearly understood. One particularly poorly understood component of snake venom is aminopeptidases. These are exo-metalloproteases, which, in mammals, are involved in important physiological functions such as the maintenance of blood pressure and brain function. Although aminopeptidase activities have been reported in some snake venoms, no detailed analysis of any individual snake venom aminopeptidases has been performed so far. As is the case for mammals, snake venom aminopeptidases may also play important roles in altering the physiological functions of victims during envenomation. In order to further understand this important group of snake venom enzymes we have isolated, functionally characterised and analysed the sequence-structure relationships of an aminopeptidase from the venom of the large, highly venomous West African gaboon viper, Bitis gabonica rhinoceros.

Methodology and Principal Findings

The venom of B. g. rhinoceros was fractionated by size exclusion chromatography and fractions with aminopeptidase activities were isolated. Fractions with aminopeptidase activities showed a pure protein with a molecular weight of 150 kDa on SDS-PAGE. In the absence of calcium, this purified protein had broad aminopeptidase activities against acidic, basic and neutral amino acids but in the presence of calcium, it had only acidic aminopeptidase activity (APA). Together with the functional data, mass spectrometry analysis of the purified protein confirmed this as an aminopeptidase A and thus this has been named as rhiminopeptidase A. The complete gene sequence of rhiminopeptidase A was obtained by sequencing the PCR amplified aminopeptidase A gene from the venom gland cDNA of B. g. rhinoceros. The gene codes for a predicted protein of 955 amino acids (110 kDa), which contains the key amino acids necessary for functioning as an aminopeptidase A. A structural model of rhiminopeptidase A shows the structure to consist of 4 domains: an N-terminal saddle-shaped β domain, a mixed α and β catalytic domain, a β-sandwich domain and a C-terminal α helical domain.

Conclusions

This study describes the discovery and characterisation of a novel aminopeptidase A from the venom of B. g. rhinoceros and highlights its potential biological importance. Similar to mammalian aminopeptidases, rhiminopeptidase A might be capable of playing roles in altering the blood pressure and brain function of victims. Furthermore, it could have additional effects on the biological functions of other host proteins by cleaving their N-terminal amino acids. This study points towards the importance of complete analysis of individual components of snake venom in order to develop effective therapies for snake bites.     

Transcriptomic analysis of the venom gland of the red-headed krait (<Emphasis Type="Italic">Bungarus flaviceps</Emphasis>) using expressed sequence tags

Background  

The Red-headed krait (Bungarus flaviceps, Squamata: Serpentes: Elapidae) is a medically important venomous snake that inhabits South-East Asia. Although the venoms of most species of the snake genus Bungarus have been well characterized, a detailed compositional analysis of B. flaviceps is currently lacking.     

Species structure of the ichthyofauna of the inland waters of Phu Quoc Island, Gulf of Thailand, Vietnam

We performed the first large-scale study of the species structure of the fishes living in the inland water bodies of Phu Quoc Island and identified 97 species of bony fishes, which represent 41 families of 13 orders. Taking into account the results of earlier studies, we compiled a total checklist of the island’s ich-thyofauna, which includes 99 species; one more species recorded in a previous publication on the ichthyo-fauna of the island can be added to this number, as well as 15–20 mostly euryhaline species, which needed clarification of their states. Only one species on the checklist, Oreochromis niloticus, has been introduced to the island; 14 species, mostly of the family Gobiidae, have never been recorded in the faunas of other regions of Vietnam. Two new species, probably endemic to the island, have been found and described (in the genera Rasbora and Barbucca). Two more species of the order Siluriformes should probably also be described as new.     

Climate tracking by freshwater fishes suggests that fish diversity in temperate lakes may be increasingly threatened by climate warming

Aim

Many freshwater fishes are migrating poleward to more thermally suitable habitats in response to warming climates. In this study, we aimed to identify which freshwater fishes are most sensitive to climatic changes and asked: (i) how fast are lakes warming? (ii) how fast are fishes moving? and (iii) are freshwater fishes tracking climate?

Location

Ontario, Canada.

Methods

We assembled a database containing time series data on climate and species occurrence data from 10,732 lakes between 1986 and 2017. We calculated the rate of lake warming and climate velocity for these lakes. Climate velocities were compared with biotic velocities, specifically the rate at which the northernmost extent of each species shifted north.

Results

Lakes in Ontario warmed by 0.2°C decade⁻¹ on average, at a climate velocity of 9.4 km decade⁻¹ between 1986 and 2017. In response, some freshwater fishes have shifted their northern range boundaries with considerable interspecific variation ranging from species moving southwards at a rate of −58.9 km decade⁻¹ to species ranges moving northwards at a rate of 83.6 km decade⁻¹ over the same time period. More freshwater fish species are moving into northern lakes in Ontario than those being lost. Generally, predators are moving their range edges northwards, whereas prey fishes are being lost from northern lakes.

Main Conclusions

The concurrent loss of cooler refugia, combined with antagonistic competitive and predatory interactions with the range expanding species, has resulted in many commercially important predators moving their range edges northwards, whereas prey species have contracted their northern range edge boundaries. Trophic partitioning of range shifts highlights a previously undocumented observation of the loss of freshwater fishes from lower trophic levels in response to climate-driven migrations.     

Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua)

Background  

Marine fishes have been shown to display low levels of genetic structuring and associated high levels of gene flow, suggesting shallow evolutionary trajectories and, possibly, limited or lacking adaptive divergence among local populations. We investigated variation in 98 gene-associated single nucleotide polymorphisms (SNPs) for evidence of selection in local populations of Atlantic cod (Gadus morhua L.) across the species distribution.     

Elasmobranch qPCR reference genes: a case study of hypoxia preconditioned epaulette sharks

Background  

Elasmobranch fishes are an ancient group of vertebrates which have high potential as model species for research into evolutionary physiology and genomics. However, no comparative studies have established suitable reference genes for quantitative PCR (qPCR) in elasmobranchs for any physiological conditions. Oxygen availability has been a major force shaping the physiological evolution of vertebrates, especially fishes. Here we examined the suitability of 9 reference candidates from various functional categories after a single hypoxic insult or after hypoxia preconditioning in epaulette shark (Hemiscyllium ocellatum).     

SjAPI,the First Functionally Characterized Ascaris-Type Protease Inhibitor from Animal Venoms

Background

Serine protease inhibitors act as modulators of serine proteases, playing important roles in protecting animal toxin peptides from degradation. However, all known serine protease inhibitors discovered thus far from animal venom belong to the Kunitz-type subfamily, and whether there are other novel types of protease inhibitors in animal venom remains unclear.

Principal Findings

Here, by screening scorpion venom gland cDNA libraries, we identified the first Ascaris-type animal toxin family, which contains four members: Scorpiops jendeki Ascaris-type protease inhibitor (SjAPI), Scorpiops jendeki Ascaris-type protease inhibitor 2 (SjAPI-2), Chaerilus tricostatus Ascaris-type protease inhibitor (CtAPI), and Buthus martensii Ascaris-type protease inhibitor (BmAPI). The detailed characterization of Ascaris-type peptide SjAPI from the venom gland of scorpion Scorpiops jendeki was carried out. The mature peptide of SjAPI contains 64 residues and possesses a classical Ascaris-type cysteine framework reticulated by five disulfide bridges, different from all known protease inhibitors from venomous animals. Enzyme and inhibitor reaction kinetics experiments showed that recombinant SjAPI was a dual function peptide with α-chymotrypsin- and elastase-inhibiting properties. Recombinant SjAPI inhibited α-chymotrypsin with a Ki of 97.1 nM and elastase with a Ki of 3.7 μM, respectively. Bioinformatics analyses and chimera experiments indicated that SjAPI contained the unique short side chain functional residues “AAV” and might be a useful template to produce new serine protease inhibitors.

Conclusions/Significance

To our knowledge, SjAPI is the first functionally characterized animal toxin peptide with an Ascaris-type fold. The structural and functional diversity of animal toxins with protease-inhibiting properties suggested that bioactive peptides from animal venom glands might be a new source of protease inhibitors, which will accelerate the development of diagnostic and therapeutic agents for human diseases that target diverse proteases.     

瓦氏黄颡鱼线粒体全基因组序列分析及系统进化

鲿科鱼类种类繁多, 外形相似, 形态学分类较为困难。为了给鲿科鱼类乃至鲇形目鱼类的系统进化研究积累基础资料, 文章采用参照近缘物种线粒体基因组设计覆盖全基因组引物的方法, 利用16对引物对瓦氏黄颡鱼(Pelteobagrus vachelli)线粒体全基因组进行扩增, PCR产物转化到质粒后测序, 最终获得线粒体基因组全序列, 其全长为16 527 bp, 包括2个rRNA基因、22个tRNA基因、13个编码蛋白质基因和一个非编码控制区。瓦氏黄颡鱼(P. vachelli)线粒体基因组结构和基因排列顺序与现已公布的鲇形目鱼类完全一致, 序列分析表明, 与鲇形目其他种属间具有较高的同源性, 与拟鲿属的同源性最高(91%)。利用鲇形目共4科6属9种及3个外群的线粒体全基因组序列, 从线粒体基因组水平探讨了鲿科鱼类及其在鲇形目的系统进化地位, 结果表明: 鲿科鱼类的瓦氏黄颡鱼(P. vachelli)、黄颡鱼(Pelteobagrus fulvidraco)、光泽黄颡鱼(Pelteobagrus nitidus)及越南拟鲿(Pseudobagrus tokiensis)构成一单系群; 拟鲿属与黄颡鱼属的关系较近; 黄颡鱼属中瓦氏黄颡鱼(P. vachelli)与光泽黄颡鱼(P.nitidus)的关系近于黄颡鱼(P. fulvidraco)。     

Mitochondrial phylogeny and phylogeography of East African squeaker catfishes (Siluriformes: Synodontis)

Background  

Squeaker catfishes (Pisces, Mochokidae, Synodontis) are widely distributed throughout Africa and inhabit a biogeographic range similar to that of the exceptionally diverse cichlid fishes, including the three East African Great Lakes and their surrounding rivers. Since squeaker catfishes also prefer the same types of habitats as many of the cichlid species, we hypothesized that the East African Synodontis species provide an excellent model group for comparative evolutionary and phylogeographic analyses.     

Genomic organization of duplicated short wave-sensitive and long wave-sensitive opsin genes in the green swordtail, <Emphasis Type="Italic">Xiphophorus helleri</Emphasis>

Background  

Long wave-sensitive (LWS) opsin genes have undergone multiple lineage-specific duplication events throughout the evolution of teleost fishes. LWS repertoire expansions in live-bearing fishes (family Poeciliidae) have equipped multiple species in this family with up to four LWS genes. Given that color vision, especially attraction to orange male coloration, is important to mate choice within poeciliids, LWS opsins have been proposed as candidate genes driving sexual selection in this family. To date the genomic organization of these genes has not been described in the family Poeciliidae, and little is known about the mechanisms regulating the expression of LWS opsins in any teleost.     

Incorporating asymmetric movement costs into measures of habitat connectivity to assess impacts of hydrologic alteration to stream fishes

Aim

Hydrologic alterations are widespread in freshwater ecosystems worldwide and often detrimentally impact fish populations. Habitat suitability models are commonly used to assess these impacts, but these models frequently rely upon observed fish–habitat relationships rather than more mechanistic underpinnings. The aim of this study was to demonstrate how to incorporate swim performance into a measure of habitat connectivity at a fine scale, providing a method for assessing the availability of suitable habitat for stream fishes.

Location

We applied this technique to an endangered species, the Western Silvery Minnow Hybognathus argyritis, in the Milk River of southern Alberta, Canada. The Milk River is an augmented system, where a diversion in nearby St. Mary River augments flow by a factor >3 × (from 1–5 m³/s to 15–20 m³/s).

Methods

We used laboratory measured swim performance of Western Silvery Minnow to develop a movement cost function that was used in conjunction with a habitat suitability model to assess habitat availability via a recently developed graph‐theoretic metric, equivalent connected area (ECA).

Results

Stream augmentation altered not only habitat suitability but also habitat connectivity for this species. During augmentation, suitable habitat area declined by 81.3%. Changes in habitat connectivity were site dependent. Movement costs between habitat patches were lower during augmentation due to current‐assisted dispersal and increased distance to patches during natural flows from dried streambeds. When movement costs were incorporated into ECA, ECA decreased by 78.0% during augmentation.

Main conclusions

With changing climate and increasing anthropogenic impacts on aquatic ecosystems, understanding how freshwater fishes relate to their habitat is critical for appropriate management. In many cases, such as the Western Silvery Minnow, mitigating habitat suitability may not be sufficient, as species are unable to reach suitable habitat. The incorporation of swim performance into habitat connectivity assessments, as carried out here, can be easily adapted to other species and situations and can improve the understanding of impacts to stream fishes and increase the effectiveness of mitigation efforts.