Diversity and Evolution of Frog Visual Opsins: Spectral Tuning and Adaptation to Distinct Light Environments.

Visual systems adapt to different light environments through several avenues including optical changes to the eye and neurological changes in how light signals are processed and interpreted. Spectral sensitivity can evolve via changes to visual pigments housed in the retinal photoreceptors through gene duplication and loss, differential and coexpression, and sequence evolution. Frogs provide an excellent, yet understudied, system for visual evolution research due to their diversity of ecologies (including biphasic aquatic-terrestrial life cycles) that we hypothesize imposed different selective pressures leading to adaptive evolution of the visual system, notably the opsins that encode the protein component of the visual pigments responsible for the first step in visual perception. Here, we analyze the diversity and evolution of visual opsin genes from 93 new eye transcriptomes plus published data for a combined dataset spanning 122 frog species and 34 families. We find that most species express the four visual opsins previously identified in frogs but show evidence for gene loss in two lineages. Further, we present evidence of positive selection in three opsins and shifts in selective pressures associated with differences in habitat and life history, but not activity pattern. We identify substantial novel variation in the visual opsins and, using microspectrophotometry, find highly variable spectral sensitivities, expanding known ranges for all frog visual pigments. Mutations at spectral-tuning sites only partially account for this variation, suggesting that frogs have used tuning pathways that are unique among vertebrates. These results support the hypothesis of adaptive evolution in photoreceptor physiology across the frog tree of life in response to varying environmental and ecological factors and further our growing understanding of vertebrate visual evolution.

Diversity and Molecular Evolution of Antimicrobial Peptides in Caecilian Amphibians.

Antimicrobial peptides (AMPs) are key molecules in the innate immune defence of vertebrates with rapid action, broad antimicrobial spectrum, and ability to evade pathogen resistance mechanisms. To date, amphibians are the major group of vertebrates from which most AMPs have been characterised, but most studies have focused on the bioactive skin secretions of anurans (frogs and toads). In this study, we have analysed the complete genomes and/or transcriptomes of eight species of caecilian amphibians (order Gymnophiona) and characterised the diversity, molecular evolution, and antimicrobial potential of the AMP repertoire of this order of amphibians. We have identified 477 candidate AMPs within the studied caecilian genome and transcriptome datasets. These candidates are grouped into 29 AMP families, with four corresponding to peptides primarily exhibiting antimicrobial activity and 25 potentially serving as AMPs in a secondary function, either in their entirety or after cleavage. In silico prediction methods were used to identify 62 of those AMPs as peptides with promising antimicrobial activity potential. Signatures of directional selection were detected for five candidate AMPs, which may indicate adaptation to the different selective pressures imposed by evolutionary arms races with specific pathogens. These findings provide encouraging support for the expectation that caecilians, being one of the least-studied groups of vertebrates, and with ~300 million years of separate evolution, are an underexplored resource of great pharmaceutical potential that could help to contest antibiotic resistance and contribute to biomedical advance.

Q586B2 is a crucial virulence factor during the early stages of Trypanosoma brucei infection that is conserved amongst trypanosomatids.

Human African trypanosomiasis or sleeping sickness, caused by the protozoan parasite Trypanosoma brucei, is characterized by the manipulation of the host's immune response to ensure parasite invasion and persistence. Uncovering key molecules that support parasite establishment is a prerequisite to interfere with this process. We identified Q586B2 as a T. brucei protein that induces IL-10 in myeloid cells, which promotes parasite infection invasiveness. Q586B2 is expressed during all T. brucei life stages and is conserved in all Trypanosomatidae. Deleting the Q586B2-encoding Tb927.6.4140 gene in T. brucei results in a decreased peak parasitemia and prolonged survival, without affecting parasite fitness in vitro, yet promoting short stumpy differentiation in vivo. Accordingly, neutralization of Q586B2 with newly generated nanobodies could hamper myeloid-derived IL-10 production and reduce parasitemia. In addition, immunization with Q586B2 delays mortality upon a challenge with various trypanosomes, including Trypanosoma cruzi. Collectively, we uncovered a conserved protein playing an important regulatory role in Trypanosomatid infection establishment.

Resistance Is Not Futile: Widespread Convergent Evolution of Resistance to Alpha-Neurotoxic Snake Venoms in Caecilians (Amphibia: Gymnophiona).

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

Odorant-binding proteins in canine anal sac glands indicate an evolutionarily conserved role in mammalian chemical communication.

BACKGROUND: Chemical communication is an important aspect of the behavioural ecology of a wide range of mammals. In dogs and other carnivores, anal sac glands are thought to convey information to conspecifics by secreting a pallet of small volatile molecules produced by symbiotic bacteria. Because these glands are unique to carnivores, it is unclear how their secretions relate to those of other placental mammals that make use of different tissues and secretions for chemical communication. Here we analyse the anal sac glands of domestic dogs to verify the secretion of proteins and infer their evolutionary relationship to those involved in the chemical communication of non-carnivoran mammals. RESULTS: Proteomic analysis of anal sac gland secretions of 17 dogs revealed the consistently abundant presence of three related proteins. Homology searches against online databases indicate that these proteins are evolutionary related to 'odorant binding proteins' (OBPs) found in a wide range of mammalian secretions and known to contribute to chemical communication. Screening of the dog's genome sequence show that the newly discovered OBPs are encoded by a single cluster of three genes in the pseudoautosomal region of the X-chromosome. Comparative genomic screening indicates that the same locus is shared by a wide range of placental mammals and that it originated at least before the radiation of extant placental orders. Phylogenetic analyses suggest a dynamic evolution of gene duplication and loss, resulting in large gene clusters in some placental taxa and recurrent loss of this locus in others. The homology of OBPs in canid anal sac glands and those found in other mammalian secretions implies that these proteins maintained a function in chemical communication throughout mammalian evolutionary history by multiple shifts in expression between secretory tissues involved in signal release and nasal mucosa involved in signal reception. CONCLUSIONS: Our study elucidates a poorly understood part of the biology of a species that lives in close association with humans. In addition, it shows that the protein repertoire underlying chemical communication in mammals is more evolutionarily stable than the variation of involved glands and tissues would suggest.

A New Family of Diverse Skin Peptides from the Microhylid Frog Genus Phrynomantis.

A wide range of frogs produce skin poisons composed of bioactive peptides for defence against pathogens, parasites and predators. While several frog families have been thoroughly screened for skin-secreted peptides, others, like the Microhylidae, have remained mostly unexplored. Previous studies of microhylids found no evidence of peptide secretion, suggesting that this defence adaptation was evolutionarily lost. We conducted transcriptome analyses of the skins of Phrynomantis bifasciatus and Phrynomantis microps, two African microhylid species long suspected to be poisonous. Our analyses reveal 17 evolutionary related transcripts that diversified from to those of cytolytic peptides found in other frog families. The 19 peptides predicted to be processed from these transcripts, named phrynomantins, show a striking structural diversity that is distinct from any previously identified frog skin peptide. Functional analyses of five phrynomantins confirm the loss of a cytolytic function and the absence of insecticidal or proinflammatory activity, suggesting that they represent an evolutionary transition to a new, yet unknown function. Our study shows that peptides have been retained in the defence poison of at least one microhylid lineage and encourages research on similarly understudied taxa to further elucidate the diversity and evolution of skin defence molecules.

Multiple Independent Recruitment of Sodefrin Precursor-Like Factors in Anuran Sexually Dimorphic Glands.

Chemical signaling in animals often plays a central role in eliciting a variety of responses during reproductive interactions between males and females. One of the best-known vertebrate courtship pheromone systems is sodefrin precursor-like factors (SPFs), a family of two-domain three-finger proteins with a female-receptivity enhancing function, currently only known from salamanders. The oldest divergence between active components in a single salamander species dates back to the Late Paleozoic, indicating that these proteins potentially gained a pheromone function earlier in amphibian evolution. Here, we combined whole transcriptome sequencing, proteomics, histology, and molecular phylogenetics in a comparative approach to investigate SPF occurrence in male breeding glands across the evolutionary tree of anurans (frogs and toads). Our study shows that multiple families of both terrestrially and aquatically reproducing frogs have substantially increased expression levels of SPFs in male breeding glands. This suggests that multiple anuran lineages make use of SPFs to complement acoustic and visual sexual signaling during courtship. Comparative analyses show that anurans independently recruited these proteins each time the gland location on the male's body allowed efficient transmission of the secretion to the female's nares.

The anuran skin peptide bradykinin mediates its own absorption across epithelial barriers of the digestive tract.

When faced with a potential predator, a wide range of frog species secrete a mixture of peptide toxins from their skin to defend themselves. We have recently shown that antimicrobial peptides (AMPs) in a frog's defensive poison enhance the uptake of these peptides across epithelia, thereby speeding up the process of predator intoxication. This study provides evidence that bradykinin, a widespread peptide toxin in anurans (frogs), is capable to pass through epithelial barriers independent of this delivery system. We quantified bradykinin peptides secreted by Bombina orientalis during acute stress, and found that at biologically relevant concentrations, bradykinin passage across model epithelia occurs even in the absence of AMPs. Monitoring of transepithelial electric resistance showed that bradykinin treatment caused a subtle yet prolonged reduction in barrier function, indicating that the peptide itself is capable to increase the permeability of epithelia. Yet, bradykinin does not cause cells to leak lactate dehydrogenase, suggesting that it does not damage cell membranes. Moreover, imaging of bradykinin-treated monolayers shows no endocytosis of fluorescent propidium iodide, indicating that the peptide does not perforate cell membranes at smaller scale and therefore is unlikely to cross epithelia via a transcellular passage. Together, these observations suggest that bradykinin, unlike other amphibian neuropeptide toxins, mediates its own passage across mucosal barriers, possibly through a paracellular route. This "self-administering" property, combined with the fact that bradykinins can potently disturb multiple physiological processes, could explain why these peptides are one of the most widespread antipredator peptides in the defensive secretions of frogs.

Antimicrobial peptides in frog poisons constitute a molecular toxin delivery system against predators.

Animals using toxic peptides and proteins for predation or defense typically depend on specialized morphological structures, like fangs, spines, or a stinger, for effective intoxication. Here we show that amphibian poisons instead incorporate their own molecular system for toxin delivery to attacking predators. Skin-secreted peptides, generally considered part of the amphibian immune system, permeabilize oral epithelial tissue and enable fast access of cosecreted toxins to the predator's bloodstream and organs. This absorption-enhancing system exists in at least three distantly related frog lineages and is likely to be a widespread adaptation, determining the outcome of predator-prey encounters in hundreds of species.

Reprint of: The non-mammalian MIF superfamily.

Macrophage migration inhibitory factor (MIF) was first described as a cytokine 50 years ago, and emerged in mammals as a pleiotropic protein with pro-inflammatory, chemotactic, and growth-promoting activities. In addition, MIF has gained substantial attention as a pivotal upstream mediator of innate and adaptive immune responses and with pathologic roles in several diseases. Of less importance in mammals is an intrinsic but non-physiologic enzymatic activity that points to MIF's evolution from an ancient defense molecule. Therefore, it is not surprising that mif-like genes also have been found across a range of different organisms including bacteria, plants, protozoa, helminths, molluscs, arthropods, fish, amphibians and birds. While Genebank analysis identifying mif-like genes across species is extensive, contained herein is an overview of the non-mammalian MIF-like proteins that have been most well studied experimentally. For many of these organisms, MIF contributes to an innate defense system or plays a role in development. For parasitic organisms however, MIF appears to function as a virulence factor aiding in the establishment or persistence of infection by modulating the host immune response. Consequently, a combined targeting of both parasitic and host MIF could lead to more effective treatment strategies for parasitic diseases of socioeconomic importance.

An orthologue of the host-defense protein psoriasin (S100A7) is expressed in frog skin.

Host-defense peptides and proteins are vital for first line protection against bacteria. Most host-defense peptides and proteins common in vertebrates have been studied primarily in mammals, while their orthologues in non-mammalian vertebrates received less attention. We found that the European Common Frog Rana temporaria expresses a protein in its skin that is evolutionarily related to the host-defense protein S100A7. This prompted us to test if the encoded protein, which is an important microbicidal protein in human skin, shows similar activity in frogs. The R. temporaria protein lacks the zinc-binding sites that are key to the antimicrobial activity of human S100A7 at neutral pH. However, despite being less potent, the R. temporaria protein does compromise bacterial membranes at low pH, similar to its human counterpart. We postulate that, while amphibian S100A7 likely serves other functions, the capacity to compromise bacterial cell membranes evolved early in tetrapod evolution.

The non-mammalian MIF superfamily.

Macrophage migration inhibitory factor (MIF) was first described as a cytokine 50 years ago, and emerged in mammals as a pleiotropic protein with pro-inflammatory, chemotactic, and growth-promoting activities. In addition, MIF has gained substantial attention as a pivotal upstream mediator of innate and adaptive immune responses and with pathologic roles in several diseases. Of less importance in mammals is an intrinsic but non-physiologic enzymatic activity that points to MIF's evolution from an ancient defense molecule. Therefore, it is not surprising that mif-like genes also have been found across a range of different organisms including bacteria, plants, protozoa, helminths, molluscs, arthropods, fish, amphibians and birds. While Genebank analysis identifying mif-like genes across species is extensive, contained herein is an overview of the non-mammalian MIF-like proteins that have been most well studied experimentally. For many of these organisms, MIF contributes to an innate defense system or plays a role in development. For parasitic organisms however, MIF appears to function as a virulence factor aiding in the establishment or persistence of infection by modulating the host immune response. Consequently, a combined targeting of both parasitic and host MIF could lead to more effective treatment strategies for parasitic diseases of socioeconomic importance.

Low structural variation in the host-defense peptide repertoire of the dwarf clawed frog Hymenochirus boettgeri (Pipidae).

THE skin secretion of many amphibians contains peptides that are able to kill a broad range of microorganisms (antimicrobial peptides: AMPs) and potentially play a role in innate immune defense. Similar to the toxin arsenals of various animals, amphibian AMP repertoires typically show major structural variation, and previous studies have suggested that this may be the result of diversifying selection in adaptation to a diverse spectrum of pathogens. Here we report on transcriptome analyses that indicate a very different pattern in the dwarf clawed frog H. boettgeri. Our analyses reveal a diverse set of transcripts containing two to six tandem repeats, together encoding 14 distinct peptides. Five of these have recently been identified as AMPs, while three more are shown here to potently inhibit the growth of gram-negative bacteria, including multi-drug resistant strains of the medically important Pseudomonas aeruginosa. Although the number of predicted peptides is similar to the numbers of related AMPs in Xenopus and Silurana frog species, they show significantly lower structural variation. Selection analyses confirm that, in contrast to the AMPs of other amphibians, the H. boettgeri peptides did not evolve under diversifying selection. Instead, the low sequence variation among tandem repeats resulted from purifying selection, recent duplication and/or concerted gene evolution. Our study demonstrates that defense peptide repertoires of closely related taxa, after diverging from each other, may evolve under differential selective regimes, leading to contrasting patterns of structural diversity.

Origin and functional diversification of an amphibian defense peptide arsenal.

The skin secretion of many amphibians contains an arsenal of bioactive molecules, including hormone-like peptides (HLPs) acting as defense toxins against predators, and antimicrobial peptides (AMPs) providing protection against infectious microorganisms. Several amphibian taxa seem to have independently acquired the genes to produce skin-secreted peptide arsenals, but it remains unknown how these originated from a non-defensive ancestral gene and evolved diverse defense functions against predators and pathogens. We conducted transcriptome, genome, peptidome and phylogenetic analyses to chart the full gene repertoire underlying the defense peptide arsenal of the frog Silurana tropicalis and reconstruct its evolutionary history. Our study uncovers a cluster of 13 transcriptionally active genes, together encoding up to 19 peptides, including diverse HLP homologues and AMPs. This gene cluster arose from a duplicated gastrointestinal hormone gene that attained a HLP-like defense function after major remodeling of its promoter region. Instead, new defense functions, including antimicrobial activity, arose by mutation of the precursor proteins, resulting in the proteolytic processing of secondary peptides alongside the original ones. Although gene duplication did not trigger functional innovation, it may have subsequently facilitated the convergent loss of the original function in multiple gene lineages (subfunctionalization), completing their transformation from HLP gene to AMP gene. The processing of multiple peptides from a single precursor entails a mechanism through which peptide-encoding genes may establish new functions without the need for gene duplication to avoid adaptive conflicts with older ones.

Squeezers and leaf-cutters: differential diversification and degeneration of the venom system in toxicoferan reptiles.

Although it has been established that all toxicoferan squamates share a common venomous ancestor, it has remained unclear whether the maxillary and mandibular venom glands are evolving on separate gene expression trajectories or if they remain under shared genetic control. We show that identical transcripts are simultaneously expressed not only in the mandibular and maxillary glands, but also in the enigmatic snake rictal gland. Toxin molecular frameworks recovered in this study were three-finger toxin (3FTx), CRiSP, crotamine (beta-defensin), cobra venom factor, cystatin, epididymal secretory protein, kunitz, L-amino acid oxidase, lectin, renin aspartate protease, veficolin, and vespryn. We also discovered a novel low-molecular weight disulfide bridged peptide class in pythonid snake glands. In the iguanian lizards, the most highly expressed are potentially antimicrobial in nature (crotamine (beta-defensin) and cystatin), with crotamine (beta-defensin) also the most diverse. However, a number of proteins characterized from anguimorph lizards and caenophidian snakes with hemotoxic or neurotoxic activities were recruited in the common toxicoferan ancestor and remain expressed, albeit in low levels, even in the iguanian lizards. In contrast, the henophidian snakes express 3FTx and lectin toxins as the dominant transcripts. Even in the constricting pythonid and boid snakes, where the glands are predominantly mucous-secreting, low-levels of toxin transcripts can be detected. Venom thus appears to play little role in feeding behavior of most iguanian lizards or the powerful constricting snakes, and the low levels of expression argue against a defensive role. However, clearly the incipient or secondarily atrophied venom systems of these taxa may be a source of novel compounds useful in drug design and discovery.

Alternative interactions define gyrase specificity in the CcdB family.

Toxin-antitoxin (TA) modules are small operons associated with stress response of bacteria. F-plasmid CcdB(F) was the first TA toxin for which its target, gyrase, was identified. Plasmidic and chromosomal CcdBs belong to distinct families. Conserved residues crucial for gyrase poisoning activity of plasmidic CcdBs are not conserved among these families. Here we show that the chromosomal CcdB(Vfi) from Vibrio fischeri is an active gyrase poison that interacts with its target via an alternative energetic mechanism. Changes in the GyrA14-binding surface of the Vibrio and F-plasmid CcdB family members illustrate neutral drift where alternative interactions can be used to achieve the same functionality. Differences in affinity between V. fischeri and F-plasmid CcdB for gyrase and their corresponding CcdA antitoxin possibly reflect distinct roles for TA modules located on plasmids and chromosomes.

Anuran radiations and the evolution of tadpole morphospace.

Anurans (frogs and toads) are unique among land vertebrates in possessing a free-living larval stage that, parallel to adult frogs, diversified into an impressive range of ecomorphs. The tempo and mode at which tadpole morphology evolved through anuran history as well as its relationship to lineage diversification remain elusive. We used a molecular phylogenetic framework to examine patterns of morphological evolution in tadpoles in light of observed episodes of accelerated lineage diversification. Our reconstructions show that the expansion of tadpole morphospace during the basal anuran radiation in the Triassic/Early Jurassic was unparalleled by the basal neobatrachian radiation in the Late Jurassic/Early Cretaceous or any subsequent radiation in the Late Cretaceous/Early Tertiary. Comparative analyses of radiation episodes indicate that the slowdown of morphospace expansion was caused not only by a drop in evolutionary rate after the basal anuran radiation but also by an overall increase in homoplasy in the characters that did evolve during later radiations. The overlapping sets of evolving characters among more recent radiations may have enhanced tadpole diversity by creating unique combinations of homoplastic traits, but the lack of innovative character changes prevented the exploration of fundamental regions in morphospace. These complex patterns transcend the four traditionally recognized tadpole morphotypes and apply to most tissue types and body parts.

Functional and structural diversification of the Anguimorpha lizard venom system.

Venom has only been recently discovered to be a basal trait of the Anguimorpha lizards. Consequently, very little is known about the timings of toxin recruitment events, venom protein molecular evolution, or even the relative physical diversifications of the venom system itself. A multidisciplinary approach was used to examine the evolution across the full taxonomical range of this ∼130 million-year-old clade. Analysis of cDNA libraries revealed complex venom transcriptomes. Most notably, three new cardioactive peptide toxin types were discovered (celestoxin, cholecystokinin, and YY peptides). The latter two represent additional examples of convergent use of genes in toxic arsenals, both having previously been documented as components of frog skin defensive chemical secretions. Two other novel venom gland-overexpressed modified versions of other protein frameworks were also recovered from the libraries (epididymal secretory protein and ribonuclease). Lectin, hyaluronidase, and veficolin toxin types were sequenced for the first time from lizard venoms and shown to be homologous to the snake venom forms. In contrast, phylogenetic analyses demonstrated that the lizard natriuretic peptide toxins were recruited independently of the form in snake venoms. The de novo evolution of helokinestatin peptide toxin encoding domains within the lizard venom natriuretic gene was revealed to be exclusive to the helodermatid/anguid subclade. New isoforms were sequenced for cysteine-rich secretory protein, kallikrein, and phospholipase A(2) toxins. Venom gland morphological analysis revealed extensive evolutionary tinkering. Anguid glands are characterized by thin capsules and mixed glands, serous at the bottom of the lobule and mucous toward the apex. Twice, independently this arrangement was segregated into specialized serous protein-secreting glands with thick capsules with the mucous lobules now distinct (Heloderma and the Lanthanotus/Varanus clade). The results obtained highlight the importance of utilizing evolution-based search strategies for biodiscovery and emphasize the largely untapped drug design and development potential of lizard venoms.

Gradual adaptation toward a range-expansion phenotype initiated the global radiation of toads.

Recent studies have identified range expansion as a potential driver of speciation. Yet it remains poorly understood how, under identical extrinsic settings, differential tendencies for geographic movement of taxa originate and subsequently affect diversification. We identified multiple traits that predict large distributional ranges in extant species of toads (Bufonidae) and used statistical methods to define and phylogenetically reconstruct an optimal range-expansion phenotype. Our results indicate that lineage-specific range-shifting abilities increased through an accumulation of adaptive traits that culminated in such a phenotype. This initiated the episode of global colonization and triggered the major radiation of toads. Evolution toward a range-expansion phenotype might be crucial to understanding both ancient widespread radiations and the evolutionary background of contemporary invasive species such as the cane toad.