Variation and covariation among life-history traits in Rumex acetosella from a successional old-field gradient

In this study morphological variation and the potential for competition to affect biomass and seedling selection of the families of five populations of Rumex acetosella L. sampled along a successional old-field gradient have been investigated. Seeds from 25 families were submitted to four competitive regimes: no competition (one plant per pot), medium competition (two plants/ pot taking plants from the same population), high within-population competition (four individuals from the same population in a pot) and high between-population competition (four individuals from two different populations in a pot). Eight traits were analysed after 3 months of growth for variation among families within populations. A significant difference among families within the two older populations was recorded for sexual biomass and related components. High sensitivity of these traits to density was observed in all populations except the youngest, suggesting specialization to particular environmental conditions in late successional populations, and a good adaptive capacity to buffer environmental variation in the pioneer population. Little significant interaction between competitive regimes and families within populations was found, i.e. genotypes within each population showed little variation in their response to environmental variation. Genotypic variance decreased with increasing competitive conditions for the majority of the traits. However, the percentage of variance in sexual reproduction explained by family was stable among treatments. Tradeoffs between vegetative reproduction and sexual reproduction were recorded at the population level along the successional gradient, with increasing competitive conditions. As succession proceeds, we observed a decrease in sexual reproduction and an increase in vegetative reproduction. At the family level, correlation among traits were similar when plants were grown in the absence of competition and at high density, with a significant negative correlation between sexual reproduction and vegetative reproduction. For both sprout number and sexual biomass, the performance of families grown under all the treatments was positively correlated. Together these results indicate allocational constraints on the reproductive biology of R. acetosella that may be favoured by natural selection and have influenced population differentiation along the successional gradient. However, they also revealed that the potential exists for evolutionary specialization through plasticity, in response to variation in environmental conditions.     

Genetic covariances promote climatic adaptation in Australian Drosophila*

Evolutionary potential for adaptation hinges upon the orientation of genetic variation for traits under selection, captured by the additive genetic variance-covariance matrix (G), as well as the evolutionary stability of G. Yet studies that assess both the stability of G and its alignment with selection are extraordinarily rare. We evaluated the stability of G in three Drosophila melanogaster populations that have adapted to local climatic conditions along a latitudinal cline. We estimated population- and sex-specific G matrices for wing size and three climatic stress-resistance traits that diverge adaptively along the cline. To determine how G affects evolutionary potential within these populations, we used simulations to quantify how well G aligns with the direction of trait divergence along the cline (as a proxy for the direction of local selection) and how genetic covariances between traits and sexes influence this alignment. We found that G was stable across the cline, showing no significant divergence overall, or in sex-specific subcomponents, among populations. G also aligned well with the direction of clinal divergence, with genetic covariances strongly elevating evolutionary potential for adaptation to climatic extremes. These results suggest that genetic covariances between both traits and sexes should significantly boost evolutionary responses to environmental change.     

A chemical perspective on the evolution of variation in Eucalyptus globulus

It is becoming increasingly easy to generate genotypic data but much harder to gather an equivalent amount of phenotypic information, particularly for chemical traits. In this study of Eucalyptus globulus ssp. globulus, we measured about 60 chemical leaf traits of trees growing in a common garden to address the following questions: (1) how much variation is there between geographic regions, populations within regions and within populations? (2) How do chemical traits vary over the species’ geographic range? (3) If so, does E. globulus ssp. globulus exhibit distinct chemotypes – plants that are morphologically similar but which differ chemically? (4) Are the affinities between E. globulus ssp. globulus and closely related subspecies apparent in the chemical variation? Variation among trees within populations contributed most variation in leaf chemistry followed by variation between geographic regions. For many traits, variation among populations within proposed races and variation among proposed races within geographic regions explained little of the total variation. There was a cline in the concentration of secondary chemicals with the lowest concentrations in Tasmanian populations and the highest in those from eastern Victoria, with intermediate concentrations in populations from Bass Strait Islands. We identified three chemotypes, characterised by specific terpenes and formylated phloroglucinol compounds. The frequency of occurrence of these chemotypes showed a geographic pattern also, with “chemotype 1” predominating in Tasmania, while “chemotypes 2 and 3” occurred at highest frequency in eastern Victoria. We suggest that the chemotypes reflect introgression between E. globulus ssp. globulus and the three closely related subspecies – E. globulus ssp. bicostata, E. globulus ssp. maidenii and E. globulus ssp. pseudoglobulus. Although the formation of land-bridges with fluctuating sea levels has no doubt shaped the evolutionary history of all four subspecies, we propose that the migratory swift parrot (Lathamus discolor), an important pollinator and a species closely associated with E. globulus, has helped shape the evolution of the four tree subspecies.     

Nonclinality of molecular variation implicates selection in maintaining a morphological cline of Drosophila melanogaster

One general approach for assessing whether phenotypic variation is due to selection is to test its correlation with presumably neutral molecular variation. Neutral variation is determined by population history, the most likely alternative explanation of spatial genetic structure, whereas phenotypic variation may be influenced by the spatial pattern of selection pressure. Several methods for comparing the spatial apportionment of molecular and morphological variation have been used. Here, we present an analysis of variance framework that compares the magnitudes of latitudinal effects for molecular and morphological variation along a body size cline in Australian Drosophila populations. Explicit incorporation of the relevant environmental gradient can result in a simple and powerful test of selection. For the Australian cline, our analysis provides strong internal evidence that the cline is due to selection.     

Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia

Drosophila melanogaster invaded Australia around 100 years ago, most likely through a northern invasion. The wide range of climatic conditions in eastern Australia across which D. melanogaster is now found provides an opportunity for researchers to identify traits and genes that are associated with climatic adaptation. Allozyme studies indicate clinal patterns for at least four loci including a strong linear cline in Adh and a non-linear cline in alpha-Gpdh. Inversion clines were initially established from cytological studies but have now been validated with larger sample sizes using molecular markers for breakpoints. Recent collections indicate that some genetic markers (Adh and In(3R)Payne) have changed over the last 20 years reflecting continuing evolution. Heritable clines have been established for quantitative traits including wing length/area, thorax length and cold and heat resistance. A cline in egg size independent of body size and a weak cline in competitive ability have also been established. Postulated clinal patterns for resistance to desiccation and starvation have not been supported by extensive sampling. Experiments under laboratory and semi-natural conditions have suggested selective factors generating clinal patterns, particularly for reproductive patterns over winter. Attempts are being made to link clinal variation in traits to specific genes using QTL analysis and the candidate locus approach, and to identify the genetic architecture of trait variation along the cline. This is proving difficult because of inversion polymorphisms that generate disequilibrium among genes. Substantial gaps still remain in linking clines to field selection and understanding the genetic and physiological basis of the adaptive shifts. However D. melanogaster populations in eastern Australia remain an excellent resource for understanding past and future evolutionary responses to climate change.     

Distinguishing between primary and secondary intergradation among morphologically differentiated populations of Anolis marmoratus

Distinguishing between primary and secondary intergradation among differentiated populations, and the relative importance of drift and selection, are persistent problems in evolutionary biology. An historical perspective on population interactions can provide insight into the nature of contacts, and thus help resolve these questions. Continuously distributed populations of Anolis marmoratus from the island of Basse Terre in the Guadeloupean archipelago of the Lesser Antilles show a striking degree of geographic variation in morphology. Initial surveys of mtDNA variation from throughout the Guadeloupean Archipelago revealed one case where levels of sequence difference and phylogenetic relationships of alleles from morphologically differentiated populations from the east coast of Basse Terre were consistent with primary intergradation. In this paper, I examine the genetic population structure of a series of populations spanning this north-south cline in morphological variation to test the hypothesis of primary intergradation. Sequences of the mitochondrial cytochrome-b gene from 50 individuals representing five populations spanning the cline were obtained and fourteen unique haplotypes (differing by 2% or less) were detected. Patterns of nucleotide substitution among haplotypes do not deviate from neutral expectation indicating no effect of selection at the level of mtDNA sequences. Estimates of population structure and gene flow were made using both summary statistics for nucleotide diversity (N_at) and cladistic methods. The results are sensitive to the choice of gene flow model, and this is discussed in detail. Mitochondrial variation in the northern populations may not be at equilibrium, and the phylogeny of alleles is consistent with a recent increase in effective population size. Estimates of nucleotide diversity, gene flow, and the phylogenetic relationships of haplotypes indicate that the southern-most population (representing the extreme of morphological variation along this cline) has been relatively isolated from populations to the north and has experienced a reduced effective population size. The apparent clinal variation between the southern population and the others may therefore reflect secondary contact and introgression rather than primary intergradation.     

Divergence of functional traits at early stages of development in Stipa tenacissima populations distributed along an environmental gradient of the Mediterranean

Assessing differences in plant functional traits (PFTs) along climatic gradients is potentially useful for understanding variation within and across populations, and for predicting their responses to climate change. This study investigates the intraspecific variability of several PFTs in Stipa tenacissima (Alpha grass) seedlings from different populations distributed across a climatic gradient. Seven populations from Tunisia to Spain within a 100–600 mm/year rainfall range were selected. Seedlings from each population were grown in a common garden. We expected the functional characteristics to differ among seedling populations according to their climatic gradient. The response patterns were helpful to predict acclimation and fitness under future climatic conditions in these populations. The seedling development analysis showed differences in PFTs among S. tenacissima populations. The biomass traits analysis revealed that higher above-ground biomass was related to higher below-ground development. The leaf traits proved that seedlings with longer leaf length would have less sclerophyllous leaves, a trade-off between productivity and drought resistance. The root traits analysis reflects seedling strategies to maximize resource uptake efficiency. PFTs showed several significant relationships with climatic conditions. The less rainfall, the higher plant allocation to root systems exploring soil. Higher mean temperatures were related to reduced root/plant development. The PFT analysis proves that species followed the ‘optimal partitioning theory’, in that plants preferentially allocate biomass to acquire the resource that most limits their development. However, both the environmental conditions and genetic diversity in S. tenacissima populations influenced seedling growth and behaviour to face ongoing climate change.

     

Wave-induced abiotic stress shapes phenotypic diversity in a coral reef fish across a geographical cline

While morphological variation across geographical clines has been well documented, it is often unclear whether such changes enhance individual performance to local environments. We examined whether the damselfish Acanthochromis polyacanthus display functional changes in swimming phenotype across a 40-km cline in wave-driven water motion on the Great Barrier Reef, Australia. A. polyacanthus populations displayed strong intraspecific variation in swimming morphology and performance that matched local levels of water motion: individuals on reefs subject to high water motion displayed higher aspect-ratio fins and faster swimming speeds than conspecifics on sheltered reefs. Remarkably, intraspecific variation within A. polyacanthus spanned over half the diversity seen among closely related damselfish species from the same region. We find that local selection driven by wave-induced abiotic stress is an overarching ecological mechanism shaping the inter- and intraspecific locomotor diversity of coral reef fishes.     

An explicit test for the contribution of environmental maternal effects to rapid clinal differentiation in an invasive plant

Population differentiation of alien invasive plants within their non‐native range has received increasingly more attention. Common gardens are typically used to assess the levels of genotypic differentiation among populations. However, in such experiments, environmental maternal effects can influence phenotypic variation among individuals if seed sources are collected from field populations under variable environmental regimes. In the present study, we investigated the causes of an altitudinal cline in an invasive plant. Seeds were collected from Senecio inaequidens (Asteraceae) populations along an altitudinal gradient in southern France. In addition, seeds from the same populations were generated by intra‐population crossings in a climatic chamber. The two seed lots were grown in a common garden in Central Belgium to identify any evidence of environmentally induced maternal effects and/or an altitudinal cline in a suite of life‐history traits. Results failed to detect any environmental maternal effects. However, an altitudinal cline in plant height and above‐ground biomass was found to be independent of the maternal environment.     

Clinal variation, intraspecific heterochrony, and microevolution in the Permian bryozoan Tabulipora carbonaria

Clinal geographic variation across an onshore-to-offshore environmental gradient occurs in the Lower Permian trepostome bryozoan Tabulipora carbonaria collected from three widespread calcareous shales of the Wreford Megacyclothem of Kansas. Separate multivariate statistical tests dclineate significant differences between populations within each shale unit, and between onshore, intermediate, and offshore populations pooled across all three shales, suggesting that even weakly-developed gradients may initiate different intraspecific morphological responses. Growth trajectories for populations along the cline also differ significantly, indicating astogenetic (developmental) heterochrony. Populations from onshore habitats are generally paedomorphic relative to those in more offshore settings, and exhibit pre- and postdisplacement and hypermorphosis in zooecial and acanthostyle characteristics. These heterochronic processes may have increased colonial reproductive potential and the efficiency of water flow for feeding and waste disposal in colonies from onshore habitats. Variation from tightly constrained development (astogenetic plasticity) decreased monotonically in an onshore-to-offshore direction; canalized growth may characterize colonies from more stable offshore habitats, whereas greater flexibility during the growth of colonies from unstable onshore biotopes may have increased their rate of survival. Populations of colonies from stratigraphically successive calcareous shales of the Wreford display patterns of growth that are nearly identical to those found in an offshore-to-onshore direction along the cline. Both clinal and temporal patterns probably resulted from selection for more paedomorphic morphologies in onshore, perhaps unstable, habitats and represent microevolution in T. carbonaria. These local adjustments to environmental conditions may produce variation that affects the rate of macroevolutionary change. □Bryozoa, clines, heterochrony, microevolution, Permian, variability.     

SPATIAL VARIATION PATTERNS IN THE MARINE MACROALGA XIPHOPHORA GLADIATA SSP. GLADIATA (PHAEOPHYTA). II. MORPHOLOGICAL VARIATION OVER LARGE SPATIAL SCALES1

The geographic expression of inter-regional morphological variation in Xiphophora gladiata ssp. gladiata was examined using a combination of discriminant function analyses, canonical correlation analysis and linear regression analysis. The pattern of that expression is an irregularly-stepped cline over scales of hundreds of kilometers, superimposed on a weak complex of microclines driven by exposure and (perhaps) other locally-varying environmental factors. The large-scale irregular cline is probably a response to the physico-chemical characteristics of regional water masses.     

濒危树种格木天然群体自由授粉子代苗期生长变异

格木是我国热带和南亚热带地区主要的珍稀濒危树种之一,具有重要的经济和生态价值。本文应用从格木主要天然分布区收集的6个群体的种子进行育苗试验,揭示格木天然群体子代芽苗和幼苗生长变异状况及其与种子性状、亲本环境因子间的相关性。结果表明,格木天然群体子代芽苗和幼苗表型性状在群体间和群体内均存在丰富的变异,苗高的变异高于地径,苗期选择以高生长为主。格木苗期生长与种子的大小、千粒重相关显著,胚根长度与经度、地径与年均气温、高径比与海拔均存在显著的相关性(p＜0.05)。研究结果为深入开展格木良种选育奠定了基础,亦将有助于其天然资源的保护经营和人工种植业的发展。     

Geographic diversity cline of R gene homologs in wild populations of Solanum pimpinellifolium (Solanaceae)

Plant resistance (R) genes tend to be highly variable within plant species and are thought to be under natural selection; however, little is known about the geographic distribution of R gene diversity within and among plant populations. To determine the possible roles of demography and selection on R gene evolution, patterns of diversity at the multigenic Cf-2 R gene family were studied in Solanum pimpinellifolium populations along the northern coast of Peru. Population diversity levels of Cf-2 homologs follow a latitudinal cline, consistent with the species's history of gradual colonization of the Peruvian coast and population variation in outcrossing levels. Although previous evidence suggests that selection has shaped the DNA sequence content of the Cf-2 genes, current results imply that the geographic distribution of Cf-2 homolog diversity has been shaped primarily by demographic factors or by selective pressures with a clinal distribution.     

Diversity and genetic differentiation among populations of Indian and Kenyan tea (Camellia sinensis (L.) O. Kuntze) revealed by AFLP markers

AFLP markers were successfully employed to detect diversity and genetic differentiation among Indian and Kenyan populations of tea (Camellia sinensis (L.) O. Kuntze). Shannon's index of diversity was used to partition the total phenotypic variation into between and within population components. On average, most of the diversity was detected within populations, with 79% of the variation being within and 21% being between populations of Indian and Kenyan tea. A dendrogram constructed on the basis of band sharing distinctly separated the three populations of tea into China type (sinensis), Assam type (assamica) and Cambod type (assamica ssp. lasiocalyx) in a manner consistent with the present taxonomy of tea, the known pedigree of some of the genotypes and their geographical origin. Principal coordinate (PCO) analysis grouped Assam genotypes both from India and Kenya supporting the suggestion that the Kenyan clones have been derived from collections made in this region. The China types were more dispersed on the PCO plot which is a reflection of wider genetic variation. As would be expected, clones collected from the same region exhibited less overall genetic variation. AFLP analysis discriminated all of the tested genotypes from India and Kenya, even those which cannot be distinguished on the basis of morphological and phenotypic traits. Received: 2 May 1996 / Accepted: 14 June 1996     

Habitat fragmentation promotes fluctuating asymmetry but not morphological divergence in two geckos

The process of species formation is often ignored in discussions on the conservation of biodiversity. Yet the clearance of vegetation may promote divergence among populations of a species through isolation, providing conditions for rapid genetic drift and novel selection pressures. Here, stepwise discriminant function analysis and fluctuating asymmetry are used to examine variation in morphology within and among non fragmented and recently fragmented populations of two species of gecko,Oedura reticulata andGehyra variegata. High reclassification error rates using discriminant function analysis, indicate that fragmentation has had no detectable effect on morphological differentiation among populations of either species. In contrast, remnant populations of both species exhibit higher mean levels of fluctuating asymmetry than do populations in undisturbed habitat. ForOedura reticulata, levels of fluctuating asymmetry are negatively correlated with the log of adult population size. These results suggest that the changes following habitat clearance have been severe enough to cause increased developmental instability in populations of both species but not detectable morphological divergence. Given the high rate of extirpation of gecko populations in the study region and the extreme vulnerability of the remaining populations, it is unlikely that species formation will be significant in maintaining reptile diversity in that region.     

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

Lack of congruence between morphometric evolution and genetic differentiation suggests a recent dispersal and local habitat adaptation of the Madeiran lizard Lacerta dugesii

Genetic differentiation among nine populations of the endemic lizard Lacerta dugesii Milne-Edwards 1829 (Lacertidae) from four groups of islands constituting the Archipelago of Madeira, was investigated by protein electrophoresis at 23 enzyme loci. Among twenty polymorphic loci, the total genetic diversity was due primarily to intra-population variation. The allele and genotypic frequencies among populations showed some heterogeneity, allowing the species to present a structuring pattern compatible with their geographical clustering. Some evidence suggests that selection acting on some loci in different ecological conditions may be responsible for the clustering of the populations studied. There was no apparent isolation effect expected under an "island" model of population divergence, and no correlation was found between genetic and geographic distances among populations. Morphological variation of the proposed three L. dugesii subspecies is not congruent with the allozyme analysis. This most probably suggests a rapid colonization of the islands followed by a strong effect of selection operating over the morphological characters used to define the subspecies.     

SPATIAL AND TEMPORAL VARIATION IN SELECTION ON CORRELATED LIFE-HISTORY TRAITS AND PLANT SIZE IN CHAMAECRISTA FASCICULATA

Intraspecific studies of selection on multiple traits of a plant's life history provide insight as to how the composite life history of an organism evolves. Current understanding of selection on plant life-history traits is deficient in three important areas: 1) the effects of selection through correlated traits, 2) the effects of selection on a trait throughout the plant's lifetime, and 3) spatial and temporal variation in selection on plant life-history traits among populations and years. This study documents spatial and temporal variation in selection on three life-history and two morphological traits for two natural populations of Chamaecrista fasciculata, a native summer annual. Life-history and morphological traits (date of seedling emergence, size at establishment, size prior to reproduction, date of initial flowering, and date of initial fruit maturation) varied significantly between sites and/or years. Selection on traits varied either spatially, between sites and among transects within one site, or temporally, between years. In addition, life-history traits were phenotypically correlated among themselves and with morphological traits; correlations were generally constant over time and space. Indirect selection caused changes in means and variances in traits not under direct selection, but which were correlated with traits under selection. Selection on date of emergence varied in direction and magnitude among different life-cycle stages, while selection on other traits varied only in magnitude among life stages of the plant. This study documents the complexity of the selective process and the importance of considering multiple life stages and traits when studying the evolution of life-history traits.     

Geographic variation in the molecular,morphological and reproductive characters of the ayuPlecoglossus altivelis (plecoglossidae) in the Japan-Ryukyu Archipelago

In an attempt to clarify geographic variation in the ayuPlecoglossus altivelis, the molecular, morphological and reproductive characters of populations from various localities in the Japan-Ryukyu Archipelago were examined. An electrophoretic analysis of enzymes encoded by 18–28 loci and univariate and multivariate analyses of 12 meristic and 17 morphometric characters indicated that there are two distinct groups of ayu in the archipelago: that of the Japan Islands and that of the Ryukyu Islands. Within the Japan Islands group, a landlocked population of Lake Biwa was found to diverge somewhat from other amphidromous populations genetically and morphologically. The Lake Biwa population was noted as being by far the most fecund, departing from a north-south cline of fecundity observed in amphidromous populations. Though no virtual divergence was observed among amphidromous populations in the Japan Islands in molecular data, divergence among these populations was considerable in morphological and fecundity data. Differences in geographic variation patterns are attributed to qualitative differences inherent in molecular and organismal characters: organismal characters are responsive to selection and also directly altered by environmental influences in part but molecular characters are hardly so. The evolution of different characters may have proceeded at different rates among geographic populations of ayu.     

Molecular Analysis of an Allozyme Cline: Alcohol Dehydrogenase in Drosophila Melanogaster on the East Coast of North America

Clines may either be selectively maintained or be the by-product of nonadaptive processes related to population structure and history. Drosophila melanogaster populations on the east coast of North America show a latitudinal cline in the frequencies of two common electrophoretically distinguishable alleles at the alcohol dehydrogenase locus (Adh), designated Adh-S and Adh-F. This cline may either be adaptive or an artifact of a possible recent dual founding of North American D. melanogaster populations in which frequencies of Adh alleles differed between founder populations. By means of a high resolution restriction-mapping technique, we studied the distribution of 113 haplotypes derived from 44 polymorphic DNA markers within the Adh region in 1533 individuals from 25 populations throughout the cline. We found significant clinal differentiation at the polymorphisms determining the mobility-difference causing amino acid replacement between Adh-F and Adh-S alleles. Hitchhiking was limited, despite extensive linkage disequilibrium, and other sites did not vary clinally. Such a pattern of differentiation implies that selection is responsible for the cline. To investigate whether selection acts only on the Adh-F/S site, we performed a ``selective equivalence'''' test under the assumption that all variability within the specified allelic class is selectively neutral. This revealed selective equivalence among Adh-S-bearing haplotypes, whose frequencies showed no differentiation throughout the cline, implying high levels of frequency-homogenizing gene flow. Geographical heterogeneity among Adh-F-bearing haplotypes implied the action of selection on one or more additional variants in linkage disequilibrium with Adh-F. In a further study of a subset of the data (n = 1076 from 18 populations), we found a combined insertion/deletion polymorphism, designated &1, located in the 5'' adult intron and in linkage disequilibrium with Adh-F, to show more marked clinal variation than Adh-F/S. Although the unequivocal identification of the precise target(s) of selection requires further study, we suggest that clinal selection may be acting epistatically on the Adh-F/S and &1 polymorphisms.