A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function.

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

The origin and diversification of a novel protein family in venomous snakes.

The genetic origins of novelty are a central interest of evolutionary biology. Most new proteins evolve from preexisting proteins but the evolutionary path from ancestral gene to novel protein is challenging to trace, and therefore the requirements for and order of coding sequence changes, expression changes, or gene duplication are not clear. Snake venoms are important novel traits that are comprised of toxins derived from several distinct protein families, but the genomic and evolutionary origins of most venom components are not understood. Here, we have traced the origin and diversification of one prominent family, the snake venom metalloproteinases (SVMPs) that play key roles in subduing prey in many vipers. Genomic analyses of several rattlesnake (Crotalus) species revealed the SVMP family massively expanded from a single, deeply conserved adam28 disintegrin and metalloproteinase gene, to as many as 31 tandem genes in the Western Diamondback rattlesnake (Crotalus atrox) through a number of single gene and multigene duplication events. Furthermore, we identified a series of stepwise intragenic deletions that occurred at different times in the course of gene family expansion and gave rise to the three major classes of secreted SVMP toxins by sequential removal of a membrane-tethering domain, the cysteine-rich domain, and a disintegrin domain, respectively. Finally, we show that gene deletion has further shaped the SVMP complex within rattlesnakes, creating both fusion genes and substantially reduced gene complexes. These results indicate that gene duplication and intragenic deletion played essential roles in the origin and diversification of these novel biochemical weapons.

A major role for noncoding regulatory mutations in the evolution of enzyme activity.

The quantitative evolution of protein activity is a common phenomenon, yet we know little about any general mechanistic tendencies that underlie it. For example, an increase (or decrease) in enzyme activity may evolve from changes in protein sequence that alter specific activity, or from changes in gene expression that alter the amount of protein produced. The latter in turn could arise via mutations that affect gene transcription, posttranscriptional processes, or copy number. Here, to determine the types of genetic changes underlying the quantitative evolution of protein activity, we dissected the basis of ecologically relevant differences in Alcohol dehydrogenase (Adh) enzyme activity between and within several Drosophila species. By using recombinant Adh transgenes to map the functional divergence of ADH enzyme activity in vivo, we find that amino acid substitutions explain only a minority (0 to 25%) of between- and within-species differences in enzyme activity. Instead, noncoding substitutions that occur across many parts of the gene (enhancer, promoter, and 5' and 3' untranslated regions) account for the majority of activity differences. Surprisingly, one substitution in a transcriptional Initiator element has occurred in parallel in two species, indicating that core promoters can be an important natural source of the tuning of gene activity. Furthermore, we show that both regulatory and coding substitutions contribute to fitness (resistance to ethanol toxicity). Although qualitative changes in protein specificity necessarily derive from coding mutations, these results suggest that regulatory mutations may be the primary source of quantitative changes in protein activity, a possibility overlooked in most analyses of protein evolution.

Birth-and-Death Evolution of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Diversification of Cuticular Hydrocarbon Synthesis in Drosophila.

The birth-and-death evolutionary model proposes that some members of a multigene family are phylogenetically stable and persist as a single copy over time, whereas other members are phylogenetically unstable and undergo frequent duplication and loss. Functional studies suggest that stable genes are likely to encode essential functions, whereas rapidly evolving genes reflect phenotypic differences in traits that diverge rapidly among species. One such class of rapidly diverging traits are insect cuticular hydrocarbons (CHCs), which play dual roles in chemical communications as short-range recognition pheromones as well as protecting the insect from desiccation. Insect CHCs diverge rapidly between related species leading to ecological adaptation and/or reproductive isolation. Because the CHC and essential fatty acid biosynthetic pathways share common genes, we hypothesized that genes involved in the synthesis of CHCs would be evolutionary unstable, whereas those involved in fatty acid-associated essential functions would be evolutionary stable. To test this hypothesis, we investigated the evolutionary history of the fatty acyl-CoA reductases (FARs) gene family that encodes enzymes in CHC synthesis. We compiled a unique data set of 200 FAR proteins across 12 Drosophila species. We uncovered a broad diversity in FAR content which is generated by gene duplications, subsequent gene losses, and alternative splicing. We also show that FARs expressed in oenocytes and presumably involved in CHC synthesis are more unstable than FARs from other tissues. Taken together, our study provides empirical evidence that a comparative approach investigating the birth-and-death evolution of gene families can identify candidate genes involved in rapidly diverging traits between species.

Extremely Divergent Haplotypes in Two Toxin Gene Complexes Encode Alternative Venom Types within Rattlesnake Species.

Natural selection is generally expected to favor one form of a given trait within a population. The presence of multiple functional variants of traits involved in activities such as feeding, reproduction, or the defense against predators is relatively uncommon within animal species. The genetic architecture and evolutionary mechanisms underlying the origin and maintenance of such polymorphisms are of special interest. Among rattlesnakes, several instances of the production of biochemically distinct neurotoxic or hemorrhagic venom types within the same species are known. Here, we investigated the genetic basis of this phenomenon in three species and found that neurotoxic and hemorrhagic individuals of the same species possess markedly different haplotypes at two toxin gene complexes. For example, neurotoxic and hemorrhagic Crotalus scutulatus individuals differ by 5 genes at the phospholipase A2 (PLA2) toxin gene complex and by 11 genes at the metalloproteinase (MP) gene complex. A similar set of extremely divergent haplotypes also underlies alternate venom types within C. helleri and C. horridus. We further show that the MP and PLA2 haplotypes of neurotoxic C. helleri appear to have been acquired through hybridization with C. scutulatus-a rare example of the horizontal transfer of a potentially highly adaptive suite of genes. These large structural variants appear analogous to immunity gene complexes in host-pathogen arms races and may reflect the impact of balancing selection at the PLA2 and MP complexes for predation on different prey.

The Deep Origin and Recent Loss of Venom Toxin Genes in Rattlesnakes.

The genetic origin of novel traits is a central but challenging puzzle in evolutionary biology. Among snakes, phospholipase A2 (PLA2)-related toxins have evolved in different lineages to function as potent neurotoxins, myotoxins, or hemotoxins. Here, we traced the genomic origin and evolution of PLA2 toxins by examining PLA2 gene number, organization, and expression in both neurotoxic and non-neurotoxic rattlesnakes. We found that even though most North American rattlesnakes do not produce neurotoxins, the genes of a specialized heterodimeric neurotoxin predate the origin of rattlesnakes and were present in their last common ancestor (∼22 mya). The neurotoxin genes were then deleted independently in the lineages leading to the Western Diamondback (Crotalus atrox) and Eastern Diamondback (C. adamanteus) rattlesnakes (∼6 mya), while a PLA2 myotoxin gene retained in C. atrox was deleted from the neurotoxic Mojave rattlesnake (C. scutulatus; ∼4 mya). The rapid evolution of PLA2 gene number appears to be due to transposon invasion that provided a template for non-allelic homologous recombination.

Expression of tandem gene duplicates is often greater than twofold.

Tandem gene duplication is an important mutational process in evolutionary adaptation and human disease. Hypothetically, two tandem gene copies should produce twice the output of a single gene, but this expectation has not been rigorously investigated. Here, we show that tandem duplication often results in more than double the gene activity. A naturally occurring tandem duplication of the Alcohol dehydrogenase (Adh) gene exhibits 2.6-fold greater expression than the single-copy gene in transgenic Drosophila This tandem duplication also exhibits greater activity than two copies of the gene in trans, demonstrating that it is the tandem arrangement and not copy number that is the cause of overactivity. We also show that tandem duplication of an unrelated synthetic reporter gene is overactive (2.3- to 5.1-fold) at all sites in the genome that we tested, suggesting that overactivity could be a general property of tandem gene duplicates. Overactivity occurs at the level of RNA transcription, and therefore tandem duplicate overactivity appears to be a previously unidentified form of position effect. The increment of surplus gene expression observed is comparable to many regulatory mutations fixed in nature and, if typical of other genomes, would shape the fate of tandem duplicates in evolution.

Gain of cis-regulatory activities underlies novel domains of wingless gene expression in Drosophila.

Changes in gene expression during animal development are largely responsible for the evolution of morphological diversity. However, the genetic and molecular mechanisms responsible for the origins of new gene-expression domains have been difficult to elucidate. Here, we sought to identify molecular events underlying the origins of three novel features of wingless (wg) gene expression that are associated with distinct pigmentation patterns in Drosophila guttifera. We compared the activity of cis-regulatory sequences (enhancers) across the wg locus in D. guttifera and Drosophila melanogaster and found strong functional conservation among the enhancers that control similar patterns of wg expression in larval imaginal discs that are essential for appendage development. For pupal tissues, however, we found three novel wg enhancer activities in D. guttifera associated with novel domains of wg expression, including two enhancers located surprisingly far away in an intron of the distant Wnt10 gene. Detailed analysis of one enhancer (the vein-tip enhancer) revealed that it overlapped with a region controlling wg expression in wing crossveins (crossvein enhancer) in D. guttifera and other species. Our results indicate that one novel domain of wg expression in D. guttifera wings evolved by co-opting pre-existing regulatory sequences governing gene activity in the developing wing. We suggest that the modification of existing enhancers is a common path to the evolution of new gene-expression domains and enhancers.

Wax, sex and the origin of species: Dual roles of insect cuticular hydrocarbons in adaptation and mating.

Evolutionary changes in traits that affect both ecological divergence and mating signals could lead to reproductive isolation and the formation of new species. Insect cuticular hydrocarbons (CHCs) are potential examples of such dual traits. They form a waxy layer on the cuticle of the insect to maintain water balance and prevent desiccation, while also acting as signaling molecules in mate recognition and chemical communication. Because the synthesis of these hydrocarbons in insect oenocytes occurs through a common biochemical pathway, natural or sexual selection on one role may affect the other. In this review, we explore how ecological divergence in insect CHCs can lead to divergence in mating signals and reproductive isolation. We suggest that the evolution of insect CHCs may be ripe models for understanding ecological speciation.

A single gene affects both ecological divergence and mate choice in Drosophila.

Evolutionary changes in traits involved in both ecological divergence and mate choice may produce reproductive isolation and speciation. However, there are few examples of such dual traits, and the genetic and molecular bases of their evolution have not been identified. We show that methyl-branched cuticular hydrocarbons (mbCHCs) are a dual trait that affects both desiccation resistance and mate choice in Drosophila serrata. We identify a fatty acid synthase mFAS (CG3524) responsible for mbCHC production in Drosophila and find that expression of mFAS is undetectable in oenocytes (cells that produce CHCs) of a closely related, desiccation-sensitive species, D. birchii, due in part to multiple changes in cis-regulatory sequences of mFAS. We suggest that ecologically influenced changes in the production of mbCHCs have contributed to reproductive isolation between the two species.

Sexually dimorphic regulation of the Wingless morphogen controls sex-specific segment number in Drosophila.

Sexual dimorphism is widespread throughout the metazoa and plays important roles in mate recognition and preference, sex-based niche partitioning, and sex-specific coadaptation. One notable example of sex-specific differences in insect body morphology is presented by the higher diptera, such as Drosophila, in which males develop fewer abdominal segments than females. Because diversity in segment number is a distinguishing feature of major arthropod clades, it is of fundamental interest to understand how different numbers of segments can be generated within the same species. Here we show that sex-specific and segment-specific regulation of the Wingless (Wg) morphogen underlies the development of sexually dimorphic adult segment number in Drosophila. Wg expression is repressed in the developing terminal male abdominal segment by the combination of the Hox protein Abdominal-B (Abd-B) and the sex-determination regulator Doublesex (Dsx). The subsequent loss of the terminal male abdominal segment during pupation occurs through a combination of developmental processes including segment compartmental transformation, apoptosis, and suppression of cell proliferation. Furthermore, we show that ectopic expression of Wg is sufficient to rescue this loss. We propose that dimorphic Wg regulation, in concert with monomorphic segment-specific programmed cell death, are the principal mechanisms of sculpting the sexually dimorphic abdomen of Drosophila.

Evolutionary origin of a novel gene expression pattern through co-option of the latent activities of existing regulatory sequences.

Spatiotemporal changes in gene expression underlie many evolutionary novelties in nature. However, the evolutionary origins of novel expression patterns, and the transcriptional control elements ("enhancers") that govern them, remain unclear. Here, we sought to explore the molecular genetic mechanisms by which new enhancers arise. We undertook a survey of closely related Drosophila species to identify recently evolved novel gene expression patterns and traced their evolutionary history. Analyses of gene expression in a variety of developing tissues of the Drosophila melanogaster species subgroup revealed high rates of expression pattern divergence, including numerous evolutionary losses, heterochronic shifts, and expansions or contractions of expression domains. However, gains of novel expression patterns were much less frequent. One gain was observed for the Neprilysin-1 (Nep1) gene, which has evolved a unique expression pattern in optic lobe neuroblasts of Drosophila santomea. Dissection of the Nep1 cis-regulatory region localized a newly derived optic lobe enhancer activity to a region of an intron that has accumulated a small number of mutations. The Nep1 optic lobe enhancer overlaps with other enhancer activities, from which the novel activity was co-opted. We suggest that the novel optic lobe enhancer evolved by exploiting the cryptic activity of extant regulatory sequences, and this may reflect a general mechanism whereby new enhancers evolve.

Generation of a novel wing colour pattern by the Wingless morphogen.

The complex, geometric colour patterns of many animal bodies have important roles in behaviour and ecology. The generation of certain patterns has been the subject of considerable theoretical exploration, however, very little is known about the actual mechanisms underlying colour pattern formation or evolution. Here we have investigated the generation and evolution of the complex, spotted wing pattern of Drosophila guttifera. We show that wing spots are induced by the Wingless morphogen, which is expressed at many discrete sites that are specified by pre-existing positional information that governs the development of wing structures. Furthermore, we demonstrate that the elaborate spot pattern evolved from simpler schemes by co-option of Wingless expression at new sites. This example of a complex design developing and evolving by the layering of new patterns on pre-patterns is likely to be a general theme in other animals.