Identical Skin Toxins by Convergent Molecular Adaptation in Frogs

Chromosome-scale assembly and quantitative trait locus mapping for major economic traits of the Culter alburnus genome using Illumina and PacBio sequencing with Hi-C mapping information

Topmouth culter (Culter alburnus) is an economically important freshwater fish with high nutritional value. However, its potential genetic advantages have not been fully exploited. Therefore, we aimed to determine the genome sequence of C. alburnus and examine quantitative trait loci (QTLs) related to major economic traits. The results showed that 24 pseudochromosomes were anchored by 914.74 Mb of the C. alburnus genome sequence. De novo sequencing identified 31,279 protein-coding genes with an average length of 8507 bp and average coding sequ ence of 1115 bp. In addition, a high-density genetic linkage map consisting of 24 linkage groups was constructed based on 353,532 high-quality single nucleotide polymorphisms and 4,710 bin markers. A total of 28 QTLs corresponding to 11 genes, 26 QTLs corresponding to 11 genes, and 12 QTLs corresponding to 5 genes were identified for sex, intermuscular spine number and body weight traits, respectively. In this study, we assembled an accurate and nearly complete genome of C. alburnus by combining Illumina, PacBio, and high-throughput Chromosome conformation capture (Hi-C) technologies. In addition, we identified QTLs that explained variances in intermuscular spine number, body weight, and sex differences in C. alburnus. These genetic markers or candidate genes associated with growth traits provide a basis for marker-assisted selection in C. alburnus.

Widespread convergence in stream fishes

Convergent evolution, the evolution of similar phenotypes among distantly related lineages, is often attributed to adaptation in response to similar selective pressures. Here, we assess the prevalence and degree of convergence in functional traits of stream fishes at the microhabitat scale in five zoogeographical regions across the world. We categorized species by microhabitat, water velocity and preference for substrate complexity and calculated the prevalence of convergence, degree of convergence and functional diversity for each category. Among species occupying similar microhabitats of small, low-gradient streams, 34% had combinations of convergent traits. Convergence occurred at higher rates than expected by chance alone, implying that adaptation to similar environmental conditions often resulted in similar evolutionary patterns along multiple niche dimensions. Two of the microhabitat groupings had significantly convergent species represented in all zoogeographical regions. Fishes occupying microhabitats with high water velocity and low structural complexity generally occupied a restricted morphospace and exhibited greater prevalence and higher degrees of convergence. This suggests that water velocity and habitat structural complexity interact, selecting a restricted distribution of trait distributions and higher degrees of convergence in stream fish assemblages. Furthermore, these results suggest that microhabitat features in streams select for fish trait distributions in a fairly predictable and deterministic manner worldwide.

Variation in Bufadienolide Composition of Parotoid Gland Secretion From Three Taxa of Japanese Toads

Toads of the genus Bufo synthesize and accumulate bufadienolides (BDs) in their parotoid glands. BDs are cardiotonic steroids that play an important role in defense against the toads' predators. Three bufonid taxa occur in mainland Japan, Bufo japonicus formosus, B. j. japonicus, and B. torrenticola. The chemical structures of BDs isolated from B. j. formosus were studied several decades ago, but there is no further information on the toxic components of Japanese toads and their metabolism. In this study, we analyzed BDs of toads from throughout Japan and compared the BD profiles by liquid chromatography/mass spectrometry (LC/MS) and hierarchical cluster analysis (HCA). We observed BDs in three taxa of Japanese toads, and identified five of the most common BDs by nuclear magnetic resonance (NMR) analyses. Of the five BDs, only bufalin was detected in all individuals. HCA of individual BD profiles divided the three taxa into five primary clusters and several subclusters. This result indicates that BD profiles differ both among and within the taxa. The clustering pattern of BDs is generally concordant with a phylogenetic tree reconstructed from the mitochondrial cytochrome b gene of Japanese toads. Our results suggest that the BDs of Japanese toads have diversified not in response to specific selective pressures, but simply due to population structuring over evolutionary time.

Insights into the skin of caecilian amphibians from gene expression profiles

Background

Gene expression profiles can provide insights into the molecular machinery behind tissue functions and, in turn, can further our understanding of environmental responses, and developmental and evolutionary processes. During vertebrate evolution, the skin has played a crucial role, displaying a wide diversity of essential functions. To unravel the molecular basis of skin specialisations and adaptations, we compared gene expression in the skin with eight other tissues in a phylogenetically and ecologically diverse species sample of one of the most neglected vertebrate groups, the caecilian amphibians (order Gymnophiona).

Results

The skin of the five studied caecilian species showed a distinct gene expression profile reflecting its developmental origin and showing similarities to other epithelial tissues. We identified 59 sequences with conserved enhanced expression in the skin that might be associated with caecilian dermal specialisations. Some of the up-regulated genes shared expression patterns with human skin and potentially are involved in skin functions across vertebrates. Variation trends in gene expression were detected between mid and posterior body skin suggesting different functions between body regions. Several candidate biologically active peptides were also annotated.

Conclusions

Our study provides the first atlas of differentially expressed sequences in caecilian tissues and a baseline to explore the molecular basis of the skin functions in caecilian amphibians, and more broadly in vertebrates.

An Economic Dilemma Between Molecular Weapon Systems May Explain an Arachno-atypical Venom in Wasp Spiders (Argiope bruennichi)

Spiders use venom to subdue their prey, but little is known about the diversity of venoms in different spider families. Given the limited data available for orb-weaver spiders (Araneidae), we selected the wasp spider Argiope bruennichi for detailed analysis. Our strategy combined a transcriptomics pipeline based on multiple assemblies with a dual proteomics workflow involving parallel mass spectrometry techniques and electrophoretic profiling. We found that the remarkably simple venom of A. bruennichi has an atypical composition compared to other spider venoms, prominently featuring members of the cysteine-rich secretory protein, antigen 5 and pathogenesis-related protein 1 (CAP) superfamily and other, mostly high-molecular-weight proteins. We also detected a subset of potentially novel toxins similar to neuropeptides. We discuss the potential function of these proteins in the context of the unique hunting behavior of wasp spiders, which rely mostly on silk to trap their prey. We propose that the simplicity of the venom evolved to solve an economic dilemma between two competing yet metabolically expensive weapon systems. This study emphasizes the importance of cutting-edge methods to encompass the lineages of smaller venomous species that have yet to be characterized in detail, allowing us to understand the biology of their venom systems and to mine this prolific resource for translational research.

The skin microbiome of elasmobranchs follows phylosymbiosis, but in teleost fishes, the microbiomes converge

BACKGROUND

The vertebrate clade diverged into Chondrichthyes (sharks, rays, and chimeras) and Osteichthyes fishes (bony fishes) approximately 420 mya, with each group accumulating vast anatomical and physiological differences, including skin properties. The skin of Chondrichthyes fishes is covered in dermal denticles, whereas Osteichthyes fishes are covered in scales and are mucous rich. The divergence time among these two fish groups is hypothesized to result in predictable variation among symbionts. Here, using shotgun metagenomics, we test if patterns of diversity in the skin surface microbiome across the two fish clades match predictions made by phylosymbiosis theory. We hypothesize (1) the skin microbiome will be host and clade-specific, (2) evolutionary difference in elasmobranch and teleost will correspond with a concomitant increase in host-microbiome dissimilarity, and (3) the skin structure of the two groups will affect the taxonomic and functional composition of the microbiomes.

RESULTS

We show that the taxonomic and functional composition of the microbiomes is host-specific. Teleost fish had lower average microbiome within clade similarity compared to among clade comparison, but their composition is not different among clade in a null based model. Elasmobranch's average similarity within clade was not different than across clade and not different in a null based model of comparison. In the comparison of host distance with microbiome distance, we found that the taxonomic composition of the microbiome was related to host distance for the elasmobranchs, but not the teleost fishes. In comparison, the gene function composition was not related to the host-organism distance for elasmobranchs but was negatively correlated with host distance for teleost fishes.

CONCLUSION

Our results show the patterns of phylosymbiosis are not consistent across both fish clades, with the elasmobranchs showing phylosymbiosis, while the teleost fish are not. The discrepancy may be linked to alternative processes underpinning microbiome assemblage, including possible historical host-microbiome evolution of the elasmobranchs and convergent evolution in the teleost which filter specific microbial groups. Our comparison of the microbiomes among fishes represents an investigation into the microbial relationships of the oldest divergence of extant vertebrate hosts and reveals that microbial relationships are not consistent across evolutionary timescales. Video abstract.

Cholecystokinin in the central nervous system of the sea lamprey Petromyzon marinus: precursor identification and neuroanatomical relationships with other neuronal signalling systems

Cholecystokinin (CCK) is a neuropeptide that modulates processes such as digestion, satiety, and anxiety. CCK-type peptides have been characterized in jawed vertebrates and invertebrates, but little is known about CCK-type signalling in the most ancient group of vertebrates, the agnathans. Here, we have cloned and sequenced a cDNA encoding a sea lamprey (Petromyzon marinus L.) CCK-type precursor (PmCCK), which contains a CCK-type octapeptide sequence (PmCCK-8) that is highly similar to gnathostome CCKs. Using mRNA in situ hybridization, the distribution of PmCCK-expressing neurons was mapped in the CNS of P. marinus. This revealed PmCCK-expressing neurons in the hypothalamus, posterior tubercle, prethalamus, nucleus of the medial longitudinal fasciculus, midbrain tegmentum, isthmus, rhombencephalic reticular formation, and the putative nucleus of the solitary tract. Some PmCCK-expressing neuronal populations were only observed in adults, revealing important differences with larvae. We generated an antiserum to PmCCK-8 to enable immunohistochemical analysis of CCK expression, which revealed that GABA or glutamate, but not serotonin, tyrosine hydroxylase or neuropeptide Y, is co-expressed in some PmCCK-8-immunoreactive (ir) neurons. Importantly, this is the first demonstration of co-localization of GABA and CCK in neurons of a non-mammalian vertebrate. We also characterized extensive cholecystokinergic fibre systems of the CNS, including innervation of habenular subnuclei. A conspicuous PmCCK-8-ir tract ascending in the lateral rhombencephalon selectively innervates a glutamatergic population in the dorsal isthmic grey. Interestingly, this tract is reminiscent of the secondary gustatory/visceral tract of teleosts. In conclusion, this study provides important new information on the evolution of the cholecystokinergic system in vertebrates.

A Hylarana latouchii Skin Secretion-Derived Novel Bombesin-Related Pentadecapeptide (Ranatensin-HLa) Evoke Myotropic Effects on the in vitro Rat Smooth Muscles

Amphibians have developed successful defensive strategies for combating predators and invasive microorganisms encountered in their broad range of environments, which involve secretion of complex cocktails of noxious, toxic and diverse bioactive molecules from the skins. In recent years, amphibian skin secretions have been considered as an extraordinary warehouse for the discovery of therapeutic medicines. In this study, through bioactivity screening of the Hylarana latouchii skin secretion-derived fractions, a novel peptide belonging to ranatensin subfamily (ranatensin-HLa) was discovered, and structurally and pharmacologically-characterised. It consists of 15 amino acid residues, pGlu-NGDRAPQWAVGHFM-NH₂, and its synthetic replicate was found to exhibit pharmacological activities on increasing the contraction of the in vitro rat bladder and uterus smooth muscles. Corresponding characteristic sigmoidal dose-response curves with EC₅₀ values of 7.1 nM and 5.5 nM were produced, respectively, in bladder and uterus. Moreover, the precursor of ranatensin-HLa showed a high degree of similarity to those of bombesin-like peptides from Odorrana grahami and Odorrana schmackeri. Hylarana latouchii skin continues to serve as a storehouse with diverse lead compounds for the development of therapeutically effective medicines.

Multi-tissue transcriptomes of caecilian amphibians highlight incomplete knowledge of vertebrate gene families

RNA sequencing (RNA-seq) has become one of the most powerful tools to unravel the genomic basis of biological adaptation and diversity. Although challenging, RNA-seq is particularly promising for research on non-model, secretive species that cannot be observed in nature easily and therefore remain comparatively understudied. Among such animals, the caecilians (order Gymnophiona) likely constitute the least known group of vertebrates, despite being an old and remarkably distinct lineage of amphibians. Here, we characterize multi-tissue transcriptomes for five species of caecilians that represent a broad level of diversity across the order. We identified vertebrate homologous elements of caecilian functional genes of varying tissue specificity that reveal a great number of unclassified gene families, especially for the skin. We annotated several protein domains for those unknown candidate gene families to investigate their function. We also conducted supertree analyses of a phylogenomic dataset of 1,955 candidate orthologous genes among five caecilian species and other major lineages of vertebrates, with the inferred tree being in agreement with current views of vertebrate evolution and systematics. Our study provides insights into the evolution of vertebrate protein-coding genes, and a basis for future research on the molecular elements underlying the particular biology and adaptations of caecilian amphibians.

Independent evolution of complex development in animals and plants: deep homology and lateral gene transfer

The evolution of multicellularity is a premier example of phenotypic convergence: simple multicellularity evolved independently many times, and complex multicellular phenotypes are found in several distant groups. Furthermore, both animal and plant lineages have independently reached extreme levels of morphological, functional, and developmental complexity. This study explores the genetic basis for the parallel evolution of complex multicellularity and development in the animal and green plant (i.e., green algae and land plants) lineages. Specifically, the study (i) identifies the SAND domain-a DNA-binding domain with important roles in the regulation of cell proliferation and differentiation, as unique to animals, green algae, and land plants; and (ii) suggests that the parallel deployment of this ancestral domain in similar regulatory roles could have contributed to the independent evolution of complex development in these distant groups. Given the deep animal-green plant divergence, the limited distribution of the SAND domain is best explained by invoking a lateral gene transfer (LGT) event from a green alga to an early metazoan. The presence of a sequence motif specifically shared by a family of SAND-containing transcription factors involved in the evolution of complex multicellularity in volvocine algae and two types of SAND proteins that emerged early in the evolution of animals is consistent with this scenario. Overall, these findings imply that (i) in addition to be involved in the evolution of similar phenotypes, deep homologous sequences can also contribute to shaping parallel evolutionary trajectories in distant lineages, and (ii) LGT could provide an additional source of latent homologous sequences that can be deployed in analogous roles and affect the evolutionary potentials of distantly related groups.

A death in the family: Sea lamprey (Petromyzon marinus) avoidance of confamilial alarm cues diminishes with phylogenetic distance

Alarm signals released after predator attack function as reliable public information revealing areas of high risk. The utility of this information can extend beyond species boundaries, benefiting heterospecifics capable of recognizing and responding appropriately to the signal. Nonmutually exclusive hypotheses explaining the acquisition of heterospecific reactivity to cues suggest it could be conserved phylogenetically following its evolution in a common ancestor (a species‐level effect) and/or learned during periods of shared risk (a population‐level effect; e.g., shared predators). Using a laboratory‐based space‐use behavioral assay, we tested the response of sea lamprey (Petromyzon marinus) to the damage‐released alarm cues of five confamilial (sympatric and allopatric) species and two distantly related out‐groups: a sympatric teleost (white sucker Catostomus commersonii) and an allopatric agnathan (Atlantic hagfish Myxine glutinosa). We found that sea lamprey differed in their response to conspecific and heterospecific odors; exhibiting progressively weaker avoidance of cues derived from more phylogenetically distant confamilials regardless of current overlap in distribution. Odors from out‐groups elicited no response. These findings suggest that a damage‐released alarm cue is at least partially conserved within the Petromyzontidae and that sea lamprey perceives predator attacks directed to closely related taxa. These findings are consistent with similar observations from gastropod, amphibian and bony fish taxa, and we discuss this in an eco‐evo context to provide a plausible explanation for the acquisition and maintenance of the response in sea lamprey.

Mesobuthus Venom-Derived Antimicrobial Peptides Possess Intrinsic Multifunctionality and Differential Potential as Drugs

Animal venoms are a mixture of peptides and proteins that serve two basic biological functions: predation and defense against both predators and microbes. Antimicrobial peptides (AMPs) are a common component extensively present in various scorpion venoms (herein abbreviated as svAMPs). However, their roles in predation and defense against predators and potential as drugs are poorly understood. Here, we report five new venom peptides with antimicrobial activity from two Mesobuthus scorpion species. These α-helical linear peptides displayed highly bactericidal activity toward all the Gram-positive bacteria used here but differential activity against Gram-negative bacteria and fungi. In addition to the antibiotic activity, these AMPs displayed lethality to houseflies and hemotoxin-like toxicity on mice by causing hemolysis, tissue damage and inducing inflammatory pain. Unlike AMPs from other origins, these venom-derived AMPs seem to be unsuitable as anti-infective drugs due to their high hemolysis and low serum stability. However, MeuTXKβ1, a known two-domain Mesobuthus AMP, is an exception since it exhibits high activity toward antibiotic resistant Staphylococci clinical isolates with low hemolysis and high serum stability. The findings that the classical AMPs play predatory and defensive roles indicate that the multifunctionality of scorpion venom components is an intrinsic feature likely evolved by natural selection from microbes, prey and predators of scorpions. This definitely provides an excellent system in which one can study how a protein adaptively evolves novel functions in a new environment. Meantimes, new strategies are needed to remove the toxicity of svAMPs on eukaryotic cells when they are used as leads for anti-infective drugs.

Enter the Dragon: The Dynamic and Multifunctional Evolution of Anguimorpha Lizard Venoms

While snake venoms have been the subject of intense study, comparatively little work has been done on lizard venoms. In this study, we have examined the structural and functional diversification of anguimorph lizard venoms and associated toxins, and related these results to dentition and predatory ecology. Venom composition was shown to be highly variable across the 20 species of Heloderma, Lanthanotus, and Varanus included in our study. While kallikrein enzymes were ubiquitous, they were also a particularly multifunctional toxin type, with differential activities on enzyme substrates and also ability to degrade alpha or beta chains of fibrinogen that reflects structural variability. Examination of other toxin types also revealed similar variability in their presence and activity levels. The high level of venom chemistry variation in varanid lizards compared to that of helodermatid lizards suggests that venom may be subject to different selection pressures in these two families. These results not only contribute to our understanding of venom evolution but also reveal anguimorph lizard venoms to be rich sources of novel bioactive molecules with potential as drug design and development lead compounds.

Physicochemical and biological characterizations of Pxt peptides from amphibian (Xenopus tropicalis) skin

Pxt peptides (Pxt-1 through Pxt-12) have been isolated from amphibian, Xenopus tropicalis Pxt-related peptides (Pxt-2, Pxt-5, Pxt-12, reverse Pxt-2, reverse Pxt-5 and reverse Pxt-12) with significant foaming properties were further characterized. In the physicochemical experiments, all Pxt-related peptides formed significant amphiphilic α-helices in 50% 2,2,2-trifluoroethanol by circular dichroism measurements. Among Pxt-related peptides, both Pxt-5 and reverse Pxt-5 were the most effective in reducing their surface tensions. Moreover, Pxt-2, Pxt-5 and reverse Pxt-5 produced constant surface tensions above their critical association concentrations, suggesting the micelle-like assemblies. In the biological experiments, Pxt-5 possessed the most potent hemolytic activity, while reverse Pxt-5 exhibited the most remarkable gene expression of interleukin 8 and heme oxygenase 1 and the most potent cytotoxicity in HaCaT cells. In contrast, Pxt-12 and reverse Pxt-12 were much weaker in antimicrobial assays for Gram-negative bacteria, Gram-positive bacteria and yeasts, as well as in hemolytic, cell viability and cytotoxicity assays in HaCaT cells. All Pxt-related peptides exhibited about 20-50% of the total cellular histamine release at 10(-5) M, as well as mastoparan and melittin in mast cells. Real-time polymerase chain reaction analysis confirmed the gene expressions of Pxt-5 in testis and Pxt-12 in muscle, in addition to skin, while Pxt-2 was only in skin.

The diversity and evolution of anuran skin peptides

Amphibians exhibit various, characteristic adaptations related to their "incomplete" shift from the aquatic to the terrestrial habitat. In particular, the integument was subject to a number of specialized modifications during the evolution of these animals. In this review, we place special emphasis on endogenous host-defence skin peptides from the cuteanous granular glands anuran amphibians (frogs and toads). The overview on the two broad groups of neuroactive and antimicrobial peptides (AMPs) goes beyond a simple itemization in that we provide a new perspective into the evolution and function of anuran AMPs. Briefly, these cationic, amphipathic and α-helical peptides are traditionally viewed as being part of the innate immune system, protecting the moist skin against invading microorganisms through their cytolytic action. However, the complete record of anuran species investigated to date suggests that AMPs are distributed sporadically (i.e., non-universally) across Anura. Together with the intriguing observation that virtually all anurans known to produce neuropeptides in their granular glands also co-secrete cytolytic peptides, we call the traditional role for AMPs as being purely antimicrobial into question and present an alternative scenario. We hypothesize AMPs to assist neuroactive peptides in their antipredator role through their cytolytic action increasing the delivery of the latter to the endocrine and nervous system of the predator. Thus, AMPs are more accurately viewed as cytolysins and their contribution to the immune system is better regarded as an accessory benefit.

Identification of novel peptides from amphibian (Xenopus tropicalis) skin by direct tissue MALDI-MS analysis

Twelve novel peptides (Pxt-1 to Pxt-12) were isolated from the skin of Xenopus tropicalis, diploid frogs, using topological MS analysis. Among them, Pxt-8, Pxt-9, and Pxt-10 were the N terminus of Pxt-1, N terminus of Pxt-3 and C terminus of Pxt-11, respectively. The Pxt-3 and Pxt-11 peptides shared significant sequence homologies with magainins 1, -2 and levitide, respectively, which all isolated from X. laevis. Pxt-12 was identical to the X. tropicalis XT-6-like precursor previously isolated by ESI-MS/MS. None of the Pxt peptides contained any Cys, Asp, Tyr or Trp, although Leu and Lys were frequently found as typical frog-skin peptides. RT-PCR analysis confirmed the gene expressions of Pxt-2, Pxt-3, Pxt-4, Pxt-5, Pxt-7 and Pxt-11 in X. tropicalis skin. Several ion peaks corresponding to all identified Pxt peptides were observed with MALDI-MS analysis of X. tropicalis secretory fluids, collected after in vivo stimulation, which suggested that Pxt peptides were definitely secretory molecules. CD studies and Schiffer-Edmundson helical wheel projections suggested that Pxt-5, as well as mastoparan, at least, could form a typical amphiphilic α helix without a phospholipid or a membrane-mimetic solvent (trifluoroethanol). Moreover, Pxt-2 and Pxt-5 showed growth inhibitory effects on both Escherichia coli (Gram-negative) and Staphylococcus aureus (Gram-positive). Measurements of dynamic light scattering and the surface tensions of Pxt peptides solutions suggested that both Pxt-2 and Pxt-5 could form associations as micelles and behave like a general surfactant. Moreover, the remarkable foaming properties of Pxt-2 and Pxt-5 were observed, as well as those of the secretory fluids of X. tropicalis.

Identification and localization of gastrointestinal hormones in the skin of the bullfrog Rana catesbeiana during periods of activity and hibernation

Amphibian skin and its secretions contain a wide variety of biogenic amines and biologically active peptides, some of which are either identical or highly homologous to gastrointestinal hormones (GHs) of higher vertebrates. This study investigated the distribution density and immunoreactive (IR) intensity of 5-hydroxytryptamine (5-HT), gastrin (GAS), somatostatin (SS), pancreatic polypeptide (PP), neuropeptide Y (NPY) and glucagon (GLU) IR cells in the skin of the bullfrog Rana catesbeiana during periods of activity and hibernation. The results indicated that the six types of GHs were all present in the bullfrog skin and were most predominant in the epidermis and mucous glands. In dorsal skin, the density of the GHs-IR cells in mucous glands was higher than that in epidermis except for GAS-IR cells. In ventral skin, the density of 5-HT, PP and NPY-IR cells in mucous glands was also higher than that in the epidermis. During hibernation, the density of the six types of GHs-IR cells and the IR intensity of GAS, SS, NPY and GLU-IR cells in the epidermis of dorsal skin increased significantly. The IR intensity of SS, PP and NPY-IR cells in granular glands of ventral skin also increased significantly during hibernation. These results suggested that multiple types of GHs-IR cells present in the skin of R. catesbeiana, may play important roles in the regulation of the physiological functions of skin. Also, adaptive changes in the density and IR intensity of GHs-IR cells occurred during hibernation.

Latent homology and convergent regulatory evolution underlies the repeated emergence of yeasts

Convergent evolution is common throughout the tree of life, but the molecular mechanisms causing similar phenotypes to appear repeatedly are obscure. Yeasts have arisen in multiple fungal clades, but the genetic causes and consequences of their evolutionary origins are unknown. Here we show that the potential to develop yeast forms arose early in fungal evolution and became dominant independently in multiple clades, most likely via parallel diversification of Zn-cluster transcription factors, a fungal-specific family involved in regulating yeast-filamentous switches. Our results imply that convergent evolution can happen by the repeated deployment of a conserved genetic toolkit for the same function in distinct clades via regulatory evolution. We suggest that this mechanism might be a common source of evolutionary convergence even at large time scales.

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.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.

Parallel signatures of sequence evolution among hearing genes in echolocating mammals: an emerging model of genetic convergence

Recent findings of sequence convergence in the Prestin gene among some bats and cetaceans suggest that parallel adaptations for high-frequency hearing have taken place during the evolution of echolocation. To determine if this gene is an exception, or instead similar processes have occurred in other hearing genes, we have examined Tmc1 and Pjvk, both of which are associated with non-syndromic hearing loss in mammals. These genes were amplified and sequenced from a number of mammalian species, including echolocating and non-echolocating bats and whales, and were analysed together with published sequences. Sections of both genes showed phylogenetic signals that conflicted with accepted species relationships, with coding regions uniting laryngeal echolocating bats in a monophyletic clade. Bayesian estimates of posterior probabilities of convergent and divergent substitutions provided more direct evidence of sequence convergence between the two groups of laryngeal echolocating bats as well as between echolocating bats and dolphins. We found strong evidence of positive selection acting on some echolocating bat species and echolocating cetaceans, contrasting with purifying selection on non-echolocating bats. Signatures of sequence convergence and molecular adaptation in two additional hearing genes suggest that the acquisition of high-frequency hearing has involved multiple loci.

A limited role for gene duplications in the evolution of platypus venom

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the "venome" of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

Caerulein-and xenopsin-related peptides with insulin-releasing activities from skin secretions of the clawed frogs, Xenopus borealis and Xenopus amieti (Pipidae)

Caerulein-related peptides were identified in norepinephrine-stimulated skin secretions of the tetraploid frog Xenopus borealis and the octoploid frog Xenopus amieti using negative ion electrospray mass spectrometry and their primary structures determined by positive ion tandem (MS/MS) mass spectrometry. X. borealis caerulein-B1 (pGlu-Gln-Asp-Tyr(SO(3))-Gly-Thr-Gly-Trp-Met-Asp-Phe.NH2) contains an additional Gly(5) residue compared with X. laevis caerulein and caerulein-B2 (pGlu-Asp-Tyr(SO(3))-Thr-Gly-Trp-Met-Asp-Phe.NH2) contains a Gln(2) deletion. X. amieti caerulein was identical to the X. laevis peptide. In addition, xenopsin, identical to the peptide from X. laevis, together with xenopsin-AM2 (pGlu-Gly-Arg-Arg-Pro-Trp-Ile- Leu) that contains the substitution Lys(3)→Arg were isolated from X. amieti secretions. X. borealis caerulein-B1, and X. amieti xenopsin and xenopsin-AM2 produced significant (P<0.05) and concentration-dependent stimulations of insulin release from the rat BRIN-BD11 clonal β cell line at concentrations ⩾30nM. The peptides did not stimulate the release of lactate dehydrogenase at concentrations up to 3μM demonstrating that the integrity of the plasma membrane had been preserved. While their precise biological role is unclear, the caerulein- and xenopsin-related peptides may constitute a component of the animal's chemical defenses against predators.

Convergence, adaptation, and constraint

Convergent evolution of similar phenotypic features in similar environmental contexts has long been taken as evidence of adaptation. Nonetheless, recent conceptual and empirical developments in many fields have led to a proliferation of ideas about the relationship between convergence and adaptation. Despite criticism from some systematically minded biologists, I reaffirm that convergence in taxa occupying similar selective environments often is the result of natural selection. However, convergent evolution of a trait in a particular environment can occur for reasons other than selection on that trait in that environment, and species can respond to similar selective pressures by evolving nonconvergent adaptations. For these reasons, studies of convergence should be coupled with other methods-such as direct measurements of selection or investigations of the functional correlates of trait evolution-to test hypotheses of adaptation. The independent acquisition of similar phenotypes by the same genetic or developmental pathway has been suggested as evidence of constraints on adaptation, a view widely repeated as genomic studies have documented phenotypic convergence resulting from change in the same genes, sometimes even by the same mutation. Contrary to some claims, convergence by changes in the same genes is not necessarily evidence of constraint, but rather suggests hypotheses that can test the relative roles of constraint and selection in directing phenotypic evolution.