Evolutionary replacement of components in a salamander pheromone signaling complex: more evidence for phenotypic-molecular decoupling

In this article we explore the evolutionary history of a functional complex at the molecular level in plethodontid salamanders. The complex consists of a proteinaceous courtship pheromone, a pheromone-producing gland on the male's chin, and a set of behaviors for delivering the pheromone to the female. Long-term evolutionary stasis is the defining feature of this complex at both the morphological and behavioral levels. However, our previous assessment of the pheromone gene, plethodontid receptivity factor (PRF), revealed rapid evolution at the molecular level despite stasis at higher levels of organization. Analysis of a second pheromone gene, sodefrin precursor-like factor (SPF), now indicates that evolutionary decoupling in this complex is pervasive. The evolutionary profiles of SPF and PRF are remarkably similar in that: (a) both genes exhibit high levels of sequence diversity both within and across taxa, (b) genetic diversity has been driven by strong positive selection, and (c) the genes have evolved heterogeneously in different salamander lineages. The composition of the pheromone signal as a whole, however, has experienced an extraordinary evolutionary transition. Whereas SPF has been retained throughout the 100 MY radiation of salamanders, PRF has only recently been recruited to a pheromone function (27 million years ago). When SPF and PRF coexist in the same clade, they show contrasting patterns of evolution. When one shows rapid evolution driven by positive selection, the other shows neutral divergence restrained by purifying selection. In one clade, the origin and subsequent rapid evolution of PRF appear to have interfered with the evolution and persistence of SPF, leading to a pattern of evolutionary replacement. Overall, these two pheromone genes provide a revealing window on the dynamics that drive the evolution of multiple traits in a signaling complex.     

Lineage-specific differences in evolutionary mode in a salamander courtship pheromone

Functionally equivalent genes may evolve heterogeneously across closely related taxa as a consequence of lineage-specific selective pressures. Such disparate evolutionary modes are especially prevalent in genes that encode postcopulatory reproductive proteins, presumably as a result of sexual selection. We might therefore expect genes that mediate reproduction prior to insemination to evolve in a similar manner. Plethodontid receptivity factor (PRF), a proteinaceous salamander pheromone produced by the male, increases female receptivity during courtship interactions. To test for lineage-specific differences in PRF's evolution, we intensively sampled PRF genes across the eastern Plethodon phylogeny (27 spp.; 34 populations) to compare gene diversification, rates of evolution, modes of selection, and types of amino acid substitution. Our analyses indicate that PRF evolutionary dynamics vary considerably from lineage to lineage. Underlying this heterogeneity, however, are two well-defined transitions in evolutionary mode. The first mode is representative of a typical protein profile, wherein neutral divergence and purifying selection are the dominant features. The second mode is characterized by incessant, cyclical evolution driven by positive selection. In this mode, the positively selected sites are bound by a limited assortment of acceptable amino acids that appear to evolve independently of other sites, resulting in a tremendous number of unique PRF alleles. Several of these selected sites are implicated in receptor binding. These sites are apparently involved in a molecular tango in which the male signal and female receptors coevolve within a confined molecular space. PRF's lineage-specific evolutionary dynamics, in combination with evidence of a molecular tango, highlight the molecular action of sexual selection on a chemical signal that is used during courtship.     

Stabilizing selection on behavior and morphology masks positive selection on the signal in a salamander pheromone signaling complex

Natural selection maintains the integration and coordination of sets of phenotypic characters that collectively perform a task. In functional complexes in which characters span molecular to behavioral levels of organization, we might then expect similar modes of selection to produce similar patterns in evolutionary divergence at each level. To test this expectation, we diagnosed selection at behavioral, morphological, and molecular levels for courtship pheromone signaling by plethodontid salamanders. At the levels of morphology and behavior tens of millions of years of stasis (stabilizing selection) occur on each side of a transition from vaccination to olfactory delivery modes. As a proxy for the molecular level, we used plethodontid receptivity factor (PRF), a protein that is an active component of the pheromone. We cloned PRF from 12 Plethodon spp. spanning the delivery transition and obtained multiple alleles from each individual surveyed. Analyses of 61 alleles for PRF identified elevated nonsynonymous over synonymous substitution rates along lineages in a molecular phylogeny, and at 8% of sites in the protein, indicating that positive (directional) selection has acted on this vertebrate pheromone gene. Structural models showed PRF is in a family of cytokines characterized by a four-alpha-helix bundle. Positive selection in PRF was associated with receptor binding sites that are under purifying selection in other cytokines of that family. The evolutionary dynamics of the plethodontid pheromone delivery complex consists of stabilizing selection on morphological and behavioral aspects of signal delivery but positive selection on the signal mediated by receptors. Thus, different selection modes prevail at different levels in this reproductive functional complex. Evolutionary studies of integrated sets of characters therefore require separate analyses of selective action at each level.     

A recombinant courtship pheromone affects sexual receptivity in a plethodontid salamander

Pheromones are important chemical signals for many vertebrates, particularly during reproductive interactions. In the terrestrial salamander Plethodon shermani, a male delivers proteinaceous pheromones to the female as part of their ritualistic courtship behavior. These pheromones increase the female's receptivity to mating, as shown by a reduction in courtship duration. One pheromone component in particular is plethodontid receptivity factor (PRF), a 22-kDa protein with multiple isoforms. This protein alone can act as a courtship pheromone that causes the female to be more receptive. We used a bacterial expression system to synthesize a single recombinant isoform of PRF. The recombinant protein was identical to the native PRF, based on mass spectrometry, circular dichroism spectra, and a behavioral bioassay that tested the effects of recombinant PRF (rPRF) on female receptivity (21% reduction in courtship duration). The rPRF appears to mimic the activity of a mixture of PRF isoforms, as well as a mixture of multiple different proteins that comprise the male courtship gland extract. Pheromones that are peptides have been characterized for some vertebrates; to date, however, rPRF is one of only 2 synthesized vertebrate proteins to retain full biological activity.     

Plethodontid modulating factor, a hypervariable salamander courtship pheromone in the three-finger protein superfamily

The soluble members of the three-finger protein superfamily all share a relatively simple 'three-finger' structure, yet perform radically different functions. Plethodontid modulating factor (PMF), a pheromone protein produced by the lungless salamander, Plethodon shermani, is a new and unusual member of this group. It affects female receptivity when delivered to the female's nares during courtship. As with other plethodontid pheromone genes, PMF is hyperexpressed in a specialized male mental (chin) gland. Unlike other plethodontid pheromone genes, however, PMF is also expressed at low levels in the skin, liver, intestine and kidneys of both sexes. The PMF sequences obtained from all tissue types were highly variable, with 103 unique haplotypes identified which averaged 35% sequence dissimilarity (range 1-60%) at the protein level. Despite this variation, however, all PMF sequences contained a conserved approximately 20-amino-acid secretion signal sequence and a pattern of eight cysteines that is also found in cytotoxins and short neurotoxins from snake venoms, as well as xenoxins from Xenopus. Although they share a common cysteine pattern, PMF isoforms differ from other three-finger proteins in: (a) amino-acid composition outside of the conserved motif; (b) length of the three distinguishing 'fingers'; (c) net charge at neutral pH. Whereas most three-finger proteins have a net positive charge at pH 7.0, PMF has a high net negative charge at neutral pH (pI range of most PMFs 3.5-4.0). Sequence comparisons suggest that PMF belongs to a distinct multigene subfamily within the three-finger protein superfamily.     

Structural Insights into the Evolution of a Sexy Protein: Novel Topology and Restricted Backbone Flexibility in a Hypervariable Pheromone from the Red-Legged Salamander,Plethodon shermani

In response to pervasive sexual selection, protein sex pheromones often display rapid mutation and accelerated evolution of corresponding gene sequences. For proteins, the general dogma is that structure is maintained even as sequence or function may rapidly change. This phenomenon is well exemplified by the three-finger protein (TFP) superfamily: a diverse class of vertebrate proteins co-opted for many biological functions – such as components of snake venoms, regulators of the complement system, and coordinators of amphibian limb regeneration. All of the >200 structurally characterized TFPs adopt the namesake “three-finger” topology. In male red-legged salamanders, the TFP pheromone Plethodontid Modulating Factor (PMF) is a hypervariable protein such that, through extensive gene duplication and pervasive sexual selection, individual male salamanders express more than 30 unique isoforms. However, it remained unclear how this accelerated evolution affected the protein structure of PMF. Using LC/MS-MS and multidimensional NMR, we report the 3D structure of the most abundant PMF isoform, PMF-G. The high resolution structural ensemble revealed a highly modified TFP structure, including a unique disulfide bonding pattern and loss of secondary structure, that define a novel protein topology with greater backbone flexibility in the third peptide finger. Sequence comparison, models of molecular evolution, and homology modeling together support that this flexible third finger is the most rapidly evolving segment of PMF. Combined with PMF sequence hypervariability, this structural flexibility may enhance the plasticity of PMF as a chemical signal by permitting potentially thousands of structural conformers. We propose that the flexible third finger plays a critical role in PMF:receptor interactions. As female receptors co-evolve, this flexibility may allow PMF to still bind its receptor(s) without the immediate need for complementary mutations. Consequently, this unique adaptation may establish new paradigms for how receptor:ligand pairs co-evolve, in particular with respect to sexual conflict.     

Are rates of species diversification correlated with rates of morphological evolution?

Some major evolutionary theories predict a relationship between rates of proliferation of new species (species diversification) and rates of morphological divergence between them. However, this relationship has not been rigorously tested using phylogeny-based approaches. Here, we test this relationship with morphological and phylogenetic data from 190 species of plethodontid salamanders. Surprisingly, we find that rates of species diversification and morphological evolution are not significantly correlated, such that rapid diversification can occur with little morphological change, and vice versa. We also find that most clades have undergone remarkably similar patterns of morphological evolution (despite extensive sympatry) and that those relatively novel phenotypes are not associated with rapid diversification. Finally, we find a strong relationship between rates of size and shape evolution, which has not been previously tested.     

Conspecific and heterospecific pheromone effects on female receptivity

Courtship pheromones play an important role in salamander reproductive behaviour. In salamanders of the family Plethodontidae, males deliver specialized pheromones to females during courtship interactions. These courtship pheromones increase female receptivity and may be involved in mate discrimination. In order to test hypotheses related to mate discrimination, we staged courtship encounters between male-female Plethodon shermani pairs in which the female received pheromones obtained from either conspecific (P. shermani) or heterospecific (P. yonahlossee orP. montanus ) males. Both conspecific and heterospecific pheromones increased female receptivity. Moreover, pheromones from both heterospecific species were as effective as the conspecific pheromone in increasing female receptivity inP. shermani females. Our results suggest that the courtship pheromone signal and function may be conserved across related species, with mate discrimination occurring before pheromone delivery. Copyright 2003 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.     

What explains patterns of species richness? The relative importance of climatic‐niche evolution,morphological evolution,and ecological limits in salamanders

A major goal of evolutionary biology and ecology is to understand why species richness varies among clades. Previous studies have suggested that variation in richness among clades might be related to variation in rates of morphological evolution among clades (e.g., body size and shape). Other studies have suggested that richness patterns might be related to variation in rates of climatic‐niche evolution. However, few studies, if any, have tested the relative importance of these variables in explaining patterns of richness among clades. Here, we test their relative importance among major clades of Plethodontidae, the most species‐rich family of salamanders. Earlier studies have suggested that climatic‐niche evolution explains patterns of diversification among plethodontid clades, whereas rates of morphological evolution do not. A subsequent study stated that rates of morphological evolution instead explained patterns of species richness among plethodontid clades (along with “ecological limits” on richness of clades, leading to saturation of clades with species, given limited resources). However, they did not consider climatic‐niche evolution. Using phylogenetic multiple regression, we show that rates of climatic‐niche evolution explain most variation in richness among plethodontid clades, whereas rates of morphological evolution do not. We find little evidence that ecological limits explain patterns of richness among plethodontid clades. We also test whether rates of morphological and climatic‐niche evolution are correlated, and find that they are not. Overall, our results help explain richness patterns in a major amphibian group and provide possibly the first test of the relative importance of climatic niches and morphological evolution in explaining diversity patterns.     

Phylogenetic analysis of ecomorphological divergence, community structure, and diversification rates in dusky salamanders (Plethodontidae: Desmognathus)

An important dimension of adaptive radiation is the degree to which diversification rates fluctuate or remain constant through time. Focusing on plethodontid salamanders of the genus Desmognathus, we present a novel synthetic analysis of phylogeographic history, rates of ecomorphological evolution and species accumulation, and community assembly in an adaptive radiation. Dusky salamanders are highly variable in life history, body size, and ecology, with many endemic lineages in the southern Appalachian Highlands of eastern North America. Our results show that life-history evolution had important consequences for the buildup of plethodontid-salamander species richness and phenotypic disparity in eastern North America, a global hot spot of salamander biodiversity. The origin of Desmognathus species with aquatic larvae was followed by a high rate of lineage accumulation, which then gradually decreased toward the present time. The peak period of lineage accumulation in the group coincides with evolutionary partitioning of lineages with aquatic larvae into seepage, stream-edge, and stream microhabitats. Phylogenetic simulations demonstrate a strong correlation between morphology and microhabitat ecology independent of phylogenetic effects and suggest that ecomorphological changes are concentrated early in the radiation of Desmognathus. Deep phylogeographic fragmentation within many codistributed ecomorph clades suggests long-term persistence of ecomorphological features and stability of endemic lineages and communities through multiple climatic cycles. Phylogenetic analyses of community structure show that ecomorphological divergence promotes the coexistence of lineages and that repeated, independent evolution of microhabitat-associated ecomorphs has a limited role in the evolutionary assembly of Desmognathus communities. Comparing and contrasting our results to other adaptive radiations having different biogeographic histories, our results suggest that rates of diversification during adaptive radiation are intimately linked to the degree to which community structure persists over evolutionary time.     

Temperature dependence of evolutionary diversification: differences between two contrasting model taxa support the metabolic theory of ecology

Biodiversity patterns are largely determined by variation of diversification rates across clades and geographic regions. Although there are multiple reasons for this variation, it has been hypothesized that metabolic rate is the crucial driver of diversification of evolutionary lineages. According to the metabolic theory of ecology (MTE), metabolic rate – and consequently speciation – is driven mainly by body size and environmental temperature. As environmental temperature affects metabolic rate in ecto‐ and endotherms differently, its impact on diversification rate should also differ between the two types of organisms. Employing two independent approaches, we analysed correlates of speciation rates and, ultimately, net diversification rates for two contrasting taxa: plethodontid salamanders and carnivoran mammals. Whereas in the ectothermic plethodontids speciation rates positively correlated with environmental temperature, in the endothermic carnivorans a reverse, negative correlation was detected. These findings comply with predictions of the MTE and suggest that similar geographic patterns of biodiversity across taxa (e.g. ecto‐ and endotherms) might have been generated by different ecological and evolutionary processes.     

Pheromonal communication in amphibians

Pheromonal communication is widespread in salamanders and newts and may also be important in some frogs and toads. Several amphibian pheromones have been behaviorally, biochemically and molecularly identified. These pheromones are typically peptides or proteins. Study of pheromone evolution in plethodontid salamanders has revealed that courtship pheromones have been subject to continual evolutionary change, perhaps as a result of co-evolution between the pheromonal ligand and its receptor. Pheromones are detected by the vomeronasal organ and main olfactory epithelium. Chemosensory neurons express vomeronasal receptors or olfactory receptors. Frogs have relatively large numbers of vomeronasal receptors that are transcribed in both the vomeronasal organ and the main olfactory epithelium. Salamander vomeronasal receptors apparently are restricted to the vomeronasal organ. To date, no chemosensory ligands have been matched to vomeronasal receptors or olfactory receptors so it is unknown whether particular receptor types are (1) specialized for detection of pheromones versus other chemosignals, or (2) specialized for detection of volatile, nonvolatile, or water-borne chemosignals. Despite progress in understanding amphibian pheromonal communication, only a small fraction of amphibian species have been examined. Study of additional species of amphibians will indicate which traits related to pheromonal communication are evolutionarily conserved and which traits have diverged over time.     

CAN PARALLEL DIVERSIFICATION OCCUR IN SYMPATRY? REPEATED PATTERNS OF BODY‐SIZE EVOLUTION IN COEXISTING CLADES OF NORTH AMERICAN SALAMANDERS

A classic paradigm in evolutionary biology is that geographically isolated clades inhabiting similar selective regimes will diversify to create similar sets of phenotypes in different locations (e.g., similar stickleback species in different lakes, similar Anolis ecomorphs on different islands). Such parallel radiations are not generally expected to occur in sympatry because the available niche space would be filled by whichever clade is diversified first. Here, we document a very different pattern, the parallel evolution of similar body-size morphs in three sympatric clades of plethodontid salamanders ( Desmognathus, Plethodon, Spelerpinae) in eastern North America. Using a comprehensive, time-calibrated phylogeny of North American plethodontids from nuclear and mitochondrial DNA sequences, we show that these three clades have undergone replicated patterns of evolution in body size and that this parallel diversification occurred in broad-scale sympatry. At the local scale, we find that coexisting species from these clades are more similar in body size than expected under a null model in which species are randomly assembled into communities. These patterns are particularly surprising in that competition is known to be important in driving phenotypic diversification and limiting local coexistence of similar-sized species within these clades. Although parallel diversification of sympatric clades may seem counterintuitive, we discuss several ecological and evolutionary factors that may allow the phenomenon to occur.     

Integrative Studies of Amphibians: From Molecules to Mating

SYNOPSIS. Studies of reproductive behavior in amphibians havebeen especially successful in synthesizing data produced frommolecular and physiological research with evolutionary parameters,such as measures of reproductive success. In particular, tworelatively new areas of amphibian research are highly amenableto synthetic studies. One area is the nature of chemical communicationby pheromone delivery in terrestrial salamanders. Male courtshippheromones, for example, are delivered to the female duringmating. These pheromones typically are received by the accessoryolfactory system and act (presumably via the hypothalamus) toincrease female receptivity. At an evolutionary level, pheromonedelivery can increase male courtship success and thus the likelihoodthat a given male will sire offspring. Variation in pheromonecomposition and effectiveness will permit us to trace the evolutionof the male pheromone on a phylogeny of related populationsand species. At a proximal level, salamander courtship pheromonesare being chemically analyzed in order to identify specificprotein components that affect female receptivity. It also willbe possible to determine whether chemical variation among malesis related to behavioral effectiveness. Thus, courtship pheromonedelivery is a behavior that has invited scrutiny from a combinationof evolutionary and mechanistic perspectives. The second researcharea that is gaining from a synthetic approach is the investigationof relationships between hormonal mechanisms and reproductivebehavior in field populations of anurans and salamanders. Amphibiansare an understudied group, and the area of hormone-behaviorrelationships is inspired by testable hypotheses based on studiesof other terrestrial vertebrates. The ability to correlate physiologicalmeasures with estimates of reproductive success identifies areasof amphibian research that will profit from continued attention.     

Peptide and protein pheromones in amphibians

Purification, characterization and biological activity of urodele and anuran sex-pheromones were reviewed. Female-attracting pheromones obtained from the abdominal gland of Cynops pyrrhogaster and C. ensicauda males are peptides consisting of 10 amino acid residues being designated sodefrin and silefrin, respectively. Each pheromone attracted only conspecific females. Molecular cloning of cDNAs encoding sodefrin and silefrin revealed that both are generated from precursor proteins. Synthesis of these pheromones is regulated by prolactin (PRL) and androgen. Responsiveness of the female vomeronasal epithelium to sodefrin is enhanced by PRL and estrogen. The submandibular gland of the male terrestrial salamander, Plethodon jardani secretes a 22-kD proteinaceous pheromone that enhances female receptivity. It was revealed that every salamander synthesizes multiple isoforms of this pheromone, Plethodontid receptivity factor. The magnificent tree frog, Litoria splendida breed in an aquatic environment. The skin glands of the male secrete a female-attracting peptide pheromone, splendipherin, comprising 25 amino acid residues. The significance of the structure of the amphibian sex-pheromone as peptide and protein is discussed in terms of their species specificity.     

Duplication and selection on abalone sperm lysin in an allopatric population

While gene duplication is a major source of evolutionary novelty, the importance of this process in reproductive protein evolution has not been widely investigated. Here, we report the first known case of gene duplication of abalone sperm lysin in an allopatric subspecies found in the Eastern Atlantic, Haliotis tuberculata coccinea. Mass spectrometry identified both copies of the lysin protein in testis tissue, and 3-dimensional structural modeling suggests that both proteins remain functional. We also detected positive selection acting on both paralogs after duplication and found evidence of a recent selective sweep. Because H. t. coccinea occurs in geographic isolation from other abalone species, these findings suggest that the evolution of lysin is not driven to create reproductive barriers to unfit hybrid formation with an overlapping species. Instead, sexual selection or sexual conflict acting during abalone fertilization could be responsible for the recent positive selection on this protein. The presence of multiple, rapidly evolving lysin genes in H. tuberculata presents an opportunity to study the early stages of diversification of a protein whose function is well understood.