Testing multivariate patterns of within-island differentiation in Podarcis dugesii from Madeira

Podarcis dugesii is a polymorphic species found on the ecologically-heterogeneous Atlantic island of Madeira. Patterns of geographic variation in size-free body dimensions (body shape), scalation, and size were portrayed on maps with isophene contours. There was no dominant pattern of among-site variation as seen in several other island lizards, although some linear combinations of morphological variables showed evidence of clines towards the east coast. Multivariate statistical techniques were developed to evaluate whether morphological variation contained an ecology-related spatial pattern, after removal of the effects of statistical non-independence of studied populations. Specifically, a partial correlation matrix comparison method provided evidence to support a proposed elevation model of among-site divergence in scalation and body size, with both an alternative rainfall model and spatial non-independence held constant. Similar divergence in scalation has been observed in other lizard species on other island archipelagos as well on some continental areas. The rainfall model was rejected for both body shape and scalation. The degree of within-population polymorphism in body shape and scalation showed no clear geographic trends in general, although the elevation and rainfall models could not be rejected for female scalation.     

ENVIRONMENTAL AND GENETIC DETERMINANTS OF PLANT SIZE IN VIOLA SORORIA

Survivorship and fecundity in the forest herb, Viola sororia, are size-dependent. The basis of size variation among individuals of Viola sororia was investigated with a uniform environment experiment. Plants collected from natural populations were vegetatively reproduced and grown under two light regimes in a greenhouse. Analysis of quantitative variation showed: 1) significant differences between light treatments for characters related to plant shape and relative growth rate; 2) significant among-genet variation for plant size, plant shape and relative growth rate but none for physiological characters; and 3) a size threshold for cleistogamous seed production and rhizome production. Heritability estimates for the characters associated with plant size and shape ranged from 0.09 to 0.39, indicating significant genetic determination for these traits. In addition, among-genet differences in relative growth rate were substantial. The results of this study suggest that the size variation found in natural populations is not solely a function of environmental heterogeneity but is significantly influenced by the genotypes composing the population.     

EVOLUTION OF VARIATION AND VARIABILITY UNDER FLUCTUATING,STABILIZING, AND DISRUPTIVE SELECTION

How variation and variability (the capacity to vary) may respond to selection remain open questions. Indeed, effects of different selection regimes on variational properties, such as canalization and developmental stability are under debate. We analyzed the patterns of among‐ and within‐individual variation in two wing‐shape characters in populations of Drosophila melanogaster maintained under fluctuating, disruptive, and stabilizing selection for more than 20 generations. Patterns of variation in wing size, which was not a direct target of selection, were also analyzed. Disruptive selection dramatically increased phenotypic variation in the two shape characters, but left phenotypic variation in wing size unaltered. Fluctuating and stabilizing selection consistently decreased phenotypic variation in all traits. In contrast, within‐individual variation, measured by the level of fluctuating asymmetry, increased for all traits under all selection regimes. These results suggest that canalization and developmental stability are evolvable and presumably controlled by different underlying genetic mechanisms, but the evolutionary responses are not consistent with an adaptive response to selection on variation. Selection also affected patterns of directional asymmetry, although inconsistently across traits and treatments.     

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

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

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome‐scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well‐defined categories of locus types (e.g., neutral vs. highly selected). Our multi‐locus analysis reveals cross‐genomic variation in selective constraints on gene flow and locus‐specific flexibility in the permeability of the interspecies membrane.     

Rapid evolution of wing size clines in Drosophila subobscura

Parallel latitudinal clines across species and continents provide dramatic evidence of the efficacy of natural selection, however little is known about the dynamics involved in cline formation. For example, several drosophilids and other ectotherms increase in body and wing size at higher latitudes. Here we compare evolution in an ancestral European and a recently introduced (North America) cline in wing size and shape in Drosophila subobscura. We show that clinal variation in wing size, spanning more than 15 degrees of latitude, has evolved in less than two decades. In females from Europe and North America, the clines are statistically indistinguishable however the cline for North American males is significantly shallower than that for European males. We document that while overall patterns of wing size are similar on two continents, the European cline is obtained largely through changing the proximal portion of the wing, whereas the North American cline is largely in the distal portion. We use data from sites collected in 1986/1988 (Pegueroles et al. 1995) and our 1997 collections to compare synchronic (divergence between contemporary populations that share a common ancestor) and allochronic (changes over time within a population) estimates of the rates of evolution. We find that, for these populations, allochronically estimated evolutionary rates within a single population are over 0.02 haldanes (2800 darwins), a value similar in magnitude to the synchronic estimates from the extremes of the cline. This paper represents an expanded analysis of data partially presented in Huey et al. (2000).     

ADAPTATION TO SEASONALITY IN A CRICKET: PATTERNS OF PHENOTYPIC AND GENOTYPIC VARIATION IN BODY SIZE AND DIAPAUSE EXPRESSION ALONG A CLINE IN SEASON LENGTH

This paper investigates patterns of phenotypic variation in the striped ground cricket (Allonemobious fasciatus) along a cline in season length and tests the hypothesis that variation in body size and diapause propensity is the result of diversifying selection due to different local conditions. We examined 83 populations and found that A. fasciatus produced a single generation per year (univoltine) north of 37°N latitude and was bivoltine south of 35°30′N. Body size generally increased with increasing season length, with a sudden drop in the region corresponding to the transition from univoltine to bivoltine life cycles, reflecting the division of total season length in two within-bivoltine populations. We reared ten populations in our laboratory and found that much of the interpopulation variance in body size observed in the field could be attributed to genetic differences. Diapause expression also varied significantly among populations and was strongly correlated with season length. The heritability of body size did not differ between populations, but full-sib estimates greatly exceeded parent–offspring estimates suggesting that there are important nonadditive genetic effects. The heritability of diapause expression, determined from analysis of full sibs, varied significantly between populations (0.33 ± 0.10 to 1.31 ± 0.21) with an average of 0.74 ± 0.16. Body size and diapause expression were genetically correlated in transition-zone populations, but not in univoltine or bivoltine populations. Our findings support the suggestion that clinal variation in body size and diapause expression in the striped ground cricket reflect adaptation to season length.     

The association between inversion In(3R)Payne and clinally varying traits in Drosophila melanogaster

In Drosophila melanogaster, inversion In(3R)Payne increases in frequency towards low latitudes and has been putatively associated with variation in size and thermal resistance, traits that also vary clinally. To assess the association between size and inversion, we obtained isofemale lines of inverted and standard karyotype of In(3R)Payne from the ends of the Australian D. melanogaster east coast cline. In the northern population, there was a significant association between In(3R)Payne and body size, with standard lines from this population being relatively larger than inverted lines. In contrast, the inversion had no influence on development time or cold resistance. We strengthened our findings further in a separate study with flies from populations from the middle of the cline as well as from the cline ends. These flies were scored for wing size and the presence of In(3R)Payne using a molecular marker. In females, the inversion accounted for around 30% of the size difference between cline ends, while in males the equivalent figure was 60%. Adaptive shifts in size but not in the other traits are therefore likely to have involved genes closely associated with In(3R)Payne. Because the size difference between karyotypes was similar in different populations, there was no evidence for coadaptation within populations.     

PHENOTYPIC PLASTICITY AND REACTION NORMS IN TEMPERATE AND TROPICAL POPULATIONS OF DROSOPHILA MELANOGASTER: OVARIAN SIZE AND DEVELOPMENTAL TEMPERATURE

The plasticity of ovariole number relative to developmental temperature was studied in three populations of Drosophila melanogaster at both ends of the cline: a temperate French population and two equatorial Congolese. Ovary size was much greater in the French flies, in agreement with an already known latitudinal cline. Among isofemale lines, significant differences in genetic variability were observed between populations with a maximum variability at intermediate temperatures. Parameters of phenotypic variability (CV and FA) were not statistically different among lines or populations, but a significant increase at low temperature was demonstrated for both. The shapes of the response curves (i.e., the norm of reaction) were analyzed by adjusting the data to a quadratic equation. The parameters of the equation were highly variable among lines. On the other hand, the temperature for maximum value of ovarioles (TMV) was much less variable and exhibited only a slightly significant difference between temperate and tropical flies (22.2°C vs. 22.7°C). During its geographic extension toward colder places, D. melanogaster underwent a large, presumably adaptative, increase in ovariole number but very little change in the norm of reaction of that trait.     

Geographical variation in the genitalia of three species of Cicadellidae (Hemiptera)

Specimens of three species of Cicadellidae, Agallia laevis Ribaut, Eupteryx stachydearum Hardy and E.urticae Fabricius, are examined from different parts of their distribution in Europe and differences in genital characters are recorded. Variation in Agallia laevis suggests a dine or stepped cline from S.E. to N.W. in the shape of the anal appendages and aedeagus, while E. stachydearum shows a N–S cline, and a similar altitudinal cline, in the length and curvature of the apical processes of the aedeagus. In E.urticae geographical variation in the size and direction of the same processes also clearly occurs, though the exact pattern is not yet clear.     

Adaptation and constraint in the evolution of Drosophila melanogaster wing shape

SUMMARY We have taken advantage of parallel instances of natural selection on body size in Drosophila melanogaster to investigate constraints and adaptation affecting wing shape. Using recently developed techniques for statistical shape analysis, we have examined variation in wing shape in similar body size clines on three continents. Gender-related shape differences were constant among all populations, suggesting that gender differences represent a developmental constraint on wing shape. In contrast, the underlying shape varied significantly between continents and shape change within each cline (i.e., between small and large body size populations) also varied between continents. Therefore, variation at these two levels presumably results from either drift or natural selection. Functional considerations suggest that shape variation between the continents is unlikely to be adaptive. However, cline-related shape change, which we show has a significant allometric component, may be adaptive. The overall range of wing shape variation, across a large range of wing size, is extremely small, and the possibility that wing shape is subject to stabilizing selection (or canalization) is discussed.     

Within-island geographic variation of the dangerous Taiwanese snake, Trimeresurus stejnegeri, in relation to ecology

Multivariate analysis shows that geographic variation occurs in the scalation and body proportions of the dangerously venomous snake Tritneresurus stejnegeri , within the island of Taiwan. Canonical correlation analysis reveals a high correlation between constellations of morphological characters and ecological factors. Numerical hypothesis testing, using partial Mantel tests, provides evidence of a causal relationship with ecology. Head shape was found to be primarily associated with mean annual temperature and altitude, and head scalation with the annual range of temperature. Body scalation was found to be influenced by altitude, mean annual temperature and mean annual rainfall. Tail length is primarily associated with the range of temperature and rainfall and secondarily with the mean annual temperature and altitude. Geographic proximity was found to be less important in the explanation of the observed geographic pattern than some ecological factors, supporting the hypothesis of an ecogenetic origin of morphological variation.     

MOLECULAR ANALYSIS OF THE ALLELES OF ALCOHOL DEHYDROGENASE ALONG A CLINE IN DROSOPHILA MELANOGASTER. I. MAINE,NORTH CAROLINA,AND FLORIDA

Clinal variation in natural populations is often assumed to be due to the operation of natural selection. However, for many clines there exist plausible neutralist explanations which suggest that aspects of population structure maintain differences among subpopulations for particular traits. We used a restriction-mapping technique to investigate the contributions of population subdivision and selection to the maintenance of the allozyme polymorphism at the alcohol dehydrogenase (Adh) locus of Drosophila melanogaster. Digestions of genomic DNAs from 270 lines of flies by seven enzymes reveal 15–20% of all possible nucleotide substitutions and virtually all of the insertion/deletion variation in a 2.7-kilobase region containing the Adh structural locus. Analysis of large samples from each of three populations along the east coast of the United States provides evidence of founder effects in the most northerly population. Although there are signs of population differentiation among the samples, similarities between two of the populations indicate that migration among populations is extensive and strengthen the argument that natural selection plays a role in maintaining the cline.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

THE GENETIC ARCHITECTURE OF WING SIZE DIVERGENCE AT VARYING SPATIAL SCALES ALONG A BODY SIZE CLINE IN DROSOPHILA MELANOGASTER

Latitudinal clines in quantitative traits are common, but surprisingly little is known about the genetic bases of these divergences and how they vary within and between clines. Here, we use line‐cross analysis to investigate the genetic architecture of wing size divergences at varying spatial scales along a body size cline in Drosophila melanogaster. Our results revealed that divergences in wing size along the cline were due to strong additive effects. Significant nonadditive genetic effects, including epistasis and maternal effects, were also detected, but they were relatively minor in comparison to the additive effects and none were common to all crosses. There was no evidence of increased epistasis in crosses between more geographically distant populations and, unlike in previous studies, we found no significant dominance effects on wing size in any cross. Our results suggest there is little variation in the genetic control of wing size along the length of the Australian cline. They also highlight marked inconsistencies in the magnitude of dominance effects across studies, which may reflect different opportunities for mutation accumulation while lines are in laboratory culture.     

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.     

EXCLUSION OF THE ROLE OF SECONDARY CONTACT IN AN ALLELE FREQUENCY CLINE IN THE MUSSEL MYTILUS EDULIS

We examined the hypothesis that secondary contact generates an allele-frequency cline at the aminopeptidase-I locus (Lap) in the marine mussel, Mytilus edulis. It has been proposed that variation at the Lap locus is neutral and that the cline results from secondary contact between differentiated oceanic and estuarine populations (Levinton, 1980). We tested this hypothesis by comparing the genotypic distributions in samples from the cline to distributions that incorporate mixing effects. We employed a statistical model that determines the degree of contact using a maximum likelihood estimator and then incorporates the mixing estimates into an expected distribution of genotypes. Wahlund effects resulting from possible admixture are thereby incorporated into the expected distribution. Failure of the model to reconcile the observed with the expected distribution of genotypes indicates that the observed population structure does not result from admixture. The null hypothesis of mixing was unable to explain about 33% of the samples. Combined tests demonstrated the general departure from the mixing model to be highly significant. The distribution of heterozygote discrepancies across the cline was inconsistent with the expectations of a mixing model. Therefore we reject explanations for the structure of the Lap cline that involve secondary contact. Selection directed at the Lap locus appears necessary to explain the genotypic structure of clinal populations.     

Morphological evolution in the day gecko Phelsuma sundbergi in the Seychelles: a multivariate study

Geographic variation was investigated in populations of the day gecko Phelsuma sundbergi from 22 islands in the Seychelles, using multivariate ordination procedures. Multiple group principal components analysis was used to negate ontogenetic variation. Seventy-eight characters from three character systems (body proportions, scalation and colour pattern) were analysed from 349 specimens. Three phenetic aggregations of granitic island populations were detected from the northwestern, northeastern and southern island groups. A comparison of the results from the three character systems analysed separately suggests that the separation of the southern form preceded that of the northern forms. Colour pattern characters reflect this closely, while body proportions and scalation characters follow evolution of body size. The populations of Phelsuma sundbergi on the coralline islands are not well differentiated from the races on the granitic islands, and probably represent recent colonizations or introductions by man.     

SPECIES DISCRIMINATION AND EVOLUTIONARY MODE OF BUCHIA (BIVALVIA: BUCHIIDAE) FROM UPPER JURASSIC–LOWER CRETACEOUS STRATA OF GRASSY ISLAND, BRITISH COLUMBIA, CANADA

Abstract: Buchiid bivalves are geographically widespread in Upper Jurassic and Lower Cretaceous strata of the Northern Hemisphere. They are often abundant and their short stratigraphic ranges make them ideal biostratigraphic index fossils; these characteristics also render them useful for study of evolutionary patterns. We used multivariate methods to determine if we could discriminate between species of Buchia and examine how morphological characters change through time within the genus. Using ten morphological characters to describe shell shape and size, we tested for taxonomic differences and morphologic change in populations of buchiids collected from a single stratigraphic section on Grassy Island, located along the west coast of Vancouver Island, British Columbia. Morphometric analysis utilized traditional morphological metrics and techniques, including linear and angular measurements as well as Fourier (outline shape) analyses. Phenetic discrimination revealed considerable overlap among the recognized species in the morphospace, as well as a fairly low discriminatory power between species when compared as a group using a step‐wise canonical variate analysis. Step‐wise discriminant analyses between species pairs gave rise to much higher classification rates, suggesting that different characters are important for distinguishing between different species pairs. Our results also indicate that single individuals and small sample sizes of Buchia specimens are insufficient for biostratigraphic discrimination (unless other rarely preserved features such as the hinge and bysuss ear are available) and that a number of previously described species variants may not be taxonomically valid. A biolog using the multivariate axis that best discriminates between species (CV1) and a random walk‐based test using a Hurst estimate analysis indicate a gradualistic evolutionary mode for the Buchia species of Grassy Island. Shell shape and size of buchiids do not appear to be closely tied to lithofacies changes over the c. 10 myr time interval, suggesting that ecophenotypic variation (as it relates to substrate changes) probably had minimal influence on morphology.     

Latitudinal clines in body size,but not in thermal tolerance or heat‐shock cognate 70 (HSC70), in the highly‐dispersing intertidal gastropod Littorina keenae (Gastropoda: Littorinidae)

Natural populations of widely‐distributed animals often exhibit clinal variation in phenotypic traits or in allele frequencies of a particular gene over their geographical range. A planktotrophic intertidal snail, Littorina keenae is broadly distributed along the north‐eastern Pacific coast through a large latitudinal range (24°50′N–43°18′N). We tested for latitudinal clines in two complex phenotypic traits – thermal tolerance and body size – and one single locus trait – heat shock cognate 70 (HSC70) – in L. keenae along almost its entire geographical range. We found only weak evidence for a latitudinal cline in the thermal tolerance and no evidence for a cline in allele frequencies at HSC70. However, as predicted by Bergmann's rule, we detected a strong latitudinal cline that accounted for 60% of the variance in body size (R² = 0.598; P < 0.001). In contrast, body size did not significantly affect thermal tolerance. HSC70 showed no genetic differentiation among the populations, supporting our previous mitochondrial gene‐based estimate of high gene flow during this snail's free‐swimming larval stage. Given that L. keenae experiences panmixia along its species range, the observed size cline may be partially or entirely caused by a phenotypically plastic response to local thermal environments rather than by genetic divergence in body size among populations in response to locally optimizing natural selection. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 494–505.