Population structure among African and derived populations of Drosophila simulans: evidence for ancient subdivision and recent admixture

A Large Panel of Drosophila simulans Reveals an Abundance of Common Variants

The rapidly expanding availability of large NGS data sets provides an opportunity to investigate population genetics at an unprecedented scale. Drosophila simulans is the sister species of the model organism Drosophila melanogaster, and is often presumed to share similar demographic history. However, previous population genetic and ecological work suggests very different signatures of selection and demography. Here, we sequence a new panel of 170 inbred genotypes of a North American population of D. simulans, a valuable complement to the DGRP and other D. melanogaster panels. We find some unexpected signatures of demography, in the form of excess intermediate frequency polymorphisms. Simulations suggest that this is possibly due to a recent population contraction and selection. We examine the outliers in the D. simulans genome determined by a haplotype test to attempt to parse the contribution of demography and selection to the patterns observed in this population. Untangling the relative contribution of demography and selection to genomic patterns of variation is challenging, however, it is clear that although D. melanogaster was thought to share demographic history with D. simulans different forces are at work in shaping genomic variation in this population of D. simulans.

The Drosophila simulans Y chromosome interacts with the autosomes to influence male fitness

The Y chromosome should degenerate because it cannot recombine. However, male-limited transmission increases selection efficiency for male-benefit alleles on the Y, and therefore, Y chromosomes should contribute significantly to variation in male fitness. This means that although the Drosophila Y chromosome is small and gene-poor, Y-linked genes are vital for male fertility in Drosophila melanogaster and the Y chromosome has large male fitness effects. It is unclear whether the same pattern is seen in the closely related Drosophila simulans. We backcrossed Y chromosomes from three geographic locations into five genetic backgrounds and found strong Y and genetic background effects on male fertility. There was a significant Y-background interaction, indicating substantial epistasis between the Y and autosomal genes affecting male fertility. This supports accumulating evidence that interactions between the Y chromosome and the autosomes are key determinants of male fitness.

Drosophila simulans: A Species with Improved Resolution in Evolve and Resequence Studies

The combination of experimental evolution with high-throughput sequencing of pooled individuals—i.e., evolve and resequence (E&R)—is a powerful approach to study adaptation from standing genetic variation under controlled, replicated conditions. Nevertheless, E&R studies in Drosophila melanogaster have frequently resulted in inordinate numbers of candidate SNPs, particularly for complex traits. Here, we contrast the genomic signature of adaptation following ∼60 generations in a novel hot environment for D. melanogaster and D. simulans. For D. simulans, the regions carrying putatively selected loci were far more distinct, and thus harbored fewer false positives, than those in D. melanogaster. We propose that species without segregating inversions and higher recombination rates, such as D. simulans, are better suited for E&R studies that aim to characterize the genetic variants underlying the adaptive response.

Patterns of Genome-Wide Variation in Glossina fuscipes fuscipes Tsetse Flies from Uganda

The tsetse fly Glossina fuscipes fuscipes (Gff) is the insect vector of the two forms of Human African Trypanosomiasis (HAT) that exist in Uganda. Understanding Gff population dynamics, and the underlying genetics of epidemiologically relevant phenotypes is key to reducing disease transmission. Using ddRAD sequence technology, complemented with whole-genome sequencing, we developed a panel of ∼73,000 single-nucleotide polymorphisms (SNPs) distributed across the Gff genome that can be used for population genomics and to perform genome-wide-association studies. We used these markers to estimate genomic patterns of linkage disequilibrium (LD) in Gff, and used the information, in combination with outlier-locus detection tests, to identify candidate regions of the genome under selection. LD in individual populations decays to half of its maximum value (r(2) max/2) between 1359 and 2429 bp. The overall LD estimated for the species reaches r(2) max/2 at 708 bp, an order of magnitude slower than in Drosophila Using 53 infected (Trypanosoma spp.) and uninfected flies from four genetically distinct Ugandan populations adapted to different environmental conditions, we were able to identify SNPs associated with the infection status of the fly and local environmental adaptation. The extent of LD in Gff likely facilitated the detection of loci under selection, despite the small sample size. Furthermore, it is probable that LD in the regions identified is much higher than the average genomic LD due to strong selection. Our results show that even modest sample sizes can reveal significant genetic associations in this species, which has implications for future studies given the difficulties of collecting field specimens with contrasting phenotypes for association analysis.

Genome assembly and annotation of a Drosophila simulans strain from Madagascar

Drosophila simulans is a close relative of the genetic model D. melanogaster. Its worldwide distribution in combination with the absence of segregating chromosomal inversions makes this species an increasingly attractive model to study the molecular signatures of adaptation in natural and experimental populations. In an effort to improve the genomic resources for D. simulans, we assembled and annotated the genome of a strain originating from Madagascar (M252), the ancestral range of D. simulans. The comparison of the M252 genome to other available D. simulans assemblies confirmed its high quality, but also highlighted genomic regions that are difficult to assemble with NGS data. The annotation of M252 provides a clear improvement with alternative splicing for 52% of the multiple-exon genes, UTRs for 70% of the genes, 225 novel genes and 781 pseudogenes being reported. We anticipate that the M252 genome will be a valuable resource for many research questions.

Variation of gene expression associated with colonisation of an anthropized environment: comparison between African and European populations of Drosophila simulans

The comparison of transcriptome profiles among populations is a powerful tool for investigating the role of gene expression change in adaptation to new environments. In this study, we use massively parallel sequencing of 3' cDNAs obtained from large samples of adult males, to compare a population of Drosophila simulans from a natural reserve within its ancestral range (eastern Africa) with a derived population collected in the strongly anthropized Rhône valley (France). The goal was to scan for adaptation linked to the invasion of new environments by the species. Among 15,090 genes retained for the analysis, 794 were found to be differentially expressed between the two populations. We observed an increase in expression of reproduction-related genes in eastern Africa, and an even stronger increase in expression of Cytochrome P450, Glutathione transferase and Glucuronosyl transferase genes in the derived population. These three gene families are involved in detoxification processes, which suggests that pesticides are a major environmental pressure for the species in this area. The survey of the Cyp6g1 upstream region revealed the insertion of a transposable element, Juan, in the regulatory sequence that is almost fixed in the Rhône Valley, but barely present in Mayotte. This shows that Cyp6g1 has undergone parallel evolution in derived populations of D. simulans as previously shown for D. melanogaster. The increasing amount of data produced by comparative population genomics and transcriptomics should permit the identification of additional genes associated with functional divergence among those differentially expressed.

Local dynamics of a fast-evolving sex-ratio system in Drosophila simulans

By distorting Mendelian transmission to their own advantage, X-linked meiotic drive elements can rapidly spread in natural populations, generating a sex-ratio bias. One expected consequence is the triggering of a co-evolutionary arms race between the sex chromosome that carries the distorter and suppressors counteracting its effect. Such an arms race has been theoretically and experimentally established and can have many evolutionary consequences. However, its dynamics in contemporary populations is still poorly documented. Here, we investigate the fate of the young X-linked Paris driver in Drosophila simulans from sub-Saharan Africa to the Middle East. We provide the first example of the early dynamics of distorters and suppressors: we find consistent evidence that the driving chromosomes have been rising in the Middle East during the last decade. In addition, identical haplotypes are at high frequencies around the two co-evolving drive loci in remote populations, implying that the driving X chromosomes share a recent common ancestor and suggesting that East Africa could be the cradle of the Paris driver. The segmental duplication associated with drive presents an unusual structure in West Africa, which could reflect a secondary state of the driver. Together with our previous demonstration of driver decline in the Indian Ocean where suppression is complete, these data provide a unique picture of the complex dynamics of a co-evolutionary arms race currently taking place in natural populations of D. simulans.

Reevaluating the infection status by the Wolbachia endosymbiont in Drosophila Neotropical species from the willistoni subgroup

Infections by the endosymbiotic bacterium Wolbachia developed a rapid global expansion within Old World Drosophila species, ultimately infecting also Neotropical species. In this sense, screenings are necessary to characterize new variants of Wolbachia or new hosts, and also in order to map the dynamics of already known infections. In this paper, we performed a double screening approach that combined Dot-blot and PCR techniques in order to reevaluate the infection status by Wolbachia in species from the willistoni subgroup of Drosophila. Genomic DNA from isofemale lines descendent from females collected in the Amazonian Rainforest (n=91) were submitted to Dot-blot, and were positive for Wolbachia, producing a gradient of hybridization signals, suggesting different infection levels, which was further confirmed through quantitative PCR. Samples with a strong signal in the Dot-blot easily amplified in the wsp-PCR, unlike most of the samples with a medium to weak signal. It was possible to molecularly characterize three Drosophila equinoxialis isofemale lines that were found to be infected in a low density by a wMel-like Wolbachia strain, which was also verified in a laboratory line of Drosophila paulistorum Amazonian. We also found Drosophila tropicalis to be infected with the wAu strain and a Drosophila paulistorum Andean-Brazilian semispecies laboratory line to be infected with a wAu-like Wolbachia. Moreover, we observed that all Drosophila willistoni samples tested with the VNTR-141 marker harbor the same Wolbachia variant, wWil, either in populations from the South or the North of Brazil. Horizontal transfer events involving species of Old World immigrants and Neotropical species of the willistoni subgroup are discussed.

Genome-wide patterns of natural variation reveal strong selective sweeps and ongoing genomic conflict in Drosophila mauritiana

Although it is well understood that selection shapes the polymorphism pattern in Drosophila, signatures of classic selective sweeps are scarce. Here, we focus on Drosophila mauritiana, an island endemic, which is closely related to Drosophila melanogaster. Based on a new, annotated genome sequence, we characterized the genome-wide polymorphism by sequencing pooled individuals (Pool-seq). We show that the interplay between selection and recombination results in a genome-wide polymorphism pattern characteristic for D. mauritiana. Two large genomic regions (>500 kb) showed the signature of almost complete selective sweeps. We propose that the absence of population structure and limited geographic distribution could explain why such pronounced sweep patterns are restricted to D. mauritiana. Further evidence for strong adaptive evolution was detected for several nucleoporin genes, some of which were not previously identified as genes involved in genomic conflict. Since this adaptive evolution is continuing after the split of D. mauritiana and Drosophila simulans, we conclude that genomic conflict is not restricted to short episodes, but rather an ongoing process in Drosophila.

Taxonomic and evolutionary analysis of Zaprionus indianus and its colonization of Palearctic and Neotropical regions

Zaprionus indianus is a dipteran (Drosophilidae) with a wide distribution throughout the tropics and temperate Palearctic and Nearctic regions. There have been proposals to reclassify the genus Zaprionus as a subgenus or group of the genus Drosophila because various molecular markers have indicated a close relationship between Zaprionus species and the immigrans-Hirtodrosophila radiation within Drosophila. These markers, together with alloenzymes and quantitative traits, have been used to describe the probable scenario for the expansion of Zaprionus indianus from its center of dispersal (Africa) to regions of Asia (ancient dispersal) and the Americas (recent dispersal). The introduction of Z. indianus into Brazil was first reported in 1999 and the current consensus is that the introduced flies came from high-latitude African populations through the importation of fruit. Once in Brazil, Z. indianus spread rapidly throughout the Southeast and then to the rest of the country, in association with highway-based fruit commerce. These and other aspects of the evolutionary biology of Z. indianus are addressed in this review, including a description of a probable route for this species’ dispersal during its recent expansion.

Genetic diversity, population structure and Wolbachia infection status in a worldwide sample of Drosophila melanogaster and D. simulans populations

Drosophila melanogaster and its close relatives have been extremely important model species in the development of population genetic models that serve to explain patterns of diversity in natural populations, a major goal of evolutionary biology. A detailed picture of the evolutionary history of these species is beginning to emerge, as the relative importance of forces including demographic changes and natural selection is established. A continuing aim is to characterise levels of genetic diversity in a large number of populations of these species, covering a wide geographic area. We have used collections from five previously un-sampled wild populations of D. melanogaster and two of D. simulans, across three continents. We estimated levels of genetic diversity within, and divergence between, these populations, and looked for evidence of genetic structure both between ancestral and derived populations, and amongst derived populations. We also investigated the prevalence of infection with the bacterial endosymbiont Wolbachia. We found that D. melanogaster populations from Sub-Saharan Africa are the most diverse, and that divergence is highest between these and non-Sub-Saharan populations. There is strong evidence for structuring of populations between Sub-Saharan Africa and the rest of the world, and some evidence for weak structure amongst derived populations. Populations from Sub-Saharan Africa also differ in the prevalence of Wolbachia infection, with very low levels of infection compared to populations from the rest of the world.

Population transcriptomics: insights from Drosophila simulans, Drosophila sechellia and their hybrids

Sequence differentiation has been widely studied between populations and species, whereas interest in expression divergence is relatively recent. Using microarrays, we compared four geographically distinct populations of Drosophila simulans and a population of Drosophila sechellia, and interspecific hybrids. We observed few differences between populations, suggesting a slight population structure in D. simulans. This structure was observed in direct population comparisons, as well as in interspecific comparisons (hybrids vs. parents, D. sechellia vs. D. simulans). Expression variance is higher in the French and Zimbabwean populations than in the populations from the ancestral range of D. simulans (Kenya and Seychelles). This suggests a large scale phenomenon of decanalization following the invasion of a new environment. Comparing D. simulans and D. sechellia, we revealed 304 consistently differentially expressed genes, with striking overrepresentation of genes of the cytochrome P450 family, which could be related to their role in detoxification as well as in hormone regulation. We also revealed differences in genes involved in Juvenile hormone and Dopamine differentiation. We finally observed very few differentially expressed genes between hybrids and parental populations, with an overrepresentation of X-linked genes.

Recurrent selection on the Winters sex-ratio genes in Drosophila simulans

Selfish genes, such as meiotic drive elements, propagate themselves through a population without increasing the fitness of host organisms. X-linked (or Y-linked) meiotic drive elements reduce the transmission of the Y (X) chromosome and skew progeny and population sex ratios, leading to intense conflict among genomic compartments. Drosophila simulans is unusual in having a least three distinct systems of X chromosome meiotic drive. Here, we characterize naturally occurring genetic variation at the Winters sex-ratio driver (Distorter on the X or Dox), its progenitor gene (Mother of Dox or MDox), and its suppressor gene (Not Much Yang or Nmy), which have been previously mapped and characterized. We survey three North American populations as well as 13 globally distributed strains and present molecular polymorphism data at the three loci. We find that all three genes show signatures of selection in North America, judging from levels of polymorphism and skews in the site-frequency spectrum. These signatures likely result from the biased transmission of the driver and selection on the suppressor for the maintenance of equal sex ratios. Coalescent modeling indicates that the timing of selection is more recent than the age of the alleles, suggesting that the driver and suppressor are coevolving under an evolutionary "arms race." None of the Winters sex-ratio genes are fixed in D. simulans, and at all loci we find ancestral alleles, which lack the gene insertions and exhibit high levels of nucleotide polymorphism compared to the derived alleles. In addition, we find several "null" alleles that have mutations on the derived Dox background, which result in loss of drive function. We discuss the possible causes of the maintenance of presence-absence polymorphism in the Winters sex-ratio genes.

A recent adaptive transposable element insertion near highly conserved developmental loci in Drosophila melanogaster

A recent genomewide screen identified 13 transposable elements that are likely to have been adaptive during or after the spread of Drosophila melanogaster out of Africa. One of these insertions, Bari-Juvenile hormone epoxy hydrolase (Bari-Jheh), was associated with the selective sweep of its flanking neutral variation and with reduction of expression of one of its neighboring genes: Jheh3. Here, we provide further evidence that Bari-Jheh insertion is adaptive. We delimit the extent of the selective sweep and show that Bari-Jheh is the only mutation linked to the sweep. Bari-Jheh also lowers the expression of its other flanking gene, Jheh2. Subtle consequences of Bari-Jheh insertion on life-history traits are consistent with the effects of reduced expression of the Jheh genes. Finally, we analyze molecular evolution of Jheh genes in both the long- and the short-term and conclude that Bari-Jheh appears to be a very rare adaptive event in the history of these genes. We discuss the implications of these findings for the detection and understanding of adaptation.

An unusual haplotype structure on human chromosome 8p23 derived from the inversion polymorphism

Chromosomal inversion is an important type of genomic variations involved in both evolution and disease pathogenesis. Here, we describe the refined genetic structure of a 3.8-Mb inversion polymorphism at chromosome 8p23. Using HapMap data of 1,073 SNPs generated from 209 unrelated samples from CEPH-Utah residents with ancestry from northern and western Europe (CEU); Yoruba in Ibadan, Nigeria (YRI); and Asian (ASN) samples, which were comprised of Han Chinese from Beijing, China (CHB) and Japanese from Tokyo, Japan (JPT)-we successfully deduced the inversion orientations of all their 418 haplotypes. In particular, distinct haplotype subgroups were identified based on principal component analysis (PCA). Such genetic substructures were consistent with clustering patterns based on neighbor-joining tree reconstruction, which revealed a total of four haplotype clades across all samples. Metaphase fluorescence in situ hybridization (FISH) in a subset of 10 HapMap samples verified their inversion orientations predicted by PCA or phylogenetic tree reconstruction. Positioning of the outgroup haplotype within one of YRI clades suggested that Human NCBI Build 36-inverted order is most likely the ancestral orientation. Furthermore, the population differentiation test and the relative extended haplotype homozygosity (REHH) analysis in this region discovered multiple selection signals, also in a population-specific manner. A positive selection signal was detected at XKR6 in the ASN population. These results revealed the correlation of inversion polymorphisms to population-specific genetic structures, and various selection patterns as possible mechanisms for the maintenance of a large chromosomal rearrangement at 8p23 region during evolution. In addition, our study also showed that haplotype-based clustering methods, such as PCA, can be applied in scanning for cryptic inversion polymorphisms at a genome-wide scale.

Positive and negative selection on noncoding DNA in Drosophila simulans

There is now a wealth of evidence that some of the most important regions of the genome are found outside those that encode proteins, and noncoding regions of the genome have been shown to be subject to substantial levels of selective constraint, particularly in Drosophila. Recent work has suggested that these regions may also have been subject to the action of positive selection, with large fractions of noncoding divergence having been driven to fixation by adaptive evolution. However, this work has focused on Drosophila melanogaster, which is thought to have experienced a reduction in effective population size (N(e)), and thus a reduction in the efficacy of selection, compared with its closest relative Drosophila simulans. Here, we examine patterns of evolution at several classes of noncoding DNA in D. simulans and find that all noncoding DNA is subject to the action of negative selection, indicated by reduced levels of polymorphism and divergence and a skew in the frequency spectrum toward rare variants. We find that the signature of negative selection on noncoding DNA and nonsynonymous sites is obscured to some extent by purifying selection acting on preferred to unpreferred synonymous codon mutations. We investigate the extent to which divergence in noncoding DNA is inferred to be the product of positive selection and to what extent these inferences depend on selection on synonymous sites and demography. Based on patterns of polymorphism and divergence for different classes of synonymous substitution, we find the divergence excess inferred in noncoding DNA and nonsynonymous sites in the D. simulans lineage difficult to reconcile with demographic explanations.

The mariner transposable element in natural populations of Drosophila simulans

In cosmopolitan species, geographical variations in copy number and/or level of transposition activity have been observed for several transposable elements (TEs). Environment, history and population structure can contribute to such variation in ways that are difficult to tease apart. For the mariner element, previous studies of the geographic variation of its somatic activity in natural populations of Drosophila simulans have shown contradictory results (latitudinal clines of divergent orientations or no apparent structure). To try and resolve these inconsistencies, we gathered all available data on the mariner somatic activity of worldwide natural populations. This includes previously published results by different groups and also new data. The correlations between the level of activity and several geoclimatic factors were tested. Although no general effect of temperature was found, a relationship with the invasion history was detected. It was also shown that recent invasive populations have a higher level of activity than the putative ancestral ones. Our results strongly suggest that variability of the mariner somatic activity among natural populations of D. simulans is mainly due to populational and historical factors probably related to the recent world colonization of this species. Indeed, this activity is correlated to the main route out of Africa (the Nile route) and the recent colonization of continents such as Australia and South America.

African Drosophila melanogaster and D. simulans populations have similar levels of sequence variability, suggesting comparable effective population sizes

Drosophila melanogaster and D. simulans are two closely related species with a similar distribution range. Many studies suggested that D. melanogaster has a smaller effective population size than D. simulans. As most evidence was derived from non-African populations, we readdressed this question by sequencing 10 X-linked loci in five African D. simulans and six African D. melanogaster populations. Contrary to previous results, we found no evidence for higher variability, and thus larger effective population size, in D. simulans. Our observation of similar levels of variability of both species will have important implications for the interpretation of patterns of molecular evolution.

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.

Population genomics, the study of genome-wide patterns of sequence variation within and between closely related species, can provide a comprehensive view of the relative importance of mutation, recombination, natural selection, and genetic drift in evolution. It can also provide fundamental insights into the biological attributes of organisms that are specifically shaped by adaptive evolution. One approach for generating population genomic datasets is to align DNA sequences from whole-genome shotgun projects to a standard reference sequence. We used this approach to carry out whole-genome analysis of polymorphism and divergence in Drosophila simulans, a close relative of the model system, D. melanogaster. We find that polymorphism and divergence fluctuate on a large scale across the genome and that these fluctuations are probably explained by natural selection rather than by variation in mutation rates. Our analysis suggests that adaptive protein evolution is common and is often related to biological processes that may be associated with gene expression, chromosome biology, and reproduction. The approaches presented here will have broad applicability to future analysis of population genomic variation in other systems, including humans.

Low-coverage genome sequences from multiple Drosophila simulans strains provide the first comprehensive view of polymorphism and divergence in the fruit fly.

Unusual pattern of nucleotide sequence variation at the OS-E and OS-F genomic regions of Drosophila simulans

Nucleotide variation at the genomic region encompassing the odorant-binding protein genes OS-E and OS-F (OS region) was surveyed in two populations of Drosophila simulans, one from Europe and the other from Africa. We found that the European population shows an atypical and large haplotype structure, which extends throughout the approximately 5-kb surveyed genomic region. This structure is depicted by two major haplotype groups segregating at intermediate frequency in the sample, one haplogroup with nearly no variation, and the other at levels more typical for this species. This pattern of variation was incompatible with neutral predictions for a population at a stationary equilibrium. Nevertheless, neutrality tests contrasting polymorphism and divergence data fail to detect any departure from the standard neutral model in this species, whereas they confirm the non-neutral behavior previously observed at the OS-E gene in D. melanogaster. Although positive Darwinian selection may have been responsible for the observed unusual nucleotide variation structure, coalescent simulation results do not allow rejecting the hypothesis that the pattern was generated by a recent bottleneck in the history of European populations of D. simulans.

Microsatellite variation and differentiation in African and non-African populations of Drosophila simulans

Drosophila simulans originated in sub-Saharan Africa or Madagascar and colonized the rest of the world after the last glaciation about 10 000 years ago. Consistent with this demographic history, sub-Saharan African populations have been shown to harbour higher levels of microsatellite and sequence variation than cosmopolitan populations. Nevertheless, only limited information is available on the population structure of D. simulans. Here, we analysed X-linked and autosomal microsatellite loci in four sub-Saharan African, one North African, one Israeli, and two European D. simulans populations. Bayesian clustering algorithms combined the North African, Israeli, and European populations into a single cosmopolitan group. The four sub-Saharan populations were split into two separate groups. Pairwise F(ST) analysis, however, indicated significant population differentiation between all eight populations surveyed. A significant signal for population reduction in cosmopolitan populations was found only for X-linked loci.

Recurrent positive selection at bgcn, a key determinant of germ line differentiation, does not appear to be driven by simple coevolution with its partner protein bam

Surveys of nucleotide sequence polymorphism in Drosophila melanogaster and Drosophila simulans were performed at 2 interacting loci crucial for gametogenesis: bag-of-marbles (bam) and benign gonial cell neoplasm (bgcn). At the polymorphism level, both loci appear to be evolving under the expectations of the neutral theory. However, ratios of polymorphism and divergence for synonymous and nonsynonymous mutations depart significantly from neutral expectations for both loci consistent with a previous observation of positive selection at bam. The deviations suggest either an excess of synonymous polymorphisms or an excess of nonsynonymous fixations at both loci. Synonymous evolution appears to conform to neutrality at bam. At bgcn, there is evidence of positive selection affecting preferred synonymous mutations along the D. simulans lineage. However, there is also a significantly higher rate of nonsynonymous fixations at bgcn within D. simulans. Thus, the deviation from neutrality detected by the McDonald-Kreitman test at these 2 loci is likely due to the selective acceleration of nonsynonymous fixations. Differences in the pattern of amino acid fixations between these 2 interacting proteins suggest that the detected positive selection is not due to a simple model of coevolution.

Positive selection versus demography: evolutionary inferences based on an unusual haplotype structure in Drosophila simulans

Coalescent simulations were used to investigate the possible role of population subdivision and history in shaping nucleotide variation in a recombining 88-kb genomic fragment of Drosophila simulans displaying an unusual large-scale haplotype structure. The multilocus analysis, based on summary statistics using specific demographic null models under recombination, indicates that the observed levels of linkage disequilibrium differed significantly from the values expected under different bottleneck and population admixture scenarios. These results indicate that demography alone may not account for the observed pattern of variation and support the previous claim that the data are better described by a model in which an adaptive mutation has not yet gone to fixation.

Estimating the genomewide rate of adaptive protein evolution in Drosophila

When polymorphism and divergence data are available for multiple loci, extended forms of the McDonald-Kreitman test can be used to estimate the average proportion of the amino acid divergence due to adaptive evolution--a statistic denoted alpha. But such tests are subject to many biases. Most serious is the possibility that high estimates of alpha reflect demographic changes rather than adaptive substitution. Testing for between-locus variation in alpha is one possible way of distinguishing between demography and selection. However, such tests have yielded contradictory results, and their efficacy is unclear. Estimates of alpha from the same model organisms have also varied widely. This study clarifies the reasons for these discrepancies, identifying several method-specific biases in widely used estimators and assessing the power of the methods. As part of this process, a new maximum-likelihood estimator is introduced. This estimator is applied to a newly compiled data set of 115 genes from Drosophila simulans, each with each orthologs from D. melanogaster and D. yakuba. In this way, it is estimated that alpha approximately 0.4+/-0.1, a value that does not vary substantially between different loci or over different periods of divergence. The implications of these results are discussed.

Multiple introductions promote range expansion of the mollusc Cyclope neritea (Nassariidae) in France: evidence from mitochondrial sequence data

Since the 1970s, the nassariid gastropod Cyclope neritea has been extending its range north along the French Atlantic coasts from the Iberian Peninsula. This may be due to natural spread because of the recent warming of the northeastern Atlantic. However, human-mediated introductions related to shellfish culture may also be a probable explanation for this sudden range expansion. To examine these two hypotheses, we carried out a comprehensive study based on mitochondrial gene sequences (cytochrome oxidase I) of the five recently colonized French bays as well as 14 populations located in the recognized native range of the species. From a total of 594 individuals, we observed 29 haplotypes to split into three divergent clades. In the native range, we observed a low molecular diversity, strong genetic structure and agreement between geography and gene genealogies. Along the French coasts, we observed the opposite: high genetic diversity and low genetic structure. Our results show that recurrent human-mediated introductions from several geographical areas in the native range may be a source for the French Atlantic populations. However, despite the low dispersal ability of C. neritea, the isolation-by-distance pattern in France suggested that this gastropod may have been present (although unnoticed) on the French Atlantic coasts before the 1970s. As C. neritea shows characteristics of a cryptogenic species, the classification of Atlantic populations as either native or introduced is not straightforward. Cryptogenic species should be studied further to determine the status of new populations close to their recognized native range.

Contrasted polymorphism patterns in a large sample of populations from the evolutionary genetics model Drosophila simulans

African populations of Drosophila simulans are thought to be ancestral in this model species and are increasingly used for testing general hypotheses in evolutionary genetics. It is often assumed that African populations are more likely to be at a neutral mutation drift equilibrium than other populations. Here we examine population structuring and the demographic profile in nine populations of D. simulans. We surveyed sequence variation in four X-linked genes (runt, sevenless, Sex-lethal, and vermilion) that have been used in a parallel study in the closely related species D. melanogaster. We found that an eastern group of populations from continental Africa and Indian Ocean islands (Kenya, Tanzania, Madagascar, and Mayotte Island) is widespread, shows little differentiation, and has probably undergone demographic expansion. The other two African populations surveyed (Cameroon and Zimbabwe) show no evidence of population expansion and are markedly differentiated from each other as well as from the populations from the eastern group. Two other populations, Europe and Antilles, are probably recent invaders to these areas. The Antilles population is probably derived from Europe through a substantial bottleneck. The history of these populations should be taken into account when drawing general conclusions from variation patterns.

The evolution of antifungal peptides in Drosophila

An essential component of the immune system of animals is the production of antimicrobial peptides (AMPs). In vertebrates and termites the protein sequence of some AMPs evolves rapidly under positive selection, suggesting that they may be coevolving with pathogens. However, antibacterial peptides in Drosophila tend to be highly conserved. We have inferred the selection pressures acting on Drosophila antifungal peptides (drosomycins) from both the divergence of drosomycin genes within and between five species of Drosophila and polymorphism data from Drosophila simulans and D. melanogaster. In common with Drosophila antibacterial peptides, there is no evidence of adaptive protein evolution in any of the drosomycin genes, suggesting that they do not coevolve with pathogens. It is possible that this reflects a lack of specific fungal and bacterial parasites in Drosophila populations. The polymorphism data from both species differed from neutrality at one locus, but this was not associated with changes in the protein sequence. The synonymous site diversity was greater in D. simulans than in D. melanogaster, but the diversity both upstream of the genes and at nonsynonymous sites was similar. This can be explained if both upstream and nonsynonymous mutations are slightly deleterious and are removed more effectively from D. simulans due to its larger effective population size.

Historical biogeography of Drosophila simulans based on Y-chromosomal sequences

Y-chromosomal sequences have been used for phylogeographic studies in humans and other mammals, but so far have been ignored as a source of historical information in Drosophila and other insects with X/Y sex determination. Here, we present the first phylogeographic study of Drosophila simulans based on the Y chromosome. Geographic distribution of Y-chromosomal haplotypes suggests a high degree of population subdivision within Africa, as well as between the African and cosmopolitan groups of populations. Consistent with earlier studies based on autosomal and X-linked loci, our results suggest that D. simulans originated in Madagascar or East Africa, and that the South and West African populations of this species are derived.

Inferring the population structure and demography of Drosophila ananassae from multilocus data

Inferring the origin, population structure, and demographic history of a species is a major objective of population genetics. Although many organisms have been analyzed, the genetic structures of subdivided populations are not well understood. Here we analyze Drosophila ananassae, a highly substructured, cosmopolitan, and human-commensal species distributed in the tropical, subtropical, and mildly temperate regions of the world. We adopt a multilocus approach (with 10 neutral loci) using 16 population samples covering almost the entire species range (Asia, Australia, and America). Analyzed with our recently developed Bayesian method, 5 populations in Southeast Asia are found to be central, while the other 11 are peripheral. These 5 central populations were sampled from localities that belonged to a single landmass ("Sundaland") during the late Pleistocene ( approximately 18,000 years ago), when sea level was approximately 120 m below the present level. The inferred migration routes of D. ananassae out of Sundaland seem to parallel those of humans in this region. Strong evidence for a population size expansion is seen particularly in the ancestral populations.

DNA sequence polymorphism and divergence at the erect wing and suppressor of sable loci of Drosophila melanogaster and D. simulans

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (approximately 50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high FST values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.

Linkage disequilibrium and recent selection at three immunity receptor loci in Drosophila simulans

Immune system genes in a California population sample of Drosophila simulans were shown to bear several hallmarks of the effects of past directional selection. One potential effect of directional selection is an increase in linkage disequilibrium among the polymorphic sites that are linked to the site under selection. In this study, we focus on three D. simulans immunity loci, Hmu, Sr-CI/Sr-CIII, and Tehao, for which the polymorphic sites are in nearly perfect linkage disequilibrium, an unusual finding even with respect to other immunity genes sampled from the same lines. The most likely explanation for this finding is that, at each locus, two divergent alleles have been selected to intermediate frequencies in the recent past. The extent to which the linkage disequilibrium extends to the flanks of each of the immunity genes is minimal, suggesting that the favored mutations actually occurred within the immunity genes themselves. Furthermore, the excess linkage disequilibrium found in the California population is not found in an African D. simulans population sample and may be a result of novel pathogen-mediated selection pressures encountered during establishment of non-African populations.

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.

Is genome size influenced by colonization of new environments in dipteran species?

Genome size differences are usually attributed to the amplification and deletion of various repeated DNA sequences, including transposable elements (TEs). Because environmental changes may promote modifications in the amount of these repeated sequences, it has been postulated that when a species colonizes new environments this could be followed by an increase in its genome size. We tested this hypothesis by estimating the genome size of geographically distinct populations of Drosophila ananassae, Drosophila malerkotliana, Drosophila melanogaster, Drosophila simulans, Drosophila subobscura, and Zaprionus indianus, all of which have known colonization capacities. There was no strong statistical differences between continents for most species. However, we found that populations of D. melanogaster from east Africa have smaller genomes than more recent populations. For species in which colonization is a recent event, the differences between genome sizes do not thus seem to be related to colonization history. These findings suggest either that genome size is seldom modified in a significant way during colonization or that it takes time for genome size of invading species to change significantly.

DNA variability and divergence at the notch locus in Drosophila melanogaster and D. simulans: a case of accelerated synonymous site divergence

DNA diversity in two segments of the Notch locus was surveyed in four populations of Drosophila melanogaster and two of D. simulans. In both species we observed evidence of non-steady-state evolution. In D. simulans we observed a significant excess of intermediate frequency variants in a non-African population. In D. melanogaster we observed a disparity between levels of sequence polymorphism and divergence between one of the Notch regions sequenced and other neutral X chromosome loci. The striking feature of the data is the high level of synonymous site divergence at Notch, which is the highest reported to date. To more thoroughly investigate the pattern of synonymous site evolution between these species, we developed a method for calibrating preferred, unpreferred, and equal synonymous substitutions by the effective (potential) number of such changes. In D. simulans, we find that preferred changes per "site" are evolving significantly faster than unpreferred changes at Notch. In contrast we observe a significantly faster per site substitution rate of unpreferred changes in D. melanogaster at this locus. These results suggest that positive selection, and not simply relaxation of constraint on codon bias, has contributed to the higher levels of unpreferred divergence along the D. melanogaster lineage at Notch.

Signature of selective sweep associated with the evolution of sex-ratio drive in Drosophila simulans

In several Drosophila species, the XY Mendelian ratio is disturbed by X-linked segregation distorters (sex-ratio drive). We used a collection of recombinants between a nondistorting chromosome and a distorting X chromosome originating from the Seychelles to map a candidate sex-ratio region in Drosophila simulans using molecular biallelic markers. Our data were compatible with the presence of a sex-ratio locus in the 7F cytological region. Using sequence polymorphism at the Nrg locus, we showed that sex-ratio has induced a strong selective sweep in populations from Madagascar and Réunion, where distorting chromosomes are close to a 50% frequency. The complete association between the marker and the sex-ratio phenotype and the near absence of mutations and recombination in the studied fragment after the sweep event indicate that this event is recent. Examples of selective sweeps are increasingly reported in a number of genomes. This case identifies the causal selective force. It illustrates that all selective sweeps are not necessarily indicative of an increase in the average fitness of populations.

Effect of breeding structure on population genetic parameters in Drosophila

The breeding structure of populations has been neglected in studies of Drosophila, even though Wright and Dobzhansky's pioneering work on the genetics of natural populations was an attempt to tackle what they regarded as an essential factor in evolution. We compared the breeding structure of sympatric populations of D. melanogaster and D. simulans, two sibling species that are widely used in evolutionary studies. We recorded changes in population density and microsatellite variation patterns for 3 years in a temperate environment of southwestern France. Results were distinctively different in the two species. Maximum population levels in summer and in autumn were similar and fluctuated greatly over years, each species being in turn the most abundant. However, genetic data showed that D. melanogaster made up a continuous breeding population in time and space of practically infinite effective size. D. simulans was fragmented into isolates with a local effective size of between 50 and 350 individuals. A consequence of this was that, while a local sample provided a reliable estimate of regional genetic variability in D. melanogaster, a sample from the same area provided an underestimate of this parameter in D. simulans. In practical terms, this means that variations in breeding structure should be accounted for in sampling schemes and in designing evolutionary genetic models. More generally, this suggests the existence of differential reactions to local environments that might contribute to several genomic differences observed between these species.

Linking phylogenetics with population genetics to reconstruct the geographic origin of a species

Reconstructing ancestral geographic origins is critical for understanding the long-term evolution of a species. Bayesian methods have been proposed to test biogeographic hypotheses while accommodating uncertainty in phylogenetic reconstruction. However, the problem that certain taxa may have a disproportionate influence on conclusions has not been addressed. Here, we infer the geographic origin of Drosophila simulans using 2,014 bp of the period locus from 63 lines collected from 18 countries. We also analyze two previously published datasets, alcohol dehydrogenase related and NADH:ubiquinone reductase 75 kDa subunit precursor. Phylogenetic inferences of all three loci support Madagascar as the geographic origin of D. simulans. Our phylogenetic conclusions are robust to taxon resampling and to the potentially confounding effects of recombination. To test our phylogenetically derived hypothesis we develop a randomization test of the population genetics prediction that sequences from the geographic origin should contain more genetic polymorphism than those from derived populations. We find that the Madagascar population has elevated genetic polymorphism relative to non-Madagascar sequences. These data are corroborated by mitochondrial DNA sequence data.

Patterns of Microsatellite Variability Among X Chromosomes and Autosomes Indicate a High Frequency of Beneficial Mutations in Non-African D. simulans

We analyzed microsatellite variability at 42 X-linked and 39 autosomal loci from African and European populations of Drosophila simulans. The African D. simulans harbored significantly more microsatellite variability than the European flies. In the European population, X-linked polymorphism was more reduced than autosomal variation, whereas there was no significant difference between chromosomes in the African population. Previous studies also observed a similar pattern but failed to distinguish between a demographic event and a selection scenario. We performed extensive computer simulations using a wide range of demographic scenarios to distinguish between the two hypotheses. Approximate summary likelihood estimates differed dramatically among X chromosomes and autosomes. Furthermore, our experimental data showed a surplus of X-linked microsatellites with a significantly reduced variability in non-African D. simulans. We conclude that our data are not compatible with a neutral scenario. Thus, the reduced variability at X-linked loci is most likely caused by selective sweeps associated with the out-of-Africa habitat expansion of D. simulans.

Strong selective sweep associated with a transposon insertion in Drosophila simulans

We know little about several important properties of beneficial mutations, including their mutational origin, their phenotypic effects (e.g., protein structure changes vs. regulatory changes), and the frequency and rapidity with which they become fixed in a population. One signature of the spread of beneficial mutations is the reduction of heterozygosity at linked sites. Here, we present population genetic data from several loci across chromosome arm 2R in Drosophila simulans. A 100-kb segment from a freely recombining region of this chromosome shows extremely reduced heterozygosity in a California population sample, yet typical levels of divergence between species, suggesting that at least one episode of strong directional selection has occurred in the region. The 5' flanking sequence of one gene in this region, Cyp6g1 (a cytochrome P450), is nearly fixed for a Doc transposable element insertion. Presence of the insertion is correlated with increased transcript abundance of Cyp6g1, a phenotype previously shown to be associated with insecticide resistance in Drosophila melanogaster. Surveys of nucleotide variation in the same genomic region in an African D. simulans population revealed no evidence for a high-frequency Doc element and no evidence for reduced polymorphism. These data are consistent with the notion that the Doc element is a geographically restricted beneficial mutation. Data from D. simulans Cyp6g1 are paralleled in many respects by data from its sister species D. melanogaster.

Influence of Two Wolbachia Strains on Population Structure of East AfricanDrosophila simulans

Drosophila simulans is hypothesized to have originated in continental East Africa or Madagascar. In this study, we investigated evolutionary forces operating on mitochondrial DNA (mtDNA) in populations of D. simulans from Zimbabwe, Malawi, Tanzania, and Kenya. Variation in mtDNA may be affected by positive selection, background selection, demographic history, and/or any maternally inherited factor such as the bacterial symbiont Wolbachia. In East Africa, the wRi and wMa Wolbachia strains associate with the siII or siIII mitochondrial haplogroups, respectively. To ask how polymorphism relates to Wolbachia infection status, we sequenced 1776 bp of mitochondrial DNA and 1029 bp of the X-linked per locus from 79 lines. The two southern populations were infected with wRi and exhibited significantly reduced mtDNA variation, while Wolbachia-uninfected siII flies from Tanzania and Kenya showed high levels of mtDNA polymorphism. These are the first known populations of D. simulans that do not exhibit reduced mtDNA variation. We observed no mitochondrial variation in the siIII haplogroup regardless of Wolbachia infection status, suggesting positive or background selection. These populations offer a unique opportunity to monitor evolutionary dynamics in ancestral populations that harbor multiple strains of Wolbachia.

Linkage Disequilibrium Patterns Across a Recombination Gradient in African Drosophila melanogaster

Previous multilocus surveys of nucleotide polymorphism have documented a genome-wide excess of intralocus linkage disequilibrium (LD) in Drosophila melanogaster and D. simulans relative to expectations based on estimated mutation and recombination rates and observed levels of diversity. These studies examined patterns of variation from predominantly non-African populations that are thought to have recently expanded their ranges from central Africa. Here, we analyze polymorphism data from a Zimbabwean population of D. melanogaster, which is likely to be closer to the standard population model assumptions of a large population with constant size. Unlike previous studies, we find that levels of LD are roughly compatible with expectations based on estimated rates of crossing over. Further, a detailed examination of genes in different recombination environments suggests that markers near the telomere of the X chromosome show considerably less linkage disequilibrium than predicted by rates of crossing over, suggesting appreciable levels of exchange due to gene conversion. Assuming that these populations are near mutation-drift equilibrium, our results are most consistent with a model that posits heterogeneity in levels of exchange due to gene conversion across the X chromosome, with gene conversion being a minor determinant of LD levels in regions of high crossing over. Alternatively, if levels of exchange due to gene conversion are not negligible in regions of high crossing over, our results suggest a marked departure from mutation-drift equilibrium (i.e., toward an excess of LD) in this Zimbabwean population. Our results also have implications for the dynamics of weakly selected mutations in regions of reduced crossing over.

Population Subdivision and Molecular Sequence Variation: Theory and Analysis of Drosophila ananassae Data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter 0398; to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of 0398;, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.

Characterization of non-cytoplasmic incompatibility inducing Wolbachia in two continental African populations of Drosophila simulans

Wolbachia is an endocellular bacterium infecting arthropods and nematodes. In arthropods, it invades host populations through various mechanisms, affecting host reproduction, the most common of which being cytoplasmic incompatibility (CI). CI is an embryonic mortality occurring when infected males mate with uninfected females or females infected by a different Wolbachia strain. This phenomenon is observed in Drosophila simulans, an intensively studied Wolbachia host, harbouring at least five distinct bacterial strains. In this study, we investigate various aspects of the Wolbachia infections occurring in two continental African populations of D. simulans: CI phenotype, phylogenetic position based on the wsp gene and associated mitochondrial haplotype. From the East African population (Tanzania), we show that (i) the siIII mitochondrial haplotype occurs in continental populations, which was unexpected based on the current views of D. simulans biogeography, (ii) the wKi strain (that rescues from CI while being unable to induce it) is very closely related to the CI-inducing strain wNo, (iii) wKi and wNo might not derive from a unique infection event, and (iv) wKi is likely to represent the same entity as the previously described wMa variant. In the West African population (Cameroon), the Wolbachia infection was found identical to the previously described wAu, which does not induce CI. This finding supports the view that wAu might be an ancient infection in D. simulans.

Worldwide distribution of transposable element copy number in natural populations of Drosophila simulans

Transposable elements (TEs), which promote various kinds of mutations, constitute a large fraction of the genome. How they invade natural populations and species is therefore of fundamental importance for understanding the dynamics of genetic diversity and genome composition. On the basis of 85 samples of natural populations of Drosophila simulans, we report the distributions of the genome insertion site numbers of nine TEs that were chosen because they have a low average number of sites. Most populations were found to have 0-3 insertion sites, but some of them had a significantly higher number of sites for a given TE. The populations located in regions outside Africa had the highest number of sites for all elements except HMS Beagle and Coral, suggesting a recent increase in the activity of some TEs associated with the colonization patterns of Drosophila simulans. The element Tirant had a very distinctive pattern of distribution: it was identified mainly in populations from East Africa and some islands in the Indian Ocean, and its insertion site number was low in all these populations. The data suggest that the genome of the entire species of Drosophila simulans may be being invaded by TEs from populations in which they are present in high copy number.

Testing Models of Selection and Demography inDrosophila simulans

We analyze patterns of nucleotide variability at 15 X-linked loci and 14 autosomal loci from a North American population of Drosophila simulans. We show that there is significantly more linkage disequilibrium on the X chromosome than on chromosome arm 3R and much more linkage disequilibrium on both chromosomes than expected from estimates of recombination rates, mutation rates, and levels of diversity. To explore what types of evolutionary models might explain this observation, we examine a model of recurrent, nonoverlapping selective sweeps and a model of a recent drastic bottleneck (e.g., founder event) in the demographic history of North American populations of D. simulans. The simple sweep model is not consistent with the observed patterns of linkage disequilibrium nor with the observed frequencies of segregating mutations. Under a restricted range of parameter values, a simple bottleneck model is consistent with multiple facets of the data. While our results do not exclude some influence of selection on X vs. autosome variability levels, they suggest that demography alone may account for patterns of linkage disequilibrium and the frequency spectrum of segregating mutations in this population of D. simulans.

Adaptive hitchhiking effects on genome variability

The continuing deluge of nucleotide polymorphism data is providing insights into the role of adaptation in shaping genome-wide patterns of variability and molecular evolution. Population genetic models in which linkage and selection interact (i.e. hitchhiking) predict that selection can leave 'footprints' in closely linked genomic regions. New analytical approaches show promise for distinguishing the signature of adaptation from that of several non-adaptive alternatives. Accounting for the effects of population structure and history poses a challenge for future investigations.

Patterns of DNA sequence variation suggest the recent action of positive selection in the janus-ocnus region of Drosophila simulans

Levels of nucleotide polymorphism in three paralogous Drosophila simulans genes, janusA (janA), janusB (janB), and ocnus (ocn), were surveyed by DNA sequencing. The three genes lie in tandem within a 2.5-kb region of chromosome arm 3R. In a sample of eight alleles from a worldwide distribution we found a significant departure from neutrality by several statistical tests. The most striking feature of this sample was that in a 1.7-kb region containing the janA and janB genes, 30 out of 31 segregating sites contained variants present only once in the sample, and 29 of these unique variants were found in the same allele. A restriction survey of an additional 28 lines of D. simulans revealed strong linkage disequilibrium over the janA-janB region and identified six more alleles matching the rare haplotype. Among the rare alleles, the level of DNA sequence variation was typical for D. simulans autosomal genes and showed no departure from neutrality. In addition, the rare haplotype was more similar to the D. melanogaster sequence, indicating that it was the ancestral form. These results suggest that the derived haplotype has risen to high worldwide frequency relatively recently, most likely as a result of natural selection.