Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans.

Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ∼ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.

Rapid evolution of plethodontid modulating factor, a hypervariable salamander courtship pheromone, is driven by positive selection.

Sexual communication in plethodontid salamanders is mediated by a proteinaceous pheromone that a male delivers to a female during courtship, boosting her receptivity. The pheromone consists of three proteins from three unrelated protein families. These proteins are among a small group of pheromones known to affect female receptivity in vertebrates. Previously, we showed that the genes of two of these proteins (PRF and SPF) are prone to incessant evolution driven by positive selection, presumably as a consequence of coevolution with female receptors. In this report, we focus on the evolution of the third pheromone protein gene family, plethodontid modulating factor (PMF), to determine whether it shows the same pattern of diversification. We used RT-PCR in mental gland cDNA to survey PMF sequences from three genera of plethodontid salamanders (27 spp.) to measure rates of evolution, level of gene diversification, modes of selection, and types of amino acid substitution. Like PRF and SPF, PMF is produced by a multigene family characterized by gene duplication and high levels of polymorphism. PMF evolution is rapid, incessant, and driven by positive selection. PMF is more extreme in these dimensions than both PRF and SPF. Nestled within this extraordinary variation, however, is a signature of purifying selection, acting to preserve important structural and biochemical features of the PMF protein (i.e., secretion signal, cysteine residues, and pI). Although a pattern of persistent diversification exists at the molecular level, the morphological and behavioral aspects of the pheromone delivery system show evolutionary stasis over millions of years.

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.

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.