DOES PHENOTYPIC PLASTICITY FOR ADULT SIZE VERSUS FOOD LEVEL IN DROSOPHILA MELANOGASTER EVOLVE IN RESPONSE TO ADAPTATION TO DIFFERENT REARING DENSITIES?

Recent studies with Drosophila have suggested that there is extensive genetic variability for phenotypic plasticity of body size versus food level. If true, we expect that the outcome of evolution at very different food levels should yield genotypes whose adult size show different patterns of phenotypic plasticity. We have tested this prediction with six independent populations of Drosophila melanogaster kept at extreme densities for 125 generations. We found that the phenotypic plasticity of body size versus food level is not affected by selection or the presence of competitors of a different genotype. However, we document increasing among population variation in phenotypic plasticity due to random genetic drift. Several reasons are explored to explain these results including the possibility that the use of highly inbred lines to make inferences about the evolution of genetically variable populations may be misleading.     

Spatial predictors of genomic and phenotypic variation differ in a lowland Middle American bird (Icterus gularis)

Spatial patterns of intraspecific variation are shaped by geographical distance among populations, historical changes in gene flow and interactions with local environments. Although these factors are not mutually exclusive and operate on both genomic and phenotypic variation, it is unclear how they affect these two axes of variation. We address this question by exploring the predictors of genomic and phenotypic divergence in Icterus gularis, a broadly distributed Middle American bird that exhibits marked geographical variation in body size across its range. We combined a comprehensive single nucleotide polymorphism and phenotypic data set to test whether genome‐wide genetic and phenotypic differentiation are best explained by (i) isolation by distance, (ii) isolation by history or (iii) isolation by environment. We find that the pronounced genetic and phenotypic variation in I. gularis are only partially correlated and differ regarding spatial predictors. Whereas genomic variation is largely explained by historical barriers to gene flow, phenotypic diversity can be best predicted by contemporary environmental heterogeneity. Our genomic analyses reveal strong phylogeographical structure coinciding with the Chivela Pass at the Isthmus of Tehuantepec that was formed during the Pleistocene, when populations were isolated in north–south refugia. In contrast, we found a strong association between body size and environmental variables, such as temperature and precipitation. The relationship between body size and local climate is consistent with a pattern produced by either natural selection or environmental plasticity. Overall, these results provide empirical evidence for why phenotypic and genomic data are often in conflict in taxonomic and phylogeographical studies.     

Quantitative genetic analysis of the body size and shape of Drosophila buzzatii

Summary Body size in Drosophila is known to be closely related to a number of traits with important life history consequences, such as fecundity, dispersal ability and mating success. We examine the quantitative genetic basis of body size in three populations of the cactophilic species Drosophila buzzatii, which inhabit climatically different areas of Australia. Flies were reared individually to eliminate any common environmental component in a full-sib design with families split between two temperatures (18° and 25 °C). The means of several size measures differ significantly among populations while the genetic correlations among these traits generally do not differ, either among populations from different natural environments or between the different laboratory temperatures. This stability of correlation structure is necessary if laboratory estimates of genetic correlations are to have any connection with the expression of genetic variation in the field. The amount of variance due to genotype-by-environment interactions (family x temperature of development) varied among populations, apparently in parallel with the magnitudes of seasonal and diurnal variation in temperature experienced by the different populations. A coastal population, inhabiting a relatively thermally benign environment, showed no interaction, while two inland populations, inhabiting thermally more extreme areas, showed interaction. This interaction term is a measure of the amount of genetic variation in the degree of phenotypic plasticity of body size in response to temperature of development. Thus the inland flies vary in their ability to attain a given body size at a particular temperature while the coastal flies do not. This phenotypic plasticity is shown to be due primarily to differences among genotypes in the amount of response to the change in temperature. A possible selective basis for the maintenance of genetic variation for the levels of phenotypic plasticity is proposed.     

Ontogenetic changes in genetic variances of age-dependent plasticity along a latitudinal gradient

The expression of phenotypic plasticity may differ among life stages of the same organism. Age-dependent plasticity can be important for adaptation to heterogeneous environments, but this has only recently been recognized. Whether age-dependent plasticity is a common outcome of local adaptation and whether populations harbor genetic variation in this respect remains largely unknown. To answer these questions, we estimated levels of additive genetic variation in age-dependent plasticity in six species of damselflies sampled from 18 populations along a latitudinal gradient spanning 3600 km. We reared full sib larvae at three temperatures and estimated genetic variances in the height and slope of thermal reaction norms of body size at three points in time during ontogeny using random regression. Our data show that most populations harbor genetic variation in growth rate (reaction norm height) in all ontogenetic stages, but only some populations and ontogenetic stages were found to harbor genetic variation in thermal plasticity (reaction norm slope). Genetic variances in reaction norm height differed among species, while genetic variances in reaction norm slope differed among populations. The slope of the ontogenetic trend in genetic variances of both reaction norm height and slope increased with latitude. We propose that differences in genetic variances reflect temporal and spatial variation in the strength and direction of natural selection on growth trajectories and age-dependent plasticity. Selection on age-dependent plasticity may depend on the interaction between temperature seasonality and time constraints associated with variation in life history traits such as generation length.     

Within-population variation in body size plasticity in response to combined nutritional and thermal stress is partially independent from variation in development time

Ongoing climate change has forced animals to face changing thermal and nutritional environments. Animals can adjust to such combinations of stressors via plasticity. Body size is a key trait influencing organismal fitness, and plasticity in this trait in response to nutritional and thermal conditions varies among genetically diverse, locally adapted populations. The standing genetic variation within a population can also influence the extent of body size plasticity. We generated near-isogenic lines from a newly collected population of Drosophila melanogaster at the mid-point of east coast Australia and assayed body size for all lines in combinations of thermal and nutritional stress. We found that isogenic lines showed distinct underlying patterns of body size plasticity in response to temperature and nutrition that were often different from the overall population response. We then tested whether plasticity in development time could explain, and therefore regulate, variation in body size to these combinations of environmental conditions. We selected five genotypes that showed the greatest variation in response to combined thermal and nutritional stress and assessed the correlation between response of developmental time and body size. While we found significant genetic variation in development time plasticity, it was a poor predictor of body size among genotypes. Our results therefore suggest that multiple developmental pathways could generate genetic variation in body size plasticity. Our study emphasizes the need to better understand genetic variation in plasticity within a population, which will help determine the potential for populations to adapt to ongoing environmental change.     

Are latitudinal clines in body size adaptive?

Body size of animals often increases with increasing latitude. These latitudinal clines in body size have interested biologists for over 150 years. However, the mechanisms that generate these clines in size are still unclear, though latitudinal gradients in temperature appear to play an important role. More importantly, many studies that examine latitudinal clines in body size and the mechanisms responsible for these clines use phenotypic data, confounding genetic (adaptive) and non‐genetic (plasticity) sources of variation. Yet, most of these studies make adaptive conclusions based on phenotypic measures of size. Here I show the dangers of making adaptive inferences from phenotypic measures of size. In addition, I use a specific form of plasticity in body size of ectotherms, called the temperature–size rule, to illustrate how confusion about genetic and non‐genetic contributions to phenotypic variation has hampered progress in understanding the evolution of latitudinal clines in size. Field‐based measurements of body size can no doubt be influenced by plasticity, but demonstrating that latitudinal clines have a genetic basis is necessary to show that these patterns are adaptive.     

Adaptive phenotypic plasticity in a clonal invader

Organisms featuring wide trait variability and occurring in a wide range of habitats, such as the ovoviviparous New Zealand freshwater snail Potamopyrgus antipodarum, are ideal models to study adaptation. Since the mid‐19th century, P. antipodarum, characterized by extremely variable shell morphology, has successfully invaded aquatic areas on four continents. Because these obligately and wholly asexual invasive populations harbor low genetic diversity compared to mixed sexual/asexual populations in the native range, we hypothesized that (1) this phenotypic variation in the invasive range might be adaptive with respect to colonization of novel habitats, and (2) that at least some of the variation might be caused by phenotypic plasticity. We surveyed 425 snails from 21 localities across northwest Europe to attempt to disentangle genetic and environmental effects on shell morphology. We analyzed brood size as proxy for fitness and shell geometric morphometrics, while controlling for genetic background. Our survey revealed 10 SNP genotypes nested into two mtDNA haplotypes and indicated that mainly lineage drove variation in shell shape but not size. Physicochemical parameters affected both shell shape and size and the interaction of these traits with brood size. In particular, stronger stream flow rates were associated with larger shells. Our measurements of brood size suggested that relatively larger slender snails with relatively large apertures were better adapted to strong flow than counterparts with broader shells and relatively small apertures. In conclusion, the apparent potential to modify shell morphology plays likely a key role in the invasive success of P. antipodarum; the two main components of shell morphology, namely shape and size, being differentially controlled, the former mainly genetically and the latter predominantly by phenotypic plasticity.     

Co-Gradient Variation in Growth Rate and Development Time of a Broadly Distributed Butterfly

Widespread species often show geographic variation in thermally-sensitive traits, providing insight into how species respond to shifts in temperature through time. Such patterns may arise from phenotypic plasticity, genetic adaptation, or their interaction. In some cases, the effects of genotype and temperature may act together to reduce, or to exacerbate, phenotypic variation in fitness-related traits across varying thermal environments. We find evidence for such interactions in life-history traits of Heteronympha merope, a butterfly distributed across a broad latitudinal gradient in south-eastern Australia. We show that body size in this butterfly is negatively related to developmental temperature in the laboratory, in accordance with the temperature-size rule, but not in the field, despite very strong temperature gradients. A common garden experiment on larval thermal responses, spanning the environmental extremes of H. merope''s distribution, revealed that butterflies from low latitude (warmer climate) populations have relatively fast intrinsic growth and development rates compared to those from cooler climates. These synergistic effects of genotype and temperature across the landscape (co-gradient variation) are likely to accentuate phenotypic variation in these traits, and this interaction must be accounted for when predicting how H. merope will respond to temperature change through time. These results highlight the importance of understanding how variation in life-history traits may arise in response to environmental change. Without this knowledge, we may fail to detect whether organisms are tracking environmental change, and if they are, whether it is by plasticity, adaptation or both.     

Integrating patterns of thermal tolerance and phenotypic plasticity with population genetics to improve understanding of vulnerability to warming in a widespread copepod

Differences in population vulnerability to warming are defined by spatial patterns in thermal adaptation. These patterns may be driven by natural selection over spatial environmental gradients, but can also be shaped by gene flow, especially in marine taxa with high dispersal potential. Understanding and predicting organismal responses to warming requires disentangling the opposing effects of selection and gene flow. We begin by documenting genetic divergence of thermal tolerance and developmental phenotypic plasticity. Ten populations of the widespread copepod Acartia tonsa were collected from sites across a large thermal gradient, ranging from the Florida Keys to Northern New Brunswick, Canada (spanning over 20° latitude). Thermal performance curves (TPCs) from common garden experiments revealed local adaptation at the sampling range extremes, with thermal tolerance increasing at low latitudes and decreasing at high latitudes. The opposite pattern was observed in phenotypic plasticity, which was strongest at high latitudes. No relationship was observed between phenotypic plasticity and environmental variables. Instead, the results are consistent with the hypothesis of a trade‐off between thermal tolerance and the strength of phenotypic plasticity. Over a large portion of the sampled range, however, we observed a remarkable lack of differentiation of TPCs. To examine whether this lack of divergence is the result of selection for a generalist performance curve or constraint by gene flow, we analyzed cytochrome oxidase I mtDNA sequences, which revealed four distinct genetic clades, abundant genetic diversity, and widely distributed haplotypes. Strong divergence in thermal performance within genetic clades, however, suggests that the pace of thermal adaptation can be relatively rapid. The combined insight from the laboratory physiological experiments and genetic data indicate that gene flow constrains differentiation of TPCs. This balance between gene flow and selection has implications for patterns of vulnerability to warming. Taking both genetic differentiation and phenotypic plasticity into account, our results suggest that local adaptation does not increase vulnerability to warming, and that low‐latitude populations in general may be more vulnerable to predicted temperature change over the next century.     

Allometric influence on phenotypic variation in the Song Sparrow (Melospiza melodia)

Song Sparrow ( Melospiza melodia ) populations found along the Pacific Coast of North America, from Baja California to the islands off the coast of Alaska, exhibit extensive morphological variation. With a multivariate analysis of size and shape, I describe a portion of this pattern and examine how it could be maintained despite gene flow among the populations. Because shape differences fall along geographic barriers, I suggest that similarities among Song Sparrow populations in multivariate shape reflect their pattern of genetic relatedness. A general pattern of Song Sparrow post-nestling growth allometry has been discovered: bill characteristics are positively allometric and all other characteristics are negatively allometric. In contrast to shape, multivariate patterns of body size variation do not correspond to geographic relationships. In combination with evidence of Song Sparrow phenotypic plasticity, it is proposed that multivariate body size is an environmentally plastic trait and that specific traits exhibit levels of phenotypic plasticity in proportion to their rate of growth with respect to body size. In this way local environmental factors which alter body size may change an entire suite of allometrically related traits and thus create striking patterns of morphological variation.     

The contribution of phenotypic plasticity to adaptation in Lacerta vivipara

Correlation between intraspecific phenotypic variability and variation of environmental conditions could reflect adaptation. Different phenotypes may result from differential expression of a genotype in different environments (phenotypic plasticity) or from expression of different genotypes (genetic diversity). Populations of Lacerta vivipara exhibit larger adult body length, lower age at maturity, higher fecundity, and smaller neonatal size in humid habitats compared to dry habitats. We conducted reciprocal transplants of juvenile L. vivipara to test for the genetic or plastic origin of this variation. We captured gravid females from four populations that differed in the relative humidity of their habitats, and during the last 2 to 4 weeks of gestation, we manipulated heat and water availability under laboratory conditions. Juveniles were released into the different populations and families were divided to compare growth rate and survival of half-sibs in two environments. Growth rate and survival were assessed using capture-recapture techniques. Growth rate was plastic in response to postnatal conditions and did not differ between populations of origin. Survival differed between populations of origin, partially because of differences in neonatal body length. The response of juvenile body length and body condition to selection in the different habitats was affected by the population of origin. This result cannot be simply interpreted in terms of adaptation; however, phenotypic plasticity of fecundity or juvenile size most probably resulted in adaptive reproductive strategies. Adaptation to the habitat by means of genetic specialization was not detected. Further investigation is needed to discriminate between genetic and long-term maternal effects.     

Development time and size-related traits in the oriental blowfly, Chrysomya megacephala along a latitudinal gradient from China

Organisms can respond to variation in temperature through the direct effect of temperature on phenotypes (phenotypic plasticity), or through long-term adaptation to temperature (and thus evolution of either mean size or thermal reaction norm). We examined the effects of various temperatures (of 20 and 30 °C) on development time, adult body size (body length and body width) and pre-adult survivorship in six populations of Chrysomya megacephala, collected at different latitudes. We found that temperature changes induced substantial plasticity in terms of development time, body size and pre-adult survivorship, indicating that developmental temperature significantly affects growth and life history traits of C. megacephala. We also detected genetic differences among populations for body size and development time, and these two traits exhibited highly significant variations in the responses of different populations to various temperature conditions, indicating genetic differences among populations in terms of thermal reaction norms. The latitude of origin of the different populations (and hence mean temperature regimes in the environments from where the populations originated) did not appear to fully explain these genetic differences. In short, changes in development time and body size in C. megacephala can be regarded as adaptations to changing thermal regimes.     

Standing variation boosted by multiple sources of introduction contributes to the success of the introduced species, Lotus corniculatus

Although ecological differences between native and introduced ranges have been considered to drive rapid expansion of invasive species, recent studies suggest that rapid evolutionary responses of invasive species to local environments may also be common. Such expansion across heterogeneous environments by adaptation to local habitats requires genetic variation. In this study, we investigated the source and role of standing variation in successful invasion of heterogeneous abiotic environments in a self-incompatible species, Lotus corniculatus. We compared phenotypic and genetic variation among cultivars, natives, and introduced genotypes, and found substantial genetic variation within both native and introduced populations. Introduced populations possessed genotypes derived from both cultivars and native populations, and had lower population differentiation, indicating multiple sources of introduction and population admixture among the sources in the introduced range. Both cultivars and introduced populations had similarly outperforming phenotypes on average, with increased biomass and earlier flowering compared with native populations, but those phenotypes were within the range of the variation in phenotypes of the native populations. In addition, clinal variation within introduced populations was detected along a climatic gradient. Multiple introductions from different sources, including cultivars, may have contributed to pre-adaptive standing variation in the current introduced populations. We conclude that both introduction of cultivar genotypes and natural selection in local environments contributed to current patterns of genetic and phenotypic variation observed in the introduced populations.     

Geographic variation in body size, sexual size dimorphism and fitness components of a seed beetle: local adaptation versus phenotypic plasticity

Variation in body size, growth and life history traits of ectotherms along latitudinal and altitudinal clines is generally assumed to represent adaptation to local environmental conditions, especially adaptation to temperature. However, the degree to which variation along these clines is due to adaptation vs plasticity remains poorly understood. In addition, geographic patterns often differ between females and males – e.g. sexual dimorphism varies along latitudinal clines, but the extent to which these sex differences are due to genetic differences between sexes vs sex differences in plasticity is poorly understood. We use common garden experiments (beetles reared at 24, 30 and 36°C) to quantify the relative contribution of genetically‐based differentiation among populations vs phenotypic plasticity to variation in body size and other traits among six populations of the seed‐feeding beetle Stator limbatus collected from various altitudes in Arizona, USA. We found that temperature induces substantial plasticity in survivorship, body size and female lifetime fecundity, indicating that developmental temperature significantly affects growth and life history traits of S. limbatus. We also detected genetic differences among populations for body size and fecundity, and genetic differences among populations in thermal reaction norms, but the altitude of origin (and hence mean temperature) does not appear to explain these genetic differences. This and other recent studies suggest that temperature is not the major environmental factor that generates geographic variation in traits of this species. In addition, though there was no overall difference in plasticity of body size between males and females (when averaged across populations), we did find that the degree to which dimorphism changed with temperature varied among populations. Consequently, future studies should be extremely cautious when using only a few study populations to examine environmental effects on sexual dimorphism.     

Experimental alteration of DNA methylation affects the phenotypic plasticity of ecologically relevant traits in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Heritable phenotypic variation in plants can be caused not only by underlying genetic differences, but also by variation in epigenetic modifications such as DNA methylation. However, we still know very little about how relevant such epigenetic variation is to the ecology and evolution of natural populations. We conducted a greenhouse experiment in which we treated a set of natural genotypes of Arabidopsis thaliana with the demethylating agent 5-azacytidine and examined the consequences of this treatment for plant traits and their phenotypic plasticity. Experimental demethylation strongly reduced the growth and fitness of plants and delayed their flowering, but the degree of this response varied significantly among genotypes. Differences in genotypes’ responses to demethylation were only weakly related to their genetic relatedness, which is consistent with the idea that natural epigenetic variation is independent of genetic variation. Demethylation also altered patterns of phenotypic plasticity, as well as the amount of phenotypic variation observed among plant individuals and genotype means. We have demonstrated that epigenetic variation can have a dramatic impact on ecologically important plant traits and their variability, as well as on the fitness of plants and their ecological interactions. Epigenetic variation may thus be an overlooked factor in the evolutionary ecology of plant populations.     

Altitudinal life history variation in the dung flies Scathophaga stercoraria and Sepsis cynipsea

 Field phenologies of high- (ca. 1500 m) and low- (ca. 500 m) altitude populations of the two most common European species of dung flies, Scathophaga stercoraria and Sepsis cynipsea, differ quite markedly due to differences in climate. To differentiate genetic adaptation due to natural selection and phenotypic plasticity, I compared standard life history characters of pairs of high- and low-altitude populations from three disjunctive sites in Switzerland in a laboratory experiment. The F1 rearing environment did not affect any of the variables of the F2 generation with which all experiments were conducted; hence, there were no carry-over or maternal effects. In Sc. stercoraria, high-altitude individuals were smaller but laid larger eggs; the latter may be advantageous in the more extreme (i.e. more variable and less predictable) high-altitude climate. Higher rearing temperature strongly decreased development time, body size and the size difference between males and females (males are larger), produced female-biased sex ratios and led to suboptimal adult emergence rates. Several of these variables also varied among the three sites, producing some interactions complicating the patterns. In Se. cynipsea, high-altitude females were marginally smaller, less long-lived and laid fewer clutches. Higher rearing temperature strongly decreased development time and body size but tended to increase the size difference between males and females (males are smaller); it also increased clutch size but decreased physiological longevity. Again, interpretation is complicated by variation across sites and some significant interactions. Overall, genetic adaptation to high-altitude conditions appears weak, probably prevented by substantial gene flow, and may be swamped by the effects of other geographic variables among populations. In contrast, phenotypic plasticity is extensive. This may be due to selection of flexible, multi-purpose genotypes. The results suggest that differences in season length between high- and low-altitude locations alone do not explain well the patterns of variation in phenology and body size. Received: 21 March 1996 / Accepted: 1 September 1996     

SIMILARITIES AND DIFFERENCES IN ALTITUDINAL VERSUS LATITUDINAL VARIATION FOR MORPHOLOGICAL TRAITS IN DROSOPHILA MELANOGASTER

Understanding how natural environments shape phenotypic variation is a major aim in evolutionary biology. Here, we have examined clinal, likely genetically based variation in morphology among 19 populations of the fruit fly (Drosophila melanogaster) from Africa and Europe, spanning a range from sea level to 3000 m altitude and including locations approximating the southern and northern range limit. We were interested in testing whether latitude and altitude have similar phenotypic effects, as has often been postulated. Both latitude and altitude were positively correlated with wing area, ovariole number, and cell number. In contrast, latitude and altitude had opposite effects on the ratio between ovariole number and body size, which was negatively correlated with egg production rate per ovariole. We also used transgenic manipulation to examine how increased cell number affects morphology and found that larger transgenic flies, due to a higher number of cells, had more ovarioles, larger wings, and, unlike flies from natural populations, increased wing loading. Clinal patterns in morphology are thus not a simple function of changes in body size; instead, each trait might be subject to different selection pressures. Together, our results provide compelling evidence for profound similarities as well as differences between phenotypic effects of latitude and altitude.     

Genetic differences and phenotypic plasticity in body size between high- and low-altitude populations of the ground beetle Carabus tosanus

The body size of a univoltine carabid beetle Carabus tosanus on Shikoku Island, Japan, was clearly smaller in higher‐altitude populations (subspecies), which possibly represents incipient speciation. To explore the determinants of altitudinal differences in body size in this species, we studied the degree of phenotypic plasticity by conducting rearing experiments at two constant temperatures and examined genetic differences through interpopulation crosses. At 15 °C, C. tosanus had a longer developmental period and a shorter adult body than at 20 °C. Nevertheless, variation in body size due to temperature effects (phenotypic plasticity) was small compared to the interpopulation differences, which suggests substantial genetic differences between populations (subspecies) at different altitudes. In F₁ offspring from crosses between a low‐altitude (subspecies tosanus) and a high‐altitude population (subspecies ishizuchianus), adult body length was affected by the genotypes of both parents, with an interaction effect of parental genotype and offspring sex. Further analyses revealed that adult body length was affected by sex‐linked factors in addition to autosomal factors. These genetic differences in body size may have resulted from adaptations to different altitudes and may be important for the process of incipient speciation because body size differences could contribute to premating reproductive isolation.     

Testing Local Host Adaptation and Phenotypic Plasticity in a Herbivore When Alternative Related Host Plants Occur Sympatrically

Host race formation in phytophagous insects can be an early stage of adaptive speciation. However, the evolution of phenotypic plasticity in host use is another possible outcome. Using a reciprocal transplant experiment we tested the hypothesis of local adaptation in the aphid Brevicoryne brassicae. Aphid genotypes derived from two sympatric host plants, Brassica oleracea and B. campestris, were assessed in order to measure the extent of phenotypic plasticity in morphological and life history traits in relation to the host plants. We obtained an index of phenotypic plasticity for each genotype. Morphological variation of aphids was summarized by principal components analysis. Significant effects of recipient host on morphological variation and life history traits (establishment, age at first reproduction, number of nymphs, and intrinsic growth rate) were detected. We did not detected genotype × host plant interaction; in general the genotypes developed better on B. campestris, independent of the host plant species from which they were collected. Therefore, there was no evidence to suggest local adaptation. Regarding plasticity, significant differences among genotypes in the index of plasticity were detected. Furthermore, significant selection on PC1 (general aphid body size) on B. campestris, and on PC1 and PC2 (body length relative to body size) on B. oleracea was detected. The elevation of the reaction norm of PC1 and the slope of the reaction norm for PC2 (i.e., plasticity) were under directional selection. Thus, host plant species constitute distinct selective environments for B. brassicae. Aphid genotypes expressed different phenotypes in response to the host plant with low or nil fitness costs. Phenotypic plasticity and gene flow limits natural selection for host specialization promoting the maintenance of genetic variation in host exploitation.