Evolutionary pace of chromosomal polymorphism in colonizing populations of Drosophila subobscura: an evolutionary time series

Abstract. Biologists have long debated the speed, uniformity, and predictability of evolutionary change. However, evaluating such patterns on a geographic scale requires time-series data on replicate sets of natural populations. Drosophila subobscura has proven an ideal model system for such studies. This fly is broadly distributed in the Old World, but was introduced into both North and South America just over two decades ago and then spread rapidly. Rapid, uniform, and predictable evolution would be demonstrated if the invading flies evolved latitudinal clines that progressively converged on those of the native populations. Evolutionary geneticists quickly capitalized on this opportunity to monitor evolutionary dynamics. Just a few years after the introduction, they surveyed chromosomal inversion frequencies in both North and South America. On both continents they detected incipient latitudinal clines in chromosome inversion frequencies that almost always had the same sign with latitude as in the Old World. Thus the initial evolution of chromosomal polymorphisms on a continental scale was remarkably rapid and consistent. Here we report newer samples of inversion frequencies for the colonizing populations: the time series now spans almost one decade for North America and almost two decades for South America. Almost all inversions in the New World continue to show the same sign of frequency with latitude as in the Old World. Nevertheless, inversion clines have not consistently increased in steepness over time; nor have they consistently continued to converge on the Old World baseline. However, five arrangements in South America show directional, continentwide shifts in frequency. Overall, the initial consistency of clinal evolutionary trajectories seen in the first surveys seems not to have been maintained.     

Strong founder effect in Drosophila pseudoobscura colonizing New Zealand from North America

The North American native species Drosophila pseudoobscura was first identified in New Zealand in the last few decades. Here, we have studied the genetic consequences of its spread across the Pacific Ocean. Using 10 microsatellites that are highly variable in North American populations, we found that the New Zealand population has substantially fewer alleles, a much lower average heterozygosity, and significantly different allele frequencies at these loci. We have discussed the relative sensitivity of these parameters for detecting the founding event. X-linked loci were more strongly differentiated between continents than autosomal loci, as reflected by larger changes in allele frequencies and greater reductions in numbers of alleles and average heterozygosity. The severity of the genetic diversity loss supports a scenario of a few D. pseudoobscura females being introduced to New Zealand from North America.     

Swift laboratory thermal evolution of wing shape (but not size) in Drosophila subobscura and its relationship with chromosomal inversion polymorphism

Latitudinal clinal variation in wing size and shape has evolved in North American populations of Drosophila subobscura within about 20 years since colonization. While the size cline is consistent to that found in original European populations (and globally in other Drosophila species), different parts of the wing have evolved on the two continents. This clearly suggests that 'chance and necessity' are simultaneously playing their roles in the process of adaptation. We report here rapid and consistent thermal evolution of wing shape (but not size) that apparently is at odds with that suggestion. Three replicated populations of D. subobscura derived from an outbred stock at Puerto Montt (Chile) were kept at each of three temperatures (13, 18 and 22 degrees C) for 1 year and have diverged for 27 generations at most. We used the methods of geometric morphometrics to study wing shape variation in both females and males from the thermal stocks, and rates of genetic divergence for wing shape were found to be as fast or even faster than those previously estimated for wing size on a continental scale. These shape changes did not follow a neat linear trend with temperature, and are associated with localized shifts of particular landmarks with some differences between sexes. Wing shape variables were found to differ in response to male genetic constitution for polymorphic chromosomal inversions, which strongly suggests that changes in gene arrangement frequencies as a response to temperature underlie the correlated changes in wing shape because of gene-inversion linkage disequilibria. In fact, we also suggest that the shape cline in North America likely predated the size cline and is consistent with the quite different evolutionary rates between inversion and size clines. These findings cast strong doubts on the supposed 'unpredictability' of the geographical cline for wing traits in D. subobscura North American colonizing populations.     

Tracking the origin of the American colonization by Drosophila subobscura: genetic comparison between Eastern and Western Mediterranean populations

Several pieces of evidence indicate a Mediterranean origin of the colonization of America by Drosophila subobscura . To ascertain whether the origin was from the Eastern or the Western Mediterranean region, samples from Barcelona (Spain) and Mt Parnes (Greece) were collected and O chromosomal inversion polymorphism and lethal genes were analysed. The frequencies of lethal chromosomes were 0.244 ± 0.039 in Barcelona and 0.336 ± 0.043 in Mt Parnes, consistent with the expectations for large populations located in the central area of the species distribution. Lethal genes seem to be distributed at random along the O chromosome in both populations. The intra-populational allelism frequencies of Barcelona and Mt Parnes were 0.016 ± 0.007 and 0.012 ± 0.005 respectively. Thus, the estimates of the effective population size were high in both populations (between 6964 and 13 004 in Barcelona and 11 874 to 26 828 in Mt Parnes). The cases of allelism in Mt Parnes were observed only between individual lethal genes, but in Barcelona some concatenated clusters of allelism were detected. This pattern of allelism can be explained by synthetic lethality, hybrid dysgenesis, the induction of recurrent lethal mutations by different factors or an effect of microdifferentiation in subpopulations. In both populations, a reduction in fitness in the heterozygotes for lethal genes has been detected. Furthermore, the estimates of the migration coefficient (between 0.0085 and 0.0120 in Barcelona, and 0.0057 and 0.0087 in Mt Parnes) confirm the existence of gene flow between Palearctic populations of D. subobscura . Our lethal genes and chromosomal inversion results are consistent with a Mediterranean origin of the colonization of America by D. subobscura , but are inconclusive with regard to the identification of the population from which these colonizers came.     

Introduction history of Drosophila subobscura in the New World: a microsatellite-based survey using ABC methods

Drosophila subobscura is a Palearctic species that was first observed in South and North America in the early 1980s, and that rapidly invaded broad latitudinal ranges on both continents. To trace the source and history of this invasion, we obtained genotypic data on nine microsatellite loci from two South American, two North American and five European populations of D. subobscura. We analysed these data with traditional statistics as well as with an approximate Bayesian computation (ABC) framework. ABC methods yielded the strongest support for the scenario involving a serial introduction with founder events from Europe into South America, and then from South America into North America. Stable effective population size of the source population was very large (around one million individuals), and the propagule size was notably smaller for the introduction into South America (i.e. high bottleneck severity index with only a few effective founders) but considerably larger for the subsequent introduction into North America (i.e. low bottleneck severity index with around 100-150 effective founders). Finally, the Mediterranean region of Europe (and most likely Barcelona from the localities so far analysed) is proposed as the source of the New World flies, based on mean individual assignment statistics.     

Rapid deployment of the five founding Amerind mtDNA haplogroups via coastal and riverine colonization

Numerous studies of variation in mtDNA in Amerindian populations established that four haplogroups are present throughout both North and South America. These four haplogroups (A, B, C, and D) and perhaps a fifth (X) in North America are postulated to be present in the initial founding migration to the Americas. Furthermore, studies of ancient mtDNA in North America suggested long-term regional continuity of the frequencies of these founding haplogroups. Present-day tribal groups possess high frequencies of private mtDNA haplotypes (variants within the major haplogroups), consistent with early establishment of local isolation of regional populations. Clearly these patterns have implications for the mode of colonization of the hemisphere. Recently, the earlier consensus among archaeologists for an initial colonization by Clovis hunters arriving through an ice-free corridor and expanding in a "blitzkrieg " wave was shown to be inconsistent with extensive genetic variability in Native Americans; a coastal migration route avoids this problem. The present paper demonstrates through a computer simulation model how colonization along coasts and rivers could have rapidly spread the founding lineages widely through North America.     

A time series of evolution in action: a latitudinal cline in wing size in South American Drosophila subobscura

Drosophila subobscura is geographically widespread in the Old World. Around the late 1970s, it was accidentally introduced into both South and North America, where it spread rapidly over broad latitudinal ranges. This invading species offers opportunities to study the speed and predictability of trait evolution on a geographic scale. One trait of special interest is body size, which shows a strong and positive latitudinal cline in many Drosophila species, including Old World D. subobscura. Surveys made about a decade after the invasion found no evidence of a size cline in either North or South America. However, a survey made in North America about two decades after the invasion showed that a conspicuous size cline had evolved and (for females) was coincident with that for Old World flies. We have now conducted parallel studies on 10 populations (13 degrees of latitude) of flies, collected in Chile in spring 1999. After rearing flies in the laboratory for several generations, we measured wing sizes and compared geographic patterns (versus latitude or temperature) for flies on all three continents. South American females have now evolved a significant latitudinal size cline that is similar in slope to that of Old World and of North American flies. Rates of evolution (haldanes) for females are among the highest ever measured for quantitative traits. In contrast, the size cline is positive but not significant for South or North American males. At any given latitude, South American flies of both sexes are relatively large; this in part reflects the relatively cool climate of coastal Chile. Interestingly, the sections of the wing that generate the size cline for females differ among all three continents. Thus, although the evolution of overall wing size is predictable on a geographic scale (at least for females), the evolution of size of particular wing components is decidedly not.     

Colonization of America by Drosophila subobscura: spatial and temporal lethal-gene allelism

About twenty years ago Drosophila subobscura, a western Palearctic species, colonized both North and South America. Lethal genes in the O chromosome has been subject to much research. Lethal gene allelisms between American populations far away have been studied. These allelisms were not negligible, but all cases were due to the lethal gene completely associated to the O5 chromosomal inversion. Here we analyze the lethal genes in a new American population of D. subobscura (Centralia, Washington), located fairly close to a previously studied population (Bellingham, Washington) and separated in space and time with other American populations (Gilroy I and II in California and Santiago de Chile). The frequencies of lethal and semilethal genes of Centralia were 16.9+/-4.6 and 6.2+/-3.0, respectively. The intrapopulational allelism of Centralia was 0.122+/-0.036. Interpopulational allelisms were studied using the lethal genes from the populations separated in space and time from Centralia. The interpopulational allelisms between Centralia and Gilroy I (California) and between Centralia and Bellingham (Washington) were higher than the intrapopulational allelism (0.155+/-0.032 and 0.153+/-0.024, respectively). In all these cases allelism was due to a complete association between a lethal gene and the O5 chromosomal inversion. Accordingly, no other lethal genes are shared in these populations.     

The ruddy duck Oxyura jamaicensis in Europe: natural colonization or human introduction?

Native to North America, ruddy ducks Oxyura jamaicensis now occur in 21 countries in the western Palaearctic (including Iceland) and their expanding population threatens the native white-headed duck, Oxyura leucocephala, through hybridization and possibly competition for food and nest sites. We used mitochondrial DNA sequences and nuclear microsatellites to test whether the European ruddy duck population is descended solely from the captive population in the UK, which traces to seven individuals imported from the USA in 1948, or, alternatively, has been augmented by natural dispersal of birds from North America. Limited genetic diversity in the European population is consistent with a founder population as small as seven birds. In addition, shifts in allele frequencies at several loci, presumably due to genetic drift in the founding population, result in significant differentiation between the European and North American populations. Despite the recent separation of these populations, almost all individuals could be unambiguously assigned based on their composite genotypes, to one of two distinct populations, one comprising all of the European ruddy ducks we sampled (including those from Iceland and captive birds in the UK) and the other comprising all North American samples. Our results confirm that the European ruddy duck population is likely to derive solely from the captive population in the UK and we find no evidence of recent arrivals from North America or of admixture between ruddy ducks from Europe and North America.     

Colonization of America by Drosophila subobscura: association between Odh gene haplotypes, lethal genes and chromosomal arrangements

The colonization of America by Drosophila subobscura has been a unique exper iment in nature that has allowed us to explore the effects of evolution on a continental scale. To analyze this evolutionary event, nucleotide sequences of the Odh (Octanol dehydrogenase) gene were obtained for 43 lethal chromosomal lines from colonizing populations of North America and 5 from South America, in addition to 5 chromosomal lines from Europe with different viabilities and 2 from laboratory marker stocks. Since 10 different Odh haplotypes were found in America, the minimum number of colonizers would be 5 (or 3 mated females). Only one Odh haplotype was found in American O(5) inversions confirming that only one copy of this inversion was included among the sample of colonizers. The same Odh haplotypes were detected in association with the same chromosomal arrangements and with identical lethal genes in both North and South America indicating that exactly the same chromosome types reached both hemispheres. These observations indicate that the two continental colonizations are not independent. They are derived from the same colonization event. The population from which the colonization started should contain the O(5) inversion, a non-negligible frequency of the O(3+4+7) arrangement and all other arrangements found in America. So far the only populations that fulfill all these requirements are those from Greece, indicating that these populations can be considered good candidates as a starting point for an in depth analysis of the origin of the American colonization by D. subobscura.     

Spatial and temporal stability of mtDNA haplogroup frequencies in native North America

Mitochondrial DNA lineage frequencies in prehistoric Aleut, eastern Utah Fremont, Southwestern Anasazi, Pyramid Lake, and Stillwater Marsh skeletal samples from northwest Nevada and the Oneota of western Illinois are compared with those in 41 contemporary aboriginal populations of North America. The ancient samples range in age from 300 years to over 6,000 years. The results indicate that the prehistoric inhabitants of North America exhibit the same level of mtDNA variability as contemporary populations of the continent. Variation in modern mtDNA haplogroup frequencies is highly geographically structured, and the prehistoric samples exhibit the same geographic pattern of variation. This indicates that differentiation of regional patterns of mtDNA lineage variation occurred early in North American prehistory (much more than 2,000 years B.P.), has remained relatively stable since its origin, and was little influenced by the disruptions hypothesized for other genetic systems as a result of population declines and relocations at contact.     

Founding population size of an aquatic invasive species

Non-native species of fish threaten native fishes throughout North America, and in the Rocky Mountains, introduced populations of lake trout threaten native populations of bull trout. Effective management of lake trout and other exotic species require understanding the dynamics of invasion in order to either suppress non-native populations or to prevent their spread. In this study, we used microsatellite genetic data to estimate the number of lake trout that invaded a population of bull trout in Swan Lake, MT. Examination of genetic diversity and allele frequencies within the Swan Lake populations showed that most of the genes in the lake trout population are descended from two founders. This emphasizes the importance of preventing even a few lake trout from colonizing new territory.     

O-chromosome lethal frequencies in Serbian and Montenegrin Drosophila subobscura populations

Lethal chromosomal frequencies were obtained from three Drosophila subobscura samples from the Mt. Avala (Serbia) population in September 2003 (0.218), June 2004 (0.204) and September 2004 (0.250). These values and those from other Balkan populations studied previously (Petnica, Kamariste, Zanjic and Djerdap) were used to analyze the possible effect of population, year, month and altitude above sea level on lethal chromosomal frequencies. According to ANOVAS no effect were observed. Furthermore, the lethal frequencies of the Balkan populations did not vary according to latitude. This is probably due to the relative proximity and high gene flow between these populations. From a joint study of all the Palearctic D. subobscura populations so far analyzed, it can be deduced that the Balkan populations are located in the central area of the species distribution. Finally, it seems that lethal chromosomal frequencies are a consequence of the genetic structure of the populations.     

Patterns of transposable element variation and clinality in Drosophila

Natural populations often exist in spatially diverse environments and may experience variation in the strength and targets of natural selection across their ranges. Drosophila provides an excellent opportunity to study the effects of spatially varying selection in natural populations, as both Drosophila melanogaster and Drosophila simulans live across a wide range of environments in North America. Here, we characterize patterns of variation in transposable elements (TEs) from six populations of D. melanogaster and nine populations of D. simulans sampled from multiple latitudes across North America. We find a nearly twofold excess of TEs in D. melanogaster relative to D. simulans, with this difference largely driven by TEs segregating at the lowest and highest allele frequencies. We find no effect of latitude on either total TE abundance or average TE allele frequencies in either species. Moreover, we show that, as a class of mutations, the most common patterns of TE variation do not coincide with the sampled latitudinal gradient, nor are they consistent with local adaptation acting on environmental differences found in the most extreme latitudes. We also do not find a cline in ancestry for North American D. melanogaster—for either TEs or single nucleotide polymorphisms—suggesting a limited role for demography in shaping patterns of TE variation. Though we find little evidence for widespread clinality among TEs in Drosophila, this does not necessarily imply a limited role for TEs in adaptation. We discuss the need for improved models of adaptation to large‐scale environmental heterogeneity, and how these might be applied to TEs.     

Diagnosis of sibling species of Drosophila involved in the colonization of North America by D. subobscura

To determine the effects of the recent colonization of the west coast of North America by the Palaearctic species Drosophila subobscura on the dynamics of the Drosophila populations, the sibling species D. athabasca and D. azteca must be classified unambiguously. We have characterized these two species using three molecular techniques: allozymes, mtDNA and RAPDs. All three techniques allow the classification of any individual as belonging to either species. The study of five localities in northern California and southern Oregon show that the area of overlap is larger than previously described.     

Phylogeographic patterns in Drosophila montana

The Drosophila virilis species group offers valuable opportunities for studying the roles of chromosomal re-arrangements and mating signals in speciation. The 13 species are divided into two subgroups, the montana and virilis 'phylads'. There is greater differentiation among species within the montana phylad in both karyotype and acoustic signals than exists among members of the virilis phylad. Drosophila montana is a divergent species which is included in the montana phylad. Here, we analyse the phylogeography of D. montana to provide a framework for understanding divergence of acoustic signals among populations. We analysed mitochondrial sequences corresponding to the cytochrome oxidase I and cytochrome oxidase II genes, as well as 16 microsatellite loci, from 108 lines of D. montana covering most of the species' range. The species shows a clear genetic differentiation between North American and Scandinavian populations. Microsatellite allele frequencies and mitochondrial DNA haplotypes gave significant FST values between populations from Canada, USA and Finland. A Bayesian analysis of population structure based on the microsatellite frequencies showed four genetically distinct groups, corresponding to these three populations plus a small sample from Japan. A network based on mitochondrial haplotypes showed two Finnish clades of very different shape and variability, and another clade with all sequences from North America and Japan. All D. montana populations showed evidence of demographic expansion but the patterns inferred by coalescent analysis differed between populations. The divergence times between Scandinavian and North American clades were estimated to range from 450,000 to 900,000 years with populations in Canada and the USA possibly representing descendants of different refugial populations. Long-term separation of D. montana populations could have provided the opportunity for differentiation observed in male signal traits, especially carrier frequency of the song, but relaxation of sexual selection during population expansion may have been necessary.     

Temperature-related genetic changes in laboratory populations of Drosophila subobscura: evidence against simple climatic-based explanations for latitudinal clines

Parallel latitudinal clines to the long-standing ones in the original Palearctic populations have independently evolved at different rates for chromosomal polymorphism and body size in South and North American populations of Drosophila subobscura since colonization around 25 years ago. This strongly suggests that (micro) evolutionary changes are largely predictable, but the underlying mechanisms are unknown. The putative role of temperature per se was investigated by using three sets of populations at each of three temperatures (13 degrees , 18 degrees , and 22 degrees C) spanning much of the tolerable range for this species. We found a lower chromosomal diversity at the warmest temperature; a quick and consistent shift in gene arrangement frequencies in response to temperature; an evolutionary decrease in wing size, mediated by both cell area and cell number, at 18 degrees C; no relationship between wing size and those inversions involved in latitudinal clines; and a shortening of the basal length of longitudinal vein IV relative to its total length with increasing standard dose. The trends for chromosomal polymorphism and body size were generally inconsistent from simple climatic-based explanations of worldwide latitudinal patterns. The findings are discussed in the light of available information on D. subobscura and results from earlier thermal selection experiments with various Drosophila species.     

Microsatellite variation within and among North American lineages of Phragmites australis

Over the past century, the spread of the common reed (Phragmites australis) has had a dramatic impact on wetland communities across North America. Although native populations of Phragmites persist, introduced invasive populations have dominated many sites and it is not clear if the two types can interbreed. This study compares patterns of differentiation in 10 microsatellite loci among North American and European Phragmites individuals with results obtained from sequencing of noncoding chloroplast DNA. Three population lineages (native, introduced and Gulf Coast) were previously identified in North America from chloroplast DNA and similar structuring was found in the nuclear genome. Each lineage was distinguished by unique alleles and allele combinations and the introduced lineage was closely related to its hypothesized source population in Europe. Size homoplasy and diagnostic base substitutions distinguishing lineages were evident at several loci, further emphasizing that native, introduced and Gulf Coast North American Phragmites lineages are genetically distinct. Gene flow between lineages was low and invasive introduced populations do not represent a hybrid population type.     

An experimental study of evolution in progress: Clines for quantitative traits in colonizing and Palearctic populations ofDrosophila

Summary Drosophila subobscura has recently colonized the American continent and is an excellent model for studying evolution in action. Previous analyses have shown that these colonizing populations have significant latitudinal clines for the frequencies of some chromosomal arrangements that parallel those clines found in the Old World. These results strongly suggest that this polymorphism is adaptive. In the present study, significant latitudinal clines for continuous morphometric variables (flies are larger in the north) have been detected in Old World populations ofD. subobscura. The adaptive nature of these clines is reinforced by the fact that parallel latitudinal clines for body size have also been detected inDrosophila obscura, a closely related sympatric species, as well as previously in otherDrosophila. On the other hand, no significant latitudinal clines for continuous morphometric traits, not even when using an overall size index, have been detected in colonizing populations ofD. subobscura. This is a rather surprising result given the number of generations that have elapsed since the species was detected in America and given that significant clines in chromosomal inversions are already established. Thus, the adaptive response of quantitative morphometric variables is not as rapid as that found for chromosomal inversions. Nevertheless, canonical correlation analysis suggests that significant latitudinal clines for body size might soon be detected in the American continent. The results obtained here are included in a projected time series with the aim of documenting size evolution in action.     

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.