Convergent evolution of novel protein function in shrew and lizard venom

Functional convergence in gastric lysozymes of foregut-fermenting rodents, ruminants, and primates is not attributed to convergent molecular evolution

Convergent evolution is a widespread phenomenon. While there are many examples of convergent evolution at the phenotypic scale, convergence at the molecular level has been more difficult to identify. A classic example of convergent evolution across scales is that of the digestive lysozyme found in ruminants and Colobine monkeys. These herbivorous species rely on foregut fermentation, which has evolved to function more optimally under acidic conditions. Here, we explored if rodents with similar dietary strategies and digestive morphologies have convergently evolved a lysozyme with digestive functions. At the phenotypic level, we find that rodents with bilocular stomach morphologies exhibited a lysozyme that maintained higher relative activities at low pH values, similar to the lysozymes of ruminants and Colobine monkeys. Additionally, the lysozyme of Peromyscus leucopus shared a similar predicted protonation state as that observed in previously identified digestive lysozymes. However, we found limited evidence of positive selection acting on the lysozyme gene in foregut-fermenting species and did not identify patterns of convergent molecular evolution in this gene. This study emphasizes that phenotypic convergence need not be the result of convergent genetic modifications, and we encourage further exploration into the mechanisms regulating convergence across biological scales.

Mammalian neurotoxins, Blarina paralytic peptides, cause hyperpolarization of human T-type Ca channel hCav3.2 activation

Among the rare venomous mammals, the short-tailed shrew Blarina brevicauda has been suggested to produce potent neurotoxins in its saliva to effectively capture prey. Several kallikrein-like lethal proteases have been identified, but the active substances of B. brevicauda remained unclear. Here, we report Blarina paralytic peptides (BPPs) 1 and 2 isolated from its submaxillary glands. Synthetic BPP2 showed mealworm paralysis and a hyperpolarization shift (-11 mV) of a human T-type Ca2+ channel (hCav3.2) activation. The amino acid sequences of BPPs were similar to those of synenkephalins, which are precursors of brain opioid peptide hormones that are highly conserved among mammals. However, BPPs rather resembled centipede neurotoxic peptides SLPTXs in terms of disulfide bond connectivity and stereostructure. Our results suggested that the neurotoxin BPPs were the result of convergent evolution as homologs of nontoxic endogenous peptides that are widely conserved in mammals. This finding is of great interest from the viewpoint of the chemical evolution of vertebrate venoms.

Systems biology as a framework to understand the physiological and endocrine bases of behavior and its evolution-From concepts to a case study in birds

Organismal behavior, with its tremendous complexity and diversity, is generated by numerous physiological systems acting in coordination. Understanding how these systems evolve to support differences in behavior within and among species is a longstanding goal in biology that has captured the imagination of researchers who work on a multitude of taxa, including humans. Of particular importance are the physiological determinants of behavioral evolution, which are sometimes overlooked because we lack a robust conceptual framework to study mechanisms underlying adaptation and diversification of behavior. Here, we discuss a framework for such an analysis that applies a "systems view" to our understanding of behavioral control. This approach involves linking separate models that consider behavior and physiology as their own networks into a singular vertically integrated behavioral control system. In doing so, hormones commonly stand out as the links, or edges, among nodes within this system. To ground our discussion, we focus on studies of manakins (Pipridae), a family of Neotropical birds. These species have numerous physiological and endocrine specializations that support their elaborate reproductive displays. As a result, manakins provide a useful example to help imagine and visualize the way systems concepts can inform our appreciation of behavioral evolution. In particular, manakins help clarify how connectedness among physiological systems-which is maintained through endocrine signaling-potentiate and/or constrain the evolution of complex behavior to yield behavioral differences across taxa. Ultimately, we hope this review will continue to stimulate thought, discussion, and the emergence of research focused on integrated phenotypes in behavioral ecology and endocrinology.

Mutation-driven parallel evolution in emergence of ACE2-utilizing sarbecoviruses

Mutation and recombination are two major genetic mechanisms that drive the evolution of viruses. They both exert an interplay during virus evolution, in which mutations provide a first ancestral source of genetic diversity for subsequent recombination. Sarbecoviruses are a group of evolutionarily related β-coronaviruses including human severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 and a trove of related animal viruses called SARS-like CoVs (SL-CoVs). This group of members either use or not use angiotensin-converting enzyme 2 (ACE2) as their entry receptor, which has been linked to the properties of their spike protein receptor binding domains (RBDs). This raises an outstanding question regarding how ACE2 binding originated within sarbecoviruses. Using a combination of analyses of phylogenies, ancestral sequences, structures, functions and molecular dynamics, we provide evidence in favor of an evolutionary scenario, in which three distinct ancestral RBDs independently developed the ACE2 binding trait via parallel amino acid mutations. In this process, evolutionary intermediate RBDs might be firstly formed through loop extensions to offer key functional residues accompanying point mutations to remove energetically unfavorable interactions and to change the dynamics of the functional loops, all required for ACE2 binding. Subsequent optimization in the context of evolutionary intermediates led to the independent emergence of ACE2-binding RBDs in the SARS-CoV and SARS-CoV-2 clades of Asian origin and the clade comprising SL-CoVs of European and African descent. These findings will help enhance our understanding of mutation-driven evolution of sarbecoviruses in their early history.

The Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

Co-option of the same ancestral gene family gave rise to mammalian and reptilian toxins

Background

Evolution can occur with surprising predictability when organisms face similar ecological challenges. For most traits, it is difficult to ascertain whether this occurs due to constraints imposed by the number of possible phenotypic solutions or because of parallel responses by shared genetic and regulatory architecture. Exceptionally, oral venoms are a tractable model of trait evolution, being largely composed of proteinaceous toxins that have evolved in many tetrapods, ranging from reptiles to mammals. Given the diversity of venomous lineages, they are believed to have evolved convergently, even though biochemically similar toxins occur in all taxa.

Results

Here, we investigate whether ancestral genes harbouring similar biochemical activity may have primed venom evolution, focusing on the origins of kallikrein-like serine proteases that form the core of most vertebrate oral venoms. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Additionally, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are genetically distinct from non-toxin ones.

Conclusions

Given the shared regulatory and genetic machinery, these findings suggest that tetrapod venoms evolved by co-option of proteins that were likely already present in saliva. We term such genes ‘toxipotent’—in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01191-1.

Venom Use in Eulipotyphlans: An Evolutionary and Ecological Approach

Venomousness is a complex functional trait that has evolved independently many times in the animal kingdom, although it is rare among mammals. Intriguingly, most venomous mammal species belong to Eulipotyphla (solenodons, shrews). This fact may be linked to their high metabolic rate and a nearly continuous demand of nutritious food, and thus it relates the venom functions to facilitation of their efficient foraging. While mammalian venoms have been investigated using biochemical and molecular assays, studies of their ecological functions have been neglected for a long time. Therefore, we provide here an overview of what is currently known about eulipotyphlan venoms, followed by a discussion of how these venoms might have evolved under ecological pressures related to food acquisition, ecological interactions, and defense and protection. We delineate six mutually nonexclusive functions of venom (prey hunting, food hoarding, food digestion, reducing intra- and interspecific conflicts, avoidance of predation risk, weapons in intraspecific competition) and a number of different subfunctions for eulipotyphlans, among which some are so far only hypothetical while others have some empirical confirmation. The functions resulting from the need for food acquisition seem to be the most important for solenodons and especially for shrews. We also present several hypotheses explaining why, despite so many potentially beneficial functions, venomousness is rare even among eulipotyphlans. The tentativeness of many of the arguments presented in this review highlights our main conclusion, i.e., insights regarding the functions of eulipotyphlan venoms merit additional study.

A Comprehensive Multi-Omic Approach Reveals a Relatively Simple Venom in a Diet Generalist, the Northern Short-Tailed Shrew, Blarina brevicauda

Animals that use venom to feed on a wide diversity of prey may evolve a complex mixture of toxins to target a variety of physiological processes and prey-defense mechanisms. Blarina brevicauda, the northern short-tailed shrew, is one of few venomous mammals, and is also known to eat evolutionarily divergent prey. Despite their complex diet, earlier proteomic and transcriptomic studies of this shrew's venom have only identified two venom proteins. Here, we investigated with comprehensive molecular approaches whether B. brevicauda venom is more complex than previously understood. We generated de novo assemblies of a B. brevicauda genome and submaxillary-gland transcriptome, as well as sequenced the salivary proteome. Our findings show that B. brevicauda's venom composition is simple relative to their broad diet and is likely limited to seven proteins from six gene families. Additionally, we explored expression levels and rate of evolution of these venom genes and the origins of key duplications that led to toxin neofunctionalization. We also found three proteins that may be involved in endogenous self-defense. The possible synergism of the toxins suggests that vertebrate prey may be the main target of the venom. Further functional assays for all venom proteins on both vertebrate and invertebrate prey would provide further insight into the ecological relevance of venom in this species.

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

Surface loops of trypsin-like serine proteases as determinants of function

Trypsin and chymotrypsin-like serine proteases from family S1 (clan PA) constitute the largest protease group in humans and more generally in vertebrates. The prototypes chymotrypsin, trypsin and elastase represent simple digestive proteases in the gut, where they cleave nearly any protein. Multidomain trypsin-like proteases are key players in the tightly controlled blood coagulation and complement systems, as well as related proteases that are secreted from diverse immune cells. Some serine proteases are expressed in nearly all tissues and fluids of the human body, such as the human kallikreins and kallikrein-related peptidases with specialization for often unique substrates and accurate timing of activity. HtrA and membrane-anchored serine proteases fulfill important physiological tasks with emerging roles in cancer. The high diversity of all family members, which share the tandem β-barrel architecture of the chymotrypsin-fold in the catalytic domain, is conferred by the large differences of eight surface loops, surrounding the active site. The length of these loops alters with insertions and deletions, resulting in remarkably different three-dimensional arrangements. In addition, metal binding sites for Na+, Ca2+ and Zn2+ serve as regulatory elements, as do N-glycosylation sites. Depending on the individual tasks of the protease, the surface loops determine substrate specificity, control the turnover and allow regulation of activation, activity and degradation by other proteins, which are often serine proteases themselves. Most intriguingly, in some serine proteases, the surface loops interact as allosteric network, partially tuned by protein co-factors. Knowledge of these subtle and complicated molecular motions may allow nowadays for new and specific pharmaceutical or medical approaches.

Arthropod venoms: Biochemistry, ecology and evolution

Comprising of over a million described species of highly diverse invertebrates, Arthropoda is amongst the most successful animal lineages to have colonized aerial, terrestrial, and aquatic domains. Venom, one of the many fascinating traits to have evolved in various members of this phylum, has underpinned their adaptation to diverse habitats. Over millions of years of evolution, arthropods have evolved ingenious ways of delivering venom in their targets for self-defence and predation. The morphological diversity of venom delivery apparatus in arthropods is astounding, and includes extensively modified pedipalps, tail (telson), mouth parts (hypostome), fangs, appendages (maxillulae), proboscis, ovipositor (stinger), and hair (urticating bristles). Recent investigations have also unravelled an astonishing venom biocomplexity with molecular scaffolds being recruited from a multitude of protein families. Venoms are a remarkable bioresource for discovering lead compounds in targeted therapeutics. Several components with prospective applications in the development of advanced lifesaving drugs and environment friendly bio-insecticides have been discovered from arthropod venoms. Despite these fascinating features, the composition, bioactivity, and molecular evolution of venom in several arthropod lineages remains largely understudied. This review highlights the prevalence of venom, its mode of toxic action, and the evolutionary dynamics of venom in Arthropoda, the most speciose phylum in the animal kingdom.

Characterization of bi-domain drosomycin-type antifungal peptides in nematodes: An example of convergent evolution

Drosomycin-type antifungal peptides (DTAFPs) are natural effectors of the innate immune system, which are restrictedly distributed in plants and ecdysozoans. Mehamycin is a bi-domain DTAFP (abbreviated as bDTAFP) firstly found in the Northern root-knot nematode Meloidogyne hapla. Here, we report its structural and functional features and the evolution of bDTAFPs in nematodes. Different from classical DTAFPs, mehamycin contains an insertion, called single Disulfide Bridge-linked Domain (abbreviated as sDBD), located in a loop region of the drosomycin scaffold. Despite this, recombinant mehamycin likely adopts a similar fold to drosomycin, as revealed by the circular dichroism spectral analysis. Functionally, it showed some weak activity against three species of fungi but relatively stronger activity against seven species of Gram-positive bacteria, indicative of functional diversification between mehamycin and classical DTAFPs. By computational data mining of the nematode databases, we identified polymorphic genes encoding mehamycin and a new multigene family of bDTAFPs (named roremycins) from Rotylenchulus reniformis. A combination of data suggests that the origination of sDBDs from M. hapla and R. reniformis is a consequence of convergent evolution, in which some probably suffered positive selection during evolution. Our study may be valuable in understanding the role of these unique antimicrobial peptides in the innate immunity of nematodes.

Engineered factor Xa variants retain procoagulant activity independent of direct factor Xa inhibitors

The absence of an adequate reversal strategy to prevent and stop potential life-threatening bleeding complications is a major drawback to the clinical use of the direct oral inhibitors of blood coagulation factor Xa. Here we show that specific modifications of the substrate-binding aromatic S4 subpocket within the factor Xa active site disrupt high-affinity engagement of the direct factor Xa inhibitors. These modifications either entail amino-acid substitution of S4 subsite residues Tyr99 and/or Phe174 (chymotrypsinogen numbering), or extension of the 99-loop that borders the S4 subsite. The latter modifications led to the engineering of a factor Xa variant that is able to support coagulation in human plasma spiked with (supra-)physiological concentrations of direct factor Xa inhibitors. As such, this factor Xa variant has the potential to be employed to bypass the direct factor Xa inhibitor-mediated anticoagulation in patients that require restoration of blood coagulation.

A major drawback in the clinical use of the oral anticoagulants that directly inhibit factor Xa in order to prevent blood clot formation is the potential for life threatening bleeding events. Here the authors describe factor Xa variants that are refractory to inhibition by these anticoagulants and could serve as rescue agents in treated patients.

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.

Target-Driven Positive Selection at Hot Spots of Scorpion Toxins Uncovers Their Potential in Design of Insecticides

Positive selection sites (PSSs), a class of amino acid sites with an excess of nonsynonymous to synonymous substitutions, are indicators of adaptive molecular evolution and have been detected in many protein families involved in a diversity of biological processes by statistical approaches. However, few studies are conducted to evaluate their functional significance and the driving force behind the evolution (i.e., agent of selection). Scorpion α-toxins are a class of multigene family of peptide neurotoxins affecting voltage-gated Na(+ )(Nav) channels, whose members exhibit differential potency and preference for insect and mammalian Nav channels. In this study, we undertook a systematical molecular dissection of nearly all the PSSs newly characterized in the Mesobuthus α-toxin family and a two-residue insertion ((19)AlaPhe(20)) located within a positively selected loop via mutational analysis of α-like MeuNaTxα-5, one member affecting both insect and mammalian Nav channels. This allows to identify hot-spot residues on its functional face involved in interaction with the receptor site of Nav channels, which comprises two PSSs (Ile(40) and Leu(41)) and the small insertion, both located on two spatially separated functional loops. Mutations at these hot-spots resulted in a remarkably decreased anti-mammalian activity in MeuNaTxα-5 with partially impaired or enhanced insecticide activity, suggesting the potential of PSSs in designing promising candidate insecticides from scorpion α-like toxins. Based on an experiment-guided toxin-channel complex model and high evolutionary variability in the receptor site of predators and prey of scorpions, we provide new evidence for target-driven adaptive evolution of scorpion toxins to deal with their targets' diversity.

A Single Glycan at the 99-Loop of Human Kallikrein-related Peptidase 2 Regulates Activation and Enzymatic Activity

Human kallikrein-related peptidase 2 (KLK2) is a key serine protease in semen liquefaction and prostate cancer together with KLK3/prostate-specific antigen. In order to decipher the function of its potential N-glycosylation site, we produced pro-KLK2 in Leishmania tarentolae cells and compared it with its non-glycosylated counterpart from Escherichia coli expression. Mass spectrometry revealed that Asn-95 carries a core glycan, consisting of two GlcNAc and three hexoses. Autocatalytic activation was retarded in glyco-pro-KLK2, whereas the activated glyco-form exhibited an increased proteolytic resistance. The specificity patterns obtained by the PICS (proteomic identification of protease cleavage sites) method are similar for both KLK2 variants, with a major preference for P1-Arg. However, glycosylation changes the enzymatic activity of KLK2 in a drastically substrate-dependent manner. Although glyco-KLK2 has a considerably lower catalytic efficiency than glycan-free KLK2 toward peptidic substrates with P2-Phe, the situation was reverted toward protein substrates, such as glyco-pro-KLK2 itself. These findings can be rationalized by the glycan-carrying 99-loop that prefers to cover the active site like a lid. By contrast, the non-glycosylated 99-loop seems to favor a wide open conformation, which mostly increases the apparent affinity for the substrates (i.e. by a reduction of K_m). Also, the cleavage pattern and kinetics in autolytic inactivation of both KLK2 variants can be explained by a shift of the target sites due to the presence of the glycan. These striking effects of glycosylation pave the way to a deeper understanding of kallikrein-related peptidase biology and pathology.

Why do we study animal toxins?

Venom (toxins) is an important trait evolved along the evolutionary tree of animals. Our knowledges on venoms, such as their origins and loss, the biological relevance and the coevolutionary patterns with other organisms are greatly helpful in understanding many fundamental biological questions, i.e., the environmental adaptation and survival competition, the evolution shaped development and balance of venoms, and the sophisticated correlations among venom, immunity, body power, intelligence, their genetic basis, inherent association, as well as the cost-benefit and trade-offs of biological economy. Lethal animal envenomation can be found worldwide. However, from foe to friend, toxin studies have led lots of important discoveries and exciting avenues in deciphering and fighting human diseases, including the works awarded the Nobel Prize and lots of key clinic therapeutics. According to our survey, so far, only less than 0.1% of the toxins of the venomous animals in China have been explored. We emphasize on the similarities shared by venom and immune systems, as well as the studies of toxin knowledge-based physiological toxin-like proteins/peptides (TLPs). We propose the natural pairing hypothesis. Evolution links toxins with humans. Our mission is to find out the right natural pairings and interactions of our body elements with toxins, and with endogenous toxin-like molecules. Although, in nature, toxins may endanger human lives, but from a philosophical point of view, knowing them well is an effective way to better understand ourselves. So, this is why we study toxins.

The first venomous crustacean revealed by transcriptomics and functional morphology: remipede venom glands express a unique toxin cocktail dominated by enzymes and a neurotoxin

Animal venoms have evolved many times. Venomous species are especially common in three of the four main groups of arthropods (Chelicerata, Myriapoda, and Hexapoda), which together represent tens of thousands of species of venomous spiders, scorpions, centipedes, and hymenopterans. Surprisingly, despite their great diversity of body plans, there is no unambiguous evidence that any crustacean is venomous. We provide the first conclusive evidence that the aquatic, blind, and cave-dwelling remipede crustaceans are venomous and that venoms evolved in all four major arthropod groups. We produced a three-dimensional reconstruction of the venom delivery apparatus of the remipede Speleonectes tulumensis, showing that remipedes can inject venom in a controlled manner. A transcriptomic profile of its venom glands shows that they express a unique cocktail of transcripts coding for known venom toxins, including a diversity of enzymes and a probable paralytic neurotoxin very similar to one described from spider venom. We screened a transcriptomic library obtained from whole animals and identified a nontoxin paralog of the remipede neurotoxin that is not expressed in the venom glands. This allowed us to reconstruct its probable evolutionary origin and underlines the importance of incorporating data derived from nonvenom gland tissue to elucidate the evolution of candidate venom proteins. This first glimpse into the venom of a crustacean and primitively aquatic arthropod reveals conspicuous differences from the venoms of other predatory arthropods such as centipedes, scorpions, and spiders and contributes valuable information for ultimately disentangling the many factors shaping the biology and evolution of venoms and venomous species.

The convergent evolution of blue iris pigmentation in primates took distinct molecular paths

How many distinct molecular paths lead to the same phenotype? One approach to this question has been to examine the genetic basis of convergent traits, which likely evolved repeatedly under a shared selective pressure. We investigated the convergent phenotype of blue iris pigmentation, which has arisen independently in four primate lineages: humans, blue-eyed black lemurs, Japanese macaques, and spider monkeys. Characterizing the phenotype across these species, we found that the variation within the blue-eyed subsets of each species occupies strongly overlapping regions of CIE L*a*b* color space. Yet whereas Japanese macaques and humans display continuous variation, the phenotypes of blue-eyed black lemurs and their sister species (whose irises are brown) occupy more clustered subspaces. Variation in an enhancer of OCA2 is primarily responsible for the phenotypic difference between humans with blue and brown irises. In the orthologous region, we found no variant that distinguishes the two lemur species or associates with quantitative phenotypic variation in Japanese macaques. Given the high similarity between the blue iris phenotypes in these species and that in humans, this finding implies that evolution has used different molecular paths to reach the same end. Am J Phys Anthropol 151:398–407, 2013.© 2013 Wiley Periodicals, Inc.

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.

Molecular evolution of α-latrotoxin, the exceptionally potent vertebrate neurotoxin in black widow spider venom

Black widow spiders (members of the genus Latrodectus) are widely feared because of their potent neurotoxic venom. α-Latrotoxin is the vertebrate-specific toxin responsible for the dramatic effects of black widow envenomation. The evolution of this toxin is enigmatic because only two α-latrotoxin sequences are known. In this study, ~4 kb α-latrotoxin sequences and their homologs were characterized from a diversity of Latrodectus species, and representatives of Steatoda and Parasteatoda, establishing the wide distribution of latrotoxins across the mega-diverse spider family Theridiidae. Across black widow species, α-latrotoxin shows ≥ 94% nucleotide identity and variability consistent with purifying selection. Multiple codon and branch-specific estimates of the nonsynonymous/synonymous substitution rate ratio also suggest a long history of purifying selection has acted on α-latrotoxin across Latrodectus and Steatoda. However, α-latrotoxin is highly divergent in amino acid sequence between these genera, with 68.7% of protein differences involving non-conservative substitutions, evidence for positive selection on its physiochemical properties and particular codons, and an elevated rate of nonsynonymous substitutions along α-latrotoxin's Latrodectus branch. Such variation likely explains the efficacy of red-back spider, L. hasselti, antivenom in treating bites from other Latrodectus species, and the weaker neurotoxic symptoms associated with Steatoda and Parasteatoda bites. Long-term purifying selection on α-latrotoxin indicates its functional importance in black widow venom, even though vertebrates are a small fraction of their diet. The greater differences between Latrodectus and Steatoda α-latrotoxin, and their relationships to invertebrate-specific latrotoxins, suggest a shift in α-latrotoxin toward increased vertebrate toxicity coincident with the evolution of widow spiders.

Convergent evolution of a reproductive trait through distinct developmental mechanisms in Drosophila

Convergent morphologies often arise due to similar selective pressures in independent lineages. It is poorly understood whether the same or different developmental genetic mechanisms underlie such convergence. Here we show that independent evolution of a reproductive trait, ovariole number, has resulted from changes in distinct developmental mechanisms, each of which may have a different underlying genetic basis in Drosophila. Ovariole number in Drosophila is species-specific, highly variable, and largely under genetic control. Convergent changes in Drosophila ovariole number have evolved independently within and between species. We previously showed that the number of a specific ovarian cell type, terminal filament (TF) cells, determines ovariole number. Here we examine TF cell development in different Drosophila lineages that independently evolved a significantly lower ovariole number than the D. melanogaster Oregon R strain. We show that in these Drosophila lineages, reduction in ovariole number occurs primarily through variations in one of two different developmental mechanisms: (1) reduced number of somatic gonad precursors (SGP cells) specified during embryogenesis; or (2) alterations of somatic gonad cell morphogenesis and differentiation in larval life. Mutations in the D. melanogaster Insulin Receptor (InR) alter SGP cell number but not ovarian morphogenesis, while targeted loss of function of bric-à-brac 2 (bab2) affects morphogenesis without changing SGP cell number. Thus, evolution can produce similar ovariole numbers through distinct developmental mechanisms, likely controlled by different genetic mechanisms.

Expression of kallikrein-related peptidase 4 in dental and non-dental tissues

Kallikrein-related peptidase 4 (KLK4) is critical for proper dental enamel formation. Klk4 null mice, and humans with two defective KLK4 alleles have obvious enamel defects, with no other apparent phenotype. KLK4 mRNA or protein is reported to be present in tissues besides teeth, including prostate, ovary, kidney, liver, and salivary gland. In this study we used the Klk4 knockout/NLS-lacZ knockin mouse to assay Klk4 expression using β-galactosidase histochemistry. Incubations for 5 h were used to detect KLK4 expression with minimal endogenous background, while overnight incubations susceptible to false positives were used to look for trace KLK4 expression. Developing maxillary molars at postnatal days 5, 6, 7, 8, and 14, developing mandibular incisors at postnatal day 14, and selected non-dental tissues from adult wild-type and Klk4(lacZ/lacZ) mice were examined by X-gal histochemistry. After 5 h of incubation, X-gal staining was observed specifically in the nuclei of maturation-stage ameloblasts in molars and incisors from Klk4(lacZ/lacZ) mice and was detected weakly in the nuclei of salivary gland ducts and in patches of prostate epithelia. We conclude that KLK4 is predominantly a tooth-specific protease with low expression in submandibular salivary gland and prostate, and with no detectable expression in liver, kidney, testis, ovary, oviduct, epididymis, and vas deferens.

Protein-protein interaction network-based detection of functionally similar proteins within species

Although functionally similar proteins across species have been widely studied, functionally similar proteins within species showing low sequence similarity have not been examined in detail. Identification of these proteins is of significant importance for understanding biological functions, evolution of protein families, progression of co-evolution, and convergent evolution and others which cannot be obtained by detection of functionally similar proteins across species. Here, we explored a method of detecting functionally similar proteins within species based on graph theory. After denoting protein-protein interaction networks using graphs, we split the graphs into subgraphs using the 1-hop method. Proteins with functional similarities in a species were detected using a method of modified shortest path to compare these subgraphs and to find the eligible optimal results. Using seven protein-protein interaction networks and this method, some functionally similar proteins with low sequence similarity that cannot detected by sequence alignment were identified. By analyzing the results, we found that, sometimes, it is difficult to separate homologous from convergent evolution. Evaluation of the performance of our method by gene ontology term overlap showed that the precision of our method was excellent.

Venomous mammals: a review

The occurrence of venom in mammals has long been considered of minor importance, but recent fossil discoveries and advances in experimental techniques have cast new light into this subject. Mammalian venoms form a heterogeneous group having different compositions and modes of action and are present in three classes of mammals, Insectivora, Monotremata, and Chiroptera. A fourth order, Primates, is proposed to have venomous representatives. In this review we highlight recent advances in the field while summarizing biochemical characteristics of these secretions and their effects upon humans and other animals. Historical aspects of venom discovery and evolutionary hypothesis regarding their origin are also discussed.

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.

Structures of Phytophthora RXLR effector proteins: a conserved but adaptable fold underpins functional diversity

Phytopathogens deliver effector proteins inside host plant cells to promote infection. These proteins can also be sensed by the plant immune system, leading to restriction of pathogen growth. Effector genes can display signatures of positive selection and rapid evolution, presumably a consequence of their co-evolutionary arms race with plants. The molecular mechanisms underlying how effectors evolve to gain new virulence functions and/or evade the plant immune system are poorly understood. Here, we report the crystal structures of the effector domains from two oomycete RXLR proteins, Phytophthora capsici AVR3a11 and Phytophthora infestans PexRD2. Despite sharing <20% sequence identity in their effector domains, they display a conserved core α-helical fold. Bioinformatic analyses suggest that the core fold occurs in ∼44% of annotated Phytophthora RXLR effectors, both as a single domain and in tandem repeats of up to 11 units. Functionally important and polymorphic residues map to the surface of the structures, and PexRD2, but not AVR3a11, oligomerizes in planta. We conclude that the core α-helical fold enables functional adaptation of these fast evolving effectors through (i) insertion/deletions in loop regions between α-helices, (ii) extensions to the N and C termini, (iii) amino acid replacements in surface residues, (iv) tandem domain duplications, and (v) oligomerization. We hypothesize that the molecular stability provided by this core fold, combined with considerable potential for plasticity, underlies the evolution of effectors that maintain their virulence activities while evading recognition by the plant immune system.

Evolutionary diversification of Mesobuthus α-scorpion toxins affecting sodium channels

α-Scorpion toxins constitute a family of peptide modulators that induce a prolongation of the action potential of excitable cells by inhibiting voltage-gated sodium channel inactivation. Although they all adopt a conserved structural scaffold, the potency and phylogentic preference of these toxins largely vary, which render them an intriguing model for studying evolutionary diversification among family members. Here, we report molecular characterization of a new multigene family of α-toxins comprising 13 members (named MeuNaTxα-1 to MeuNaTxα-13) from the scorpion Mesobuthus eupeus. Of them, five native toxins (MeuNaTxα-1 to -5) were purified to homogeneity from the venom and the solution structure of MeuNaTxα-5 was solved by nuclear magnetic resonance. A systematic functional evaluation of MeuNaTxα-1, -2, -4, and -5 was conducted by two-electrode voltage-clamp recordings on seven cloned mammalian voltage-gated sodium channels (Na(v)1.2 to Na(v)1.8) and the insect counterpart DmNa(v)1 expressed in Xenopus oocytes. Results show that all these four peptides slow inactivation of DmNa(v)1 and are inactive on Na(v)1.8 at micromolar concentrations. However, they exhibit differential specificity for the other six channel isoforms (Na(v)1.2 to Na(v)1.7), in which MeuNaTxα-4 shows no activity on these isoforms and thus represents the first Mesobuthus-derived insect-selective α-toxin identified so far with a half maximal effective concentration of 130 ± 2 nm on DmNa(v)1 and a half maximal lethal dose of about 200 pmol g(-1) on the insect Musca domestica; MeuNaTxα-2 only affects Na(v)1.4; MeuNaTxα-1 and MeuNaTxα-5 have a wider range of channel spectrum, the former active on Na(v)1.2, Na(v)1.3, Na(v)1.6, and Na(v)1.7, whereas the latter acting on Na(v)1.3-Na(v)1.7. Remarkably, MeuNaTxα-4 and MeuNaTxα-5 are two nearly identical peptides differing by only one point mutation at site 50 (A50V) but exhibit rather different channel subtype selectivity, highlighting a switch role of this site in altering the target specificity. By the maximum likelihood models of codon substitution, we detected nine positively selected sites (PSSs) that could be involved in functional diversification of Mesobuthus α-toxins. The PSSs include site 50 and other seven sites located in functional surfaces of α-toxins. This work represents the first thorough investigation of evolutionary diversification of α-toxins derived from a specific scorpion lineage from the perspectives of sequence, structure, function, and evolution.

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.

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.

Thermostability in endoglucanases is fold-specific

Background

Endoglucanases are usually considered to be synergistically involved in the initial stages of cellulose breakdown-an essential step in the bioprocessing of lignocellulosic plant materials into bioethanol. Despite their economic importance, we currently lack a basic understanding of how some endoglucanases can sustain their ability to function at elevated temperatures required for bioprocessing, while others cannot. In this study, we present a detailed comparative analysis of both thermophilic and mesophilic endoglucanases in order to gain insights into origins of thermostability. We analyzed the sequences and structures for sets of endoglucanase proteins drawn from the Carbohydrate-Active enZymes (CAZy) database.

Results

Our results demonstrate that thermophilic endoglucanases and their mesophilic counterparts differ significantly in their amino acid compositions. Strikingly, these compositional differences are specific to protein folds and enzyme families, and lead to differences in intramolecular interactions in a fold-dependent fashion.

Conclusions

Here, we provide fold-specific guidelines to control thermostability in endoglucanases that will aid in making production of biofuels from plant biomass more efficient.

Evolutionary history of tissue kallikreins

The gene family of human kallikrein-related peptidases (KLKs) encodes proteins with diverse and pleiotropic functions in normal physiology as well as in disease states. Currently, the most widely known KLK is KLK3 or prostate-specific antigen (PSA) that has applications in clinical diagnosis and monitoring of prostate cancer. The KLK gene family encompasses the largest contiguous cluster of serine proteases in humans which is not interrupted by non-KLK genes. This exceptional and unique characteristic of KLKs makes them ideal for evolutionary studies aiming to infer the direction and timing of gene duplication events. Previous studies on the evolution of KLKs were restricted to mammals and the emergence of KLKs was suggested about 150 million years ago (mya). In order to elucidate the evolutionary history of KLKs, we performed comprehensive phylogenetic analyses of KLK homologous proteins in multiple genomes including those that have been completed recently. Interestingly, we were able to identify novel reptilian, avian and amphibian KLK members which allowed us to trace the emergence of KLKs 330 mya. We suggest that a series of duplication and mutation events gave rise to the KLK gene family. The prominent feature of the KLK family is that it consists of tandemly and uninterruptedly arrayed genes in all species under investigation. The chromosomal co-localization in a single cluster distinguishes KLKs from trypsin and other trypsin-like proteases which are spread in different genetic loci. All the defining features of the KLKs were further found to be conserved in the novel KLK protein sequences. The study of this unique family will further assist in selecting new model organisms for functional studies of proteolytic pathways involving KLKs.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

Kallikreins on steroids: structure, function, and hormonal regulation of prostate-specific antigen and the extended kallikrein locus

The 15 members of the kallikrein-related serine peptidase (KLK) family have diverse tissue-specific expression profiles and putative proteolytic functions. The kallikrein family is also emerging as a rich source of disease biomarkers with KLK3, commonly known as prostate-specific antigen, being the current serum biomarker for prostate cancer. The kallikrein locus is also notable because it is extraordinarily responsive to steroids and other hormones. Indeed, at least 14 functional hormone response elements have been identified in the kallikrein locus. A more comprehensive understanding of the transcriptional regulation of kallikreins may help the field make more informed hypotheses about the physiological functions of kallikreins and their effectiveness as biomarkers. In this review, we describe the organization of the kallikrein locus and the structure of kallikrein genes and proteins. We also focus on the transcriptional regulation of kallikreins by androgens, progestins, glucocorticoids, mineralocorticoids, estrogens, and other hormones in animal models and human prostate, breast, and reproductive tract tissues. The interaction of the androgen receptor with androgen response elements in the promoter and enhancer of KLK2 and KLK3 is also summarized in detail. There is evidence that all kallikreins are regulated by multiple nuclear receptors. Yet, apart from KLK2 and KLK3, it is not clear whether all kallikreins are direct transcriptional targets. Therefore, we argue that gaining more detailed information about the mechanisms that regulate kallikrein expression should be a priority of future studies and that the kallikrein locus will continue to be an important model in the era of genome-wide analyses.

Drosotoxin, a selective inhibitor of tetrodotoxin-resistant sodium channels

The design of animal toxins with high target selectivity has long been a goal in protein engineering. Based on evolutionary relationship between the Drosophila antifungal defensin (drosomycin) and scorpion depressant Na(+) channel toxins, we exploited a strategy to create a novel chimeric molecule (named drosotoxin) with high selectivity for channel subtypes, which was achieved by using drosomycin to substitute the structural core of BmKITc, a depressant toxin acting on both insect and mammalian Na(+) channels. Recombinant drosotoxin selectively inhibited tetrodotoxin-resistant (TTX-R) Na(+) channels in rat dorsal root ganglion (DRG) neurons with a 50% inhibitory concentration (IC(50)) of 2.6+/-0.5muM. This chimeric peptide showed no activity on K(+), Ca(2+) and TTX-sensitive (TTX-S) Na(+) channels in rat DRG neurons and Drosophila para/tipE channels at micromolar concentrations. Drosotoxin represents the first chimeric toxin and example of a non-toxic core scaffold with high selectivity on mammalian TTX-R Na(+) channels.

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.