Selective sweep at the Drosophila melanogaster Suppressor of Hairless locus and its association with the In(2L)t inversion polymorphism.

The hitchhiking model of population genetics predicts that an allele favored by Darwinian selection can replace haplotypes from the same locus previously established at a neutral mutation-drift equilibrium. This process, known as "selective sweep," was studied by comparing molecular variation between the polymorphic In(2L)t inversion and the standard chromosome. Sequence variation was recorded at the Suppressor of Hairless (Su[H]) gene in an African population of Drosophila melanogaster. We found 47 nucleotide polymorphisms among 20 sequences of 1.2 kb. Neutrality tests were nonsignificant at the nucleotide level. However, these sites were strongly associated, because 290 out of 741 observed pairwise combinations between them were in significant linkage disequilibrium. We found only seven haplotypes, two occurring in the 9 In(2L)t chromosomes, and five in the 11 standard chromosomes, with no shared haplotype. Two haplotypes, one in each chromosome arrangement, made up two-thirds of the sample. This low haplotype diversity departed from neutrality in a haplotype test. This pattern supports a selective sweep hypothesis for the Su(H) chromosome region.     

High nucleotide sequence variation in a region of low recombination in Drosophila simulans is consistent with the background selection model

We surveyed nucleotide sequence variation at glucose dehydrogenase (Gld), in a region of low recombination on chromosome 3R, from a population sample of Drosophila simulans. The levels of nucleotide variation were surprisingly high. There was no departure from the expectation of a neutral model for the level of polymorphism, indicating no evidence of a selective sweep in this region. There was a significant deficiency of singleton polymorphisms according to the Fu and Li test, although Tajima and Hudson, Kreitman, and Aguade (HKA) tests do not provide evidence of a significant elevation of variation due to balancing selection. Genetic map data for the D. simulans third chromosome were used to calculate expected values of pi for Gld under a current model of background selection, varying the values for the parameter sh (selection coefficient against deleterious mutations). We show that the recombinational landscape of D. simulans is sufficiently different from that of D. melanogaster that we expect higher variation under the background selection model, even when effective population sizes are assumed to be equal. The data for Gld were tested against the predictions using computer simulations of the distribution of the number of segregating sites conditioned on pi. Background selection alone can explain our observations as long as sh is larger than 0.005 and species-level effective population size is assumed to be several- fold larger than in D. melanogaster. Alternatively, the deleterious mutation rate may be smaller in D. simulans, or balancing selection may be acting nearby, thereby reducing the effect of background selection.      

Regulation of genome stability by TEL1 and MEC1, yeast homologs of the mammalian ATM and ATR genes

DNA sequence surveys of Drosophila melanogaster populations show a strong positive correlation between the recombination rate experienced by a locus and its level of nucleotide polymorphism. In particular, surveys of the fourth chromosome gene ci(D) show greatly reduced levels of nucleotide variation; this observation was originally interpreted in terms of selective sweeps occurring on the nonrecombining fourth chromosome. Subsequent theoretical work has, however, uncovered several other selective processes that can reduce variation. In this study, we revisit the Drosophila fourth chromosome, investigating variation in 5-6 kb of the gene ankyrin in D. melanogaster and D. simulans. Silent nucleotide site diversity is approximately 5 x 10(-4) for both species, consistent with the previous observations of low variation at ci(D). Given the observed frequency spectra at ankyrin, coalescent simulations indicate that reduced diversity in the region is unlikely to be due to a selective sweep alone. We find evidence for recombinational exchange at this locus, and both species appear to be fixed for an insertion of the transposable element HB in an intron of ankyrin.     

Mechanisms regulating the copy numbers of six LTR retrotransposons in the genome of Drosophila melanogaster

There has been debate over the mechanisms that control the copy number of transposable elements in the genome of Drosophila melanogaster. Target sites in D. melanogaster populations are occupied at low frequencies, suggesting that there is some form of selection acting against transposable elements. Three main theories have been proposed to explain how selection acts against transposable elements: insertions of a copy of a transposable element are selected against; chromosomal rearrangements caused by ectopic exchange between element copies are selected against; or the process of transposition itself is selected against. The three theories give different predictions for the pattern of transposable element insertions in the chromosomes of D. melanogaster. We analysed the abundance of six LTR (long terminal repeat) retrotransposons on the X and fourth chromosomes of multiple strains of D. melanogaster, which we compare with the predictions of each theory. The data suggest that no one theory can account for the insertion patterns of all six retrotransposons. Comparing our results with earlier work using these transposable element families, we find a significant correlation between studies in the particular model of copy number regulation supported by the proportion of elements on the X for the different transposable element families. This suggests that different retrotransposon families are regulated by different mechanisms.     

Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster

Restriction map variation in 64 X chromosome lines extracted from three different populations of Drosophila melanogaster was investigated with seven six-nucleotide-recognizing restriction enzymes for a 106-kb region encompassing the yellow gene and the achaete-scute complex that is located in the region of reduced crossing over close to the telomere. Nine restriction site polymorphisms (out of 176 sites scored) and 19 length polymorphisms (15 insertions and 4 deletions) were detected. The estimated level of heterozygosity per nucleotide, H = 0.0003, is much lower than that reported for autosomal and sex-linked loci located in regions with normal levels of crossing over. The overall frequency of polymorphic restriction sites is reduced. Six out of nine restriction site polymorphisms are unique and the other three have frequencies less than 0.17. Some large insertions have reached relatively high frequencies, 0.08 to 0.17. Consistent with the theoretically predicted negative relationship between crossing over and the magnitude of linkage disequilibrium, an increase in the relative number of nonrandom associations was observed in the y-ac-sc region.     

Restriction-map variation at the zeste-tko region in natural populations of Drosophila melanogaster

Restriction-map variation in 64 X chromosome lines extracted from three different natural populations of Drosophila melanogaster was investigated with seven six-nucleotide-recognizing enzymes for a 20-kb region including the zeste and tko genes. Ten restriction-site and four length polymorphisms (two insertions and two deletions) were detected. Contrary to the predicted lower level of variation for genes on the X chromosome, the level of variation attributable to nucleotide substitution (estimated heterozygosity/nucleotide = 0.004) was similar to that previously reported for autosomal loci. The amount of insertion/deletion variation in the studied region was within the range observed in autosomal regions and thus not explainable by a simple selection model against the effects of insertional mutations. A general lack of linkage disequilibrium between polymorphic sites was observed.      

Chromosomal patterns of microsatellite variability contrast sharply in African and non-African populations of Drosophila melanogaster.

Levels of neutral variation are influenced by background selection and hitchhiking. The relative contribution of these evolutionary forces to the distribution of neutral variation is still the subject of ongoing debates. Using 133 microsatellites, we determined levels of variability on X chromosomes and autosomes in African and non-African D. melanogaster populations. In the ancestral African populations microsatellite variability was higher on X chromosomes than on autosomes. In non-African populations X-linked polymorphism is significantly more reduced than autosomal variation. In non-African populations we observed a significant positive correlation between X chromosomal polymorphism and recombination rate. These results are consistent with the interpretation that background selection shapes levels of neutral variability in the ancestral populations, while the pattern in derived populations is determined by multiple selective sweeps during the colonization process. Further research, however, is required to investigate the influence of inversion polymorphisms and unequal sex ratios.     

Selective sweep near the In(2L)t inversion breakpoint in an African population of Drosophila melanogaster

Chromosomal inversions largely inhibit recombination and may be associated with selective forces, such as hitch-hiking effects: the effect of positive selection on linked loci. A West African population of Drosophila melanogaster showed a high frequency (0.61) of the In(2L)t inversion. Departure from neutrality statistically associated with the inversion polymorphism was previously recorded at Su(H), a locus distant from the proximal breakpoint of the inversion. These results were consistent with hitch-hiking effects with recombination. The present sequence polymorphism survey involves a 1 kb fragment of the Vha68-1 locus located closer to the proximal breakpoint of the inversion. It shows a significant deficit of polymorphism with respect to divergence when compared with other loci studied in the same population, thus suggesting selective effects. Only 11 polymorphic sites are present in a sample of 20 chromosomes and these sites present a significant excess of rare-frequency variants. The major haplotype shows an unexpectedly high frequency. Our estimate of the background selection effect is not sufficient to account for the observed reduction of polymorphism. Intraspecific variation is structured between inverted and standard chromosomes; there are no shared polymorphisms but also no fixed differences between them. This pattern, together with that found on other loci previously studied near this inversion breakpoint, suggests hitch-hiking effects enhanced by the inversion.     

DNA sequence variation and the recombinational landscape in Drosophila pseudoobscura: a study of the second chromosome.

The relationship between rates of recombination and DNA sequence polymorphism was analyzed for the second chromosome of Drosophila pseudoobscura. We constructed integrated genetic and physical maps of this chromosome using molecular markers at 10 loci spanning most of its physical length. The total length of the map was 128.2 cM, almost twice that of the homologous chromosome arm (3R) in D. melanogaster. There appears to be very little centromeric suppression of recombination, and rates of recombination are quite uniform across most of the chromosome. Levels of sequence variation (theta(W), based on the number of segregating sites) at seven loci (tropomyosin 1, Rhodopsin 3, Rhodopsin 1, bicoid, Xanthine dehydrogenase, Myosin light chain 1, and ribosomal protein 49) varied from 0.0036 to 0.0167. Generally consistent with earlier studies, the average estimate of theta(W) at total sites is 1.5-fold higher than that in D. melanogaster, while average theta(W) at silent sites is almost 3-fold higher. These estimates of variation were analyzed in the context of a background selection model under the same parameters of mutation rate and selection as have been proposed for D. melanogaster. It is likely that a significant fraction of the higher level of sequence variation in D. pseudoobscura can be explained by differences in regional rates of recombination rather than a larger species-level effective population size. However, the distribution of variation among synonymous, nonsynonymous, and noncoding sites appears to be quite different between the species, making direct comparisons of neutral variation, and hence inferences about effective population size, difficult. Tajima's D statistics for 6 out of the 7 loci surveyed are negative, suggesting that D. pseudoobscura may have experienced a rapid population expansion in the recent past or, alternatively, that slightly deleterious mutations constitute an important component of standing variation in this species.     

Nucleotide variation and recombination along the fourth chromosome in Drosophila simulans

The fourth chromosome of Drosophila melanogaster and its sister species are believed to be nonrecombining and have been a model system for testing predictions of the effects of selection on linked, neutral variation. We recently examined nucleotide variation along the chromosome of D. melanogaster and revealed that a low average level of recombination could be associated with considerably high levels of nucleotide variation. In this report, we further investigate the variation along the fourth chromosome of D. simulans. We sequenced 12 gene regions evenly distributed along the fourth chromosome for a worldwide collection of 11 isofemale lines and 5 gene regions in a local population of 10 isofemale lines from South America. In contrast to predictions for regions of very low recombination, these data reveal that the variation levels in many gene regions, including an intron region of the ci gene, vary considerably along the fourth chromosome. Nucleotide diversity ranged from 0.0010 to 0.0074 in 9 gene regions interspersed with several regions of greatly reduced variation. Tests of recombination indicate that the recombination level is not as low as previously thought, likely an order of magnitude higher than that in D. melanogaster. Finally, estimates of the recombination parameters are shown to support a crossover-plus-conversion model.     

The lack of recombination drives the fixation of transposable elements on the fourth chromosome of Drosophila melanogaster

In regions of suppressed recombination, where selection is expected to be less efficient in removing slightly deleterious mutations, transposable element (TE) insertions should be more likely to drift to higher frequencies, and even to reach fixation. In the absence of excision events, once a TE is fixed it cannot be eliminated from the population, and accumulation of elements thus should become an irreversible process. In the long term, this can drive the degeneration of large non-recombining fractions of the genomes. Chromosome 4 of Drosophila melanogaster has very low levels of recombination, if any, and this could be causing its degeneration. Here we report the results of a PCR-based analysis of the population frequencies of TE insertions in a sample from three African natural populations. We investigated 27 insertions from 12 TE families, located in regions of either suppressed or free recombination. Our results suggest that TE insertions tend to be fixed in the non-recombining regions, particularly on the fourth chromosome. We have also found that this involves all types of elements, and that fixed insertions are significantly shorter and more divergent from the canonical sequence than those segregating in the sample (28.1% vs 86.3% of the canonical length, and average nucleotide divergence (D(XY)) = 0.082 vs 0.008, respectively). Finally, DNA-based elements seem to show a greater tendency to reach fixation than retrotransposons. Implications of these findings for the population dynamics of TEs, and the evolutionary forces that shape the patterns of genetic variation in regions of reduced recombination, are discussed.     

Linkage disequilibrium mapping of molecular polymorphisms at the scabrous locus associated with naturally occurring variation in bristle number in Drosophila melanogaster

We evaluated the hypothesis that the Drosophila melanogaster second chromosome gene scabrous (sca), a candidate sensory bristle number quantitative trait locus (QTL), contributes to naturally occurring variation in bristle number. Variation in abdominal and sternopleural bristle number was quantified for wild-derived sca alleles in seven genetic backgrounds: as homozygous second chromosomes (C2) in an isogenic background, homozygous lines in which approximately 20 cM including the sca locus had been introgressed into the isogenic background (sca BC), as C2 and sca BC heterozygotes and hemizygotes against a P element insertional sca allele and a P-induced sca deficiency in the same isogenic background, and as sca BC heterozygotes against the wild-type sca allele of isogenic strain. Molecular restriction map variation was determined for a 45 kb region including the sca locus, and single-stranded conformational polymorphism (SSCP) was examined for the third intron and parts of the third and fourth exons. Associations between each of the 27 molecular polymorphisms and bristle number were evaluated within each genotype and on the first principal component score determined from all seven genotypes, separately for each sex and bristle trait. Permutation tests were used to assess the empirical significance thresholds, accounting for multiple, correlated tests, and correlated markers. Three sites in regulatory regions were associated with female-specific variation in abdominal bristle number, one of which was an SSCP site in the region of the gene associated with regulation of sca in embryonic abdominal segments.     

Reduced sequence variability on the Neo-Y chromosome of Drosophila americana americana.

Sex chromosomes are generally morphologically and functionally distinct, but the evolutionary forces that cause this differentiation are poorly understood. Drosophila americana americana was used in this study to examine one aspect of sex chromosome evolution, the degeneration of nonrecombining Y chromosomes. The primary X chromosome of D. a. americana is fused with a chromosomal element that was ancestrally an autosome, causing this homologous chromosomal pair to segregate with the sex chromosomes. Sequence variation at the Alcohol Dehydrogenase (Adh) gene was used to determine the pattern of nucleotide variation on the neo-sex chromosomes in natural populations. Sequences of Adh were obtained for neo-X and neo-Y chromosomes of D. a. americana, and for Adh of D. a. texana, in which it is autosomal. No significant sequence differentiation is present between the neo-X and neo-Y chromosomes of D. a. americana or the autosomes of D. a. texana. There is a significantly lower level of sequence diversity on the neo-Y chromosome relative to the neo-X in D. a. americana. This reduction in variability on the neo-Y does not appear to have resulted from a selective sweep. Coalescent simulations of the evolutionary transition of an autosome into a Y chromosome indicate there may be a low level of recombination between the neo-X and neo-Y alleles of Adh and that the effective population size of this chromosome may have been reduced below the expected value of 25% of the autosomal effective size, possibly because of the effects of background selection or sexual selection.     

A scan of molecular variation leads to the narrow localization of a selective sweep affecting both Afrotropical and cosmopolitan populations of Drosophila melanogaster

Drosophila melanogaster originated in tropical Africa but has achieved a cosmopolitan distribution in association with human habitation. Cosmopolitan populations of D. melanogaster are known to have reduced genetic variation, particularly on the X chromosome. However, the relative importance of population bottlenecks and selective sweeps in explaining this reduction is uncertain. We surveyed variation at 31 microsatellites across a 330-kb section of the X chromosome located between the white and kirre genes. Two linked clusters of loci were observed with reduced variation and a skew toward rare alleles in both an Ecuador and a Zimbabwe population sample. Examining Zimbabwe DNA sequence polymorphism within one of these regions allowed us to localize a selective sweep to a 361-bp window within the 5' regulatory region of the roughest gene, with one nucleotide substitution representing the best candidate for the target of selection. Estimates of sweep age suggested that this fixation event occurred prior to the expansion of D. melanogaster from sub-Saharan Africa. For both putative sweep regions in our data set, cosmopolitan populations showed wider footprints of selection compared to those in Zimbabwe. This pattern appears consistent with the demographic amplification of preexisting sweep signals due to one or more population bottlenecks.     

Genetic recombination and adaptation to fluctuating environments: selection for geotaxis in Drosophila melanogaster

Heritable variation in fitness is the fuel of adaptive evolution, and sex can generate new adaptive combinations of alleles. If the generation of beneficial combinations drives the evolution of recombination, then the level of recombination should result in changes in the response to selection. Three types of lines of Drosophila melanogaster varying in their level of genetic recombination were selected over 38 generations for geotaxis. The within-chromosome recombination level of these lines was controlled for 60% of the genome: chromosome X and chromosome II. The full recombination lines had normal, unmanipulated levels of recombination on these two chromosomes. Conversely, nonrecombination lines had recombination effectively eliminated within the X and second chromosomes. Finally, partial recombination lines had the effective rate of within-chromosome recombination lowered to 10% of natural levels for these two chromosomes. The rate of response to selection was measured for continuous negative geotaxis and for a fluctuating environment (alternating selection for negative and positive geotaxis). All selected Drosophila lines responded to selection and approximately 36% of the response to selection was because of the X and second chromosomes. However, recombination did not accelerate adaptation during either directional or fluctuating selection for geotaxis.     

The evolutionary history of Drosophila buzzatti. XXVI. Macrogeographic patterns of inversion polymorphism in New World populations

Inversion polymorphisms in the second and fourth chromosomes of the cactophilic Drosophila buzzatti in the native distribution range of the species are described. Over 5,000 flies from 26 localities were scored revealing interesting geographic structuring of arrangement frequencies. Multiple regression and partial correlation approaches showed that the frequencies of second and fourth chromosome arrangements vary clinically along latitudinal and altitudinal gradients and to a lesser extent with longitude. Although many non selective explanations can account for this pattern, its resemblance to the clinal pattern described in recently established Australian populations of Drosophila buzzatii, strongly suggests a selective explanation. Additionally, the correlated variation observed between the frequencies of arrangements 2St on the second chromosome and 4St on the fourth suggests a pattern of interchromosomal association, which, when considering the vast area surveyed, might be explained as the result of epistatic interactions. The analysis of population structure revealed a significant regional pattern, concordant with previously described phytogeographic regions. F-statistics showed that the patterns of variation were different not only between the second and fourth chromosomes, but also between second chromosome arrangements, suggesting that selective differentiation might have contributed to population structure. Since D. buzzatii breeds and feeds on the decaying tissues of diverse cactus species present in different phytogeographic regions, it is difficult to distinguish the underlying causes of the geographic patterns observed. However, inversion heterozygosity is not correlated with the diversity of potential cactus hosts. The evidence presented suggests that differential selection may be the main cause for the population structure. It is also possible to conclude that the inversion polymorphism of D. buzzatti is flexible rather than rigid.     

What is the impact of transposable elements on host genome variability?

The spread of a transposable element family through a wild population may be of astonishing rapidity. At least three families of transposable genetic elements have recently invaded Drosophila melanogaster worldwide, including the P element. The mechanism has been a process of effectively replicative transposition, and, for the P element, has occurred notwithstanding the sterility induced by unrestricted movement. This element's invasion into D. melanogaster has been accompanied by the development of heterogeneity between P sequences, most of which now have internal deletions. Increasing evidence suggests that some deleted elements can repress P transposition, thereby protecting the host from the harmful effects of complete elements. Such repressing elements may rise to high frequencies in populations as a result of selection at the level of the host. We here investigate selective sweeps invoked by the spread of P sequences in D. melanogaster populations. Numerous high-frequency sites have been identified on the X chromosome, which differ in frequency between populations, and which are associated with repression of P-element transposition. Unexpectedly, sequences adjacent to high-frequency P-element sites do not show reduced levels of genetic diversity, and DNA variability is in linkage equilibrium with the presence or absence of a P element at the adjacent selected site. This might be explained by multiple insertions or through a selection for recombination analogous to that seen in 'hitchhiking'.     

Palliating the impact of fixation of a major gene on the genetic variation of artificially selected polygenes

Selective sweeps of variation caused by fixation of major genes may have a dramatic impact on the genetic gain from background polygenic variation, particularly in the genome regions closely linked to the major gene. The response to selection can be restrained because of the reduced selection intensity and the reduced effective population size caused by the increase in frequency of the major gene. In the context of a selected population where fixation of a known major gene is desired, the question arises as to which is the optimal path of increase in frequency of the gene so that the selective sweep of variation resulting from its fixation is minimized. Using basic theoretical arguments we propose a frequency path that maximizes simultaneously the effective population size applicable to the selected background and the selection intensity on the polygenic variation by minimizing the average squared selection intensity on the major gene over generations up to a given fixation time. We also propose the use of mating between carriers and non-carriers of the major gene, in order to promote the effective recombination between the major gene and its linked polygenic background. Using a locus-based computer simulation assuming different degrees of linkage, we show that the path proposed is more effective than a similar path recently published, and that the combination of the selection and mating methods provides an efficient way to palliate the negative effects of a selective sweep.     

S-element insertions are associated with the evolution of the Hsp70 genes in Drosophila melanogaster

The "selfish DNA" theory postulates that transposable elements (TEs) are intragenomic parasites, and that natural selection against deleterious effects associated with their presence is the main force preventing their genomic spread in natural populations. In agreement with this model, TEs in Drosophila melanogaster populations are usually found at low frequencies in most genomic locations. Only a few cases of fixation of TE insertions have been reported, usually in regions of low recombination, where selection is expected to be less effective. Here, we report a population genetics study on the apparent fixation of an S-element in a highly recombining region in two natural populations of D. melanogaster. Three similar fragments of an S-element are inserted into the 5' regions of three members of a heat shock gene family, Hsp70 (Hsp70Aa and Hsp70Ab in polytene chromosome band 87A, and Hsp70Bb in 87C). A PCR-based analysis suggests that the insertions are fixed or at high frequencies in the entire species. A population survey of the levels of nucleotide sequence variation at the insertion site in 87C in two natural populations of D. melanogaster provided evidence for reduced levels of variation in the region, normal levels of recombination, and selection, reflected in a significant departure from neutrality of the variant frequency spectrum. This was particularly strong for the S-element inverted repeats (IRs) and suggests that these are of functional significance for the host.