Molecular evidence to suggest the origin of a colonization: Drosophila subobscura in America

The recent colonization of America by Drosophila subobscura represents a great opportunity for evolutionary biology studies. Knowledge of the populations from which the colonization started would provide an understanding of how genetic composition changed during adaptation to the new environment. Thus, a 793 nucleotide fragment of the Odh (Octanol dehydrogenase) gene was sequenced in 66 chromosomal lines from Barcelona (western Mediterranean) and in 66 from Mt. Parnes (Greece, eastern Mediterranean). No sequence of Odh fragment in Barcelona or Mt. Parnes was identical to any of those previously detected in America. However, an Odh sequence from Barcelona differed in only one nucleotide from another found in American populations. In both cases, the chromosomal lines presented the same inversion: O(7), and the Odh gene was located within this inversion. This evidence suggests a possible western Mediterranean origin for the colonization. Finally, the molecular and inversion data indicate that the colonization was not characterized by multiple reintroductions.     

Colonization of the Americas by Drosophila Subobscura: Lethal-Gene Allelism and Association with Chromosomal Arrangements

Drosophila subobscura is a Palearctic species that has recently colonized the Americas. It was first found in 1978 in Puerto Montt, Chile, and in 1982 in Port Townsend, WA. The colonization and rapid expansion of the species in western South and North America provides distinctive opportunities for investigating the process of evolution in action. The inversion polymorphism in the O chromosome from populations of central California and northern Washington, separated by 1300 km, corresponds to a previously observed latitudinal cline, also observed in Europe. Recessive lethal genes are not randomly distributed among the chromosomal arrangements. The incidence of lethal allelism is high, yielding unrealistically low estimates of the effective size of these populations (on the order of 1000 individuals). The high incidence of lethal allelism is likely to be a consequence of the low number of the American colonizers (on the order of 10-100 individuals), but the persistence of the allelism over several years suggests that some lethal-carrying chromosomes may be heterotic owing to shared associations between lethal and other genes.     

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.     

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.     

COLONIZATION OF AMERICA BY DROSOPHILA SUBOBSCURA: LETHAL GENES AND THE PROBLEM OF THE O5 INVERSION

In this work, the process of colonization of North and South America by the species Drosophila subobscura has been studied by analyzing the variability of lethal genes. The genetic structures of a Palearctic natural central population (Bordils, Spain) and a colonizer population from America (Gilroy, California) have been compared. The frequencies of lethal chromosomes and their allelism are 29.007% and 0.0069 in the first population and 14.414% and 0.0526 in the American population. A founder effect is detected after the computation of some population parameters (N_e, h, he and the lethal load). Furthermore, the allelism of lethal chromosomes has revealed a strong association between a lethal gene and the O₅ inversion both in Gilroy and in the population of Puerto Montt (Chile). The interpopulation allelism shows that the O₅ arrangement from the USA and Chile is the same, confirming that the colonizing processes of North and South America are correlated. The O₅ arrangement can also be useful as a genetic marker to trace the origin of the colonization. The frequency of the O₅ arrangement in the original population of the colonization could be used to estimate the number of colonizers. This population is still unknown, but taking the extreme values of the frequency of the O₅ inversion in natural Palearctic populations (1–15%), the number of colonizers could vary between 9 and 149 individuals.     

Lethal genes in O5, chromosomes of Drosophila subobscura from Europe and America

In populations of D. subobscura , a species that is know for its high chromosomal polymorphism, the O₅ inversion has a rather erratic frequency distribution in the Palearctic region. An O₅ lethal chromosomal line obtained from a Balkan population near Zanjic (South Adriatic, Montenegro, Yugoslavia) was tested for lethal allelism with other O₅ lethal chromosomal lines derived from American (USA and Chile) colonizing populations, and from the French population of Taulé. No allelism was found between the Balkan lethal gene and those from America and France. Thus, the lethal genes of the O₅ inversions are not of the same origin and it is most probable that the American colonizations did not start from the Zanjic population. The general difference in the chromosomal inversion polymorphism corroborates this conclusion. The cytological analysis confirms the assumption that all O₅ chromosomes studied are identical with respect to breakage points.     

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.     

Tracking the origin of an invasive species: Drosophila subobscura in Argentina

Biological invasions are excellent opportunities to study the evolutionary forces leading to the adaptation of a species to a new habitat. Knowledge of the introduction history of colonizing species helps tracking colonizing routes and assists in defining management strategies for invasive species. The Palearctic species Drosophila subobscura is a good model organism for tracking colonizations since it was detected in Chile and western North America three decades ago and later on in the Atlantic coast of Argentina. To unravel the origin of the Argentinean colonizers two populations have been analysed with several genetic markers. Chromosomal arrangements and microsatellite alleles found in Argentina are almost similar to those observed in Chile and USA. The lethal allelism test demonstrates that the lethal gene associated with the O₅ inversions in Argentina is identical to that found in Chile and USA, strongly supporting the hypothesis that all the American colonizing populations originated from the same colonization event. A secondary bottleneck is detected in the Argentinean populations and the genetic markers suggest that these populations originated from the invasion of 80–150 founding individuals from Chile.     

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.     

Seasonal changes of chromosomal inversion polymorphism in a Drosophila subobscura natural population from southeastern European continental refugium of the last glaciation period

Chromosomal inversion polymorphism was studied in a natural population of Drosophila subobscura from Djerdap gorge of Danube river, that represents one of southeastern European continental refugiums during last glaciation period. The existence for changes of chromosomal inversion polymorphism during different periods of time was tested. Observed very high abundance of chromosomal arrangements and lower IFR value grouped this population to central European populations. It was characterized with O7 and O(3+4+5) arrangements, not found in other populations from the region analyzed. Except in this population from warmer habitat, observed E(1+2+9+12) chromosomal arrangement was found, in the region analyzed, only in one population from Mediterranean refugial area. Many similarities have been observed between population from Djerdap gorge and other populations from the region analyzed, indicating that process of post-glacial migration may have started from this southeastern European refugium area. Seasonal changes in the frequencies for Est and O(3+4+7) chromosomal arrangements where Est was more frequent during the spring period and O(3+4+7) during the summer period were found. The changes for J-chromosome (Jst/Jst), E-chromosome (Est/E(1+2+9); E(8)/E(1+2+9) and O-chromosome (Ost/O(3+4+22); O(3+4)/O(3+4); O(3+4+1)/O(3+4+7) whose karyotypic combination frequencies were associated with seasonal changes were also found. These changes most probably represent associations with seasonal climatic changes in the population analyzed, supporting adaptivity of chromosomal inversions in D. subobscura.     

O chromosome inversion polymorphism in Northern and Atlantic Europe and its implications in the American colonization by Drosophila subobscura

The Western Mediterranean populations of D. subobscura have a very similar chromosomal inversion polymorphism to the American colonizing populations, with the exception of the O₅ inversion. This inversion, which is found in appreciable frequencies in colonizing populations and is distributed according to a significant latitudinal cline, both in North and South America, has never been reported in that Palearctic region. Therefore, an analysis of the distribution of this inversion, along a latitudinal cline, has been carried out in Europe. The O Chromosome inversion polymorphism has been studied in seven populations: Gävle and Lilla Edet (Sweden), Gesten (Denmark), Ter-Apel (Holland) and Crézy, Aizenay and Taulé (France). The O₅ inversion was only detected in three populations and in low frequencies: Gävle (3.7%), Lilla Edet (1.8 %) and Taulé (1.2 %). However, none of these three populations has all the inversions found in the American populations. The clinal distribution of some O chromosomal arrangements has also been studied according to latitude and the temporal changes of this polymorphism.     

Genetic loads and coadaptivity of chromosomal inversions I. Inversion polymorphism and genetic load on chromosome O in a Drosophila subobscura population from Petnica

Inversion polymorphism on chromosome O and polymorphism for the viability of determining genes have been studied in a natural population of Drosophila subobscura from Petnica (Serbia). The range of inversion polymorphism and the abundance of particular gene arrangements in the study population agree with a general pattern of inversion polymorphism of D. subobscura in Europe. The data obtained on the amount of genetic loads show that the D. subobscura population from Petnica displays a moderate degree of that polymorphism, compared to the other studied populations of these species. Therefore, the D. subobscura population from Petnica could be tentatively classified as an ecologically central population. Examination association of chromosomal, thus, inversion polymorphism with gene polymorphism, in the form of genetic loads show that differences exist in the mean viability among certain gene arrangements. The distribution of deleterious genes among chromosome O gene arrangements were non-random.     

PHYLOGEOGRAPHY OF BATRACHOSPERMUM MACROSPORUM (BATRACHOSPERMALES,RHODOPHYTA) FROM NORTH AND SOUTH AMERICA1

Phylogeographic trends in Batrachospermum macrosporum Mont. were investigated using the mitochondrial intergenic spacer between the cytochrome oxidase subunit 2 and 3 genes (cox2‐3). A total of 11 stream segments were sampled with seven in the coastal plain of North America and four in tropical areas of South America. Fifteen thalli were sampled from seven streams, 14 thalli from two streams, and eight thalli from two streams. There were 16 haplotypes detected using 149 individuals. Of the eight haplotypes from locations in North America, all were 334 base pairs (bp) in length, and of those from South America, five were 344 bp, and three were 348 bp. Two individual networks were produced: one for the haplotypes from North America and another for those from South America, and these could not be joined due to the large number of base pair differences. This split between haplotypes from North and South America was confirmed with sequence data of the rbcL gene. There was very little genetic variation among the haplotypes from the North American locations, leading us to hypothesize that these are fairly recent colonization events along the coastal plain. In contrast, there was high variation among haplotypes from South America, and it would appear that the Amazon serves as a center of diversity. We detected considerable variation in haplotypes among streams, but frequently, a single haplotype in an individual stream segment, which is consistent with data from previous studies of other batrachospermalean taxa, may suggest a single colonization event per stream.     

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.     

Chromosomal inversion polymorphism leads to extensive genetic structure: a multilocus survey in Drosophila subobscura

The adaptive character of inversion polymorphism in Drosophila subobscura is well established. The O(ST) and O(3+4) chromosomal arrangements of this species differ by two overlapping inversions that arose independently on O(3) chromosomes. Nucleotide variation in eight gene regions distributed along inversion O(3) was analyzed in 14 O(ST) and 14 O(3+4) lines. Levels of variation within arrangements were quite similar along the inversion. In addition, we detected (i) extensive genetic differentiation between arrangements in all regions, regardless of their distance to the inversion breakpoints; (ii) strong association between nucleotide variants and chromosomal arrangements; and (iii) high levels of linkage disequilibrium in intralocus and also in interlocus comparisons, extending over distances as great as approximately 4 Mb. These results are not consistent with the higher genetic exchange between chromosomal arrangements expected in the central part of an inversion from double-crossover events. Hence, double crossovers were not produced or, alternatively, recombinant chromosomes were eliminated by natural selection to maintain coadapted gene complexes. If the strong genetic differentiation detected along O(3) extends to other inversions, nucleotide variation would be highly structured not only in D. subobscura, but also in the genome of other species with a rich chromosomal polymorphism.     

Microsatellite variation in colonizing and palearctic populations of Drosophila subobscura

The recent colonization of North America by Drosophila subobscura has provided a great opportunity to analyze a colonization process from the beginning. A comparative study using 10 microsatellite loci was conducted for five European and two North American populations. No genetic differentiation between European populations was detected, indicating that gene flow is high among them and that the microsatellites used in the present work represent neutral markers not subject to differentiation due to selection. Extensive reduction in the number of alleles and a significant decrease in heterozygosity in colonizing populations were detected that could be explained by the founder effect and a subsequent quick but not infinite expansion. Assuming that all alleles present in the colonized area were carried by the sample of colonizers, we estimated that most probably 4-11 individuals expanded in the new area. F(ST) and the chord distance measures reflect the colonization process more accurately, since drift has been the major force in differentiating the Old and New World populations, and thus other measures considering allele size differences, such as Rho(ST) and deltamu2, are less reliable for studying nonequilibrium populations. Finally, our results were consistent with the two-phase microsatellite mutational model, indicating that most alleles are generated by gain or loss of a repeat unit, while some alleles originate by more complex mutations.     

Changes in chromosomal polymorphism and global warming: The case of Drosophila subobscura from Apatin (Serbia)

In this study, chromosomal inversion polymorphism data for a natural population of Drosophila subobscura from a swampy region near the town of Apatin (Serbia) were compared with data for the same population collected approximately 15 years earlier. The pattern of chromosomal inversion polymorphism changed over time. There were significant increases in the frequency of characteristic southern latitude ("warm" adapted) chromosomal arrangements and significant decreases in the frequency of characteristic northern latitude ("cold" adapted) chromosomal arrangements in the O and U chromosomes. The chromosomal arrangements O(3+4) and O(3+4) (+) (22) (derived from the O(3+4) arrangement) showed significant increases in 2008 and 2009 with regard to the 1994 sample. There was also a significant increase (～50%) in the U(1) (+) (2) arrangement, while U(1+8) (+) (2) (a typical southern arrangement) was detected for the first time. Since the Apatin swampy population of D. subobscura has existed for a long time in a stable habitat with high humidity that has not been changed by man our results indicate that natural selection has produced chromosomal changes in response to the increase in temperature that has occurred in the Balkan Peninsula of central southeastern European.     

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.     

Viabilities of originally natural O-chromosomal inversion homo- and heterokaryotypes in Drosophila subobscura

The influence of epistatic interactions of lethal and non-lethal genes upon viability of Drosophila inversion karyotypes is poorly known. In this paper we present comparative results for viabilities of 21 originally natural O-inversion homo- and 38 heterokaryotypes in a D. subobscura population. We observed strong heterotic effect in viability of O-lethal heterozygotes irrespective of different inversion backgrounds, which indicates a mechanism for protection of a great number of lethal genes. In O-non-lethal heterozygotes the heterotic effect in viability was absent irrespective of different inversion backgrounds. Our results showed a great extent of genetic load and high abundance of O-chromosomal arrangements in the population analyzed. It belongs to the set of central European populations. An epistasis of lethal genes present in O-inversion hetero- and not present in O-inversion homokaryotypes of moderate or low frequencies could be good example for co-adaptation of chromosomal inversions with regard to the genetic load. This represents a more efficient mechanism for limitation of genetic load than alternative mechanisms for protection of lethals. Except for lethal genes, possible epistatic interactions of mildly deleterious (subvital) genes, could also be responsible for limiting the extent of genetic load in natural D. subobscura populations. We demonstrated a non-random distribution of several combinations of viability classes among different O-inversion homo- and heterokaryotypes. As a consequence of that, the viabilities of the O-inversion homokaryotypes compared to heterokaryotypes were significantly higher at low frequencies than in moderate or high frequencies. This shows frequency-dependence as a mechanism of balancing selection for protection of O-chromosomal inversions in natural D. subobscura populations. In addition, the viabilities of the O-inversion homokaryotypes of lower frequency, compared to homokaryotypes of moderate or high frequency, were significantly higher. This again indicates the existence of supergene selection as another mechanism for protection of chromosomal inversions, as co-adapted complexes in natural D. subobscura populations.     

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

Abstract.— 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.