Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

Similar Levels of X-linked and Autosomal Nucleotide Variation in African and non-African populations of Drosophila melanogaster

Nadia D Singh1*, J Michael Macpherson1, Jeffrey D Jensen2 and Dmitri A Petrov1

Author Affiliations

1 Department of Biological Sciences, Stanford University, Stanford, CA 94305 USA

2 Department of Molecular Biology and Genetics, Cornell University, Ithaca NY 14853 USA

For all author emails, please log on.

BMC Evolutionary Biology 2007, 7:202  doi:10.1186/1471-2148-7-202

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/7/202


Received:26 April 2007
Accepted:25 October 2007
Published:25 October 2007

© 2007 Singh et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Levels of molecular diversity in Drosophila have repeatedly been shown to be higher in ancestral, African populations than in derived, non-African populations. This pattern holds for both coding and noncoding regions for a variety of molecular markers including single nucleotide polymorphisms and microsatellites. Comparisons of X-linked and autosomal diversity have yielded results largely dependent on population of origin.

Results

In an attempt to further elucidate patterns of sequence diversity in Drosophila melanogaster, we studied nucleotide variation at putatively nonfunctional X-linked and autosomal loci in sub-Saharan African and North American strains of D. melanogaster. We combine our experimental results with data from previous studies of molecular polymorphism in this species. We confirm that levels of diversity are consistently higher in African versus North American strains. The relative reduction of diversity for X-linked and autosomal loci in the derived, North American strains depends heavily on the studied loci. While the compiled dataset, comprised primarily of regions within or in close proximity to genes, shows a much more severe reduction of diversity on the X chromosome compared to autosomes in derived strains, the dataset consisting of intergenic loci located far from genes shows very similar reductions of diversities for X-linked and autosomal loci in derived strains. In addition, levels of diversity at X-linked and autosomal loci in the presumably ancestral African population are more similar than expected under an assumption of neutrality and equal numbers of breeding males and females.

Conclusion

We show that simple demographic scenarios under assumptions of neutral theory cannot explain all of the observed patterns of molecular diversity. We suggest that the simplest model is a population bottleneck that retains an ancestral female-biased sex ratio, coupled with higher rates of positive selection at X-linked loci in close proximity to genes specifically in derived, non-African populations.

Background

Drosophila melanogaster has been extensively used to investigate the role of various evolutionary forces in the generation and maintenance of sequence variation in natural populations. D. melanogaster is a cosmopolitan species believed to be African in origin [1] and has only recently colonized the rest of the world. This demographic history makes D. melanogaster an attractive system for assessing the impact of factors such as demography and natural selection on levels of extant sequence variation.

Studies of molecular polymorphism in D. melanogaster abound (for review see [2]). However, despite plentiful polymorphism data, the history of selection and demography of D. melanogaster remains relatively obscure. Some of the difficulty stems from the fact that the current data are biased in favor of genic (synonymous sites, UTRs, and intronic) and mostly X-linked loci (for review see [2]). While there have been some studies of intergenic sequence polymorphism, most of the studied regions are in close proximity to genes [3-5]. In addition, features such as local mutation rate and local recombination rate, both of which are likely to influence patterns of polymorphism, have not been controlled for in previous studies. As a result, patterns of polymorphism at these loci are likely to reflect not only genome-wide processes such as demographic history but also local, locus-specific effects of selection, mutation rate variability and Hill-Robertson [6,7] effects.

All studies to date do support a general reduction in levels of variation in derived versus ancestral strains of D. melanogaster. This was initially reported in studies relying on electrophoretic variants of proteins as well as variations in restriction maps (e.g. [8-11]); more recently single nucleotide polymorphism and microsatellite data have been used to show similar patterns (for review see [2]). Within these studies of nucleotide or microsatellite variability, however, the magnitude of this reduction appears to be heavily dependent on the choice of loci. Initial reports investigating X-linked loci consistently showed reductions in sequence variation in non-African, derived strains [9,12,13], while autosomal loci appeared to have more variable reductions in sequence polymorphism in non-African strains relative to their African progenitors [14-21]. More recent, large-scale studies show similar trends; sequence variation is markedly reduced on the X chromosome in derived strains while reductions in diversity at autosomal regions are less uniform across loci [22,23]. These results were generally interpreted as evidence that purely demographic processes could not fully account for patterns of nucleotide variability in non-African populations of D. melanogaster. Multilocus studies evaluating the frequency distributions of nucleotide variants also seem to implicate selective forces in the generation of patterns of nucleotide variability in these derived populations [3,4,24], as do studies using microsatellites [25-29].

Because patterns of variation show locus-specific effects, direct comparisons of diversity levels between the X chromosome and the autosomes have also shown variable results. Some studies suggest that X-linked diversity is indeed lower than autosomal diversity in D. melanogaster and its cosmopolitan sister species D. simulans [27,29]. However, the suggestion has also been made that while this is the case for non-African strains of these species, X-linked diversity is equal to or even greater than autosomal diversity in African strains [22,28]. There are also some data to support increased levels of variation of autosomal loci in African and some non-African strains but not in others [23].

It is thus clear that sampling plays an important role in the studies of molecular polymorphism in Drosophila. To shed more light on X-linked and autosomal diversity in ancestral and derived strains of D. melanogaster, we conducted a careful survey of single nucleotide polymorphism in D. melanogaster. We chose loci so as to minimize the potential effects of locus-specific forces that affect nucleotide variation such as selective constraint on protein function, Hill-Robertson effects, and variation in mutation rates. The chosen loci are intergenic, located in regions of high recombination, have a similar GC content, are as far as possible from genes, and have a similar level of divergence with D. simulans. By using stringent and uniform criteria for candidate locus selection, we hoped to minimize locus-specific effects such that we might be able to compare autosomal and X-linked diversity between ancestral and derived strains in a consistent manner.

Our findings indicate that levels of nucleotide diversity in North American strains of D. melanogaster are reduced compared to those in African strains, with similar reductions in diversity at the X-linked and autosomal loci. Moreover, levels of X-linked and autosomal diversity are more similar than expected under an assumption of strict neutrality and an equal sex ratio in both ancestral and derived populations of D. melanogaster. If we assume that levels of sequence diversity in extant African populations are indicative of levels of diversity in ancestral populations, these data are consistent with an historical excess of breeding females over the number of breeding males and a population bottleneck where this excess is retained. These findings are particularly relevant to the observation of consistently higher codon bias of X-linked genes in D. melanogaster [30-32] and may imply that natural selection might be more efficient at X-linked loci in general.

We also combined our data with much available published single nucleotide polymorphism data and conducted a meta-analysis. Results from this analysis suggest that while levels of diversity in African strains are greater than diversity levels in non-African strains, this effect is much more pronounced on the X chromosome than it is on the autosomes. In addition, the compiled data from ancestral populations suggest that levels of nucleotide polymorphism are similar between X-linked and autosomal loci. If we again use extant sequence diversity in African populations as indicators of ancestral polymorphism, these data cannot be fully explained by a population bottleneck that retains the presumed unequal ancestral sex ratio.

We discuss the differences between the two datasets and consider several hypotheses aimed at reconciling the observed patterns. While the data from the intergenic loci presented here are consistent with a female-biased sex ratio in ancestral and derived populations coupled with a population bottleneck, the results from the meta-analysis are qualitatively consistent with either a model wherein the ancestral population has a female-biased sex ratio that is not retained during and after the bottleneck, or where the ratio of the effective sizes of the X and autosomes is less than one in the ancestral population and this ratio is preserved during and after a bottleneck. Presently, it appears that a demographic model incorporating a female-biased sex ratio in the ancestral population and a population bottleneck is likely to be the most appropriate. However, our confidence in this model would benefit from further consistent and careful sampling of polymorphism in the D. melanogaster genome, particularly on the autosomes.

Results

Candidate locus selection

We identified eight autosomal and eight X-linked loci in the D. melanogaster genome based on several criteria designed to minimize the effects of natural selection on levels of sequence variation. We chose regions of the genome that are nongenic and located as far from genes as possible. The chosen loci are at least 5 kb and on average 21 kb from the nearest annotated gene as measured from the midpoint of the candidate locus to the 5' or 3' end (depending on the orientation) of the nearest coding sequence. They ranged in size from 322 to 487 bp (Table 1). All the studied loci are located in high recombination regions of the genome [33] to minimize linkage to functionally important sites. We computed divergence with D. simulans and picked regions with similar rates of substitution to control for both selective constraint and mutation rate variation as much as possible. Using these criteria, we selected loci least likely to be affected by natural selection acting either on the loci themselves or on nearby linked loci, and we controlled for variation in recombination rate and mutation rate as much as possible.

Table 1. Summary statistics of X-linked and autosomal loci

Polymorphism and divergence at intergenic loci

We surveyed these sixteen noncoding loci for levels of nucleotide diversity using twelve African strains and twelve North American strains of D. melanogaster (see Methods). Sequence data were collected from between nine and twelve strains per locus for each population. Levels of nucleotide diversity estimated based on the number segregating sites (θs) are presented in Table 1 and Figure 1. To test for selective neutrality at each locus, Tajima's D [34] was computed, and in all cases was not significantly different from zero (P > 0.05, all loci). Tajima's D for these data grouped by population or chromosome was also not significantly different from zero (P > 0.2, sign test, both cases).

thumbnailFigure 1. Frequency distribution of observed nucleotide diversity across loci in a) Singh et al. 16 X-linked and autosomal loci from African strains, b) Singh et al. 16 X-linked and autosomal loci from North American strains c) compiled data from X-linked and autosomal loci for African strains and d) compiled data from X-linked and autosomal loci for African strains.

Average levels of nucleotide diversity per site (estimated by θs) were 17.52 × 10-3 (15.73 × 10-3–19.17 × 10-3) and 9.00 × 10-3 (7.70 × 10-3–10.20 × 10-3) for African and American strains (Table 1), respectively, and this increase in the level of DNA polymorphism in African strains was statistically significant (P = 0.0005, paired sign test). When partitioned by chromosomal location, X-linked loci also show significantly elevated diversity in African strains (17.30 × 10-3 versus 8.15 × 10-3, P = 0.008, paired sign test), while the increase in African diversity at autosomal loci was only marginally statistically significant (17.79 × 10-3 versus 9.39 × 10-3, P = 0.07, paired sign test). Based on these estimates of diversity we estimate that the ratio of polymorphism between African and non-African strains is consistently greater than one. When both X-linked and autosomal loci are considered, θ A f r i c a n θ A m e r i c a n MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGbbqqcqWGMbGzcqWGYbGCcqWGPbqAcqWGJbWycqWGHbqycqWGUbGBaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabd2gaTjabdwgaLjabdkhaYjabdMgaPjabdogaJjabdggaHjabd6gaUbqabaaaaaaa@443B@ = 1.96 (95% confidence interval, 1.65–2.33; see Methods), when only X-linked loci are considered, θ A f r i c a n θ A m e r i c a n MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGbbqqcqWGMbGzcqWGYbGCcqWGPbqAcqWGJbWycqWGHbqycqWGUbGBaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabd2gaTjabdwgaLjabdkhaYjabdMgaPjabdogaJjabdggaHjabd6gaUbqabaaaaaaa@443B@ = 2.04 (1.58–2.63), and when only autosomal loci are considered, θ A f r i c a n θ A m e r i c a n MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGbbqqcqWGMbGzcqWGYbGCcqWGPbqAcqWGJbWycqWGHbqycqWGUbGBaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabd2gaTjabdwgaLjabdkhaYjabdMgaPjabdogaJjabdggaHjabd6gaUbqabaaaaaaa@443B@ = 1.90 (1.46–2.39) (Table 2).

Table 2. Estimates of Ratios of Diversity

Estimates of nucleotide diversity at orthologous pairs of noncoding loci between African and American strains were significantly correlated (Kendall's τ = 0.40; P = 0.03). This correlation is at least partly due to mutation rate heterogeneity across loci; if local mutation rate is taken into account by dividing diversity by divergence with D. simulans, levels of sequence variation are only marginally correlated between African and American strains (Kendall's τ = 0.31; P = 0.07).

Thus, while choosing candidate loci with similar divergence with D. simulans accounts for mutation rate heterogeneity to some degree, there remains some residual effect of local mutation rate on locus-specific patterns of polymorphism. To take into account the role of this mutation rate variation in nucleotide polymorphism levels, we divided our diversity estimates by divergence with D. simulans. These normalized estimates of diversity are also significantly higher in African versus North American strains overall (3.34 × 10-1 versus 1.72 × 10-1, P = 0.0005, paired sign test) as well as for X-linked loci (3.38 × 10-1 versus 1.67 × 10-1, P = 0.02, paired sign test), and marginally significantly higher for autosomal loci (3.35 × 10-1 versus 1.79 × 10-1, P = 0.07, paired sign test) (Table 1). Using these scaled levels of diversity to estimate the ratio of sequence polymorphism between ancestral and derived populations yields values significantly exceeding one in all cases; the ratios of African to American corrected diversities are 1.95 (1.60–2.36), 1.90 (1.42–2.51), and 2.06 (1.46–2.78) for all loci, autosomal loci, and X-linked loci, respectively (Table 2). The slight differences between the results of these two statistical comparisons of diversity in ancestral and derived populations may reflect the conservative nature of the paired-sign test.

Nucleotide polymorphism can also be compared between X-linked and autosomal loci. Average levels of nucleotide polymorphism uncorrected for sequence divergence were 13.28 × 10-3 (95% confidence interval 11.73 × 10-3–14.82 × 10-3) and 13.10 × 10-3 (11.59 × 10-3–14.62 × 10-3) (Table 1) for autosomal and X-linked loci when all strains were grouped together. When African and non-African populations were considered separately, X-linked loci and autosomal loci also show similar levels of diversity per site (Figure 2, Table 1). There was no significant difference in estimated polymorphism between X-linked and autosomal loci, either when both populations were considered together or separately (P = 0.94, P = 0.83, and P = 0.46 for both populations, African strains only, and American strains only, respectively, Mann-Whitney U-test). Using these estimates of nucleotide diversity to infer the ratio of X-linked to autosomal variation yielded mean ratios greater than the expected 3/4 under a model of equal numbers of breeding males and females. These ratios were estimated as 0.99 (0.83–1.17), 0.99 (0.81–1.20), and 0.92 (0.68–1.21) for both populations together, African strains, and North American strains, respectively (Table 2). These confidence intervals indicate that we can reject the expected 3/4 ratio of X to autosomal levels of variation for the African population data and for the combined populations.

thumbnailFigure 2. Mean diversities for X-linked and autosomal loci in African and non-African populations of D. melanogaster, comparing the Singh et al. dataset with the compiled dataset.

To again take into account the effects of heterogeneity in local mutation rates, we normalized our estimates of nucleotide diversity by sequence divergence with D. simulans. Average levels of normalized nucleotide diversity are 2.53 × 10-1 (2.11 × 10-1–3.01 × 10-1) and 2.55 × 10-1 (2.12 × 10-1–3.08 × 10-1) for autosomal and X-linked loci, respectively, when both African and American populations are considered, and this similarity between X-linked and autosomal loci was echoed when African and North American populations were treated separately (Table 1). These differences between the X and the autosomes were not significant in any comparison (P = 0.85, P = 0.92, and P = 0.53 for both populations, African strains, and American strains, respectively, Mann-Whitney U-test). The ratios of normalized diversity levels between X-linked and autosomal loci are 1.03 (0.79–1.32), 1.02 (0.77–1.32), and 0.95 (0.64–1.31) for both populations, African strains and American strains, respectively (Table 2). Much like the uncorrected estimates of nucleotide diversity, all samples show a higher ratio of X-linked to autosomal variation than expected, and for African populations and for the combined dataset we can reject (with 95% confidence) a ratio of 3/4.

Meta-analysis

To place our results within the context of previous results, we compiled polymorphism data from several published studies (see Methods) (Table 3). When all the data were included, levels of DNA polymorphism were significantly higher in African populations than in non-African populations (Figure 1). Mean diversity in African strains was 13.60 × 10-3 (12.43 × 10-3–14.96 × 10-3) while mean diversity in non-African strains was 5.96 × 10-3 (4.33 × 10-3–6.18 × 10-3) (Table 4). Autosomal loci show the same pattern of increased variation in African strains; average polymorphism per site was 13.92 × 10-3 (9.43 × 10-3–18.04 × 10-3) versus 13.25 × 10-3 (3.90 × 10-3–13.52 × 10-3) for African and non-African strains, respectively. Sequence variation is significantly higher in African populations for X-linked loci as well, with mean diversities of 13.55 × 10-3 (12.38 × 10-3–14.96 × 10-3) and 4.73 × 10-3 (3.97 × 10-3–5.36 × 10-3) for African and non-African strains, respectively (Table 4).

Table 3. Data used in the meta analysis published data

Table 4. Bootstrapped mean diversities and bootstrapped mean diversities normalized by divergence for the meta-analysis dataset

These differences in mean diversities between populations are recapitulated in our estimates of the ratio of nucleotide diversities in African and non-African strains (Table 4). When all loci are considered, the African/non-African ratio of diversity is significantly greater than one (mean = 2.30, 95% confidence interval = 2.11–2.95). This is also the case when only X-linked loci are included (mean = 2.88, 95% confidence interval = 2.35–3.38). When only autosomal loci are considered the ratio is greater than one but not significantly (mean = 1.09, 95% confidence interval = 0.77–1.87). The lack of difference for autosomal loci might be due to the relative paucity of autosomal loci in this dataset, as there are only 21 autosomal loci contained in this dataset (Table 3). Normalizing diversity levels by divergence further supports higher levels of diversity in African strains, the mean ratio of normalized African to non-African diversity is 2.18 (1.64–2.67) when all loci are included, 1.12 (0.60–1.90) for autosomal loci and 2.68 (1.81–3.23) for X-linked loci (Table 4).

However, factors contributing to sequence variation such as mutation rate variation and recombination were not taken into account when these loci were selected for study. To investigate the relationship between population of origin and sequence variation more thoroughly, we calculated partial correlation coefficients between nucleotide diversity and chromosomal location controlling for the combined effects of local recombination rate and sequence divergence with D. simulans. Partial correlation analysis among population (coded as a binary metric with African origin = 1 and non-African origin = 0), nucleotide diversity, recombination rate and sequence divergence revealed an overall strong, positive correlation between population and diversity (Kendall's partial τ = 0.53, P << 0.0001). This trend was echoed when this analysis was restricted to only autosomal or X-linked loci, (Kendall's partial τ = 0.31, τ = 0.61, P = 0.02, P <<0.0001 for autosomal and X-linked loci, respectively).

This compiled dataset can also be used to specifically compare X-linked and autosomal variation. Overall, levels of diversity (estimated through θs) were 13.59 × 10-3 (8.17 × 10-3–14.61 × 10-3) and 9.27 × 10-3 (8.56 × 10-3–10.06 × 10-3) for autosomal and X-linked loci when both ancestral and derived populations were considered. Levels of X-linked and autosomal sequence variation were comparable in African strains, though X-linked variation was markedly depressed relative to autosomal diversity in non-African strains (Figure 2, Table 4). These estimates correspond to ratios of effective population sizes of the X to the autosomes of 0.69 (0.59–0.94) for both populations considered together, 1.00 (0.56–1.23) for African strains and 0.37 (0.296–0.60) for non-African strains. Normalizing diversity values by divergence generates a ratio of X-linked to autosomal diversity of 0.57 (0.34–0.71) for both populations, 0.80 (0.21–0.92) for African strains, and 0.32 (0.17–0.48) for the non-African strains (Table 4).

Accounting for sequence divergence and recombination rate revealed an overall weak yet statistically significant association between autosomal linkage and increased polymorphism (Kendall's partial τ = -0.17, P = 0.02, where chromosomal location is coded as a binary metric with autosomal linkage = 0 and X-linkage = 1). This may be due entirely to the sequence data from non-African strains, where a significant association is also found (Kendall's partial τ = -0.28, P = 0.0002), as no such relationship found within the data from African strains (Kendall's partial τ = -0.02, P = 0.79).

Simulations

In order to better understand the patterns of variation revealed by these two datasets and to determine if our results are consistent with neutral demographic expectations, we explored the effects of three alternative demographic scenarios on patterns of variability at X-linked and autosomal loci. In the first scenario, we assume that the ancestral population of D. melanogaster had similar levels of X-linked and autosomal variation. Under this scenario, a population bottleneck during and after which the ratio of X-linked to autosomal variation and selective patterns remain unchanged should lead to reductions in diversity of similar magnitude on the X chromosome and autosomes. This is because in this model the effective population sizes of the X chromosome and autosomes are equal before, during and after the bottleneck, which leads to identical expected coalescence histories.

In the second scenario, we considered the possibility that patterns of selection or mating structure change during and/or after the population bottleneck. To model the effects of such a change, simulations were conducted using the Thornton and Andolfatto [35] estimated D. melanogaster bottleneck parameters (see Methods). Consistent with previous results [36], these simulations suggest that levels of sequence variation on the X chromosome and autosomes can be affected differently by the same demographic event. For instance, if the ratio of the effective number of breeding females to males decreases from 7 (leading to the expected θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ of 1 in Africa) in the pre-bottleneck population to 1 during and after the bottleneck, the expected θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ ratio in the derived population is reduced to 0.75.

Finally, we considered the possibility that the θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ ratio in the ancestral population of the European and North American populations is substantially different from one. This requires either that the ancestral population has not in fact been sampled in any of the studies of African populations, or that extant sequence variation in African strains are not reflective of ancestral levels of diversity. If θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ was less than 1 in the ancestral population, a bottleneck should reduce X-linked variation to a greater extent. Simulation results indicate that if, for instance, θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ was 0.75 in the ancestral population, a bottleneck should reduce the expected θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ ratio to 0.58 in the derived populations.

Discussion

Drosophila melanogaster is a human commensal, and has a cosmopolitan distribution [1]. This species is thought to have originated in Africa, and the demographic history of D. melanogaster is marked by an expansion out of Africa into Europe and Asia, followed by the colonization of the Americas [1]. Colonization events are generally comprised of a population bottleneck during migration followed by population expansion in a new habitat, and these phases of the colonization process may have profound consequences for within-population variability in the derived population relative to the ancestral population. Population bottlenecks are expected to depress levels of standing genetic variation [22,36]. The numbers of breeding males and females involved in the founding event can also cause the levels of diversity to differ between the X chromosome and the autosomes [36]. In addition, adaptation to a novel habitat can reduce levels of diversity at selected sites as well as at nearby linked sites [37]; this effect may manifest differentially on the X versus the autosomes if rates of adaptation differ systematically between these chromosome sets.

Because both demographic and selective models predict a reduction in population variability in derived populations, it is very challenging to disentangle the roles of these two forces (for review see [38]). Distinguishing between purely selective and purely demographic models has often relied on the inference that because population bottlenecks affect the entire genome, this demographic event should systematically reduce variation at a genomic scale. In contrast, because adaptive evolution is a locus-specific force, repeated bouts of positive selection due to local adaptation should result in extreme reductions in polymorphism only at selected loci. However, the inherent stochasticity in population bottlenecks can lead to substantial variation across loci in the reduction of nucleotide diversity due to the bottleneck, further complicating the distinction between the effects of demography and natural selection [39]. With a 1:1 sex ratio, the effective size of the X chromosome is 3/4 that of the autosomes, which results in differences in scaling time between the X and the autosomes. When this difference in scaling is properly accounted for, the population bottleneck, although less recent on the X chromosome, becomes more severe (in that it has a smaller minimal size and is effectively longer) on the X chromosome than it is on the autosomes, which leads to a greater reduction of diversity at X-linked loci relative to that found at autosomal loci.

Comparing patterns of variability between the two datasets

We can directly compare estimates of diversity between the intergenic loci presented here and the data contained in the meta-analysis to identify those populations and/or regions which appear to differ the most dramatically between the two datasets. Given mutation rate variation, the most appropriate metric is polymorphism normalized by divergence. There are no significant differences between the ratio of ancestral to derived (normalized) diversity between the two datasets for either X-linked or autosomal loci (Table 2, Table 4), and both datasets do support a reduction in variation in the non-African strains, which is consistent with a population bottleneck. In African populations, overlapping confidence intervals of the normalized estimates of diversity between the two datasets for both autosomal and X-linked loci (Table 1, Table 4) suggest that these estimates of diversity are similar between the two sets of data. Though the mean estimates of normalized diversity for autosomal regions in non-African populations do differ considerably between the two datasets (3.35 versus 1.79 for the compiled dataset and intergenic dataset, respectively), these estimates also show overlapping confidence intervals between the two datasets, suggesting that the intergenic loci presented here have levels of diversity comparable to those primarily coding regions in the compiled dataset. However, this lack of statistical significance could also be due to lack of power, given that there are only 8 loci in the intergenic dataset and 21 loci in the compiled dataset. In contrast, the two datasets do differ significantly with respect to levels of variation at X-linked loci in non-African populations, with greatly reduced normalized estimates of diversity in the compiled dataset relative to the intergenic dataset. Thus, while the majority of the estimates of diversity are consistent across both datasets, X-linked diversity in derived populations show marked differences between the intergenic loci and the compiled loci.

X-linked and autosomal variation in African populations

Both the newly collected data and the meta-analysis indicate that θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ in African populations is elevated above the expected 0.75 (under an assumption of equal numbers of breeding males and females); the data presented here suggest that the ratio is 0.99 and the estimate of the ratio is nearly identical in the meta-analysis dataset, at 1.00 (Table 2, Table 4). Given that we now have information from a variety of loci including sequences that are both coding and noncoding, located far and close to genes, and represent different types of molecular markers, we believe these results to be robust. Moreover, these surveys utilized strains collected in different parts of Africa (Table 3), which hints at the possibility that these patterns are not restricted to any one particular African population, although given population structure in Africa [40,41] these patterns may be restricted to certain geographic locales within the continent. Finally, our analyses suggest that these results are also robust to variation in recombination rate, levels of selective constraint, and mutation rates.

There are several possible explanations for the increase in X-linked sequence variation relative to expectation. One possibility is that the variance of reproductive success in D. melanogaster males is larger than that in females. Recall that the effective population size of an autosomal locus ( N e A = 4 N M N F N M + N F MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGobGtdaWgaaWcbaGaemyzau2aaSbaaWqaaiabdgeabbqabaaaleqaaOGaeyypa0ZaaSaaaeaacqaI0aancqWGobGtdaWgaaWcbaGaemyta0eabeaakiabd6eaonaaBaaaleaacqWGgbGraeqaaaGcbaGaemOta40aaSbaaSqaaiabd2eanbqabaGccqGHRaWkcqWGobGtdaWgaaWcbaGaemOrayeabeaaaaaaaa@3D5D@ ) and the effective size of an X-linked locus ( N e X = 9 N M N F 4 N M + 2 N F MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGobGtdaWgaaWcbaGaemyzau2aaSbaaWqaaiabdIfaybqabaaaleqaaOGaeyypa0ZaaSaaaeaacqaI5aqocqWGobGtdaWgaaWcbaGaemyta0eabeaakiabd6eaonaaBaaaleaacqWGgbGraeqaaaGcbaGaeGinaqJaemOta40aaSbaaSqaaiabd2eanbqabaGccqGHRaWkcqaIYaGmcqWGobGtdaWgaaWcbaGaemOrayeabeaaaaaaaa@3F7D@ ) depends on the effective numbers of breeding males (NM) and breeding females (NF) in the system [42]. When there are equal numbers of breeding males and females, the expected ratio of effective sizes of the X and the autosomes is 3/4. To achieve a X/A ratio of 1, we can calculate from these equations that N M N F MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWccaqaaiabd6eaonaaBaaaleaacqWGnbqtaeqaaaGcbaGaemOta40aaSbaaSqaaiabdAeagbqabaaaaaaa@31A2@ .

Consequently, as the ratio of effectively reproducing females and males (F/M) increases from 1 to 7, the expected θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ ratio increases 0.75 to 1, and can increase to as much as 1.125 as the F/M ratio is elevated even further.

It is also possible that the relative elevation of variability on the X-chromosome is due to different patterns of natural selection acting on the X-linked and autosomal genes. Background selection could in principle generate this pattern as well, as this model predicts higher levels of neutral variation on the X chromosome [43-46]. This is a consequence of the greater efficacy of background selection on X chromosome due to the hemizygosity of the X chromosome in males, which results in a larger proportion of X chromosomes that are free of deleterious mutations. Whether background selection, especially in regions of high recombination, could be sufficiently strong to elevate the θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ substantially and, more specifically, from 0.75 to 1 is unclear.

Finally, it is possible that the ratio of X-linked to autosomal diversity is greater than expected because of a reduction in levels of polymorphism on the autosomes rather than an increase in polymorphism on the X chromosome. It is possible that adaptation is more frequent and/or involves adaptive events of greater selective strength on the autosomes than on the X chromosome. This would reduce autosomal variation relative to that on the X chromosome, and is theoretically possible if adaptive evolution operates primarily on standing variation as opposed to novel beneficial mutations [47].

It has also been suggested that a reduction in autosomal variation could result from the presence of polymorphic inversions [22]; because inversion heterozygosity could in theory suppress recombination and thus increase linkage among sites, the reduction in diversity due to selection may be more pronounced in genomic regions with inversions. Moreover, levels of variability within inversions of young to intermediate age are lower than expected [48,49]. Inversions do appear to be more common on the autosomes than on the X chromosome in African populations of D. melanogaster [50] which may contribute to a reduction in autosomal diversity relative to expectation in these populations. Consistent with this model, four of the eight autosomal intergenic regions presented here are contained within the breakpoints of five of the large inversions that are frequent in Africa (In(2L)t, In(2R)NS, In(3R)K, In(3R)P and In(3L)P) [51]; levels of polymorphism at these regions is approximately 1/2 that found in the four intergenic regions outside of these known inversions in the African population (data not shown). However, it is important to note that if an inversion is sufficiently old and is maintained as a balanced polymorphism, it may also serve to increase standing diversity, acting as a balancer and thereby preventing the coalescence of inverted and non-inverted chromosomes. Thus, the effects of inversions on levels of molecular variation in Drosophila merits further investigation.

X-linked and autosomal variation in non-African populations

The results from the non-African populations collected for the present study and from the meta-analysis are more difficult to interpret. Both sets of results support a reduction in diversity in these populations, although the magnitude of the reduction on the X and the autosomes differs substantially between the two datasets. The regions within or nearby genes collated and reanalyzed here argue for a much more severe reduction on the X chromosome than on autosomes (Table 4), which, in addition to being largely consistent with the individual studies in the meta-analysis [4,5,22,23,52], are also consistent with studies of microsatellite variability [29]. In contrast, the analysis of sixteen intergenic regions presented in this paper shows very similar reductions in diversity on the X chromosome and the autosomes. This is primarily due to a greater reduction in X-linked diversity in the compiled dataset, as reductions in diversity at autosomal loci are similar between the two datasets.

To reconcile the differences between these two datasets, we considered several alternative demographic scenarios, and assessed which aspects of the polymorphism data were consistent with each model. First, we considered a model with equal effective population sizes of the X and the autosomes in the ancestral population in combination with a population bottleneck during and after which the ratio of effective sizes of the X to the autosomes is unchanged. Under such a scenario, we expect similar levels of X-linked and autosomal variation in ancestral populations as well as in derived populations. Thus, this model is entirely consistent with the patterns of X-linked and autosomal variability from both African and non-African strains at the intergenic loci presented in the current study, and is similarly consistent with patterns of variation at X-linked and autosomal loci in African populations of the compiled dataset. In this case, the reduced diversity at X-linked loci seen in derived populations in the compiled dataset must be explained. One possibility for the reduction in diversity on the X chromosome in derived strains is genetic hitchhiking; under this model nucleotide diversity on the X is lower than autosomal diversity due to more frequent or more recent selective sweeps on the X chromosome relative to the autosomes, perhaps due to adaptation to temperate climates in these derived strains.

A second plausible demographic model is one in which the ancestral population has similar levels of autosomal and X-linked variation, but that this ratio is shifted towards a higher effective size of the autosomes relative to the X during and after the bottleneck. Under this model, we expect that levels of X-linked diversity will be more severely reduced than levels of diversity on the autosomes; this demographic model is thus consistent with all of the African polymorphism data presented here as well as the non-African polymorphism data presented in the combined dataset.

Therefore, the increased X-linked diversity (or decreased autosomal diversity) in North American populations at the intergenic regions presented here requires explanation. One possibility is that patterns of background selection at intergenic loci have changed between putatively ancestral and derived populations. As described earlier, the greater efficacy of background selection at X-linked loci can increase standing variation in these regions relative to the autosomes, which suggests that if background selection plays a larger role at intergenic loci in derived populations, then the reduction of nucleotide diversity due to background selection would be more pronounced on the autosomes than on the X chromosome.

Finally, we considered a demographic model in which the θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ ratio in the ancestral population is less than 1, coupled with a population bottleneck. This scenario also leads to more extreme reductions of diversity at X-linked loci, which is what is observed in the non-African strains at the loci in the combined dataset. It is not clear whether biologically plausible scenarios under this model can generate a θ X c h r o m o s o m e θ A u t o s o m e s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaGGaciab=H7aXnaaBaaaleaacqWGybawcqGHsislcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBcqWGTbqBcqWGVbWBcqWGZbWCcqWGVbWBcqWGTbqBcqWGLbqzaeqaaaGcbaGae8hUde3aaSbaaSqaaiabdgeabjabdwha1jabdsha0jabd+gaVjabdohaZjabd+gaVjabd2gaTjabdwgaLjabdohaZbqabaaaaaaa@4CD5@ ratio as low as 0.37, as observed in the meta-analysis dataset, but these simulation results suggest that an exaggerated depression of diversity at X-linked loci can be at least partly explained by demographic scenarios. This model thus requires that either extant levels of sequence variation in African strains are not reflective of ancestral polymorphism levels or that the ancestral populations of this species have not been sampled in the datasets presented here. Moreover, this model cannot explain the similar levels of sequence variation at X-linked and autosomal loci found at the intergenic regions collected for the current study; the excess of X-linked polymorphism or dearth of autosomal variation in these regions could result from a change in the patterns of background selection during and after the population bottleneck.

It is indeed difficult to distinguish between these alternative explanations for the differences between the datasets presented here. Based on the currently available data, however, we believe the demographic model of a female-biased ancestral sex ratio coupled with a population bottleneck that retains this ratio to be most appropriate for two reasons. First, the increased levels of diversity at the intergenic loci chosen in the current study as compared to the levels of diversity at loci in the compiled dataset (Figure 2) are consistent with reduced effects of background selection and genetic hitchhiking in these regions. This suggests that patterns of variation at these regions are more likely to reflect demographic rather than selective processes. Second, to the extent that extant levels of sequence variation in African strains are indicative of ancestral levels of polymorphism, the comparable levels of X-linked and autosomal nucleotide diversity in sampled African populations of D. melanogaster are suggestive of an unequal sex ratio in ancestral populations of this species.

Conclusion

On balance, therefore, we believe that given the currently available data, a demographic model with unequal numbers of breeding males and females may be the most appropriate explanation for the current data. This model is suggestive of a ratio of breeding females to breeding males of 7:1, and while it is not beyond reason that males could have a larger variance in reproductive success than females in Drosophila, it is not clear whether a ratio of 7:1 is realistic. Indeed, it is challenging to assess the plausibility of this ratio, given the lack of relevant parameter estimates from natural populations of this species.

It is important to note that since this analysis was undertaken, additional polymorphism datasets have been made available [3,53]. We believe our results are likely robust to the inclusion of the new X-linked loci presented in [3], as our current meta-analysis includes almost one-half of the loci in the more recent study, and no major differences between the initial [4] and later [3] X-linked datasets were reported. In addition, the recent survey of autosomal polymorphism, while supportive of a female-biased sex ratio in Africa, seem to suggest a male-biased sex ratio in Europe [53]. Thus, it remains to be seen whether the patterns we observe are characteristic of the genome and moreover, whether the patterns in non-African strains vary with population of origin; as more polymorphism data becomes available, we will have increasing power to test this explicitly.

The demographic model we present assumes that the effective numbers of breeding males and females remains unchanged throughout the demographic history of this species, and can fully account for the patterns of diversity at X-linked and autosomal loci in African and North American strains of D. melanogaster from the intergenic regions presented here. The exaggerated reduction of diversity of X-linked regions in derived populations contained within the compiled dataset unfortunately remains difficult to explain. One possibility is that this reduction is due to increased hitchhiking on the X chromosome in these populations; that divergence at X-linked loci is significantly higher than autosomal divergence in the compiled dataset (P = 0.002, Mann-Whitney U-test) supports this hypothesis.

Importantly, the hitchhiking aspect of this model must be discussed in the context of several previous observations. First, X-linkage does not systematically increase rates of molecular evolution; while some previous studies do support a faster-X model [54,55], the most recent and most comprehensive test of this model shows no evidence for faster-X evolution [56]. The polymorphism data presented here, however, hint at the possibility that the effects of genetic hitchhiking may be more pronounced on the D. melanogaster X chromosome than on the autosomes in association with population expansion to new habitats. These results could be reconciled either by suggesting that the increase in rates of positive selection is too episodic to elevate the rate of evolution overall on the X chromosome, or that X-linked loci are also subject to stronger purifying selection. More frequent fixation of allelic variants by positive selection could be balanced by a lower rate of fixation of slightly deleterious alleles on the X chromosome.

Second, if our results are indeed reflective of recent adaptive events on the X chromosome, we do expect to observe the signature of selective sweeps on this chromosome. While previous studies have identified regions of the X-chromosome that have patterns of polymorphism consistent with selective sweeps, patterns of variation at these regions are also largely consistent with demographic effects (for review see [57]). Indeed, identifying loci whose patterns of sequence variation cannot be explained by demography alone depends on the demographic model employed and the population bottleneck parameters (for review see[57]). Thus, further investigation of the effects of population bottlenecks under a variety of bottleneck parameters on X-linked and autosomal variation, as well as additional examination of what bottleneck parameters may be appropriate for D. melanogaster may be warranted.

Another observation that merits discussion is the estimation of selective constraint in the D. melanogaster genome. Recent studies of constraint in Drosophila appear to provide evidence that in addition to genic regions such as introns and synonymous sites, intergenic regions also appear to be subject to selective constraint [58-60]. This bears particular relevance given that our simulations are based on the evolution of neutral sequences; to the extent that none of the regions presented here, intergenic or otherwise, are truly neutral, the results from our simulations may not prove accurate. Consequently, modeling the joint effects of demography and selection on X-linked and autosomal variation should prove illuminating.

Importantly, this model of higher rates of adaptation at X-linked regions in close proximity to genes in derived populations makes several predictions that can be explicitly tested in the future. For instance, a model of recurrent selection will only produce decreased nucleotide variation on the X chromosome in a species such as D. melanogaster (with recombination in females only) if natural selection acts on new mutations and these new mutations are at least partially recessive [61]. This is of particular importance given that the majority of polymorphic sites in non-African strains are polymorphic in African strains as well [4]. However, there are certainly mutations that are novel to derived strains of D. melanogaster, and given a specific scenario of the population expansion out of Africa, it may be possible to theoretically estimate the beneficial mutation rate required to produce such a depression of X-linked polymorphism and assess whether such estimates are biologically plausible.

In addition, this model predicts systematic differences in the site-frequency spectrum between the X and the autosomes; more frequent or more recent episodes of adaptive evolution should produce a skew in the frequency spectrum towards an excess of rare variants on the X relative to the autosomes [36]. Comprehensive analysis of the frequency spectrum of segregating variants at a large collection of X-linked and autosomal loci will thus shed light on the potential relevance of this explanation.

Further, this model predicts that rates of adaptation are higher in derived populations than in ancestral populations. Given the currently available data, this does not appear to be the case, at least for X-linked loci[57], but additional data would certainly make the estimates of the rate of adaptive substitution more accurate, and may reveal differences between ancestral and derived populations, or between X-linked and autosomal loci. In addition, this model predicts that there should be marked divergence between African and non-African strains, particularly for X-linked loci, due to increased rates of adaptation in derived strains. Further investigation of this question will also shed light on the viability of this model.

Moreover, given the effects of genetic hitchhiking on sequence diversity, this model predicts that bottleneck parameters estimated using X chromosome data will differ from those parameters estimated using autosomal data. The most recent estimates of these parameters from D. melanogaster [35,62] are indeed based on a dataset which, while wonderfully rich, is comprised entirely of X-linked loci. There is little question that this dataset of X-linked loci is the most comprehensive dataset available for such analyses, but if there is a greater incidence of hitchhiking on the X in derived populations, the severity of the bottleneck may have been overestimated. A perhaps fruitful avenue of future inquiry would be to reevaluate the estimation of bottleneck parameters using only autosomal data, should such data become available.

Finally, if this model is correct, it will have a number of important implications. First, patterns of variation of X-linked and autosomal diversity depend heavily on sequence type, which suggests that the evolutionary forces responsible for the generation and maintenance of sequence diversity are heterogeneous across the genome, with marked differences between sequences that are either genic or in close proximity to genic sequence and regions that are distant to genes. Second, the similarity of the effective population sizes of the X chromosome and autosomes in the African populations should make natural selection more effective for X-linked loci in general, at least in these populations. For example, for identical co-dominant mutations on the X chromosome and the autosomes we expect that the absolute value of the effective strength of selection (Nes) to be 33% larger on the X chromosome [30]. It is possible that this effect is sufficient to account for the higher codon bias of the X-linked genes in D. melanogaster [30-32]. The higher codon bias of the X-linked genes in D. pseudoobscura [30] might imply that effective population size of the X chromosomes in general is higher than 3/4 of the autosomal effective population size throughout the Drosophila genus; this effect could be due to a generally higher variance of reproductive success in males than females throughout Drosophila.

Methods

Selection of Loci

We selected eight autosomal and eight X-linked loci from noncoding regions of the D. melanogaster genome. We began with all noncoding loci (not contained within the coordinates of coding sequences as reported in version 3.1 of the D. melanogaster genome) that were at least 5 kb from annotated genes in both directions, as measured from the midpoint of the noncoding locus to the edge (5' or 3' end, depending on the genic orientation) of the nearest coding sequence. We restricted these candidates to those with no significant BLASTX hits with an E-value cutoff of 0.001 in the hopes of restricting ourselves to loci that were truly noncoding. We further limited ourselves to loci located in regions of high recombination, with estimates of local recombination rate between 3.25 and 3.5 cM/Mb [33] (see also [63]). We allowed for a maximum GC content of 45%, and used only regions for which there was an Exelixis deletion strain encompassing the entire locus available (see below). From this list of candidate loci, we selected two regions from each of the four major autosomal chromosome arms (2L, 2R, 3L, and 3R) as well as eight regions from the X chromosome, maximizing the physical distance among loci for each chromosome arm. Finally, we computed sequence divergence between each locus and the orthologous locus in D. simulans (based on the consensus sequence) and used these values to identify regions that were evolving as close as possible to expectation under neutrality. We chose X-linked and autosomal loci such that the distributions of divergence values for our autosomal and X-linked loci were not significantly different (P = 0.59, Mann-Whitney U-test) to control for variation in mutation rate.

Fly Strains and Genomic DNA Extraction

We have used 12 North American strains (6 from North Carolina (gift from G. Gibson) and 6 from Davis, California (gift from S. Nuzhdin)) and 12 African strains from Malawi. Because the American strains had been passed through > 30 generations of brother-sister mating (S. Nuzhdin and G. Gibson, personal communications), there was very little nucleotide variability within in each strain. We therefore assumed and then confirmed that individuals from these strains were nearly homozygous at each of our studied loci.

In contrast, the strains from Malawi had a great deal of nucleotide diversity and measures were taken to effectively homozygose these strains at each of our eight autosomal loci. We took advantage of the D. melanogaster deficiency collection created by Exelixis [64] and purchased eight deficiency strains, each of which had a deletion at one of the studied autosomal loci. We then crossed each of our twelve Malawi strains with each of the 8 deletion strains and used the appropriate 96 F1's for our analyses. These F1's had only a single copy of one of the parental Malawi haplotypes at the studied locus.

We extracted genomic DNA from a single male from each of the strains according to a protocol described by G. Gloor and W. Engels (personal communication). Each fly was crushed with the end of a pipette tip and subsequently immersed in a buffered solution (10 mM Tris-Cl pH 8.2, 1 mM EDTA, 25 mM NaCl, 200 μg/mL proteinase K). This solution was incubated at 37°C for 30 minutes and then at 95°C to inactivate the proteinase K.

Primer Design, PCR Amplification and Sequencing

We designed primers for our sixteen loci using Primer3 [65] and ensured that each primer pair would produce a unique amplicon using Virtual PCR [66]. Primer sequences are available upon request. Our PCR reactions were 20 μL and contained 10 μL of RedTaq Ready Mix (Qiagen), 1 μL of each 20 μM primer, 1 μL of genomic DNA and 7 μL water. We used a touchdown PCR program to amplify our loci to facilitate direct sequencing. The amplifying conditions were as follows: 94°/30 seconds (s), Tm/30s, 72°/30s, where Tm began at 62° and was decreased by 0.5° each cycle, followed by 20 cycles of 94°/30s, 55°/30s, 72°/30s, with a final 10 minute extension period at 72°. These PCR products were sequenced directly at Genaissance Pharmaceuticals using our amplifying primers.

Data Analysis

Sequences were manually aligned using Sequencher version 4.2.2, and values of Tajima's D were calculated using DnaSP version 3.1. Statistical analyses were performed in Statview version 5.0. Mean levels of nucleotide diversity (θs) and confidence intervals for these estimates were obtained by bootstrapping over all sites in each respective dataset in order to compare diversity across loci and populations; bootstrapping was done with scripts that were developed in-house for the purposes of this project. The values reported are the means of their respective bootstrap distributions. Each confidence interval derives from 1000 bootstrap replicates in the case of the intergenic dataset and 10,000 replicates in the case of the compiled dataset.

Meta-Analysis

We also compiled polymorphism data from several previously published results. Included in our meta-analysis are seven intergenic loci from Kern and Begun [5], 20 coding regions from Andolfatto [22], 10 coding regions from Haddrill et al. [52], 105 noncoding regions from Glinka et al. [4] and 8 coding regions from Sheldahl et al. [23] (Table 3). These data include both X-linked and autosomal loci. For the purpose of calculating θs and θs /divergence, loci which were available in more than one strain were given equal weights, such that total weight of the locus was 1. Both the observed values and bootstrap values of θs and θs /divergence were calculated using these weights. Estimates of polymorphism were taken directly from the literature, as were estimates of divergence with D. simulans if reported. If divergence estimates were not presented, we retrieved the studied sequences and calculated divergence against the consensus sequence of the D. simulans genome. We also calculated local recombination rate at each locus following a procedure detailed previously [33].

Simulations

We conducted coalescent simulations to estimate the expected effect of a bottleneck under different demographic scenarios. Simulations were carried out using ms [67]. The bottleneck parameters for the X chromosome are taken from Thornton and Andolfatto [35]. Specifically we take the population bottleneck on the X chromosome to last from 0.022 to 0.0042 in units of the expected time to coalescence between two alleles. During the bottleneck the effective population size is reduced to 0.029 of that before and after the bottleneck. The simulation command line for the scenario in which F/M ratio goes from 7 in the ancestral population to 1 during and after the bottleneck is:

X: ms 20 100000 -t 40.0 -r 40.0 10000 -eN 0.0042 0.029 -eN 0.022 1.0

A: ms 20 100000 -t 53.4 -r 53.4 10000 -eN 0.00315 0.029 -eN 0.0165 0.75

The simulation code for the scenario in which F/M ratio stays at 1 throughout the simulations is:

X: ms 20 10000 -t 30.0 -r 30.0 10000 -eN 0.0042 0.029 -eN 0.022 1.0

A: ms 20 10000 -t 40.0 -r 40.0 10000 -eN 0.00315 0.029 -eN 0.0165 1.0

Authors' contributions

NDS and DAP conceived the idea for this project. JMM and DAP carried out coalescent simulations. NDS, JMM and DAP participated in the analysis and contributed to the interpretation of the simulation results. NDS, JMM, JDJ and DAP contributed to the interpretation of the empirical results. NDS, JMM, and DAP were responsible for the writing of this manuscript, and NDS, JMM, JDJ and DAP were responsible for the editing of this manuscript. All authors have given their final approval for its publication.

Acknowledgements

The authors gratefully acknowledge K. Thornton and P. Andolfatto for helpful discussions on this manuscript. Comments from the five reviewers significantly improved this manuscript as well. The authors also extend thanks to J. Hatzidakis for her assistance with the flywork required for this experiment. This work was supported in part by the NSF (DEB-031717) and NIH (GM077368) grants (DAP), Stanford Genome Training Program funded by 5 T32 HG00044 from the NHGRI (NDS). JMM was supported by a Howard Hughes predoctoral fellowship. JDJ was supported by National Science Foundation grant DMS-0201037 to R. Durrett, C.F. Aquadro and R. Nielsen.

References

  1. Lachaise D, Cariou ML, David JR, Lemeunier F, Tsacas L, Ashburner M: Historical biogeography of the Drosophila melanogaster species subgroup.

    Evol Biol 1988, 22:159-225. OpenURL

  2. Mousset S, Derome N: Molecular polymorphism in Drosophila melanogaster and D. simulans: what have we learned from recent studies?

    Genetica 2004, 120(1–3):79-86. PubMed Abstract | Publisher Full Text OpenURL

  3. Ometto L, Glinka S, De Lorenzo D, Stephan W: Inferring the effects of demography and selection on Drosophila melanogaster populations from a chromosome-wide scan of DNA variation.

    Molec Biol Evol 2005, 22:2119-2130. PubMed Abstract | Publisher Full Text OpenURL

  4. Glinka S, Ometto L, Mousset S, Stephan W, De Lorenzo D: Demography and natural selection have shaped genetic variation in Drosophila melanogaster: A Multi-locus approach.

    Genetics 2003, 165:1269-1278. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Kern AD, Begun DJ: Patterns of polymorphism and divergence from noncoding sequences of Drosophila melanogaster and D. simulans: Evidence for nonequilibrium processes.

    Molec Biol Evol 2005, 22(1):51-62. PubMed Abstract | Publisher Full Text OpenURL

  6. Hill WG, Robertson A: The effect of linkage on limits to artificial selection.

    Genet Res 1966, 8:269-294. PubMed Abstract OpenURL

  7. Felsenstein J: The evolutionary advantage of recombination.

    Genetics 1974, 78:737-756. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Choudhary M, Singh RS: A Comprehensive Study of Genic Variation in Natural Populations of Drosophila melanogaster III. Variations in Genetic Strucutre and Their Causes Between Drosophila melanogaster and Its Sibling Species Drosophila simlans.

    Genetics 1987, 117:697-710. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Begun DJ, Aquadro CF: African and North American populations of Drosophila melanogaster are very different at the DNA level.

    Nature 1993, 365:548-550. PubMed Abstract | Publisher Full Text OpenURL

  10. Singh RS, Choudhary M, David JR: Contrasting patterns of geographic variation in cosmopolitan sibling species D. melanogaster and D. simulans.

    Biochem Genet 1987, 25:27-40. PubMed Abstract | Publisher Full Text OpenURL

  11. Singh RS, Hickey DA, David JR: Genetic differentiation between geographically distinct populations of Drosophila melanogaster.

    Genetics 1982, 101:235-256. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Begun DJ, Aquadro CF: Molecular variation at the vermilion locus in geographically diverse populations of Drosophila melanogaster and D. simulans.

    Genetics 1995, 140(3):1019-1032. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Langley CH, Lazzaro BP, Phillips W, Heikkinen E, Braverman JM: Linkage disequilibria and the site frequency spectra in the su(s) and su(wa) regions of the Drosophila melanogaster X chromosome.

    Genetics 2000, 156(4):1837-1852. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Begun DJ, Betancourt AJ, Langley CH, Stephan W: Is the fast/slow allozyme variation at the Adh locus of Drosophila melanogaster an ancient balanced polymorphism?

    Mol Biol Evol 1999, 16(12):1816-1819. PubMed Abstract | Publisher Full Text OpenURL

  15. Aguade M: Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

    Genetics 1999, 152(2):543-551. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Aguade M: Different forces drive the evolution of the Acp26Aa and Acp26Ab accessory gland genes in the Drosophila melanogaster species complex.

    Genetics 1998, 150(3):1079-1089. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Clark AG, Wang L: Molecular population genetics of Drosophila immune system genes.

    Genetics 1997, 147(2):713-724. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Andolfatto P, Kreitman M: Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans.

    Genetics 2000, 154(4):1681-1691. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Tsaur SC, Ting CT, Wu CI: Positive selection driving the evolution of a gene of male reproduction, Acp26Aa, of Drosophila: II. Divergence versus polymorphism.

    Mol Biol Evol 1998, 15(8):1040-1046. PubMed Abstract | Publisher Full Text OpenURL

  20. Balakirev ES, Ayala FJ: Molecular population genetics of the beta-esterase gene cluster of Drosophila melanogaster.

    Journal of Genetics 2003, 82(3):115-131. PubMed Abstract | Publisher Full Text OpenURL

  21. Balakirev ES, Ayala FJ: Nucleotide variation of the Est-6 gene region in natural populations of Drosophila melanogaster.

    Genetics 2003, 165(4):1901-1914. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Andolfatto P: Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans.

    Mol Biol Evol 2001, 18(3):279-290. PubMed Abstract | Publisher Full Text OpenURL

  23. Sheldahl LA, Weinreich DM, Rand DM: Recombination, dominance and selection on amino acid polymorphism in the Drosophila genome: Contrasting patterns on the X and fourth chromosomes.

    Genetics 2003, 165(3):1195-1208. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Hale LR, Singh RS: A comprehensive study of genic variation in natural populations of Drosophila melanogaster. IV. Mitochondrial DNA variation and the role of history vs. selection in the genetic structure of geographic populations.

    Genetics 1991, 129(1):103-118. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Harr B, Zangerl B, Brem G, Schlotterer C: Conservation of locus-specific microsatellite variability across species: A comparison of two Drosophila sibling species, D. melanogaster and D. simulans.

    Mol Biol Evol 1998, 15(2):176-184. PubMed Abstract | Publisher Full Text OpenURL

  26. Schlotterer C, Vogl C, Tautz D: Polymorphism and locus-specific effects on polymorphism at microsatellite loci in natural Drosophila melanogaster populations.

    Genetics 1997, 146(1):309-320. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Irvin SD, Wetterstrand KA, Hutter CM, Aquadro CF: Genetic variation and differentiation at microsatellite loci in Drosophila simulans: Evidence for founder effects in new world populations.

    Genetics 1998, 150(2):777-790. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Schoefl G, Schloetterer C: Patterns of microsatellite variability among X chromosomes and autosomes indicate a high frequency of beneficial mutations in non-african D. simulans.

    Mol Biol Evol 2004, 21(7):1384-1390. PubMed Abstract | Publisher Full Text OpenURL

  29. Kauer M, Zangerl B, Dieringer D, Schlotterer C: Chromosomal patterns of microsatellite variability contrast sharply in African and non-African populations of Drosophila melanogaster.

    Genetics 2002, 160(1):247-256. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Singh ND, Davis JC, Petrov DA: X-linked genes evolve higher codon bias in Drosophila and Caenorhabditis.

    Genetics 2005, 171:145-155. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Comeron JM, Kreitman M, Aguade M: Natural Selection on Synonymous Sites is Correlated with Gene Length and Recombination in Drosophila.

    Genetics 1999, 151:239-249. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Hambuch TM, Parsch J: Patterns of Synonymous Codon Usage in Drosophila melanogaster Genes with Sex-biased Expression.

    Genetics 2005, 170:1691-1700. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Singh ND, Arndt PF, Petrov DA: Genomic Heterogeneity of Background Substitutional Patterns in Drosophila melanogaster.

    Genetics 2004, 169:709-722. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.

    Genetics 1989, 123(3):585-596. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Thornton K, Andolfatto P: Approximate Bayesian inference reveals evidence for a recent, severe bottleneck in a Netherlands population of Drosophila melanogaster.

    Genetics 2006, 172:1607-1619. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Wall JD, Andolfatto P, Przeworski M: Testing models of selection and demography in Drosophila simulans.

    Genetics 2002, 162(1):203-216. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Maynard Smith J, Haigh J: The hitch-hiking effect of a favourable gene.

    Genet Res 1974, 23:23-35. PubMed Abstract OpenURL

  38. Thornton KR, Jensen JD, Bequet C, Andolfatto P: Progress and prospects in mapping recent selection in the genome.

    Heredity 2007, in press. PubMed Abstract | Publisher Full Text OpenURL

  39. Thornton KR, Jensen JD: Controlling the false positive rate in multilocus genome scans for selection.

    Genetics 2006, 175(2):737-750.

    genetics.106.064642

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Pool JE, Aquadro CF: History and structure of sub-saharan populations of Drosophila melanogaster.

    Genetics 2006, 174(2):915-929. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Baudry E, Virginier B, Veuille M: Non-African populations of Drosophila melanogaster have a unique origin.

    Mol Biol Evol 2006, 21(8):1482-1491. Publisher Full Text OpenURL

  42. Hartl DL, Clark AG: Principles of Population Genetics. Sunderland, Massachusetts: Sinauer Associates, Inc; 1989.

  43. Aquadro CF, Begun DJ, Kindahl EC: Selection, recombination, and DNA polymorphism in Drosophila. In Non-Neutral Evolution. Edited by Golding B. London: Chapman & Hall; 1994:46-56. OpenURL

  44. Begun DJ, Whitley P: Reduced X-linked nucleotide polymorphism in Drosophila simulans.

    Proceedings of the National Academy of Sciences of the United States of America 2000, 97(11):5960-5965. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Charlesworth B: Background selection and patterns of genetic diversity in Drosophila melanogaster.

    Genet Res 1996, 68(2):131-149. PubMed Abstract OpenURL

  46. Charlesworth B, Morgan MT, Charlesworth D: The effect of deleterious mutations on neutral molecular variation.

    Genetics 1993, 134:1289-1303. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Orr HA, Betancourt AJ: Haldane's Sieve and Adaptation From the Standing Genetic Variation.

    Genetics 2001, 157:875-884. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Andolfatto P, Depaulis F, Navarro A: Inversion polymorphisms and nucleotide variability in Drosophila.

    Genet Res 2001, 77(1):1-8. PubMed Abstract | Publisher Full Text OpenURL

  49. Navarro A, Barbadilla A, Ruiz A: Effect of Inversion Polymorphism on the Neutral Nucleotide Variability of Linked Chromosomal Regions in Drosophila.

    Genetics 2000, 155:685-698. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Lemeunier F, Aulard S: Inversion polymorphism frequency in Drosophila melanogaster. In Drosphila inversion polymorphism. Edited by Krimbas CB, Powell JR. Boca Raton, FL: CRC Press; 1992:309-405. OpenURL

  51. Aulard S, David JR, Lemeunier F: Chromsomal inversion polymorphism in Afrotropical populations of Drosophila melanogaster.

    Genet Res 2002, 79(1):49-63. PubMed Abstract | Publisher Full Text OpenURL

  52. Haddrill PR, Thornton KR, Charlesworth B, Andolfatto P: Multilocus patterns of nucleotide variability and the demographic and selection history of Drosophila melanogaster populations.

    Genome Res 2005, 15:790-799. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Hutter S, Li H, Beisswanger S, De Lorenzo D, Stephan W: Distinctly Different Sex Ratios in African and European Populations of Drosophila melanogaster Inferred From Chromosomewide Single Nucleotide Polymorphism Data.

    Genetics 2007, 177(1):469-480. PubMed Abstract | Publisher Full Text OpenURL

  54. Counterman BA, Ortiz-Barrientos C, Noor MAF: Using Comparative Genomic Data to Test for Fast-X Evolution.

    Evolution 2004, 58(3):656-660. PubMed Abstract OpenURL

  55. Thornton K, Long M: Rapid divergence of gene duplicates on the Drosophila melanogaster X chromosome.

    Mol Biol Evol 2002, 19(6):918-925. PubMed Abstract | Publisher Full Text OpenURL

  56. Thornton K, Bachtrog D, Andolfatto P: X chromosomes and autosomes evolve at similar rates in Drosophila: no evidence for faster-X protein evolution.

    Genome Res 2006, 16(4):498-504.

    gr.4447906

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Li H, Stephan W: The recent demographic and adaptive history of Drosophila melanogaster.

    Heredity 2007, 98:65-68. PubMed Abstract | Publisher Full Text OpenURL

  58. Andolfatto P: Adaptive evolution of non-coding DNA in Drosophila.

    Nature 2005, 437:1149-1152. PubMed Abstract | Publisher Full Text OpenURL

  59. Bachtrog D, Andolfatto P: Selection, recombination and demographic history in Drosophla miranda.

    Genetics 2006, in press. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Halligan DL, Keightley PD: Ubiquitous selective constraints in the Drosophila genome revealed by genome-wide interspecies comparison.

    Genome Res 2006, 16:875-884. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Betancourt AJ, Kim Y, Orr HA: A pseudohitchhiking model of × vs. autosomal diversity.

    Genetics 2004, 168(4):2261-2269. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Li H, Stephan W: Inferring the demographic history and rate of adaptive substitution in Drosophila.

    PLoS Genet 2006, 2:e166. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. The Drosophila melanogaster recombination rate calculator [http://cgi.stanford.edu/~lipatov/recombination/recombination-rates.txt] webcite

  64. Parks AL, Cook KR, Belvin M, Dompe NA, Fawcett R, Huppert K, Tan LR, Winter CG, Bogart KP, Deal JE, et al.: Systematic generation of high-resolution deletion coverage of the Drosophila melanogaster genome.

    Nat Genet 2004, 36(3):288-292. PubMed Abstract | Publisher Full Text OpenURL

  65. Primer3 [http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi] webcite

  66. Virtual PCR [http://grup.cribi.unipd.it/cgi-bin/mateo/vpcr2.cgi] webcite

  67. Hudson RR: Generating samples under a Wright-Fisher neutral model.

    Bioinformatics 2002, 18:337-338. PubMed Abstract | Publisher Full Text OpenURL