Climate change and recent genetic flux in populations of Drosophila robusta
BMC Evolutionary Biology volume 5, Article number: 4 (2005)
Studied since the early 1940's, chromosomal polymorphisms in the deciduous woods species Drosophila robusta have been characterized by well-defined latitudinal, longitudinal, and elevational clines, but – until at least ten years ago – stable, local population frequencies. Recent biogeographical analyses indicate that D. robusta invaded North America from southeast Asia and has persisted in eastern temperate forests for at least 20–25 my without speciating. The abundant chromosome polymorphisms found across the range of D. robusta are thus likely to be relatively ancient, having accumulated over many well known climatic cycles in North America. Sufficient long-term data are now available such that we can now gauge the rate of these evolutionary changes in natural populations due to environmental change.
Recent local collections have revealed significant changes in the frequencies of several chromosomal forms. New data presented here extend the range of these changes to six states, three in the northeastern United States and three west of the Mississippi River. These data reinforce recent directional changes in which the frequencies of three gene arrangements have reached percentage levels typical of distant southern populations consistent with regional climatic changes. Another gene arrangement has been steadily decreasing in frequency at a number of the sites studied. Meteorological records from 1945 to 2003 indicate temperature increases at all study sites, particularly average minimum air temperatures.
Observation of parallel genetic flux suggests that these long-term temporal frequency shifts in widely disparate populations of D. robusta are evolutionary responses to environmental change. Since these chromosomes are known to be sensitive to ambient temperature, regional climatic shifts associated with global warming are likely to be responsible.
In recent years, numerous publications – over 50 pages of them in current listings of the bibliographic search engine Medline – have detailed changes in biological systems and organisms that appear associated with the climatic changes commonly referred to as global warming (e.g., [1–5], dealing with observed or impending shifts in the habitats of various organisms). Few studies, however, have documented long-term changes in the genetic structure of species populations on a regional scale, necessary for understanding the microevolutionary consequences of global change . In the genus Drosophila, comprised of over 1500 species , only four species have received such attention. While some long term increases have been documented in the frequencies of chromosomal gene arrangements Standard (ST), Tree Line (TL), and Pike's Peak (PP) in Drosophila pseudoobscura in the western United States and Canada over ca 40 years [8–10] and of several arrangements in Drosophila melanogaster in Japan , little is known of the causes for these changes.
Orengo and Prevosti  first suggested that climatic warming was responsible for long-term temporal changes in the frequencies of certain gene arrangements in European populations of D. subobscura. Later, Rodriguez-Trelles and co-workers [13–15] and Sole', et al.  convincingly demonstrated that climatic changes, including global warming, were likely the driving force of microevolutionary changes in these populations.
Recently, significant changes have also been documented in the chromosomal variation of D. robusta Sturtevant, some of them possibly attributable to global warming [17–19]. Here we describe additional data that underscore the variety of historical changes being experienced by populations of this species. Clearly, the chromosomal polymorphisms in D. robusta are also dynamic, and when compared with the considerable geographical and experimental data available for this species collected over the last 60 years, strongly implicate regional climatic changes as a cause for these temporal frequency shifts.
Regional patterns of climate change were revealed from ANCOVA analysis of temperature and precipitation data from 1945 – 2003. No long term tends were detected for precipitation except for a significant increase in Central Park, NY (data available from the authors). Significant long-term temperature changes were apparent at all six 2003 study sites (Fig. 1; plus Philadelphia, last studied in 2002) for the three available temperature indicators: average monthly minimum temperature (MINTMP), average monthly temperature (AVETMP), and average maximum monthly temperature (MAXTMP; all 3 ANCOVA models, P < 0.0001). AVTEMP and MAXTMP varied significantly from site to site in different years, but MINTMP significantly increased over the 58-year period (P < 0.0001) with no heterogeneity among sites. Thus, the most consistent change in climate across sites in this study was due to significant increases in minimum monthly temperatures.
Evidence for site-specific temperature increases from 1945 – 2003 was evident for these three temperature indicators (Fig. 2). Five of the seven sites showed linear or curvilinear increases in MINTMP, usually starting in the 1970's and extending until 2000. In contrast, both Fayetteville and Englewood, NJ have experienced significant temporal decreases in MAXTMP (with a concurrent significant increase in MINTMP in Fayetteville). Causes for these decreases are obscure, but in New Jersey may be due to the exaggerated temperature fluctuations in the early 1990's (as well as missing data for 1992 and 1993, Fig. 2).
Concurrent with these temperature shifts, chromosome elements in populations of D. robusta show systematic temporal frequency changes (all frequency data are available from the authors). X chromosome combinations 1S and S1 show substantial frequency changes with time, with S1 doing so in a site-specific manner (Fig. 3). Frequencies of combination 1S have decreased over time, and are negatively correlated with increasing temperatures (Tables 1, 2). Increases in S1 are only marginally correlated with temperature, but this may be due to other factors, including its low frequency in eastern populations (Fig. 3). Gene arrangements 2L-1 and 3R-1 have also increased in frequency since 1945 (Fig. 4, 5; Tables 1, 2), and the changes are significantly site-specific as indicated by the significant year × site interactions. While temporal changes in 2L-1 and 3R-1 are correlated (r = 0.419, P = 0.0016, n = 54), only those of 3R-1 are significantly correlated with increasing temperatures (Fig. 5; Table 2).
The most consistent historical changes have been gains in the frequency of arrangement 3R-1 with concomitant decreases of its allelic form, 3R (Figure 5). Coupled with the clinal tendency of 3R-1, varying between 100% in the southernmost latitudes and zero in the northernmost , its similar directional increases of recent years in at least four states, some far apart, and concordance with long term temperature increases indicates that 3R-1 is responding to regional climatic shifts.
The frequency of arrangement 2L-1 increased directionally in most of the same localities as 3R-1 (Figure 4). The frequency of 2L-2, which is also more common in the south than in the north, has been decreasing steadily, almost to the point of extinction, in the four northern 2003 populations. This is probably not a temperature effect, however, as in the past 2L-2 has been irregularly distributed across the southern states, rarely rising above frequencies of 20 percent, without evidence of a latitudinal cline [20–22].
Like 2L-1 and 3R-1, the frequency of gene arrangement XL (the sum of data for SS, S1, and S2), tends to be higher in southern populations than in northern ones . Except for a few deviations in small samples, it has increased steadily – with concomitant decreases in northern arrangement XL-1 – in all the localities studied except Olivette (where it is already about 95%) and Iowa. At Englewood, for example, it rose from about 40% in 1948 to an average of 78% in 2000–2003, that of northern arrangement XL-1 falling from about 60% to about 22% in the same period. Frequencies of X chromosome combination 1S (XL-1.XR) are clearly consistent with this pattern (Tables 1, 2).
Studies of Drosophila inversion polymorphisms have now provided historical insights into environmental change. On different continents, temporal genetic changes correlated with increasing temperatures in populations of D. subobscura [12–16] and D. robusta [[17, 18], this study] bear the imprint of large scale climate change. For D. robusta, not all temporal changes were consistent with this hypothesis. In the course of comparing plant growth in rural and urban areas in relation to ozone exposure, rural areas tended to have lower average temperatures than urban ones ; this may explain the smaller (or lack of) change in 3R-1 frequencies in relatively rural Fayetteville and Iowa City as compared to the more urban locales. Also, there are fewer data available for these sites collected over time spans long enough to document temporal trends. Another possibility is that the climatic changes have not been of equal intensity in all localities, as indicated by regional increases and decreases in local temperatures (Fig. 2) suggested by the site by year interactions for AVTMP and MAXTMP in the ANCOVA analyses.
Carson  noted that the chromosomal variation in populations of Drosophila robusta at Olivette, Missouri over a 10-year period was characterized by "extraordinary stability," wherein "certain frequencies may shift significantly as compared with the previous year, but in every case the observed frequencies approximate some previously observed level." This pattern continued until at least 1967 (Fig. 3, 4, 5). Temperature records from this period seem rather stable before significantly increasing through 2000 (AVTMP and MINTMP second order regression slopes are both positive; Fig 2). Similar stability, frequency changes of less than 10 percent, has also been recorded near Blacksburg, Virginia from 1950 to 1962 (M. Levitan, unpublished data) and in northeastern New Jersey from 1948 until at least 1975 (Fig. 3, 4, 5).
By contrast, every population sampled in 2003 evidenced at least one significant change of chromosomal polymorphism frequency compared to the numbers in the same area ten or more years previously, with indication that the changes in at least three are part of similar, if not identical, directional historical processes. Therefore, documentation of temporal genetic changes in these populations requires at least 10 to 20 years of comparative data.
Influences of natural selection due to ambient temperature variation on frequencies of these gene arrangements and X chromosome associations have been demonstrated in laboratory experiments [24, 25], and inferences from latitudinal and multiple elevational clines [26–31]. Frequencies of arrangements XL-1 and 2L-3 are clearly associated with cooler temperatures at higher elevations and latitudes, and 2L-1 increases in frequency in the laboratory under warmer temperatures. Carriers of 2L-3 have shorter egg to adult development times expressed in cooler temperatures, explaining increases in this gene arrangement with elevation and latitude .
Evidence of strong natural selection maintaining earlier observed genetic stability has also come from perturbation experiments in the wild . At four time periods, large numbers of flies carrying X-chromosome combinations S2 and 22, autosomal combination 2L-1.2R-1 (11), and 3R-1 from South Carolina, Alabama, and Mississippi were released in midsummer in the Englewood, NJ woods. In other years the released flies carried XL-1, 2L-3, and 3R from Minnesota and Michigan. Despite evidence of hybridization of the introduced flies with the local population, in each case the following spring saw return to frequencies of the previous year.
According to Wright , "The alternative to natural selection in a changing environment is, as noted by Dobzhansky , the emergence of superior genetic systems" in the gene arrangements that have been increasing in frequency. He envisioned this happening in one population by recombination in inversion homozygotes, with the new adaptive gene complex in one or a few flies spreading to other localities by "occasional very long dispersion by the wind." He conceded that such a hypothesis would be most satisfactory for cases where a very rare inversion suddenly started increasing, such as the rises of TL and PP of D. pseudoobscura mentioned above, since high frequencies of inversion homozygotes would undermine the structural stability of the new adaptive complex. The sudden rises of D. robusta X-chromosome combination S1 in New York and New Jersey and of S2 in Missouri could be cases in point, but it would not explain the near doubling of S1 in Arkansas from a 35% base between 2001 and 2003 nor the three most consistent directional changes, those of 3R-1, 2L-1, and XL (Figures 3, 4, 5). Indeed, the nearly simultaneous changes in places as far apart as New York and Allentown, Pennsylvania, let alone New York and Missouri, would be very unlikely to depend on chance wind-blown dispersions.
Here, causes for systematic frequency changes in widespread populations can be inferred to be a result of common environmental causes given the large amount of background data available for D. robusta. When compared to the inversion frequency data collected over 60 years from more than 150 natural populations [20, 21] and results from numerous laboratory experiments, temporal frequency shifts across the range of D. robusta strongly suggest temperature variation as a likely mechanism driving microevolutionary change [17, 18]. Such regional frequency shifts suggest a common response to temperature fluctuations or causes correlated with them, whether due to regional climatic changes or global warming.
Although direct evidence for climatic change has been accumulating for many years, its consequences for causing evolutionary changes have only recently been observed. Chromosome polymorphisms in Drosophila species have been historically important genetic systems for understanding mechanisms of evolutionary change, and have now been studied long enough to begin revealing widespread, systematic temporal frequency shifts in response to environmental change. These polymorphisms thus represent excellent indicators of future climatic shifts.
There are 14 commonly encountered gene arrangements segregating in natural populations of D. robusta located on five of the six arms of the 3 major chromosomes. The "Standard" arrangements were labeled for the respective chromosome arms: XL, XR, 2L, etc. Others were named in order of their discovery, e.g., XL-1, XL-2, XR-1, 2L-1 [described and configured in [20, 24]].
The Standard arrangement of each arm was dubbed "S," and the other arrangements are referred to by the Arabic numerals in their names . A fly with karyotype XL/XL-1, XR/XR-2, for example, would be S/1, S/2 in this notation. Depending on the linkage combination of the arrangements, it is also either SS/12 or S2/1S.
Linkage relationships are inferred from karyotypic analyses of adult males and females . Female D. robusta, unlike many other drosophilids, quickly deplete stored sperm in the absence of remating so that wild-caught females can be despermed by repeated transfers to fresh food vials and then crossed to homokaryotypic males in a controlled fashion. Karyotypes of at least 6 larvae from these test crosses were prepared in order to infer the linkage combination of X chromosome gene arrangements. Salivary gland smears from larvae derived from matings in the wild, so-called "egg sample" data, were included when collected females did not survive the desperming transfers.
This report compares data obtained in 2003 to earlier summer collections at a number of geographically isolated populations going back in some cases to 1946: Olivette, a suburb of St. Louis, Missouri; woods alongside Route 4 in Englewood, New Jersey; Trexler Memorial Park on the western outskirts of Allentown, Pennsylvania; the North Woods of Central Park in New York City; Fayetteville, Arkansas; and woods along the Iowa River at Iowa City, Iowa.
The significance of year-to-year differences in frequency of each chromosome or chromosome arm obtained prior to 2003 was determined by G-tests . Chronologically contiguous results that proved statistically homogeneous are combined in the figures. Individual X- and second-chromosome arrangements could not be tested in this way due to problems of independence.
Analysis of meteorological data
Temperature and precipitation data were obtained online from the National Climatic Data Center http://www.ncdc.noaa.gov/oa/ncdc.html for stations with complete records extending from 1945–2003 nearest to each site. Continuous data for Fayetteville and Englewood were not available for the closest weather stations, so records from those stations with complete records were used. Correlations between temperatures from these stations were highly significant: Drake Field and the Agricultural Experiment Station for Fayetteville (r = 0.48 – 0.90, P < 0.01), Little Falls and Ridgefield for Englewood (r = 0.84 – 0.98, P < 0.0001); see Fig. 2. Analysis of covariance in PROC GLM  was used to test for overall significance of temperature trends and chromosome frequency changes across sites, and to evaluate regional patterns of temperature and precipitation change. Polynomial regression analyses of temperature and precipitation data with time were performed with PROC REG . In all cases, linear or second order polynomial regression explained the most variation for a particular model. Pearson product-moment correlations between arcsin transformed chromosome frequencies and the average temperature of the month for each collection were calculated with PROC CORR . Non-parametric correlation analyses produced equivalent results. We also assessed correlations with temperatures for the month prior to collection, and the average temperature of the 3 months prior to collection: none were significant.
Parmesan C, Yohe G: A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003, 421: 37-42. 10.1038/nature01286.
Parmesan C, Galbraith H: Observed Impacts of Global Climate Change in the U.S. Pew Center on Global Climate Change, Arlington, Virginia. 2004, 67-
Root TL, Price JT, Hall KR, Schneider SH, Rosenzweig C, Pounds JA: Fingerprints of global warming on wild animals and plants. Nature. 2003, 421: 57-60. 10.1038/nature01333.
Williams SE, Botho EE, Fox S: Climate change in Australian tropical rainforests: an impending environmental catastrophe. Proc Roy Soc London, B (Biol Sci). 2003, 270: 1887-1892. 10.1098/rspb.2003.2464.
Cotton PA: Avian migration phenology and global climate change. Proc Natl Acad Sci USA. 2003, 100: 12219-12222. 10.1073/pnas.1930548100.
Rice KJ, Emery NC: Managing microevolution: restoration in the face of global change. Front Ecol Environ. 2003, 1: 469-478.
Wheeler MR: Additions to the catalog of the world's Drosophilidae. In: The Genetics and Biology of Drosophila. Edited by: Ashburner M, Carson HL, Thompson JN Jr. 1986, London: Academic Press, 3e: 395-409.
Dobzhansky T, Anderson WW, Pavlovsky O: Genetics of natural populations. XXXVIII. Continuity and change in populations of Drosophila pseudoobscura in western United States. Evol. 1966, 20: 418-427.
Anderson WW, Arnold J, Baldwin DG, Beckenbach AT, Brown CJ, Bryant SH, Coyne JA, Harshman LG, Heed WB, Jeffrey DE, Klaczko LB, Moore BC, Porter JM, Powell JR, Prout T, Schaeffer SW, Stephens JC, Taylor CE, Turner ME, Williams GO, Moore JA: Four decades of inversion polymorphism in Drosophila pseudoobscura. Proc Natl Acad Sci USA. 1991, 88: 10367-10371.
Powell JR: Inversion polymorphism in Drosophila pseudoobscura and Drosophila persimilis. In Drosophila Inversion Polymorphism. Edited by: Krimbas CB, Powell JR. 1992, Boca Raton, FL: CRC Press, 73-126.
Inoue Y, Watanabe T, Watanabe TK: Evolutionary change of the chromosomal polymorphism in Drosophila melanogaster populations. Gene. 1984, 38: 753-765.
Orengo D-J, Prevosti A: Temporal changes in chromosomal polymorphism of Drosophila subobscura related to climatic changes. Evol. 1996, 50: 1346-1350.
Rodriguez-Trelles F, Alvarez G, Zapata C: Time-series analysis of seasonal changes of the O inversion polymorphism of Drosophila subobscura. Genet. 1996, 142: 179-187.
Rodriguez-Trelles F, Rodríguez M: Rapid micro-evolution and loss of chromosomal diversity in Drosophila in response to climate warming. Evol Ecol. 1998, 12: 829-838. 10.1023/A:1006546616462.
Rodriguez-Trelles F, Rodríguez M, Scheiner S: Tracking the genetic effects of global warming: Drosophila and other model systems. Cons Ecol. 1998, 2: 2-[http://www.ecologyandsociety.org/vol2/iss2/art2/]
Sole' E, Balanya J, Sperlich D, Serra L: Long-term changes in the chromosomal inversion polymorphism of Drosophila subobscura. I. Mediterranean populations from southwestern Europe. Evol. 2002, 56: 830-835.
Levitan M: Studies of linkage in populations. XIV. Historical changes in frequencies of gene arrangements and arrangement combinations in natural populations of Drosophila robusta. Evol. 2001, 55: 2359-2362.
Levitan M: Climatic factors and increased frequencies of 'southern' chromosome forms in natural populations of Drosophila robusta. Evol Ecol Res. 2003, 55: 597-604.
Schilthuizen M: Fly genome warms to climate change. 2003, [http://sciencenow.sciencemag.org/cgi/content/full/2003/521/2]
Levitan M: Chromosomal variation in Drosophila robusta Sturtevant. In Drosophila Inversion Polymorphism. Edited by: Krimbas CB, Powell JR. 1992, Boca Raton, FL: CRC Press, 221-238.
Etges WJ, Levitan M: Paleoclimatic variation, adaptation and biogeography of inversion polymorphisms in natural populations of Drosophila robusta. Biol J Linn Soc. 2004, 51: 395-411. 10.1111/j.1095-8312.2004.00306.x.
Carson HL: Genetic conditions which promote or retard the formation of species. Cold Spring Harbor Symp Quant Biol. 1959, 24: 87-105.
Gregg JW, Jones CG, Dawson TE: Urbanization effects on tree growth in the vicinity of New York City. Nat. 2003, 424: 183-187. 10.1038/nature01728.
Carson HL: The population genetics of Drosophila robusta. Adv Genet. 1958, 9: 1-40.
Levitan M: Experiments on chromosomal variability in Drosophila robusta. Genet. 1951, 36: 285-305.
Etges WJ: Chromosomal influences on life history variation along an altitudinal transect in Drosophila robusta. Am Nat. 1989, 133: 83-110. 10.1086/284903.
Stalker HD, Carson HL: Morphological variation in natural populations of Drosophila robusta. Evol. 1947, 1: 237-248.
Stalker HD, Carson HL: An altitudinal transect of Drosophila robusta Sturtevant. Evol. 1948, 2: 295-305.
Etges WJ: Genetic structure and change in natural populations of Drosophila robusta : Systematic inversion and inversion association frequency shifts in the Great Smoky Mountains. Evol. 1984, 38: 675-688.
Levitan M: Studies of linkage in populations. IX. The effects of altitude on X-chromosome arrangement combinations in Drosophila robusta. Genet. 1978, 89: 751-763.
Levitan M, Scheffer SJ: Studies of linkage in populations. X. Altitude and autosomal gene arrangements in Drosophila robusta. Genet Res Cambr. 1993, 61: 9-20.
Wright S: Evolution and the Genetics of Populations. Variability Within and Among Natural Populations. 1978, Chicago and London: University of Chicago Press, 4: 118-123.
Dobzhansky T: Genetics of the Evolutionary Process. 1970, New York: Columbia University Press, 274-276.
Carson HL: The effect of inversions on crossing over in Drosophila robusta. Gen. 1953, 38: 168-186.
Levitan M: Studies of linkage in populations. I. Associations of second chromosome inversions in Drosophila robusta. Evol. 1955, 9: 62-74.
Sokal RR, Rohlf FJ: Biometry: The Principles and Practice of Statistics in Biological Research. 1995, San Francisco: Freeman, 3
SAS Institute: SAS/STAT User's Guide, Version 6. 1989, SAS Institute Inc., Cary, NC
The authors thank the 1996 M. M. Kaplan Family Foundation for grant GCO 00-1223 to ML and NSF grant DEB-0211125 to WJE; A. R. Templeton for the 2002 and 2003 collections in Missouri; B. McAllister for the 2003 collections in Iowa; M. M. Kaplan for advice and comments on the manuscript; S. Beaupré for statistical advice; N. Maj for computer and graphics assistance; and S. Hanada, N. Khanna, S. Gmuca, and K. Chan, sponsored by the New York Academy of Sciences, for assistance in the analysis of the 2003 collections.
ML collected the Allentown, New York, and New Jersey flies and carried out the chromosomal analyses. WJE contributed the Fayetteville flies. Both authors contributed to planning the study, analyzing the data, and writing and revising the manuscript.
About this article
Cite this article
Levitan, M., Etges, W.J. Climate change and recent genetic flux in populations of Drosophila robusta. BMC Evol Biol 5, 4 (2005). https://doi.org/10.1186/1471-2148-5-4
- Regional Climatic Change
- Gene Arrangement
- Microevolutionary Change
- Average Maximum Monthly Temperature
- Average Monthly Minimum Temperature