The best of all worlds or the best possible world? Developmental constraint in the evolution of beta-tubulin and the sperm tail axoneme

Characterization of the Drosophila suzukii β2-tubulin gene and the utilization of its promoter to monitor sex separation and insemination

The Drosophila melanogaster β2-tubulin gene (Dm-β2t) controls the function of microtubules in the testis and sperm, and has been evaluated for use in biocontrol strategies based on the sterile insect technique, including sexing and the induction of male sterility. The spotted-wing Drosophila (Drosophila suzukii) is native to eastern Asia but has spread globally as an invasive pest of fruit crops, so biocontrol strategies are urgently required for this species. We therefore isolated the β2tubulin ortholog Ds-β2t from the USA laboratory strain of D. suzukii and confirmed the presence of functional motifs by aligning orthologs from multiple insects. The developmental expression profile of Ds-β2t was determined by RT-PCR using gene-specific primers and was similar to that of Dm-β2t. We then isolated the Ds-β2t promoter and used it to generate transgenic strains expressing a testis-specific fluorescent protein starting from the thirdinstar larvae. Efficient sexing was achieved based on fluorescence detection, and the transgenic males showed a similar survival rate to wild-type males. Fluorescence imaging and PCR were also used to confirm the insemination of wild-type females by transgenic males. We therefore confirm that D. suzukii strains expressing fluorescent markers under the control of the Ds-β2t promoter can be used for sexing and the confirmation of mating, and we discuss the wider potential of the Ds-β2t promoter in the context of genetic control strategies for D. suzukii.

Branchial cilia and sperm flagella recruit distinct axonemal components

Eukaryotic cilia and flagella have highly conserved 9 + 2 structures. They are functionally diverged to play cell-type-specific roles even in a multicellular organism. Although their structural components are therefore believed to be common, few studies have investigated the molecular diversity of the protein components of the cilia and flagella in a single organism. Here we carried out a proteomic analysis and compared protein components between branchial cilia and sperm flagella in a marine invertebrate chordate, Ciona intestinalis. Distinct feature of protein recruitment in branchial cilia and sperm flagella has been clarified; (1) Isoforms of α- and β-tubulins as well as those of actins are distinctly used in branchial cilia or sperm flagella. (2) Structural components, such as dynein docking complex, tektins and an outer dense fiber protein, are used differently by the cilia and flagella. (3) Sperm flagella are specialized for the cAMP- and Ca2+-dependent regulation of outer arm dynein and for energy metabolism by glycolytic enzymes. Our present study clearly demonstrates that flagellar or ciliary proteins are properly recruited according to their function and stability, despite their apparent structural resemblance and conservation.

Ultrastructure of the sperm axoneme and molecular analysis of axonemal dynein in Ephemeroptera (Insecta)

The Ephemeroptera sperm axoneme is devoid of outer dynein arms (ODA) and exhibits a pronounced modification of the central pair complex (CPC), which is substituted by the central sheath (CS): a tubular element of unknown molecular composition. We performed a detailed ultrastructural analysis of sperm axonemes in the genera Cloeon and Ecdyonurus using quick-freeze, deep-etch electron microscopy, showing that the loss of the conventional CPC is not only concomitant with the loss of ODA, but also with a substantial modification in the longitudinal distribution of both radial spokes (RS) and inner dynein arms (IDA). Such structures are no longer distributed following the alternation of different repeats as in the 9 + 2 axoneme, but instead share a 32 nm longitudinal repeat: a multiple of the 8 nm repeat observed along the CS wall. Differently from the conventional CPC, the CS and the surrounding RS possess a ninefold symmetry, coherently with the three-dimensional pattern of motility observed in Cloeon free spermatozoa. Biochemical analyses revealed that ultrastructural modifications are concomitant with a reduced complexity of the IDA heavy chain complement. We propose that these structural and molecular modifications might be related to the relief from the evolutionary constraints imposed by the CPC on the basal 9 + 9 + 2 axoneme and could also represent the minimal set compatible with flagellar beating and progressive motility mechanically regulated as suggested by the geometric clutch hypothesis. © 2014 Wiley Periodicals, Inc.

Distinct functional roles of β-tubulin isotypes in microtubule arrays of Tetrahymena thermophila, a model single-celled organism

BACKGROUND

The multi-tubulin hypothesis proposes that each tubulin isotype performs a unique role, or subset of roles, in the universe of microtubule function(s). To test this hypothesis, we are investigating the functions of the recently discovered, noncanonical β-like tubulins (BLTs) of the ciliate, Tetrahymena thermophila. Tetrahymena forms 17 distinct microtubular structures whose assembly had been thought to be based on single α- and β-isotypes. However, completion of the macronuclear genome sequence of Tetrahymena demonstrated that this ciliate possessed a β-tubulin multigene family: two synonymous genes (BTU1 and BTU2) encode the canonical β-tubulin, BTU2, and six genes (BLT1-6) yield five divergent β-tubulin isotypes. In this report, we examine the structural features and functions of two of the BLTs (BLT1 and BLT4) and compare them to those of BTU2.

METHODOLOGY/PRINCIPAL FINDINGS

With respect to BTU2, BLT1 and BLT4 had multiple sequence substitutions in their GTP-binding sites, in their interaction surfaces, and in their microtubule-targeting motifs, which together suggest that they have specialized functions. To assess the roles of these tubulins in vivo, we transformed Tetrahymena with expression vectors that direct the synthesis of GFP-tagged versions of the isotypes. We show that GFP-BLT1 and GFP-BLT4 were not detectable in somatic cilia and basal bodies, whereas GFP-BTU2 strongly labeled these structures. During cell division, GFP-BLT1 and GFP-BLT4, but not GFP-BTU2, were incorporated into the microtubule arrays of the macronucleus and into the mitotic apparatus of the micronucleus. GFP-BLT1 also participated in formation of the microtubules of the meiotic apparatus of the micronucleus during conjugation. Partitioning of the isotypes between nuclear and ciliary microtubules was confirmed biochemically.

CONCLUSION/SIGNIFICANCE

We conclude that Tetrahymena uses a family of distinct β-tubulin isotypes to construct subsets of functionally different microtubules, a result that provides strong support for the multi-tubulin hypothesis.

Drosophila Dynein intermediate chain gene, Dic61B, is required for spermatogenesis

This study reports the identification and characterization of a novel gene, Dic61B, required for male fertility in Drosophila. Complementation mapping of a novel male sterile mutation, ms21, isolated in our lab revealed it to be allelic to CG7051 at 61B1 cytogenetic region, since two piggyBac insertion alleles, CG7051(c05439) and CG7051(f07138) failed to complement. CG7051 putatively encodes a Dynein intermediate chain. All three mutants, ms21, CG7051(c05439) and CG7051(f07138), exhibited absolute recessive male sterility with abnormally coiled sperm axonemes causing faulty sperm individualization as revealed by Phalloidin staining in Don Juan-GFP background. Sequencing of PCR amplicons uncovered two point mutations in ms21 allele and confirmed the piggyBac insertions in CG7051(c05439) and CG7051(f07138) alleles to be in 5'UTR and 4(th) exon of CG7051 respectively, excision of which reverted the male sterility. In situ hybridization to polytene chromosomes demonstrated CG7051 to be a single copy gene. RT-PCR of testis RNA revealed defective splicing of the CG7051 transcripts in mutants. Interestingly, expression of cytoplasmic dynein intermediate chain, α, β, γ tubulins and α-spectrin was normal in mutants while ultra structural studies revealed defects in the assembly of sperm axonemes. Bioinformatics further highlighted the homology of CG7051 to axonemal dynein intermediate chain of various organisms, including DNAI1 of humans, mutations in which lead to male sterility due to immotile sperms. Based on these observations we conclude that CG7051 encodes a novel axonemal dynein intermediate chain essential for male fertility in Drosophila and rename it as Dic61B. This is the first axonemal Dic gene of Drosophila to be characterized at molecular level and shown to be required for spermatogenesis.

Tubulin evolution in insects: gene duplication and subfunctionalization provide specialized isoforms in a functionally constrained gene family

Background

The completion of 19 insect genome sequencing projects spanning six insect orders provides the opportunity to investigate the evolution of important gene families, here tubulins. Tubulins are a family of eukaryotic structural genes that form microtubules, fundamental components of the cytoskeleton that mediate cell division, shape, motility, and intracellular trafficking. Previous in vivo studies in Drosophila find a stringent relationship between tubulin structure and function; small, biochemically similar changes in the major alpha 1 or testis-specific beta 2 tubulin protein render each unable to generate a motile spermtail axoneme. This has evolutionary implications, not a single non-synonymous substitution is found in beta 2 among 17 species of Drosophila and Hirtodrosophila flies spanning 60 Myr of evolution. This raises an important question, How do tubulins evolve while maintaining their function? To answer, we use molecular evolutionary analyses to characterize the evolution of insect tubulins.

Results

Sixty-six alpha tubulins and eighty-six beta tubulin gene copies were retrieved and subjected to molecular evolutionary analyses. Four ancient clades of alpha and beta tubulins are found in insects, a major isoform clade (alpha 1, beta 1) and three minor, tissue-specific clades (alpha 2-4, beta 2-4). Based on a Homarus americanus (lobster) outgroup, these were generated through gene duplication events on major beta and alpha tubulin ancestors, followed by subfunctionalization in expression domain. Strong purifying selection acts on all tubulins, yet maximum pairwise amino acid distances between tubulin paralogs are large (0.464 substitutions/site beta tubulins, 0.707 alpha tubulins). Conversely orthologs, with the exception of reproductive tissue isoforms, show little sequence variation except in the last 15 carboxy terminus tail (CTT) residues, which serve as sites for post-translational modifications (PTMs) and interactions with microtubule-associated proteins. CTT residues overwhelming comprise the co-evolving residues between Drosophila alpha 2 and beta 3 tubulin proteins, indicating CTT specializations can be mediated at the level of the tubulin dimer. Gene duplications post-dating separation of the insect orders are unevenly distributed, most often appearing in major alpha 1 and minor beta 2 clades. More than 40 introns are found in tubulins. Their distribution among tubulins reveals that insertion and deletion events are common, surprising given their potential for disrupting tubulin coding sequence. Compensatory evolution is found in Drosophila beta 2 tubulin cis-regulation, and reveals selective pressures acting to maintain testis expression without the use of previously identified testis cis-regulatory elements.

Conclusion

Tubulins have stringent structure/function relationships, indicated by strong purifying selection, the loss of many gene duplication products, alpha-beta co-evolution in the tubulin dimer, and compensatory evolution in beta 2 tubulin cis-regulation. They evolve through gene duplication, subfunctionalization in expression domain and divergence of duplication products, largely in CTT residues that mediate interactions with other proteins. This has resulted in the tissue-specific minor insect isoforms, and in particular the highly diverse α3, α4, and β2 reproductive tissue-specific tubulin isoforms, illustrating that even a highly conserved protein family can participate in the adaptive process and respond to sexual selection.

Cooperativity between the beta-tubulin carboxy tail and the body of the molecule is required for microtubule function

Using Drosophila spermatogenesis as a model, we show that function of the beta-tubulin C-terminal tail (CTT) is not independent of the body of the molecule. For optimal microtubule function, the beta-tubulin CTT and body must match. beta2 is the only beta-tubulin used in meiosis and spermatid differentiation. beta1-tubulin is used in basal bodies, but beta1 cannot replace beta2. However, when beta1 is co-expressed with beta2, both beta-tubulins are equally incorporated into all microtubules, and males exhibit near wild type fertility. In contrast, co-expression of beta2beta1C and beta1beta2C, two reciprocal chimeric molecules with bodies and tails swapped, results in defects in meiosis, cytoskeletal microtubules, and axonemes; males produce few functional sperm and few or no progeny. In these experiments, all the same beta-tubulin parts are present, but unlike the co-assembled native beta-tubulins, the "trans" configuration of the co-assembled chimeras is poorly functional. Our data thus reveal essential intra-molecular interactions between the CTT and other parts of the beta-tubulin molecule, even though the CTT is a flexible surface feature of tubulin heterodimers and microtubules. In addition, we show that Drosophila sperm tail length depends on the total tubulin pool available for axoneme assembly and spermatid elongation. D. melanogaster and other Drosophila species have extraordinarily long sperm tails, the length of which is remarkably constant in wild type flies. We show that in males of experimental genotypes that express wild type tubulins but have half the amount of the normal tubulin pool size, sperm tails are substantially shorter than wild type.

Centrioles to basal bodies in the spermiogenesis of Mastotermes darwiniensis (Insecta, Isoptera)

In addition to their role in centrosome organization, the centrioles have another distinct function as basal bodies for the formation of cilia and flagella. Centriole duplication has been reported to require two alternate assembly pathways: template or de novo. Since spermiogenesis in the termite Mastotermes darwiniensis lead to the formation of multiflagellate sperm, this process represents a useful model system in which to follow basal body formation and flagella assembly. We present evidence of a possible de novo pathway for basal body formation in the differentiating germ cell. This cell also contains typical centrosomal proteins, such as centrosomin, pericentrin-like protein, gamma-tubulin, that undergo redistribution as spermatid differentiation proceeds. The spermatid centrioles are long structures formed by nine doublet rather than triplet microtubules provided with short projections extending towards the surrounding cytoplasm and with links between doublets. The sperm basal bodies are aligned in parallel beneath the nucleus. They consist of long regions close to the nucleus showing nine doublets in a cartwheel array devoid of any projections; on the contrary, the short region close to the plasma membrane, where the sperm flagella emerge, is characterized by projections similar to those observed in the centrioles linking the basal body to the plasma membrane. It is hypothesized that this appearance is in connection with the centriole elongation and further with the flagellar axonemal organization. Microtubule doublets of sperm flagellar axonemes are provided with outer dynein arms, while inner arms are rarely visible.

Axoneme-dependent tubulin modifications in singlet microtubules of the Drosophila sperm tail

Drosophila melanogaster sperm tubulins are posttranslationally glutamylated and glycylated. We show here that axonemes are the substrate for these tubulin C-terminal modifications. Axoneme architecture is required, but full length, motile axonemes are not necessary. Tubulin glutamylation occurs during or shortly after assembly into the axoneme; only glutamylated tubulins are glycylated. Tubulins in other testis microtubules are not modified. Only a small subset of total Drosophila sperm axoneme tubulins have these modifications. Biochemical fractionation of Drosophila sperm showed that central pair and accessory microtubules have the majority of poly-modified tubulins, whereas doublet microtubules have only small amounts of mono- and oligo-modified tubulins. Glutamylation patterns for different beta-tubulins experimentally assembled into axonemes were consistent with utilization of modification sites corresponding to those identified in other organisms, but surrounding sequence context was also important. We compared tubulin modifications in the 9 + 9 + 2 insect sperm tail axonemes of Drosophila with the canonical 9 + 2 axonemes of sperm of the sea urchin Lytichinus pictus and the 9 + 0 motile sperm axonemes of the eel Anguilla japonica. In contrast to Drosophila sperm, L. pictus sperm have equivalent levels of modified tubulins in both doublet and central pair microtubule fractions, whereas the doublets of A. japonica sperm exhibit little glutamylation but extensive glycylation. Tubulin C-terminal modifications are a prevalent feature of motile axonemes, but there is no conserved pattern for placement or amount of these

New insights into the cell biology of insect axonemes

Insects do not possess ciliated epithelia, and cilia/flagella are present in the sperm tail and--as modified cilia--in mechano- and chemosensory neurons. The core cytoskeletal component of these organelles, the axoneme, is a microtubule-based structure that has been conserved throughout evolution. However, in insects the sperm axoneme exhibits distinctive structural features; moreover, several insect groups are characterized by an unusual sperm axoneme variability. Besides the abundance of morphological data on insect sperm flagella, most of the available molecular information on the insect axoneme comes from genetic studies on Drosophila spermatogenesis, and only recently other insect species have been proposed as useful models. Here, we review the current knowledge on the cell biology of insect axoneme, including contributions from both Drosophila and other model insects.

Sperm structure of Limoniidae and their phylogenetic relationship with Tipulidae (Diptera, Nematocera)

The sperm ultrastructure of a few species of Limoniidae (Limonia nigropunctata; L. nubeculosa; Chionea n. sp.; C. alpina; C. lutescens) was studied. The two species of Limonia have a monolayered acrosome with crystallized material, a three-lobed nucleus in cross section, a ring of centriole adjunct material and a flagellum which consists of a 9+9+1 axoneme and a single mitochondrial derivative. The central axonemal tubule is provided with 15 protofilaments in its tubular wall, while the accessory tubules have 13 protofilaments and are flanked by the electron-dense intertubular material. The three species of Chionea share a monolayered acrosome, a nucleus with two longitudinal grooves, a centriole adjunct material which surrounds the centriole and the initial part of the axoneme. The axoneme is of conventional type, with 9+9+2 microtubular pattern, with accessory tubules provided with 13 protofilaments and intertubular material. However, in C. lutescens the accessory tubules start with 15 protofilaments and transform into a tubule with 13 protofilaments. These data are discussed in the light of the phylogenetic relationship between Limoniidae and Tipulidae. For this purpose, the sperm ultrastructure of Nephrotoma appendiculata was also considered comparatively.

Ciliary tubulin and its post-translational modifications

Tubulin, the most abundant axonemal protein, is extensively modified by several highly conserved post-translational mechanisms including acetylation, detyrosination, glutamylation, and glycylation. We discuss the pathways that contribute to the assembly and maintenance of axonemal microtubules, with emphasis on the potential functions of post-translational modifications that affect tubulin. The recent identification of a number of tubulin modifying enzymes and mutational studies of modification sites on tubulin have allowed for significant functional insights. Polymeric modifications of tubulin (glutamylation and glycylation) have emerged as important determinants of the 9 + 2 axoneme assembly and motility.

Evolutionary conservation of lampbrush-like loops in drosophilids

Background

Loopin-1 is an abundant, male germ line specific protein of Drosophila melanogaster. The polyclonal antibody T53-F1 specifically recognizes Loopin-1 and enables its visualization on the Y-chromosome lampbrush-like loop named kl-3 during primary spermatocyte development, as well as on sperm tails. In order to test lampbrush-like loop evolutionary conservation, extensive phase-contrast microscopy and immunostaining with T53-F1 antibody was performed in other drosophilids scattered along their genealogical tree.

Results

In the male germ line of all species tested there are cells showing giant nuclei and intranuclear structures similar to those of Drosophila melanogaster primary spermatocytes. Moreover, the antibody T53-F1 recognizes intranuclear structures in primary spermatocytes of all drosophilids analyzed. Interestingly, the extent and conformation of the staining pattern is species-specific. In addition, the intense staining of sperm tails in all species suggests that the terminal localization of Loopin-1 and its orthologues is conserved. A comparison of these cytological data and the data coming from the literature about sperm length, amount of sperm tail entering the egg during fertilization, shape and extent of both loops and primary spermatocyte nuclei, seems to exclude direct relationships among these parameters.

Conclusion

Taken together, the data reported strongly suggest that lampbrush-like loops are a conserved feature of primary spermatocyte nuclei in many, if not all, drosophilids. Moreover, the conserved pattern of the T53-F1 immunostaining indicates that a Loopin-1-like protein is present in all the species analyzed, whose localization on lampbrush-like loops and sperm tails during spermatogenesis is evolutionary conserved.

Cells, molecules and morphogenesis: the making of the vertebrate ear

The development and evolution of mechanosensory cells and the vertebrate ear is reviewed with an emphasis on delineating the cellular, molecular and developmental basis of these changes. Outgroup comparisons suggests that mechanosensory cells are ancient features of multicellular organisms. Molecular evidence suggests that key genes involved in mechanosensory cell function and development are also conserved among metazoans. The divergent morphology of mechanosensory cells across phyla is interpreted here as 'deep molecular homology' that was in parallel shaped into different forms in each lineage. The vertebrate mechanosensory hair cell and its associated neuron are interpreted as uniquely derived features of vertebrates. It is proposed that the vertebrate otic placode presents a unique embryonic adaptation in which the diffusely distributed ancestral mechanosensory cells became concentrated to generate a large neurosensory precursor population. Morphogenesis of the inner ear is reviewed and shown to depend on genes expressed in and around the hindbrain that interact with the otic placode to define boundaries and polarities. These patterning genes affect downstream genes needed to maintain proliferation and to execute ear morphogenesis. We propose that fibroblast growth factors (FGFs) and their receptors (FGFRs) are a crucial central node to translate patterning into the complex morphology of the vertebrate ear. Unfortunately, the FGF and FGFR genes have not been fully analyzed in the many mutants with morphogenetic ear defects described thus far. Likewise, little information exists on the ear histogenesis and neurogenesis in many mutants. Nevertheless, a molecular mechanism is now emerging for the formation of the horizontal canal, an evolutionary novelty of the gnathostome ear. The existing general module mediating vertical canal growth and morphogenesis was modified by two sets of new genes: one set responsible for horizontal canal morphogenesis and another set for neurosensory formation of the horizontal crista and associated sensory neurons. The dramatic progress in deciphering the molecular basis of ear morphogenesis offers grounds for optimism for translational research toward intervention in human morphogenetic defects of the ear.

Tubulin glycylation and glutamylation deficiencies in unconventional insect axonemes

Though the 9+2 axonemal organization has generally been conserved throughout metazoan evolution, insect spermatozoa possess a substantial variety in axoneme ultrastructure, displaying different axonemal patterns. Therefore, insects provide a wide range of models that may be useful for the study of the mechanisms of axoneme assembly. We have used antibodies specific for glutamylated, monoglycylated, and polyglycylated tubulin to investigate the tubulin isoform content expressed in the unorthodox sperm axonemes of four insect species belonging to both of the superorders Palaeoptera and Neoptera. Each one of these axonemal models exhibits distinctive structural features, either showing the typical radial organization endowed with a ninefold symmetry or consisting of an helical arrangement with up to 200 microtubular doublets, but in all cases these axonemes share the absence of a microtubule central pair. Our results showed that all these atypical patterns are characterized by a dramatic decrease in both tubulin glycylation and glutamylation levels or even lack of both polymodifications. These data provide the first examples of a simultaneous extreme reduction or even absence of both polymodifications in axonemal tubulin. Given the unrelated positions of the analyzed species in the insect phylogenetic tree, this common feature is probably not due to evolutionary relationships. Therefore, our findings support the hypothesis of the existence of a correlation between the low level of polymodifications and the lack of a microtubule central pair in these peculiar insect flagellar axonemes, similarly as was previously proposed for cilia of Tetrahymena glycylation site mutants.

The proximal region of the β-tubulin C-terminal tail is sufficient for axoneme assembly

We have used Drosophila testis-specific beta2-tubulin to determine sequence requirements for different microtubules. The beta2-tubulin C-terminal tail has unique sperm-specific functions [Dev Biol 158:267-286 (2003)] and is also important for forming stable heterodimers with alpha-tubulin, a general function common to all microtubules [Mol Biol Cell 12(7):2185-2194 (2001)]. beta-tubulins utilized in motile 9 + 2 axonemes contain a C-terminal sequence "axoneme motif" [Science 275 (1997) 70-73]. C-terminal truncated beta2-tubulin cannot form the sperm tail axoneme. Here we show that a partially truncated beta2-tubulin (beta2Delta7) containing only the proximal portion of the C-terminal tail, including the axoneme motif, can support production of functional motile sperm. We conclude that these proximal eight amino acids specify the binding site for protein(s) essential to support assembly of the motile axoneme. Males that express beta2Delta7, although they are fertile, produce fewer sperm than wild type males. Beta2Delta7 causes a slightly increased error rate in spermatogenesis attributable to loss of stabilizing properties intrinsic to the full-length C-terminal tail. Therefore, beta2Delta7 males would be at a selective disadvantage and it is likely that the full-length C-terminus would be essential in the wild and in evolution.

Common pattern of evolution of gene expression level and protein sequence in Drosophila

Sequence divergence scaled by variation within species has been used to infer the action of selection upon individual genes. Applying this approach to expression, we compared whole-genome whole-body RNA levels in 10 heterozygous Drosophila simulans genotypes and a pooled sample of 10 D. melanogaster lines using Affymetrix Genechip. For 972 genes expressed in D. melanogaster, the transcript level was below detection threshold in D. simulans, which may be explained either by sequence divergence between the primers on the chip and the mRNA transcripts or by down-regulation of these genes. Out of 6,707 genes that were expressed in both species, transcript level was significantly different between species for 534 genes (at P < 0.001). Genes whose expression is under stabilizing selection should exhibit reduced genetic variation within species and reduced divergence between species. Expression of genes under directional selection in D. simulans should be highly divergent from D. melanogaster, while showing low genetic variation in D. simulans. Finally, the genes with large variation within species but modest divergence between species are candidates for balancing selection. Rapidly diverging, low-polymorphism genes included those involved in reproduction (e.g., Mst 3Ba, 98Cb; Acps 26Aa, 63F; and sperm-specific dynein). Genes with high variation in transcript abundance within species included metallothionein and hairless, both hypothesized to be segregating in nature because of gene-by-environment interactions. Further, we compared expression divergence and DNA substitution rate in 195 genes. Synonymous substitution rate and expression divergences were uncorrelated, whereas there was a significant positive correlation between nonsynonymous substitution rate and expression divergence. We hypothesize that as a substantial fraction of nonsynonymous divergence has been shown to be adaptive, much of the observed expression divergence is likewise adaptive.