Section of Evolution and Ecology, University of CaliforniaDavis
Several studies of Drosophila melanogaster have reported patterns of DNA variation in genes located near centromeres or near the telomere of the X chromosome (Aguadé, Miyashita, and Langley 1989
; Langley et al. 1993, 2000
; Wayne and Kreitman 1996
). There is good evidence of reduced crossing over for these regions in D. melanogaster. However, there are relatively few population data for such regions in Drosophila simulans (Begun and Aquadro 1991
; Martin-Campos et al. 1992
; Hilton, Kliman, and Hey 1994
; Wayne and Kreitman 1996
). Here I report data from population samples of four genes located near the centromere of chromosome 3 in D. simulans: Hem-protein, CKII-
, Gelsolin, and Amalgam. The physical locations of these genes in D. melanogaster are 79E2, 80A, 82A13, and 84A5, respectively. Because there are no detectable chromosomal rearrangements between species in these regions (Ashburner 1989
), I assume that the locations of these genes are the same in D. simulans. The sequences used were from a set of inbred lines derived from flies collected in the Wolfskill Orchard in Winters, Calif. (Begun and Whitley 2000
). I also report sequence data from Drosophila yakuba for Hem-protein, CKII-
, and Gelsolin. For some analyses, I used previously published data from sequences of 40 D. simulans genes distributed across the X chromosome and chromosome arm 3R (Begun and Whitley 2000
). Sequences were analyzed with the DnaSP program (Rozas and Rozas 1999
). Silent mutations were classified as preferred or unpreferred (Sharp and Lloyd 1993
) by using the outgroup method as described by Akashi (1996)
. New sequences reported here can be found in GenBank under accession numbers AY052156AY052190.
The four genes located near the centromere of chromosome 3 of D. simulans (table 1
) show reduced heterozygosity () at silent sites compared with other genes on chromosome 3. The mean silent
for these four genes (0.009) is significantly lower (Mann-Whitney, P = 0.015) than the mean for 19 genes located more distally (0.035). Despite the clear difference in the levels of silent polymorphism for centromeric versus more distal genes, the mean silent divergence of four centromeric genes (0.097) is not significantly different from the mean of other genes (0.108) on chromosome 3. Mean replacement
values for noncentromeric (
= 0.0013) and centromeric genes (
= 0.0009) are not significantly different. Although the ratio of replacement
to replacement divergence is smaller for the centromeric genes (0.15) than for the noncentromeric genes (0.44), the difference between regions is not significant.
|
The different statistical conclusions on the effect of chromosomal location for silent versus replacement variation in D. simulans might be a result of reduced statistical power associated with low levels of amino acid variation (compared with those of silent variation) or might reflect a genuine difference between the dynamics of silent and replacement variation. More extensive sampling of replacement variation across the genome would be required to resolve this issue.
The two sampled genes located closest to the centromere are CKII-, about 26 bands to the left of the centromere on 3L, and Gelsolin, about 8 bands to the right of the centromere on 3R. Both have severely reduced levels of silent heterozygosity (table 1 ) relative to the average for the chromosome (0.035; Begun and Whitley 2000
). Silent
for Hem-protein is relatively high (
= 0.015), although this gene is located only
42 bands distal to the centromere on 3L. The same is true for Amalgam (about 97 bands to the right of the centromere on 3R), which is about as polymorphic as Hem-protein. These data suggest that severe reductions in silent heterozygosity are restricted to only a very small euchromatic region (perhaps 1,000 kb) on each side of the centromere of chromosome 3 in D. simulans. Few data exist on patterns of recombination near D. melanogaster centromeres, and even less information is available for D. simulans. However, the small amount of D. simulans genetic data (Sturtevant 1929
; True, Mercer, and Laurie 1996
) indicates a minimal centromere effect relative to D. melanogaster. Given the established relationship between recombination and polymorphism, the genetic and population data from D. simulans are consistent with a dramatic reduction in recombination over a small physical region near the centromere. Such regions may be ideal for detailed study of the effects of variable recombination on sequence evolution.
Regions of normal recombination in D. simulans exhibit a ratio of unpreferred to preferred fixations of 2 (Takano 1998;
Begun 2001
). For three genes located near the centromere of chromosome 3 (table 3)
, the ratio of unpreferred to preferred fixations (1:8) is significantly different (G-test, P = 0.01) from the ratio observed for more distally located genes on chromosome 3 (37:18). There should be no effect of recombination rate on the proportion of preferred versus unpreferred fixations for genes evolving at mutation-selection-drift equilibriumwe expect roughly equal numbers of fixations in the two presumptive fitness classes regardless of the recombinational environment (Akashi 1996
). A possible explanation for the observed heterogeneity is that the recombination rate in the centromere region of D. simulans has recently increased. Genes with a long history of reduced recombination are expected to have a higher proportion of unpreferred to preferred codons at equilibrium than are genes located in regions of more extensive crossing over. This is because low rates of crossing over decrease the efficacy of purifying selection against very slightly deleterious (e.g., unpreferred) mutations. If the recombination rate in a region of historically low recombination recently increased, then we would expect to observe a transient increase in the fixation rate for very slightly beneficial mutations as such regions "recover" from a history of ineffectual purifying selection. Further investigation of silent fixations in centromeric regions of D. simulans and further genetic analysis of species in the melanogaster subgroup will be required to evaluate this hypothesis.
|
|
Acknowledgements
I thank P. Whitley for lab work and C. Langley, W. Stephan, and anonymous reviewers for comments. This work was supported by the NIH and the Sloan Foundation.
Footnotes
Wolfgang Stephan, Reviewing Editor
Keywords: Drosophila simulans
protein variation
DNA polymorphism
natural selection
population
Address for correspondence and reprints: David J. Begun, Section of Evolution and Ecology, University of California, Davis, California 95616. djbegun{at}ucdavis.edu
.
References
Akashi H., 1996 Molecular evolution between Drosophila melanogaster and D. simulans: reduced codon bias, faster rates of amino acid substitution, and larger proteins in D. melanogaster Genetics 144:1297-1307
Aguad M., N. Miyashita, C. H. Langley, 1989 Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster Genetics 122:607-615
Ashburner M., 1989 Drosophila: a laboratory handbook Cold Spring Harbor Press, New York
Begun D. J., 2001 The frequency distribution of nucleotide variation in Drosophila simulans Mol. Biol. Evol 18:1343-1352
Begun D. J., C. F. Aquadro, 1991 Molecular population genetics of the distal region of the X chromosome in Drosophila: evidence for genetic hitchhiking of the yellow-achaete region Genetics 129:1147-1158
Begun D. J., P. Whitley, 2000 Reduced X-linked nucleotide variation in Drosophila simulans Proc. Natl. Acad. Sci. USA 97:5960-5965
Hilton H., R. M. Kliman, J. Hey, 1994 Using hitchhiking genes to study adaptation and divergence during speciation within the Drosophila melanogaster species complex Evolution 48:1900-1913[ISI]
Langley C. H., B. P. Lazzaro, W. Phillips, E. Heikkinen, J. M. Braverman, 2000 Linkage disequilibrium and the site frequency spectra in the su(s) and su(wa) regions of the Drosophila melanogaster X chromosome Genetics 156:1837-1852
Langley C. H., J. MacDonald, N. Miyashita, M. Aguad, 1993 Lack of correlation between interspecific divergence and intraspecific polymorphism at the suppressor of forked region in Drosophila melanogaster and Drosophila simulans Proc. Natl. Acad. Sci. USA 90:1800-1803[Abstract]
Martin-Campos J. M., J. M. Comeron, N. Miyashita, M. Aguad, 1992 Intraspecific and interspecific variation at the y-ac-sc region of Drosophila simulans and Drosophila melanogaster Genetics 130:805-816
Rozas J., R. Rozas, 1999 DnaSP 3: an integrated program for molecular population genetics and molecular evolution analysis Bioinformatics 15:174-175
Sharp P. M., A. T. Lloyd, 1993 Codon usage Pp. 378397 in G. Maroni, ed. An atlas of Drosophila genes: sequences and molecular features. Oxford University Press, Oxford, England
Sturtevant A. H., 1929 Contributions to the genetics of Drosophila simulans Carnegie Inst. Wash. Publ 399:1-62
Takano T. S., 1998 Rate variation of DNA sequence evolution in the Drosophila lineages Genetics 149:959-970
True J. R., J. M. Mercer, C. C. Laurie, 1996 Differences in crossover frequency and distribution among three sibling species of Drosophila Genetics 142:507-523
Watterson G. A., 1975 On the number of segregating sites in genetical models without recombination Theor. Popul. Biol 7:256-276[ISI][Medline]
Wayne M. L., M. Kreitman, 1996 Reduced variation at concertina, a heterochromatic locus in Drosophila Genet. Res 69:101-108[ISI]