Institut für Tierzucht und Genetik, Veterinärmedizinische Universität Wien
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
For species with a low dispersal rate, the enormous biotic differences are expected to have a significant impact on the populations living on the two slopes of the lower Nahel Oren canyon. Indeed, significant differences in mutation frequencies were reported for wild strains of the fungus Sordaris fimicola isolated from both slopes of the canyon. Most importantly, some of the between slope difference was inherited through two generations of selfing, suggesting fixed genetic differences between the two slopes (Lamb et al. 1998
). Similarly, an almost threefold difference in the number of full-length BARE-1 elements was observed for wild barley individuals collected in the lower Nahel Oren canyon, possibly representing an adaptation to the different environments (Kalendar et al. 2000
).
For a highly mobile species, such as Drosophila melanogaster, which could migrate several miles per day (Coyne and Milstead 1987
), the predictions are less clear. One recent study demonstrated significant differences between SFS and NSF D. melanogaster isofemale lines (Nevo et al. 1998
) for many traits, including survival, longevity, and behavior. Most importantly, the observed differences suggest that they represent an adaptation to the habitat. Therefore, the authors concluded that strong microclimatic natural selection is overriding migration in D. melanogaster at the lower Nahel Oren canyon (Nevo et al. 1998
). Further support for the differentiation between NFS and SFS D. melanogaster populations was provided by a study on the mating behavior of flies collected at this microsite. Using two different sets of synthetic populations for each slope which were maintained in the laboratory for 12 and 48 generations, respectively, the authors showed a significant nonrandom mating pattern (Korol et al. 2000
). Flies significantly preferred mating partners from the same slope. Hence, restriction in gene flow between the two slopes through the mating behavior of the two fly populations may be a crucial component of the adaptation to the different habitats on the NFS and SFS.
Provided D. melanogaster shows significant assortative mating and phenotypic divergence, the SFS and NSF D. melanogaster populations should form separate gene pools. If divergent selection between the two populations is not a recent phenomenon, then NFS and SFS populations are expected to be well differentiated by genetic drift and selection. In this report we use a large number (48) of microsatellite loci, which are highly informative markers, to quantify the genetic divergence of the SFS and NFS populations in the lower Nahel Oren canyon.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
Genetic distances were calculated using the proportion of shared alleles implemented in the MICROSAT software (Minch et al. 1995
). The distance matrix obtained was converted into a dendrogram using the Neighbor-Joining algorithm (Saitou and Nei 1987
) provided with the PHYLIP software package (Felsenstein 1991
) and graphically displayed with TREEVIEW (Page 1996
).
The model-based clustering method of multilocus genotypes was performed with the software STRUCTURE (Pritchard, Stephens, and Donnelly 2000
). One critical parameter for Markov Chain Monte Carlo methods, such as those used in STRUCTURE, is the number of performed iterations. We followed the recommendation of Pritchard, Stephens, and Donnelly (2000)
and used 10,000 iterations for burn-in followed by 100,000 iterations for data collection (Pritchard, Stephens, and Donnelly 2000
). To test whether lnP(X|K) had converged, we made multiple runs and performed additional runs with 500,000 iterations. STRUCTURE assumes Hardy-Weinberg equilibrium within populations and linkage equilibrium among loci.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Genetic differentiation measured by FST was determined for SFS, NFS, and D. simulans. FST values for each pairwise comparison of D. simulans with SFS and NFS D. melanogaster populations were high and statistically significant (0.415 and 0.410, respectively). The FST value between SFS and NFS populations was very small (0.0012) and not statistically significant when was adjusted for multiple testing (Sokal and Rohlf 1995
). Despite the lack of statistical significance, some differentiation between SFS and NFS groups was observed (FST = 0.0012). To further investigate this, we extended our analysis and grouped our samples not only according to site of collection but also according to shipment (table 1
). Only those groups which consisted of at least three individuals were considered. Interestingly, this analysis was not consistent with the hypothesis of a genetic differentiation between slopes. Only the pairwise comparison of SFS II and SFS III was statistically significant after Bonferroni correction for multiple testing (table 3
). Notably, these two samples were collected on different days, and SFS-III flies were second generation rather than first-generation individuals (table 1
). The identical result was obtained when the analysis was repeated with a data set, which was haploidized by randomly discarding one allele at each locus and individual (data not shown). Hence, we could rule out that the significant differentiation between SFS II and III is an artifact produced by the usage of first- and second-generation flies.
|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In the following, we will discuss several scenarios which could potentially reconcile the lack of genetic differentiation for neutral microsatellites in our data set with previous observations.
Genetic Differentiation is Limited to Genomic Regions Subjected to Selection
Microsatellites are highly informative neutral markers (Schlötterer 2000
). In addition to mutation and genetic drift, microsatellite allele frequencies are also affected by linkage to a selected genomic region. If recombination rates are high and population sizes are large, only very tightly linked microsatellites are affected by selection. Assuming that selection has shaped the phenotypic differences between SFS and NFS D. melanogaster populations, genomic regions not targeted by selection may not be strongly affected, as long as population sizes are large and recombination rates are high. Hence, if microsatellites are decoupled from selected regions, they are not expected to show differentiation between SFS and NFS.
In the presence of assortative mating, however, genetic drift will gradually result in genetic differentiation between the two slopes. To reconcile assortative mating and phenotypic differences with the absence of genetic differentiation for neutral microsatellites, one has to assume that the separation of the populations has occurred only recently. Unfortunately, in the absence of information about selection coefficients, the number of selected loci, dominance effects, epistasis, and migration rates it remains speculative whether this scenario is plausible. One prediction could be made, however, that a small number of microsatellite loci, which are linked to selected chromosomal regions, should differentiate SFS and NFS very well. A genome scan using informative markers such as microsatellites could be highly instrumental to identify such genomic regions (Schlötterer and Wiehe 1999
).
Temporal Heterogeneity
A sample of 11 flies which were collected 5 days after the initial collection were found to be significantly differentiated from a sample consisting of 15 individuals from the initial collection. Interestingly, all samples were collected at the same slope, suggesting a larger temporal differentiation than spatial differentiation between SFS and NFS flies. Given that only a difference of a few days resulted in a significant differentiation among samples, one could assume that samples collected at larger time intervals may have even more pronounced differences. On the other hand, the differentiation was modest between temporal samples (FST = 0.03) and could only be detected if individuals were grouped for the analysis. Neither a tree of individuals nor model-based clustering of multilocus genotypes supported the presence of more than a single population. A recent analysis of microsatellite variability in synthetic NFS and SFS populations, which were established from isofemale lines collected in the same year but at different days than our lines, revealed strong population differentiation between the two populations (Michalak et al. 2001
). Hence, unless the different sampling regime in both studies could account for the contrasting microsatellite data (see subsequently), temporal heterogeneity remains an important explanation for the discrepancy between our microsatellite-based results and previous results.
Sampling Regime
In this report we laid special emphasis on the use of freshly collected lines. By the exclusive analysis of first- or second-generation lines, the allele spectrum is not expected to be influenced by the propagation in the laboratory. In contrast, previous studies used special synthetic populations, which were founded by 10 females and males from each of the 25 isofemale lines established for both slopes (Korol et al. 2000
; Michalak et al. 2001
). These synthetic populations have then been propagated in the laboratory under standard conditions. Although this propagation regime is not expected to select for different phenotypes in the SFS and NFS synthetic population, genetic drift could substantially affect them. Without careful recording of the population size in the synthetic populations, it is difficult to account for accelerated drift by minor bottlenecks during the maintenance of the synthetic populations.
In addition to drift, sexual selection during the propagation of the synthetic population could lead to unpredictable results. A recent study demonstrated significant interaction between male and female genotypes determining the success of the sperm (Clark, Begun, and Prout 1999
). Hence, genetic drift could be further accelerated by sexual selection. As a consequence, each synthetic population may have a set of male and female genotypes that have coevolved, assuring the highest fitness of each sex. Although the effect of antagonistic coevolution of males and females on mating preferences has not been tested, it is conceivable that a preference for mating partners of the same line or synthetic population could be observed. For this reason, nonrandom mating would be expected between lines that have been maintained independently for several generations. This hypothesis could account for the observed difference in mating preference between the SFS and NFS synthetic populations (Korol et al. 2000
).
One way to estimate the impact of genetic drift would be a comparison of microsatellite variability in our lines and the synthetic populations. Because of the different levels of variability of D. melanogaster microsatellites (Harr et al. 1998b
), it is important to use the identical loci in such a comparison. Two microsatellite loci are shared between our study and the recent report by Michalak et al. (2001)
. Table 5
compares the gene diversity and number of alleles for these two loci. Despite that a slightly higher number of flies were analyzed for the synthetic populations, the number of alleles was always higher in the freshly collected populations. Similarly, in almost all comparisons, gene diversity was lower in the synthetic populations. Only the NFS population had a slightly higher variability at the DROGPDHA locus in the synthetic population (0.48 vs. 0.43). Hence, although the strong genetic differentiation between SFS and NFS populations reported by Michalak et al. (2001)
may be the result of genetic drift during the propagation of the lines, the important question is whether the behavioral and phenotypic differences could be also attributed to the culture conditions.
|
In summary, our results suggest that on the basis of highly informative neutral markers no significant population differentiation between freshly collected SFS and NFS D. melanogaster populations could be detected. Although several scenarios could be envisioned to reconcile the lack of genetic differentiation with previously described phenotypic and behavioral differences, the most convincing one is the scenario of a recent divergence between SFS and NFS flies. However, further studies on the genetic differentiation of the synthetic populations as well as phenotypic examination of freshly collected flies are required.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Address for correspondence and reprints: Christian Schlötterer, Institut für Tierzucht und Genetik, Veterinärmedizinische Universität Wien, Josef Baumann Gasse 1, 1210 Vienna, Austria. christian.schloetterer{at}vu-wien.ac.at
.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Clark A. G., D. J. Begun, T. Prout, 1999 Female x male interactions in Drosophila sperm competition Science 283:217-220
Coyne J. A., 1992 Genetics and speciation Nature 355:511-515[ISI][Medline]
Coyne J. A., B. Milstead, 1987 Long-distance migration of Drosophila. Dispersal of D. melanogaster alleles from a Maryland orchard Am. Nat 130:70-82[ISI]
Felsenstein J., 1991 PHYLIP (phylogeny inference package). Version 3.57c University of Washington, Seattle
Goudet J., 1995 FSTAT Version 1.2, a program for IBM PC compatibles to calculate Weir and Cockerham's estimators of F-statistics J. Hered 86:485-486[ISI]
Goudet J., M. Raymond, T. de Meeüs, F. Rousset, 1996 Testing differentiation in diploid populations Genetics 144:1933-1940
Harr B., S. Weiss, J. R. David, G. Brem, C. Schlötterer, 1998a. A microsatellite-based multilocus phylogeny of the Drosophila melanogaster species complex Curr. Biol 8:1183-1186[ISI][Medline]
Harr B., B. Zangerl, G. Brem, C. Schlötterer, 1998b. Conservation of locus specific microsatellite variability across species: a comparison of two Drosophila sibling species D. melanogaster and D. simulans Mol. Biol. Evol 15:176-184[Abstract]
Kalendar R., J. Tanskanen, S. Immonen, E. Nevo, A. H. Schulman, 2000 Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence Proc. Natl. Acad. Sci. USA 97:6603-6607
Korol A., E. Rashkovetsky, K. Iliadi, P. Michalak, Y. Ronin, E. Nevo, 2000 Nonrandom mating in Drosophila melanogaster laboratory populations derived from closely adjacent ecologically contrasting slopes at "Evolution canyon." Proc. Natl. Acad. Sci. USA 97:12637-12642
Lamb B. C., M. Saleem, W. Scott, N. Thapa, E. Nevo, 1998 Inherited and environmentally induced differences in mutation frequencies between wild strains of Sordaria fimicola from "Evolution Canyon." Genetics 149:87-99
Meyer A., T. D. Kocher, P. Basasibwaki, A. C. Wilson, 1990 Monophyletic origin of Lake Victoria cichlid fishes suggested by mitochondrial DNA sequences Nature 347:550-553[ISI][Medline]
Michalak P., I. Minkov, A. Helin, D. N. Lerman, B. R. Bettencourt, M. E. Feder, A. B. Korol, E. Nevo, 2001 Genetic evidence for adaptation-driven incipient speciation of Drosophila melanogaster along a microclimatic contrast in "Evolution Canyon," Israel Proc. Natl. Acad. Sci. USA 98:13195-13200
Miller S. A., D. D. Dykes, H. F. Polesky, 1988 A simple salting out procedure for extracting DNA from human nucleated cells Nucleic Acids Res 16:1215.[ISI][Medline]
Minch E., A. Ruiz-Linares, D. Goldstein, M. Feldman, L. L. Cavalli-Sforza, 1995 Microsat (Version 1.4d): a computer program for calculating various statistics on microsatellite allele data
Nevo E., 1995 Asian, African and European biota meet at Evolution Canyon, Israel: local tests of global biodiversity and genetic diversity patterns Proc. R. Soc. Lond. Ser. B 262:149-155[ISI]
Nevo E., E. Rashkovetsky, T. Pavlicek, A. Korol, 1998 A complex adaptive syndrome in Drosophila caused by microclimatic contrasts Heredity 80:9-16[ISI][Medline]
Page R. D. M., 1996 TREEVIEW: an application to display phylogenetic trees on personal computers Comput. Appl. Biosci 12:357-358[Medline]
Pritchard J. K., M. Stephens, P. Donnelly, 2000 Inference of population structure using multilocus genotype data Genetics 155:945-959
Saitou R. K., M. Nei, 1987 The neighbor-joining method: a new method for reconstructing phylogenetic trees Mol. Biol. Evol 4:406-425[Abstract]
Schlötterer C., 1998 Microsatellites Pp. 237261 in A. R. Hoelzel, ed. Molecular genetic analysis of populations: a practical approach. 2nd edition. Oxford University Press, Oxford
. 2000 Evolutionary dynamics of microsatellite DNA Chromosoma 109:365-371[ISI][Medline]
Schlötterer C., T. Wiehe, 1999 Microsatellites, a neutral marker to infer selective sweeps Pp. 238248 in D. Goldstein and C. Schlötterer, eds. Microsatellitesevolution and applications. Oxford University Press, Oxford
Schlötterer C., B. Zangerl, 1999 The use of imperfect microsatellites for DNA fingerprinting and population genetics Pp. 153165 in J. T. Epplen and T. Lubjuhn, eds. DNA profiling and DNA fingerprinting. Birkhäuser, Basel
Sokal R. R., F. J. Rohlf, 1995 Biometry W. H. Freeman and Company, New York