Department of Ecology and Evolutionary Biology, University of Arizona
Correspondence: E-mail: jgood{at}email.arizona.edu.
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: Mus musculus pleiotropy sexual selection sperm development X chromosome
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
First, sexual selection may be driving the evolution of male reproductive proteins, either through male-male competition or through male-female interactions (Swanson and Vacquier 2002). In mice, spermatogenesis begins at birth and consists of approximately 11 days of mitotic cell growth, 10 days of meiosis, followed by 14 days during which haploid spermatids undergo a dramatic morphological transformation into sperm (Eddy 2002). If sexual selection on male gametes acts predominantly on mutations that influence the form and function of sperm, then genes expressed late in spermatogenesis might be expected to evolve more quickly and to show a greater signature of positive selection. In particular, genes involved in the morphological transformation from spermatids into sperm are likely targets for sexual selection.
Second, variation in evolutionary rate among reproductive genes might be due to differences in functional constraint rather than (or in addition to) the action of positive selection. Genes expressed early in development (Raff 1996) or in many different tissues (Duret and Mouchiroud 2000) are expected to exhibit more pleiotropic effects. Hence, genes that are expressed late in development or are tissue specific might be expected to evolve more rapidly on average. Some studies support these predictions while others do not. For example, in mammals, tissue-specific genes evolve more rapidly and exhibit a larger variance in evolutionary rate compared to widely expressed genes (Duret and Mouchiroud 2000; Winter, Goodstadt, and Ponting 2004). On the other hand, there is little evidence for a simple relationship between rates of molecular evolution and developmental timing of expression in Caenorhabditis elegans (Castillo-Davis and Hartl 2002, 2003; Cutter and Ward 2005).
Thus, hypotheses for the rapid evolution of male reproductive genes based on sexual selection or reduced functional constraint both predict that genes expressed late in spermatogenesis will evolve more quickly than early-expressed genes. However, sexual selection may lead to a ratio of nonsynonymous to synonymous nucleotide changes that is greater than one for some late-expressed genes. This signature of positive selection is not expected for genes under relaxed functional constraint.
A third way in which developmental timing of expression may be important for understanding the molecular evolution of male reproductive proteins relates to the genomic location of these genes. In mammals, recessive mutations that benefit males are expected to accumulate on the X chromosome because selection will be more effective in the hemizygous sex (Rice 1984; Charlesworth, Coyne, and Barton 1987). However, the X chromosome is inactivated in the male germline during meiosis, and this ought to select against X linkage for genes that are expressed late in spermatogenesis (Betrán, Thornton, and Long 2002; Wu and Xu 2003). Consistent with this, there appears to be an excess of X-linked genes in the mouse that are expressed early in spermatogenesis (Wang et al. 2001), while late-expressed genes are rarely on the X chromosome (Betrán, Thornton, and Long 2002; Khil et al. 2004). The overrepresentation of X-linked genes early in spermatogenesis might be caused by (1) an evolutionary shift towards earlier expression of genes that were always X linked (Wu and Xu 2003) or (2) an evolutionary increase in the number of X-linked genes with male-specific function, either through recruitment of genes from the autosomes or through the evolution of novel function at preexisting X-linked genes. The second hypothesis predicts that early-expressed, X-linked spermatogenic genes may show a signature of positive selection.
The recent publication of patterns of gene expression in testes and other tissues in the house mouse, Mus musculus (Schultz, Hamra, and Garbers 2003), presents an opportunity to test these predictions and ask whether rates of evolution depend on when genes are expressed during development. Here, we investigate the molecular evolution of genes expressed in the male germline of mice at nine developmental time points during spermatogenesis (Schultz, Hamra, and Garbers 2003). We discovered a significant positive correlation between developmental timing of expression and the rate of protein evolution. Moreover, many late-expressed, rapidly evolving genes were testes specific in their expression, and many showed evidence of positive selection in a maximum likelihood analysis. These observations are consistent with a strong role for sexual selection in driving the evolution of male reproductive genes in mice, particularly in later developmental stages.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Two approaches were used to associate genes with specific time points or phases of spermatogenesis based on their expression profiles. First, Schultz, Hamra, and Garbers (2003) applied a Gaussian clustering algorithm to their data and identified eight groups of genes based on overall similarity of expression profiles across all time points (ARRAYMINER 4.0, Optimal Designs, Brussels, Belgium). Here, we used these clusters to identify groups of genes whose expression increased or decreased at particular times during spermatogenesis. Second, for each gene we noted the day of highest expression. Expression data collected at day 60 were not considered here because spermatogenesis is cyclic and all stages overlap in adult testis (Eddy 2002).
Tissue-specific expression profiles were calculated using data from the Genomics Institute of the Novartis Research Foundation (GNF) SymAtlas (v0.8; symatlas.gnf.org/SymAtlas/), an Affymetrix array database generated by Su et al. (2004). Expression was quantified based on the average signal difference between match and mismatch probes. GNF probe numbers were matched to Ensembl genes using the University of California, Santa CRUZ (UCSC) Mouse Gene Sorter (genome.ucsc.edu). For each gene, we measured tissue specificity as the proportion of expression in the most highly expressed tissue relative to the total level of expression summed over all tissues (Winter, Goodstadt, and Ponting 2004). Using the same approach, we also calculated the proportion of expression in the testis
We define testis- and tissue-specific genes as those with
and
respectively. This definition corresponds to approximately four times the median
or
of all the genes in the GNF SymAtlas database (N = 36,118). To control for correlated expression across tissues we only considered 31 tissues previously reported to not have significantly similar expression profiles in the mouse (Winter, Goodstadt, and Ponting 2004). Our measure of testis specificity should be broadly equivalent to male germline specificity because the genes identified by Schultz, Hamra, and Garbers (2003) were enriched for transcripts expressed in spermatogenic cells and not in the somatic cells of the testis.
Rates of Molecular Evolution
Protein-coding domains for 987 mouse and rat one-to-one orthologs were obtained from the Ensembl database. For all evolutionary analyses we only considered annotated gene pairs with orthology based on best reciprocal BLAST hits. We translated each pair into amino acid sequences, aligned them in-frame using the PileUp routine (Wisconsin Package v10.2, Genetics Computer Group, Madison, Wis.), and back translated to nucleotide sequence. We then calculated the number of nonsynonymous substitutions per nonsynonymous site (KA) and the number of synonymous substitutions per synonymous site (KS) using the method of Li (1993) as implemented in the Diverge routine (Wisconsin Package v10.2, Genetics Computer Group). All pairs with fewer than 100 aligned codons were discarded. For genes with multiple transcripts we aligned all pairwise combinations and retained those pairs with the fewest number of ambiguous codons and/or the lowest KA/KS values. In all group comparisons we used R (www.r-project.org/) to construct 95% confidence intervals about group means and medians based on 10,000 bootstrap replicates. All other statistical tests were conducted using JMP (v5.0, SAS Institute Inc., Cary, N.C.) unless otherwise noted.
To test for positive selection we examined 387 mouse-rat-human orthologs (also based on Ensembl annotation) using a maximum likelihood framework. Genes considered in the mouse-rat pairwise comparisons were excluded from this analysis if they lacked an annotated human ortholog or were represented by more than one transcript in at least one species. First, orthologous peptide sequences were aligned with CLUSTALW using default parameters (v1.8; Thompson, Higgins, and Gibson 1994). Then nucleotide sequences were aligned in-frame with RevTrans (v1.3, written by R. Wernersson) using the amino acid alignment as a guide. Several model-based tests of positive selection have been proposed (Yang 2002). We chose to conduct a relatively conservative test using the codeml program in the PAML package (Yang 1997). For each gene tree, we fit the data to two models (models M7 and M8; Yang et al. 2000). M7 assumes that KA/KS (dN/dS) values follow a beta distribution constrained between 0 and 1, while the second model (M8) allows for an additional category of sites with KA/KS > 1. The likelihood of data under each model can be directly compared using a standard likelihood ratio test because M7 is a nested subset of M8. Positive selection is inferred if M8 provides a significantly greater fit when tested using a chi-square distribution. We chose a significance threshold of P < 0.01 to be conservative. M8 has been reported to be sensitive to starting values of
(
KA/KS) and tends to get stuck in local optima. We used an initial
of 1 to identify genes with a significantly better fit to M8. We then further evaluated M7 and M8 in this subset using an initial
value of 0.1 to check for starting point dependency. Only the highest likelihood scores are reported here.
Functional Annotation of Genes
We retrieve functional annotation of genes from Ensembl (Build 32; http://www.ensembl.org/Mus_musculus/)based on the Gene Ontology database (Ashburner et al. 2000). All levels of the annotation hierarchy were retrieved from each of the three Gene Ontology categories (i.e., biological process, cellular component, and molecular function).
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
|
|
|
|
|
We tested if early-expressed genes showed higher rates of evolution when X linked. The median KA/KS value for early-expressed X-linked genes (0.142, C.I. 0.0690.290, N = 21) was found to be slightly but not significantly higher than autosomal genes (0.118, C.I. 0.1040.133, N = 329).
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Pleiotropy and Constraint
There is weak evidence for reduced constraint in some of the genes examined here. For example, widely expressed genes showed higher rates of evolution late in spermatogenesis compared to early in spermatogenesis, although this difference was not significant (fig. 2B). In addition, tissue-specific and testis-specific genes showed higher rates of evolution than widely expressed genes, a pattern that might reflect differences in level of constraint, but also likely reflects differences in the strength of positive selection (see below).
Spermatogenesis in mice appears to be exceptionally sensitive to genetic disruption. This conclusion stems primarily from the observation that knockout phenotypes disproportionately involve male sterility (Escalier 2001), including a number of genetic targets not previously implicated in male reproduction (Tourtellotte et al. 1999; Xu et al. 1999). However, there is little evidence that early events in spermatogenesis are more vulnerable. Particularly delicate steps have been identified across all phases of spermatogenesis, including primordial germ cells and spermatogonia (premeiotic), primary spermatocytes (meiotic), and elongating spermatids (postmeiotic; Escalier 2001). Consistent with this, we found only a weak trend relating rates of protein evolution to developmental timing of expression independent of testis specificity.
In mammals, tissue-specific genes evolve more rapidly than widely expressed genes (Duret and Mouchiroud 2000). However, there is a large variance in rates of evolution across tissues, and tissue-specific biology may play a larger role than general models of pleiotropic constraint in explaining this pattern (Winter, Goodstadt, and Ponting 2004). For example, genes specific to the liver evolve more than three times as fast as those specific to the dorsal root ganglion in rodents (Winter, Goodstadt, and Ponting 2004). The correspondence between increased rates of protein evolution and increased testis specificity late in spermatogenesis was striking (table 2). Moreover, variation in rates of protein evolution depended more heavily on variation in testis specificity across spermatogenesis than overall levels of tissue specificity
illustrating the importance of testis-specific biology versus overall levels of constraint in these data. We found that testis-specific genes
expressed late in spermatogenesis have a median KA/KS (0.231, N = 216) of more than twice the genomic rate of 0.11 and over 70% higher than previous estimates for testis-specific genes (Winter, Goodstadt, and Ponting 2004).
Positive Selection and Spermatogenesis
The high rate of protein evolution in the later stages of spermatogenesis appears to be driven also by a higher prevalence of positive selection. Forty-four of the 50 most rapidly evolving genes (KA/KS 1.950.56, pairwise mouse-rat comparisons) were from late-expressed clusters. Likewise, a significant proportion of the genes found to be under positive selection within the subset we analyzed using maximum likelihood were from late-expressed clusters (table 4).
Why would positive selection be more frequent late in spermatogenesis? First, development of the male germline becomes increasingly specialized following the onset of meiosis, and most testis-specific genes are expressed late in spermatogenesis (table 2; Hecht 1993; Schultz, Hamra, and Garbers 2003). In general, there appears to be strong selection for the evolution of novel reproductive function in the mouse (Eddy 2002; Waterston et al. 2002). This increase in testis specificity might facilitate adaptive evolution if mutations that are beneficial when expressed in the testis are free from pleiotropic effects in other tissues.
Second, many of the gametic phenotypes that are likely to be targets of sexual selection, either due to sperm competition or interaction with the female reproductive tract, develop during the dramatic postmeiotic morphogenesis of round spermatids into spermatozoa. During this period, the nucleus is reshaped and condensed, the acrosome and flagellum are formed, and much of the cytoplasm is removed (Meistrich 1993). Genes in cluster 7 showed on average an almost 40-fold increase in expression during this transformation (days 2129, table 1). Furthermore, this cluster had one of the highest median rates of protein evolution (fig. 3), the largest proportion of very rapidly evolving genes (table S1), and the highest incidence of positive selection in our maximum likelihood analysis (table S2).
Some of these genes are thought to be involved in the formation of the sperm cell surface and are logical candidates for playing a role in fertilization (Schultz, Hamra, and Garbers 2003). For example, we identified multiple genes involved in protein binding and sperm-egg interactions that are rapidly evolving and/or positively selected (e.g., Spam, Spaca3, Zpr3; tables S1 and S2), consistent with the findings of previous studies (Torgerson, Kulathinal, and Singh 2002; Swanson, Nielsen, and Yang 2003; Castillo-Davis et al. 2004). Given their direct interaction with the female gamete, these proteins are logical targets for sexual selection through sperm competition and/or sexual conflict (Torgerson, Kulathinal, and Singh 2002). Interestingly, several of the rapidly evolving genes we identified in cluster 7 are involved in transcription, nuclear reorganization, and condensation (tables S1 and S2). Usually proteins involved in DNA packaging are highly conserved across animal taxa (Thatcher and Gorovsky 1994). However, nuclear reorganization and DNA condensation may play a key role in sperm competition and ultimately fertilization ability, and, therefore, genes underlying these phenotypes could be major targets of selection. Consistent with this, multiple proteins involved in DNA packaging in sperm appear to be under positive selection within primates (Wyckoff, Wang, and Wu 2000). Extensive divergence in sperm head morphology has been documented within mammals, even among closely related species (Breed 2004). However, the genotypic basis of this variation remains poorly known. Given the large number of testis-specific genes expressed during spermiogenesis, the potential number of selective targets could be quite large.
Finally, other forces might lead to strong positive selection on phenotypes that develop during spermatogenesis, independent of sexual selection. In particular, natural selection related to immune defense could contribute to divergence in reproductive genes. In mammals, genes involved in immunity often evolve rapidly (Waterston et al. 2002; Castillo-Davis et al. 2004). Furthermore, multiple genes involved in mammalian pathogen defense have been identified that are expressed in reproductive tissue and appear to be under positive selection (Maxwell, Morrison, and Dorin 2003; Lynn et al. 2004). However, it is unclear if testis-specific genes involved in immune defense are common enough to contribute to the general pattern of increased rates of molecular evolution across spermatogenesis that we observed.
Adaptive Evolution and the X Chromosome
Previous studies have reported that a disproportionately high number of X-linked genes are expressed in the mouse testis during early spermatogenesis (Wang et al. 2001; Khil et al. 2004) and that these genes evolve more quickly on average than comparable autosomal genes (Torgerson and Singh 2003). Both patterns are consistent with some theoretical predictions of the evolution of sexually selected genes (Rice 1984). In mammals, this phenomenon is expected to be mostly restricted to early spermatogenesis because X inactivation during early meiosis clearly selects against X linkage (Betrán, Thornton, and Long 2002; Wu and Xu 2003). However, we found little evidence of either rapid evolution or differential accumulation of early-expressed genes on the X chromosome. There was a weak, nonsignificant trend towards increased KA/KS in premeiotic X-linked genes. Moreover, only 3 of the 50 most rapidly evolving genes and 1 of the 29 genes identified with maximum likelihoodbased tests of selection were X linked. However, none of the eight X-linked genes considered by Torgerson and Singh (2003) were found within our set of 1,420 genes. We also found only a slight trend towards increased X linkage early in spermatogenesis. A recently proposed theoretical model that incorporates pleiotropy does not predict increased X linkage of sexually selected genes involved in male reproduction (Fitzpatrick 2004). As mentioned above, we found a much lower degree of testis specificity early in spermatogenesis. If this lack of testis specificity reduces the efficacy of selection, it is not surprising that genes expressed early in spermatogenesis identified in the current study do not conform to the classic predictions that do not consider pleiotropy (Rice 1984). For this and the developmental reasons mentioned above, the X chromosome seems a less likely target of sexual selection on genes expressed in the male germline.
However, our result is in contrast to a stronger bias of early testis-expressed genes towards the X reported in a recent analysis of numerous data sets (Wang et al. 2001; Khil et al. 2004), including some of the data collected by Schultz, Hamra, and Garbers (2003). A number of differences exist between our study and the work of Khil et al. (2004). For example, we required that all genes have at least a threefold change in expression during spermatogenesis and be absent in expression profiles from somatic tissues of the testis. Consequently, there appears to be surprisingly little overlap between the set of testis-expressed genes considered by Khil et al. (2004) and those in the current study. It also remains unclear if the previously described patterns are driven by sexual selection or adaptive evolution in general. Other functional classes presumably not influenced by sexual selection, such as genes expressed in muscle and the brain, have also been reported to be overrepresented on the X (Bortoluzzi et al. 1998; Zechner et al. 2001), and many aspects of the unusual genetic composition of this chromosome have yet to be fully explained (Vallender and Lahn 2004).
Conclusions
Rates of protein evolution and tissue specificity dramatically increased with timing of expression in mouse spermatogenesis. The highest rates of protein evolution were found in testis-specific genes expressed during the morphological transformation of round spermatids into spermatozoa. Conversely, genes expressed early in spermatogenesis were found to evolve at rates similar to the genome median. Tests of natural selection using a likelihood framework confirmed that at least part of this pattern was driven by a higher incidence of positive selection later in sperm development. Overall, these data suggest that the strength of positive selection commonly associated with the evolution of male gametes varies considerably across development and acts primarily on phenotypes that develop during the postmeiotic process of spermiogenesis.
![]() |
Supplementary Materials |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Table S2. Twenty-nine genes inferred to be under positive selection using a maximum likelihood framework.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Ashburner, M., C. A. Ball, J. A. Blake et al. (20 co-authors). 2000. Gene Ontology: tool for the unification of biology. Nat. Genet. 25:2529.[CrossRef][ISI][Medline]
Begun, D. J., P. Whitley, B. L. Todd, H. M. Waldrip-Dail, and A. G. Clark. 2000. Molecular population genetics of male accessory gland proteins in Drosophila. Genetics 156:18791888.
Betrán, E., K. Thornton, and M. Long. 2002. Retroposed new genes out of the X in Drosophila. Genome Res. 12:18541859.
Bortoluzzi, S., L. Rampoldi, B. Simionati et al. (13 co-authors). 1998. A comprehensive, high-resolution genomic transcript map of human skeletal muscle. Genome Res. 8:817825.
Breed, W. G. 2004. The spermatozoon of Eurasian murine rodents: its morphological diversity and evolution. J. Morphol. 261:5269.[CrossRef][ISI][Medline]
Castillo-Davis, C. I., and D. L. Hartl. 2002. Genome evolution and developmental constraint in Caenorhabditis elegans. Mol. Biol. Evol. 19:728735.
Castillo-Davis, C. I., and D. L. Hartl. 2003. Conservation, relocation and duplication in genome evolution. Trends Genet. 19:593597.[CrossRef][ISI][Medline]
Castillo-Davis, C. I., F. A. Kondrashov, D. L. Hartl, and R. J. Kulathinal. 2004. The functional genomic distribution of protein divergence in two animal phyla: coevolution, genomic conflict, and constraint. Genome Res. 14:802811.
Charlesworth, B., J. A. Coyne, and N. H. Barton. 1987. The relative rates of evolution of sex chromosomes and autosomes. Am. Nat. 130:113146.[CrossRef][ISI]
Cutter, A. D., and S. Ward. 2005. Sexual and temporal dynamics of molecular evolution in C. elegans development. Mol. Biol. Evol. 22:178188.
Duret, L., and D. Mouchiroud. 2000. Determinants of substitution rates in mammalian genes: expression pattern affects selection intensity but not mutation rate. Mol. Biol. Evol. 17:6874.
Eddy, E. M. 2002. Male germ cell gene expression. Recent Prog. Horm. Res. 57:103128.
Escalier, D. 2001. Impact of genetic engineering on the understanding of spermatogenesis. Hum. Reprod. Update 7:191210.
Fitzpatrick, M. J. 2004. Pleiotropy and the genomic location of sexually selected genes. Am. Nat. 163:800808.[CrossRef][ISI][Medline]
Gibbs, R. A., G. M. Weinstock, M. L. Metzker et al. (230 co-authors). 2004. Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature 428:493521.[CrossRef][ISI]
Glassey, B., and A. Civetta. 2004. Positive selection at reproductive ADAM genes with potential intercellular binding activity. Mol. Biol. Evol. 21:851859.
Hecht, N. B. 1993. Gene expression during male germ cell development. Pp. 400432 in C. Desjardins and L. L. Ewing, eds. Cell and molecular biology of the testis. Oxford University Press, New York.
Khil, P. P., N. A. Smirnova, P. J. Romanienko, and R. D. Camerini-Otero. 2004. The mouse X chromosome is enriched for sex-biased genes not subject to selection by meiotic sex chromosome inactivation. Nat. Genet. 36:642646.[CrossRef][ISI][Medline]
Li, W.-H. 1993. Unbiased estimation of the rates of synonymous and nonsynonymous substitution. J. Mol. Evol. 36:9699.[ISI][Medline]
Lynn, D. J., A. T. Lloyd, M. A. Fares, and C. O' Farrelly. 2004. Evidence of positively selected sites in mammalian -defensins. Mol. Biol. Evol. 21:819827.
Maxwell, A. I., G. M. Morrison, and J. R. Dorin. 2003. Rapid sequence divergence in mammalian ß-defensins by adaptive evolution. Mol. Imunnol. 40:413421.[CrossRef]
McCarrey, J. R. 1993. Development of the germ cell. Pp. 5889 in C. Desjardins and L. L. Ewing, eds. Cell and molecular biology of the testis. Oxford University Press, New York.
Meistrich, M. L. 1993. Nuclear morphogenesis during spermiogenesis. Pp. 6797 in D. de Kretser, ed. Molecular biology of the male reproductive system. Academic Press, San Diego, Calif.
Raff, R. A. 1996. The shape of life: genes, development, and the evolution of animal form. University of Chicago Press, Chicago, Ill.
Rice, W. R. 1984. Sex chromosomes and the evolution of sexual dimorphism. Evolution 38:735742.[ISI]
Schultz, N., F. K. Hamra, and D. L. Garbers. 2003. A multitude of genes expressed solely in meiotic or postmeiotic spermatogenic cells offers a myriad of contraceptive targets. Proc. Natl. Acad. Sci. USA 100:1220112206.
Su, A.I., T. Wiltshire, S. Batalov et al. (13 co-authors). 2004. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc. Natl. Acad. Sci. USA 101:60626067.
Swanson, W. J., R. Nielsen, and Q. F. Yang. 2003. Pervasive adaptive evolution in mammalian fertilization proteins. Mol. Biol. Evol. 20:1820.
Swanson, W. J., and V. D. Vacquier. 2002. The rapid evolution of reproductive proteins. Nat. Rev. Genet. 3:137144.[CrossRef][ISI][Medline]
Thatcher, T. H., and M. A. Gorovsky. 1994. Phylogenetic analysis of the core histones H2a, H2b, H3, and H4. Nucleic Acids Res. 22:174179.[Abstract]
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:46734680.[Abstract]
Torgerson, D. G., R. J. Kulathinal, and R. S. Singh. 2002. Mammalian sperm proteins are rapidly evolving: evidence of positive selection in functionally diverse genes. Mol. Biol. Evol. 19:19731980.
Torgerson, D. G., and R. S. Singh. 2003. Sex-linked mammalian sperm proteins evolve faster than autosomal ones. Mol. Biol. Evol. 20:17051709.
Tourtellotte, W. G., R. Nagarajan, A. Auyeng, C. Mueller, and J. Milbrandt. 1999. Infertility associated with incomplete spermatogenic arrest and oligozoospermia in Egr4-deficient mice. Development 126:50615071.
Vallender, E. J., and B. T. Lahn. 2004. How mammalian sex chromosomes acquired their peculiar gene content. BioEssays 26:159169.[CrossRef][ISI][Medline]
Wang, P. J., J. R. McCarrey, F. Yang, and D. C. Page. 2001. An abundance of X-linked genes expressed in spermatogonia. Nat. Genet. 27:422426.[CrossRef][ISI][Medline]
Waterston, R. H., K. Lindblad-Toh, E. Birney et al. (222 co-authors). 2002. Initial sequencing and comparative analysis of the mouse genome. Nature 420:520562.[CrossRef][ISI][Medline]
Winter, E. E., L. Goodstadt, and C. P. Ponting. 2004. Elevated rates of protein secretion, evolution, and disease among tissue-specific genes. Genome Res. 14:5461.
Wu, C. I., and E. Y. Xu. 2003. Sexual antagonism and X inactivation: the SAXI hypothesis. Trends Genet. 19:243247.[CrossRef][ISI][Medline]
Wyckoff, G. J., W. Wang, and C. I. Wu. 2000. Rapid evolution of male reproductive genes in the descent of man. Nature 403:304309.[CrossRef][ISI][Medline]
Xu, X., P. A. Toselli, L. D. Russell, and D. C. Seldin. 1999. Globozoospermia in mice lacking the casein kinase II alpha' catalytic subunit. Nat. Genet. 23:118121.[CrossRef][ISI][Medline]
Yang, Z. H. 1997. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13:555556.[Medline]
Yang, Z. H. 2002. Inference of selection from multiple species alignments. Curr. Opin. Genet. Dev. 12:688694.[CrossRef][ISI][Medline]
Yang, Z. H., R. Nielsen, N. Goldman, and A. M. K. Pedersen. 2000. Codon substitution models for heterogeneous selection pressure at amino acid sites. Genetics 155:431449.
Zechner, U., M. Wilda, H. Kehrer-Sawatzki, W. Vogel, R. Fundele, and H. Hameister. 2001. A high density of X-linked genes for general cognitive ability: a run-away process shaping human evolution? Trends Genet. 17:697701.[CrossRef][ISI][Medline]