MBE Advance Access originally published online on January 29, 2008
Molecular Biology and Evolution 2008 25(4):617-619; doi:10.1093/molbev/msn020
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Letters |
Evidence for Male-Driven Evolution in Drosophila
Division of Biological Sciences, University of California, San Diego
E-mail: dbachtrog{at}ucsd.edu.
 |
Abstract |
|---|
 TOP Abstract AcknowledgementsReferences |
|---|
In several vertebrate taxa studied to date, mutation rates are higher in males than females (male-driven evolution). The male-to-female mutation rate (
) can be estimated by contrasting DNA divergence data at X-linked, Y-linked, and autosomal loci. Previous studies in Drosophila, comparing X-linked and autosomal divergence, have found no evidence for male-driven evolution in this genus. Here, I compare levels of nucleotide divergence between homologous X- and Y-linked loci in Drosophila miranda. Using divergence at both synonymous sites and at short introns, I estimate to be approximately 2. This study thus provides the first evidence for male-biased mutation rates outside vertebrates, supporting the view that DNA sequence evolution is male driven in a wide variety of taxa.
Key Words: molecular evolution • male-driven evolution • Drosophila
Haldane (1935)
proposed that in humans the mutation rate in males would be higher than in females because the male germ line goes through many more cell divisions per generation than the female. Thus, if mutations arise by replication error, most mutations would occur in males and DNA sequence evolution would be male driven. The male-to-female mutation rate () can be estimated by contrasting DNA divergence data at X-linked, Y-linked, and autosomal loci as these chromosomes reside differently in males and females. In particular, the Y chromosome is restricted to males, autosomes spend half their times in males and females, whereas X chromosomes are transmitted only 1/3 of the time through males. Thus, under male-driven evolution, Y chromosomes should have the highest mutation rate and X chromosomes the lowest. Assume the mutation rates for a X-linked, Y-linked, and autosomal sequence are X, Y, and A, respectively. Miyata et al. (1987) showed that Y/X = 3/(2 + ), Y/A = 2/(1 + ), and X/A = (2/3)(2 + )/(1 + ). From these formulas, one can estimate if the Y/X, Y/A, or X/A ratio is known. To estimate these ratios, divergence at a pair of homologous nonfunctional sequences on the different chromosomes can be used because in a nonfunctional sequence the rate of nucleotide substitution is equal to the rate of mutation (Kimura 1983).
In several vertebrate taxa, including fish, birds, and mammals, researchers found evidence for male-driven evolution based on comparisons of homologous X and Y or autosomal sequences (Ellegren and Fridolfsson 1997, 2003; Makova and Li 2002; Sandstedt and Tucker 2005). In different organisms, estimates of were found to vary roughly between 2 and 10 (see Li et al. [2002] for a recent review). To date, male-driven evolution has not been found outside vertebrates. Previous studies in Drosophila have contrasted X and autosomal divergence and found no evidence for male-driven evolution (Bauer and Aquadro 1997; Betancourt et al. 2002). However, there are potential problems when comparing levels of divergence at nonhomologous loci. In particular, neutral mutation rates may differ between loci due to differences in base composition or chromosomal effects. Also, the number of loci compared in previous studies was relatively small, and the power to detect male-driven evolution is maximized in contrasts involving X- versus Y-linked data (Miyata et al. 1987). However, the lack of homologous sequences between the X and Y chromosome has prevented such comparisons in the Drosophila melanogaster species group.
Here, I test for male-driven evolution in Drosophila taking advantage of the recently formed neo-sex chromosomes of Drosophila miranda. In this species, a pair of autosomes became sex linked only about 1 MYA, creating so-called neo-sex chromosomes (Bachtrog and Charlesworth 2002). The neo-sex chromosomes of D. miranda show high levels of sequence similarity over wide stretches (Bachtrog 2005), which allow me to compare rates of molecular evolution at homologous neo-X- and neo-Y-linked loci. To test for male-driven evolution, a total of 205 homologous gene pairs located on the neo-sex chromosomes were analyzed (Bachtrog et al. 2008). Changes along the neo-X and neo-Y chromosome of D. miranda are polarized using Drosophila pseudoobscura as an outgroup species, assuming parsimony. Substitution rates along the neo-X and neo-Y chromosome were calculated using Jukes–Cantor. Given the low amount of sequence divergence between the neo-sex chromosomes (
3%, see table 1), parsimony reconstruction to infer the ancestral state is accurate (Kimura 1983).
|
To estimate neutral mutation rates of the neo-sex chromosomes, I consider 2 potentially unconstrained classes of neutral DNA: synonymous sites and short introns. In Drosophila, synonymous codons are not used randomly but instead some codons are used preferentially over a different codon encoding for the same amino acid; for example, preferred (P) and unpreferred (U) codons (Bulmer 1991
; Akashi and Schaeffer 1997). Thus, synonymous sites are not necessarily free of selective constraints. Nevertheless, selection for codon bias is strongly reduced in D. miranda, and synonymous sites are evolving probably close to neutral in this species (Bartolomé and Charlesworth 2006; Bachtrog 2007). Synonymous changes accumulate on the neo-Y chromosome at a higher rate than on the neo-X (dS neo-Y = 1.88% vs. dS neo-X = 1.33%, table 1), corresponding to a ratio of Y/X divergence of roughly 1.4. A similar excess of neo-Y-linked mutations is estimated using the absolute number of mutations accumulated on the neo-X or neo-Y branch (table 2).
|
The efficacy of natural selection is reduced on the neo-Y chromosome, due to its smaller effective population size and interference selection, which could result in a faster accumulation of unpreferred substitutions on the neo-Y (Kliman and Hey 1993
; Bachtrog and Charlesworth 2002). To distinguish whether synonymous substitutions accumulate on the neo-Y either due to a higher mutation rate (i.e., male-driven evolution) or due to reduced selective constraints against unpreferred synonymous codons, I classified synonymous changes into deleterious (P 
U), beneficial (U P), and presumably neutral changes (i.e., U U and P P), using preferred codons identified in D. pseudoobscura (Bachtrog 2007). If increased mutation rates on the neo-Y account for the elevated rate in synonymous substitutions, both the neutral and the deleterious class should show a faster rate of accumulation on the neo-Y. If elevated divergence on the neo-Y solely results from a reduced efficacy of natural selection against unpreferred synonymous changes, only the deleterious class should show an increased number of changes on the neo-Y. I find that both deleterious and neutral changes accumulate similarly and significantly faster on the neo-Y (table 2). This suggests an elevated mutation rate on the neo-Y chromosome as the cause for the increased rate of fixations of synonymous substitutions. Only preferred changes show no significant acceleration in their rate of accumulation on the neo-Y (table 2). This is consistent with a larger mutation rate of the neo-Y and weak positive selection for preferred codons operating on the neo-X (Bartolomé and Charlesworth 2006; Bachtrog 2007).
Another class of putatively unconstrained sequence in Drosophila is short introns. In comparisons between D. melanogaster and Drosophila simulans or between D. miranda and D. pseudoobscura, short introns are found to evolve significantly faster than long introns at rates comparable with synonymous sites (Haddrill et al. 2005; Bachtrog and Andolfatto 2006; Halligan and Keightley 2006). Thus, I also classified mutations at short introns (<100 bp) to accumulate either on the neo-X or on the neo-Y branch. Like synonymous sites, mutations at short introns accumulate faster on the neo-Y compared with the neo-X; the ratio of Y/X divergence is roughly 1.2–1.4 (tables 1 and 2), similar to the ratio inferred using synonymous site. Thus, both synonymous sites and short introns suggest a higher mutation rate for the neo-Y chromosome compared with the neo-X.
Elevated mutation rates on the neo-Y chromosome indicate a higher mutation rate in males relative to females. One can use the amount of nucleotide divergence on the neo-X and neo-Y chromosome to calculate , the male-to-female ratio of mutation rate. Confidence intervals (CIs) were estimated following Makova and Li (2002). I infer that = 1.8 (95% CI 1.5–2.1) using synonymous changes and = 1.7 (95% CI 1.2–2.3) using short introns (see tables 1 and 2), that is, the mutation rate in male Drosophila might almost be twice as high as the female mutation rate. Previous studies in Drosophila, which compared sequence divergence at X-linked and autosomal loci found no evidence for male-driven evolution (Bauer and Aquadro 1997; Betancourt et al. 2002). This probably is due to a combination of the relatively small number of loci considered, lower power in X–autosome than X–Y comparisons (Miyata et al. 1987), and the use of nonhomologous loci. A recent study comparing X and autosomal divergence of the genome sequences of D. melanogaster, D. simulans, and Drosophila yakuba found that X-linked sequences—including intron, intergenic, and synonymous sites—generally evolve faster than autosomal sequences (Begun et al. 2007). However, no distinction was made between potentially functional and nonfunctional sites (i.e., neutral synonymous changes or short introns), preventing inferences about male-driven evolution in this species group. It is also possible that either the number of germ line cell divisions in males and females differ between the D. melanogaster and D. pseudoobscura species group or the age of reproduction, such that only the latter shows male-biased evolution (Bauer and Aquadro 1997). In particular, in D. melanogaster the ratio of germ line cell divisions appears to change over time from indicating a weak female bias to a male bias as the age of reproduction increases (Bauer and Aquadro 1997). Thus, differences in the average age of reproduction among Drosophila species could account for varying levels of male-driven evolution in this genus, but little is known about differences in life history strategies among Drosophila species in nature.
It should be noted that = 1.8 is likely an underestimate of the extent of the male mutation bias, due to different amounts of segregating polymorphism on the neo-sex chromosomes (Bartolome et al. 2005). Sampling of a single allele from the neo-X and neo-Y chromosome, as done in this study, does not allow fixed and segregating mutations to be distinguished. That is, a small fraction of the mutations assigned as a fixed difference occurring on the neo-X chromosome might in fact be segregating in the population. The severe reduction in levels of variability on the neo-Y chromosome implies that nearly all mutations on the neo-Y are fixed within the population (Bachtrog 2004; Bartolomé and Charlesworth 2006). However, a few of the variants that are assigned as a fixed difference occurring on the neo-X might be polymorphic in the population, causing us to underestimate the amount of male-driven evolution. Also note that the neo-Y chromosome of D. miranda lacks recombination, whereas the neo-X recombines in females. Thus, if recombination were mutagenic this would also result in underestimating . There is, however, no evidence for mutagenic effects of recombination in Drosophila (Begun and Aquadro 1992).
Mutation rates have been found to be higher in males in several vertebrate taxa studied to date. My study is the first to report male-biased mutation rates outside vertebrates. Taking advantage of the recently formed neo-sex chromosomes of D. miranda, I estimate that there are approximately twice as many mutations in male Drosophila compared with females. This provides evidence that DNA sequence evolution is male driven in a wide variety of taxa, even outside vertebrates.
|
Acknowledgements |
|---|
|
TOPAbstractAcknowledgementsReferences |
|---|
This work was funded by a National Institute of Health Grant (GM076007) to D.B.
|
Footnotes |
|---|
Jody Hey, Associate Editor
|
References |
|---|
|
TOPAbstractAcknowledgementsReferences |
|---|
Akashi H, Schaeffer SW. Natural selection and the frequency distributions of "silent" DNA polymorphism in Drosophila. Genetics. (1997) 146:295–307.[Abstract]
Bachtrog D. Evidence that positive selection drives Y-chromosome degeneration in Drosophila miranda. Nat Genet. (2004) 36:518–522.[CrossRef][Web of Science][Medline]
Bachtrog D. Sex chromosome evolution: molecular aspects of Y chromosome degeneration in Drosophila. Genome Res. (2005) 15:1393–1401.
Bachtrog D. Reduced selection for codon usage bias in Drosophila miranda. J Mol Evol. (2007) 64:586–590.[CrossRef][Web of Science][Medline]
Bachtrog D, Andolfatto P. Selection, recombination and demographic history in Drosophila miranda. Genetics. (2006) 174:2045–2059.
Bachtrog D, Charlesworth B. Reduced adaptation of a non-recombining neo-Y chromosome. Nature. (2002) 416:323–326.[CrossRef][Medline]
Bachtrog D, Hom E, Wong K, Maside X, de Jong P. Genomic degradation of a young Y chromosome in Drosophila miranda. Genome Biol. (2008) 9:R30.[Medline]
Bartolomé C, Charlesworth B. Evolution of amino-acid sequences and codon usage on the Drosophila miranda neo-sex chromosomes. Genetics. (2006) 174:2033–2044.
Bartolome C, Maside X, Yi S, Grant A, Charlesworth B. Patterns of selection on synonymous and nonsynonymous variants in Drosophila miranda. Genetics. (2005) 169:1495–1507.
Bauer VL, Aquadro CF. Rates of DNA sequence evolution are not sex-biased in Drosophila melanogaster and D. simulans. Mol Biol Evol. (1997) 14:1252–1257.[Abstract]
Begun D, Holloway A, Stevens K, et al, 13 co-authors. Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans. PLoS Biol. (2007) 5:e310.[CrossRef][Medline]
Begun DJ, Aquadro CF. Levels of naturally occurring DNA polymorphism correlate with recombination rates in Drosophila melanogaster. Nature. (1992) 356:519–520.[CrossRef][Medline]
Betancourt A, Presgraves D, Swanson W. A test for faster X evolution in Drosophila. Mol Biol Evol. (2002) 19:1816–1819.
Bulmer M. The selection-mutation-drift theory of synonymous codon usage. Genetics. (1991) 129:897–907.[Abstract]
Ellegren H, Fridolfsson AK. Male-driven evolution of DNA sequences in birds. Nat Genet. (1997) 17:182–184.[CrossRef][Web of Science][Medline]
Ellegren H, Fridolfsson AK. Sex-specific mutation rates in salmonid fish. J Mol Evol. (2003) 56:458–463.[CrossRef][Web of Science][Medline]
Haddrill P, Charlesworth B, Halligan D, Andolfatto P. Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content. Genome Biol. (2005) 6:R67.[CrossRef][Medline]
Haldane JBS. The rate of spontaneous mutation of a human gene. J Genet. (1935) 31:317–326.[CrossRef][Web of Science]
Halligan D, Keightley P. Ubiquitous selective constraints in the Drosophila genome revealed by a genome-wide interspecies comparison. Genome Res. (2006) 16:875–884.
Kimura M. The neutral theory of molecular evolution. (1983) Cambridge: Cambridge University Press.
Kliman RM, Hey J. Reduced natural selection associated with low recombination in Drosophila melanogaster. Mol Biol Evol. (1993) 10:1239–1258.[Abstract]
Li W, Yi S, Makova K. Male-driven evolution. Curr Opin Genet Dev. (2002) 12:650–656.[CrossRef][Web of Science][Medline]
Makova K, Li W. Strong male-driven evolution of DNA sequences in humans and apes. Nature. (2002) 416:624–626.[CrossRef][Medline]
Miyata T, Hayashida H, Kuma K, Mitsuyasu K, Yasunaga T. Male-driven molecular evolution: a model and nucleotide sequence analysis. Cold Spring Harb Symp Quant Biol. (1987) 52:863–867.
Sandstedt S, Tucker P. Male-driven evolution in closely related species of the mouse genus Mus. J Mol Evol. (2005) 61:138–144.[CrossRef][Web of Science][Medline]
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  C. J. Pink, S. K. Swaminathan, I. Dunham, J. Rogers, A. Ward, and L. D. Hurst Evidence That Replication-Associated Mutation Alone Does Not Explain Between-Chromosome Differences In Substitution Rates Gen Biol Evol, June 22, 2009; 2009(0): 13 - 22. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
V. L. Bauer DuMont, N. D. Singh, M. H. Wright, and C. F. Aquadro Locus-Specific Decoupling of Base Composition Evolution at Synonymous Sites and Introns along the Drosophila melanogaster and Drosophila sechellia Lineages Gen Biol Evol, June 22, 2009; 2009(0): 67 - 74. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Engelstadter Muller's Ratchet and the Degeneration of Y Chromosomes: A Simulation Study Genetics, October 1, 2008; 180(2): 957 - 967. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||