Genetic variation and differentiation at microsatellite loci in Drosophila simulans. Evidence for founder effects in new world populations

Fixation times of de novo and standing beneficial variants in subdivided populations

The rate at which beneficial alleles fix in a population depends on the probability of and time to fixation of such alleles. Both of these quantities can be significantly impacted by population subdivision and limited gene flow. Here, we investigate how limited dispersal influences the rate of fixation of beneficial de novo mutations, as well as fixation time from standing genetic variation. We investigate this for a population structured according to the island model of dispersal allowing us to use the diffusion approximation, which we complement with simulations. We find that fixation may take on average fewer generations under limited dispersal than under panmixia when selection is moderate. This is especially the case if adaptation occurs from de novo recessive mutations, and dispersal is not too limited (such that approximately FST<0.2). The reason is that mildly limited dispersal leads to only a moderate increase in effective population size (which slows down fixation), but is sufficient to cause a relative excess of homozygosity due to inbreeding, thereby exposing rare recessive alleles to selection (which accelerates fixation). We also explore the effect of metapopulation dynamics through local extinction followed by recolonization, finding that such dynamics always accelerate fixation from standing genetic variation, while de novo mutations show faster fixation interspersed with longer waiting times. Finally, we discuss the implications of our results for the detection of sweeps, suggesting that limited dispersal mitigates the expected differences between the genetic signatures of sweeps involving recessive and dominant alleles.

Genetic Diversity and Variation in Mitochondrial COI Gene in Wild and Hatchery Populations of Saxidomus purpuratus

To investigate the genetic diversity and genetic variations of four wild (Geoje, Jinhae, Yeosu, and Boryeong) and two hatchery (Goheung and Geoje) populations of purplish Washington clam (Saxidomus purpuratus), 421 bp sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene were analyzed. A total of 149 haplotypes were identified from 358 individuals from the four wild and two hatchery populations with 109 substitutions. The genetic diversity of the wild populations and Geoje hatchery population were high, whereas the total number of haplotypes, population-specific haplotypes, and haplotype diversity were comparatively low in the Goheung hatchery population. The fixation index (F_ST) indicated that there was no significant genetic difference between the four wild populations. However, the Goheung hatchery population was significantly different from that of the Geoje hatchery, exhibiting the most pronounced difference, and two wild populations (Jinhae and Yeosu). The low genetic diversity indices exhibited by the Goheung hatchery population might have resulted from farm propagation using a limited parental stock. Therefore, to maintain genetic diversity, a proper breeding management program using more progenitors is required in hatcheries, in addition to regular monitoring of both hatchery and wild populations.

Light dependent courtship behavior in Drosophila simulans and D. melanogaster

Differences in courtship signals and perception are well-known among Drosophila species. One such described difference is the dependency on light, and thus presumably vision, for copulation success. Many studies have described a difference in light-dependent copulation success between D. melanogaster and D. simulans, identifying D. simulans as a light-dependent species, and D. melanogaster as a light-independent one. However, many of these studies use assays of varying design and few strains to represent the entire species. Here, we attempt to better characterize this purported difference using 11 strains of each species, paired by collection location, in behavioral assays conducted at two different exposure times. We show that, while there is a species-wide difference in magnitude of light-dependent copulation success, D. melanogaster copulation success is, on average, still impaired in the dark at both exposure times we measured. Additionally, there is significant variation in strain-specific ability to copulate in the dark in both species across two different exposure times. We find that this variation correlates strongly with longitude in D. melanogaster, but not in D. simulans. We hypothesize that differences in species history and demography may explain behavioral variation. Finally, we use courtship assays to show that light-dependent copulation success in one D. simulans strain is driven in part by both males and females. We discuss potential differences in courtship signals and/or signal importance between these species and potential for further comparative studies for functional characterization.

Distinguishing migration from isolation using genes with intragenic recombination: detecting introgression in the Drosophila simulans species complex

BACKGROUND

Determining the presence or absence of gene flow between populations is the target of some statistical methods in population genetics. Until recently, these methods either avoided the use of recombining genes, or treated recombination as a nuisance parameter. However, genes with recombination contribute additional information for the detection of gene flow (i.e. through linkage disequilibrium).

METHODS

We present three summary statistics based on the spatial arrangement of fixed differences, and shared and exclusive polymorphisms that are sensitive to the presence and direction of gene flow. Power and false positive rate for tests based on these statistics are studied by simulation.

RESULTS

The application of these tests to populations from the Drosophila simulans species complex yielded results consistent with migration between D. simulans and its two endemic sister species D. mauritiana and D. sechellia, and between populations D. mauritiana on the islands of the Mauritius and Rodrigues.

CONCLUSIONS

We demonstrate the sensitivity of the developed statistics to the presence and direction of gene flow, and characterize their power as a function of differentiation level and recombination rate. The properties of these statistics make them especially suitable for analyzing high-throughput sequencing data or for their integration within the approximate Bayesian computation framework.

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.

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.

Population genomic inferences from sparse high-throughput sequencing of two populations of Drosophila melanogaster

Short-read sequencing techniques provide the opportunity to capture genome-wide sequence data in a single experiment. A current challenge is to identify questions that shallow-depth genomic data can address successfully and to develop corresponding analytical methods that are statistically sound. Here, we apply the Roche/454 platform to survey natural variation in strains of Drosophila melanogaster from an African (n = 3) and a North American (n = 6) population. Reads were aligned to the reference D. melanogaster genomic assembly, single nucleotide polymorphisms were identified, and nucleotide variation was quantified genome wide. Simulations and empirical results suggest that nucleotide diversity can be accurately estimated from sparse data with as little as 0.2x coverage per line. The unbiased genomic sampling provided by random short-read sequencing also allows insight into distributions of transposable elements and copy number polymorphisms found within populations and demonstrates that short-read sequencing methods provide an efficient means to quantify variation in genome organization and content. Continued development of methods for statistical inference of shallow-depth genome-wide sequencing data will allow such sparse, partial data sets to become the norm in the emerging field of population genomics.

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.

Similar levels of X-linked and autosomal nucleotide variation in African and non-African populations of Drosophila melanogaster

Background

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

Results

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

Conclusion

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

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.

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.

Microsatellite analysis of olive fly populations in the Mediterranean indicates a westward expansion of the species

Bactrocera oleae is the major insect pest of the olive fruit. Twelve microsatellite loci isolated from the genome of this insect were used in a Mediterranean-wide population analysis. These loci were highly polymorphic with a mean number of alleles per locus of 10.42 and a mean effective number of alleles of 2.76. The analysis was performed on a sample of 671 flies collected from nineteen locations around the European part of the Mediterranean basin. Despite the high level of gene flow across the Mediterranean, results support the notion of a differentiation of three subpopulations: one of the Iberian Peninsula, one of Greece and Italy and one of Cyprus. In addition, the gradual decrease of heterozygosity from the Eastern to the Western part of the Mediterranean indicates a westward expansion of the species.

A pseudohitchhiking model of X vs. autosomal diversity

We study levels of X-linked vs. autosomal diversity using a model developed to analyze the hitchhiking effect. Repeated bouts of hitchhiking are thought to lower X-linked diversity for two reasons: first, because sojourn times of beneficial mutations are shorter on the X, and second, because adaptive substitutions may be more frequent on the X. We investigate whether each of these effects does, in fact, cause reduced X-linked diversity under hitchhiking. We study the strength of the hitchhiking effect on the X vs. autosomes when there is no recombination and under two different recombination schemes. When recombination occurs in both sexes, X-linked vs. autosomal diversity is reduced by hitchhiking under a broad range of conditions, but when there is no recombination in males, as in Drosophila, the required conditions are considerably more restrictive.

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.

Similar mating and sperm displacement patterns in two highly divergent D. simulans populations from Africa and Europe

The frequency of remating in Drosophila melanogaster is affected by both genetic and ecological factors. We studied the remating behaviour in one European (Italy) and one African (Uganda) Drosophila simulans population using six highly polymorphic microsatellite markers. Despite that the populations were genetically distinct (F(ST) = 0.18) and originated from very dissimilar ecological settings with different population densities, we inferred a very similar mating pattern. The remating parameter alpha was similar in both populations (a = 1.3-1.4). No more than two distinct paternal genotypes per family were detected in each population.

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.

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.

Rapid recovery of genetic diversity of dogwhelk (Nucella lapillus L.) populations after local extinction and recolonization contradicts predictions from life-history characteristics

The dogwhelk Nucella lapillus is a predatory marine gastropod populating North Atlantic rocky shores. As with many other gastropod species, N. lapillus was affected by tributyltin (TBT) pollution during the 1970s and 1980s, when local populations became extinct. After a partial ban on TBT in the United Kingdom in 1987, vacant sites have been recolonized. N. lapillus lacks a planktonic larval stage and is therefore expected to have limited dispersal ability. Relatively fast recolonization of some sites, however, contradicts this assumption. We compared levels of genetic diversity and genetic structuring between recolonized sites and sites that showed continuous population at three localities across the British Isles. No significant genetic effects of extinction/recolonization events were observed in SW Scotland and NE England. In SW England we observed a decrease in genetic diversity and an increase in genetic structure in recolonized populations. This last result could be an artefact, however, due to the superposition of other local factors influencing the genetic structuring of dogwhelk populations. We conclude that recolonization of vacant sites was accomplished by a relatively high number of individuals originating from several source populations (the 'migrant-pool' model of recolonization), implying that movements are more widespread than expected on the basis of development mode alone. Comparison with published data on genetic structure of marine organisms with contrasted larval dispersal supports this hypothesis. Our results also stress the importance of local factors (geographical or ecological) in determining genetic structure of dogwhelk populations.

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.

Ecological determinants of genetic diversity in an expanding population of the shrub Myrica cerifera

The ecological mechanisms that contribute to the acquisition of genetic diversity in an expanding population of the shrub, Myrica cerifera, on an island habitat were investigated. Genealogical reconstruction was used to assess the contribution of early reproductive colonists to subsequent recruitment. In addition, through determination of parentage, the source of recruiting seedlings was identified and the contribution of seed and pollen dispersal into the colonizing sites was inferred. The relative contribution of different sources of gene flow was determined directly and an investigation was made into how variability in breeding patterns may have contributed to observed levels of genetic variability. It was expected that early colonists that could flower would contribute to subsequent recruiting cohorts, and that the limited number of such early reproductive colonists would lead to variance in mating success, inbreeding, or bottlenecks which could reduce genetic diversity and increase genetic differentiation among subsequent recruiting cohorts. Analyses of parentage (with paternity exclusion probability > 95%) for all recruiting plants demonstrated that in fact, there was little contribution by the early reproductive colonists to subsequent cohorts, and that immigration from outside the study sites in the form of seed dispersal accounted for over 94% of the recruitment in the study plots, with pollen dispersal accounting for less than 3% gene flow. No genetic bottleneck or evidence of reproductive skew in the recruiting cohorts were found, suggesting that propagule dispersal was from many source individuals in other established populations.

The population genetics of a biological control introduction: mitochondrial DNA and microsatellie variation in native and introduced populations of Aphidus ervi, a parisitoid wasp

Introductions of biological control agents may cause bottlenecks in population size despite efforts to avoid them. We examined the population genetics of Aphidius ervi (Hymenoptera: Braconidae), a parasitoid that was introduced to North America from Western Europe in 1959 to control pea aphids. To explore the phylogeographical relationships of A. ervi we sequenced 1249 bp of mitochondrial DNA (mtDNA) from 27 individuals from the native range and 51 individuals from the introduced range. Most individuals from Western Europe, the Middle East and North America shared one of two common haplotypes, consistent with the known history of the introduction. However, some A. ervi from the Pacific Northwest have a haplotype that is most similar to haplotypes found in Japan, raising the possibility of a second accidental introduction. To examine population structure and assess whether a bottleneck occurred upon introduction to North America, we assayed variation at 5 microsatellite loci in 62 individuals from 2 native populations and 230 individuals from 6 introduced populations. Introduced samples had fewer rare alleles than native samples (F1,34 = 13.5, P = 0.0008), but heterozygosity did not differ significantly. These results suggest that a mild bottleneck occurred in spite of the introduction of over 1000 individuals. Using a hierarchical Bayesian approach, the founding population size was estimated to be 245 individuals. amova showed significant genetic differentiation between the European and North American samples, and a Bayesian assignment approach clustered individuals into four groups, with most European individuals in one group and most North American individuals in the other three. These results highlight that genetic changes are associated with founder events in rapidly growing natural populations, even when the founding population size is relatively large.

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.

Nucleotide polymorphism in the Est6 promoter, which is widespread in derived populations of Drosophila melanogaster, changes the level of Esterase 6 expressed in the male ejaculatory duct

Previous analysis of an Australian population of D. melanogaster revealed two predominant Est6 promoter haplotypes, P1 and P7. These haplotypes, which differ at 14 sites over a 325-bp region, are associated with a 15-20% difference in male EST6 activity. Here we show that the P1/P7 sequence difference causes the male activity variation by recreating the activity difference among >60 independently transformed lines containing representative P1 or P7 promoter alleles fused to an identical Est6 coding region. Furthermore we find that the whole fly difference reflects about a twofold difference in EST6 activity in the anterior sperm ejaculatory duct. EST6 activity variation in this tissue is known to affect reproductive fitness. Using a combination of RFLP analysis and DNA sequencing, we show that P1 and P7 are predominant in six populations from America, Asia, and Australia, albeit less frequent in a population from the presumptively ancestral east African range of the species. The sequence data show significant departures from neutral expectations for the derived American and Australian populations but not the presumptively ancestral Zimbabwean population. Thus the P1/P7 difference could be a major source of adaptively significant EST6 activity variation through much of the now cosmopolitan range of D. melanogaster.

The neutral theory in the genomic era

A number of tests have been developed to detect positive selection at the molecular level. These tests are based on DNA polymorphism within and divergence between species. Applications of these tests have revealed a large collection of genes that have evolved under positive selection and some general insights into adaptive evolution. Recently, these tests have been applied on a genomic scale and have provided estimates of the frequency of adaptive substitutions and a critical test of the neutral theory.

Genetic evidence for adaptation-driven incipient speciation of Drosophila melanogaster along a microclimatic contrast in "Evolution Canyon," Israel

Substantial genetic differentiation, as great as among species, exists between populations of Drosophila melanogaster inhabiting opposite slopes of a small canyon. Previous work has shown that prezygotic sexual isolation and numerous differences in stress-related phenotypes have evolved between D. melanogaster populations in "Evolution Canyon," Israel, in which slopes 100-400 m apart differ dramatically in aridity, solar radiation, and associated vegetation. Because the canyon's width is well within flies' dispersal capabilities, we examined genetic changes associated with local adaptation and incipient speciation in the absence of geographical isolation. Here we report remarkable genetic differentiation of microsatellites and divergence in the regulatory region of hsp70Ba which encodes the major inducible heat shock protein of Drosophila, in the two populations. Additionally, an analysis of microsatellites suggests a limited exchange of migrants and lack of recent population bottlenecks. We hypothesize that adaptation to the contrasting microclimates overwhelms gene flow and is responsible for the genetic and phenotypic divergence between the populations.

Dynamics of microsatellite divergence under stepwise mutation and proportional slippage/point mutation models

Recently Kruglyak, Durrett, Schug, and Aquadro showed that microsatellite equilibrium distributions can result from a balance between polymerase slippage and point mutations. Here, we introduce an elaboration of their model that keeps track of all parts of a perfect repeat and a simplification that ignores point mutations. We develop a detailed mathematical theory for these models that exhibits properties of microsatellite distributions, such as positive skewness of allele lengths, that are consistent with data but are inconsistent with the predictions of the stepwise mutation model. We use our theoretical results to analyze the successes and failures of the genetic distances (delta(mu))(2) and D(SW) when used to date four divergences: African vs. non-African human populations, humans vs. chimpanzees, Drosophila melanogaster vs. D. simulans, and sheep vs. cattle. The influence of point mutations explains some of the problems with the last two examples, as does the fact that these genetic distances have large stochastic variance. However, we find that these two features are not enough to explain the problems of dating the human-chimpanzee split. One possible explanation of this phenomenon is that long microsatellites have a mutational bias that favors contractions over expansions.

The frequency distribution of nucleotide variation in Drosophila simulans

Patterns of codon bias in Drosophila suggest that silent mutations can be classified into two types: unpreferred (slightly deleterious) and preferred (slightly beneficial). Results of previous analyses of polymorphism and divergence in Drosophila simulans were interpreted as supporting a mutation-selection-drift model in which slightly deleterious, silent mutants make significantly greater contributions to polymorphism than to divergence. Frequencies of unpreferred polymorphisms were inferred to be lower than frequencies of other silent polymorphisms. Here, I analyzed additional D. simulans data to reevaluate the support for these ideas. I found that D. simulans has fixed more unpreferred than preferred mutations, suggesting that this lineage has not been at mutation-selection-drift equilibrium at silent sites. Frequencies of polarized unpreferred polymorphisms are not skewed toward rare alleles. However, frequencies of unpolarized unpreferred codons are lower in high-bias genes than in low-bias genes. This supports the idea that unpreferred codons are borderline deleterious mutations. Purifying selection on silent sites appears to be stronger at twofold-degenerate codons than at fourfold-degenerate codons. Finally, I found that X-linked polymorphisms occur at a higher average frequency than polymorphisms on chromosome arm 3R, even though an average X-linked site is significantly less likely to be polymorphic than an average site on 3R. This result supports a previous analysis of D. simulans indicating different population genetics of X-linked versus autosomal mutations.

Regions of lower crossing over harbor more rare variants in African populations of Drosophila melanogaster

A correlation between diversity levels and rates of recombination is predicted both by models of positive selection, such as hitchhiking associated with the rapid fixation of advantageous mutations, and by models of purifying selection against strongly deleterious mutations (commonly referred to as "background selection"). With parameter values appropriate for Drosophila populations, only the first class of models predicts a marked skew in the frequency spectrum of linked neutral variants, relative to a neutral model. Here, we consider 29 loci scattered throughout the Drosophila melanogaster genome. We show that, in African populations, a summary of the frequency spectrum of polymorphic mutations is positively correlated with the meiotic rate of crossing over. This pattern is demonstrated to be unlikely under a model of background selection. Models of weakly deleterious selection are not expected to produce both the observed correlation and the extent to which nucleotide diversity is reduced in regions of low (but nonzero) recombination. Thus, of existing models, hitchhiking due to the recurrent fixation of advantageous variants is the most plausible explanation for the data.

Microsatellite variation and fine-scale population structure in the wood frog (Rana sylvatica)

We investigated genetic population structure in wood frogs (Rana sylvatica) from a series of Prairie Pothole wetlands in the northern Great Plains. Amphibians are often thought to exist in demographic metapopulations, which require some movement between populations, yet genetic studies have revealed strong subdivision among populations, even at relatively fine scales (several km). Wood frogs are highly philopatric and studies of dispersal suggest that they may exhibit subdivision on a scale of approximately 1-2 km. We used microsatellites to examine population structure among 11 breeding assemblages separated by as little as 50 m up to approximately 5.5 km, plus one population separated from the others by 20 km. We found evidence for differentiation at the largest distances we examined and among a few neighbouring ponds, but most populations were strikingly similar in allele frequencies, suggesting high gene flow among all but the most distant populations. We hypothesize that the few significant differences among neighbouring populations at the finest scale may be a transient effect of extinction-recolonization founder events, driven by periodic drying of wetlands in this hydrologically dynamic landscape.

The Y chromosomes of Drosophila simulans are highly polymorphic for their ability to suppress sex-ratio drive

The sex-ratio trait, known in several species of Drosophila including D. simulans, results from meiotic drive of the X chromosome against the Y. Males that carry a sex-ratio X chromosome produce strongly female-biased progeny. In D. simulans, drive suppressors have evolved on the Y chromosome and on the autosomes. Both the frequency of sex-ratio X and the strength of the total drive suppression (Y-linked and autosomal) vary widely among geographic populations of this worldwide species. We have investigated the pattern of Y-linked drive suppression in six natural populations representative of this variability. Y-linked suppressors were found to be a regular component of the suppression, with large differences between populations in the mean level of suppression. These variations did not correspond to differences in frequency of discrete types of Y chromosomes, but to a more or less wide continuum of phenotypes, from nonsuppressor to partial or total suppressor. We concluded that a large diversity of Y-linked suppressor alleles exists in D. simulans and that some populations are highly polymorphic. Our results support the hypothesis that a Y-chromosome polymorphism can be easily maintained by a balance between meiotic drive and the cost of drive suppression.

Increase of genetic variation over time in a recently founded population of great reed warblers (Acrocephalus arundinaceus) revealed by microsatellites and DNA fingerprinting

Genetic similarity within pairs of individuals was examined using both 10 polymorphic microsatellite loci and multi-locus DNA fingerprinting profiles in a semi-isolated population of great reed warblers at Lake Kvismaren, south Central Sweden, in 1987-1993. The population was founded by a few individuals in 1978, followed by a gradual increase in numbers until 1988, since when the population has remained relatively stable with about 60 breeding birds. We have previously found that high genetic similarity between pair-mates in the population during the early part of the study period reduced egg hatching success, and hence reproductive success. The measures of pairwise genetic similarity, microsatellite allele sharing and DNA fingerprinting band sharing, were highly correlated with pedigree-based relatedness. Both microsatellite and DNA fingerprinting similarities between pair-mates declined significantly over the study period, and the pattern was most pronounced in the DNA fingerprinting data. Analyses restricted to the microsatellite data showed that the average annual microsatellite similarity between pairwise combinations of individuals, as well as individual homozygosity in males, declined significantly over the study period, and that several immigrants carrying novel alleles entered the population during the study. Hence, the temporal decline in genetic similarity of mates in the population is probably a consequence of increased immigration, facilitated by the recent expansion of the species in the region. These results suggest that the population has now recovered genetically, or is in the process of recovering, from a recent founder event.

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.

A genome-wide departure from the standard neutral model in natural populations of Drosophila

We analyze nucleotide polymorphism data for a large number of loci in areas of normal to high recombination in Drosophila melanogaster and D. simulans (24 and 16 loci, respectively). We find a genome-wide, systematic departure from the neutral expectation for a panmictic population at equilibrium in natural populations of both species. The distribution of sequence-based estimates of 2Nc across loci is inconsistent with the assumptions of the standard neutral theory, given the observed levels of nucleotide diversity and accepted values for recombination and mutation rates. Under these assumptions, most estimates of 2Nc are severalfold too low; in other words, both species exhibit greater intralocus linkage disequilibrium than expected. Variation in recombination or mutation rates is not sufficient to account for the excess of linkage disequilibrium. While an equilibrium island model does not seem to account for the data, more complicated forms of population structure may. A proper test of alternative demographic models will require loci to be sampled in a more consistent fashion.

High density of long dinucleotide microsatellites in Drosophila subobscura

We isolated 96 dinucleotide repeats with five or more tandemly repeated units from a subgenomic Drosophila subobscura library. The mean repeat unit length of microsatellite clones in D. subobscura is 15, higher than that observed in other Drosophila species. Population variation was assayed in 32-40 chromosomes from Barcelona, Spain, using 18 randomly chosen microsatellite loci. Positive correlation between measures of variation and perfect repeat length measures (mean size, most common, and longest allele) is consistent with a higher mutation rate in loci with longer repeat units. Levels of microsatellite variation measured as variance in repeat number and heterozygosity in D. subobscura were similar to those of Drosophila pseudoobscura and higher than those of Drosophila melanogaster and Drosophila simulans. Our data suggest that higher levels of microsatellite variation, and possibly density, in D. subobscura compared with D. melanogaster are due to both a higher average effective population and a higher intrinsic slippage rate in the former species.

Microsatellite polymorphism in the Mediterranean fruit fly, Ceratitis capitata

A total of forty-three simple sequence repeats (SSRs) were identified in the Mediterranean fruit fly (medfly) Ceratitis capitata. The most common SSR was the dinucleotide (TG)n/(CA)n occurring in thirty of the forty-three microsatellite loci. Polymorphism at ten dinucleotide markers was investigated in 122 flies from six natural populations sampled in the native and colonized areas. A very high level of allelic diversity was detected in the species range. An average of 13.6 alleles was found over all the ten loci indicating the informativeness of SSRs as genetic markers for the medfly. The distribution of microsatellite polymorphism in the species range reflects the medfly colonization history.

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.

Genetic variation in the spread of Drosophila subobscura from a nonequilibrium population

Drosophila subobscura was first identified in North America in the early 1980s, and a newer D. subobscura population in Utah appears to have been established more than 10 years later. In this study, we use nuclear microsatellite allele frequencies, mitochondrial restriction fragment length polymorphism (RFLP) allele frequencies, and computer simulations to investigate possible scenarios of how this species has spread across North America. Our method develops a 95% confidence interval for the maximum and minimum number of founders that could have colonized the new population given various scenarios for spread. Unlike many other methods, it may be applied to nonequilibrium source populations given certain conditions. We find that observed allele frequency differences between newer and older D. subobscura populations are consistent with very few inseminated females being transported east from the coast in a single step or with larger numbers of colonizers invading after several intermediate steps. They are not consistent with a large, panmictic population of D. subobscura colonizing Utah in a single step.

Unusual haplotype structure at the proximal breakpoint of In(2L)t in a natural population of Drosophila melanogaster

The existence of temporally stable frequency clines for In(2L)t in natural populations of Drosophila melanogaster suggests a role for selection in the maintenance of this polymorphism. We have collected nucleotide polymorphism data from the proximal breakpoint junction regions of In(2L)t to infer its evolutionary history. The finding of a novel LINE-like element near the In(2L)t breakpoint junction in sampled inverted chromosomes supports a transposable element-mediated origin for this inversion. An analysis of nucleotide variation in a Costa Rican population sample of standard and inverted chromosomes indicates a unique and relatively recent origin for In(2L)t. Additional In(2L)t alleles from three geographically diverse populations reveal no detectable geographic differentiation. Low levels of In(2L)t nucleotide polymorphism suggest a recent increase in the inversion's frequency in tropical populations. An unusual feature of our sample of standard alleles is a marked heterogeneity in levels of linkage disequilibrium among polymorphic sites across the breakpoint region. We introduce a test of neutral equilibrium haplotype structure that corrects both for multiple tests and for an arbitrarily chosen window size. It reveals that an approximately 1.4-kb region immediately spanning the breakpoint has fewer haplotypes than expected under the neutral model, given the expected level of recombination in this genomic region. Certain features of our data suggest that the unusual pattern in standard chromosomes is the product of selection rather than demography.

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

Previous studies based on allozyme variation have found little evidence for genetic differentiation in Drosophila simulans. On the basis of DNA sequence variation at two nuclear loci in four African populations of D. simulans, we show that there is significant structure to D. simulans populations within Africa. Variation at one of the loci, vermilion, appears to be neutral and supports an eastern African origin for European and American populations. Samples from the West Indies, Europe, and North America had a nucleotide diversity lower than that of African populations at vermilion and show nonequilibrium haplotype distributions at both vermilion and G6pd, consistent with a hypothesis of recent bottleneck and possibly also admixture in the history of these populations. Directional selection, previously documented at G6pd, appears to have occurred within the coalescence time of the species, obscuring deep population history.