A comparison of nuclear and mitochondrial cline shapes in a hybrid zone in the Sceloporus grammicus complex (Squamata; Phrynosomatidae)

The F5 and FM2 chromosome races of the Sceloporus grammicus complex form a hybrid zone in the Mexican state of Hidalgo. Previous studies of this zone have assessed genetic structure by averaging estimates of shape and width across three diagnostic chromosome markers. This approach is likely to mask subtle differences in cline shape among loci (e.g. selected vs. neutral), and obscure any displacement of cline centres (if present). Here we use maximum likelihood methods to construct the best fitting individual clines for three chromosomal markers, and also add two new markers; the mitochondrial DNA (mtDNA) locus, and the nuclear ribosomal DNA (rDNA) repeat. For each locus, hybrid zone models were fitted by cline shape and width, and the position and number of segments describing the centre of the zone. Pairwise comparisons between all clines revealed concordance between chromosomes 2 and 6, but significant discordance in cline structure among all other paired combinations. The concordance of chromosomes 2 and 6 suggests that these clines are maintained by genome-wide forces. The discordance of the chromosome 1 cline suggests an influence of asymmetric introgression, while the mtDNA cline is probably influenced by selection and drift. The rDNA locus reveals a pattern best explained by either extreme asymmetric introgression or gene conversion. The structure of zone indicates that genome-wide processes and locus specific selective forces as well as drift, are operating to different degrees on different loci. The locus-by-locus approach used here permits a finer discrimination among possible mechanisms responsible for the maintenance of the individual clines.     

Numerical analysis of random drift in a cline

The equilibrium state of a diffusion model for random genetic drift in a cline is analyzed numerically. The monoecious organism occupies an unbounded linear habitat with constant, uniform population density. Migration is homogeneous, symmetric and independent of genotype. A single diallelic locus with a step environment is investigated in the absence of dominance and mutation. The flattening of the expected cline due to random drift is very slight in natural populations. The ratio of the variance of either gene frequency to the product of the expected gene frequencies decreases monotonically to a nonzero constant. The correlation between the gene frequencies at two points decreases monotonically to zero as the separation is increased with the average position fixed; the decrease is asymptotically exponential. The correlation decreases monotonically to a positive constant depending on the separation as the average position increasingly deviates from the center of the cline with the separation fixed. The correlation also decreases monotonically to zero if one of the points is fixed and the other is moved outward in the habitat, the ultimate decrease again being exponential. Some asymptotic formulae are derived analytically.——The loss of an allele favored in an environmental pocket is investigated by simulating a chain of demes exchanging migrants, the other assumptions being the same as above. For most natural populations, provided the allele would be maintained in the population deterministically, this process is too slow to have evolutionary importance.     

Numerical analysis of weak random drift in a cline

The equilibrium state of a diffusion model for weak random genetic drift in a cline is analyzed numerically. The monoecious organism occupies an unbounded linear habitat with constant, uniform population density. Migration is homogeneous, symmetric, and independent of genotype. A single diallelic locus with a step environment is investigated in the absence of dominance and mutation. The ratio of the variance of either gene frequency to the product of the expected gene frequencies decreases monotonically to a non-zero constant. The correlation between the gene frequencies at two points decreases monotonically to zero as the separation is increased with the average position fixed; the decrease is asymptotically exponential. The correlation decreases monotonically to a positive constant depending on the separation as the average position increasingly deviates from the center of the cline with the separation fixed. The correlation also decreases monotonically to zero if one of the points is fixed and the other is moved outward in the habitat, the ultimate decrease again being exponential. All the results are parameter free. Some asymptotic formulae are derived analytically.     

EVOLUTIONARY IMPLICATIONS OF DIVERGENT CLINES IN AN AVIAN (MANACUS: AVES) HYBRID ZONE

Abstract A previous study of the hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus) documented the unidirectional introgression of vitellinus male secondary sexual traits across the zone. Here, we examine the hybrid zone in greater genetic and morphological detail. Statistical comparisons of clines are performed using maximum‐likelihood and nonparametric bootstrap methods. Our results demonstrate that an array of six molecular and two morphometric markers agree in cline position and width. Clines for male collar and belly color are similar in width to the first eight clines, but are shifted in position by at least five cline widths. The result is that birds in intervening populations are genetically and morphometrically very like parental candei, but males have the plumage color of parental vitellinus. Neither neutral diffusion nor nonlinearity of color scales appear to be viable explanations for the large cline shifts. Genetic dominance of vitellinus plumage traits is another potential explanation that will require breeding experiments to test. Sexual selection remains a plausible explanation for the observed introgression of vitellinus color traits in these highly dimorphic, polygynous, lek‐mating birds. Two other clines, including a nondiagnostic isozyme locus, are similar in position to the main cluster of clines, but are broader in width. Thus, introgression at some loci is greater than that detected with diagnostic markers. Assuming that narrow clines are maintained by selection, variation in cline width indicates that selection is not uniform throughout the genome and that diagnostic markers are under more intense selective pressure. The traditional focus on diagnostic markers in studies of hybrid zones may therefore lead to underestimates of average introgression. This effect may be more pronounced in organisms with low levels of genetic divergence between hybridizing taxa.     

Classic clover cline clues

Adaptive clines are striking examples of natural selection in action, yet few have been studied in depth. In this issue of Molecular Ecology, Kooyers & Olsen (2012) introduce modern analyses and thinking towards studies of a classical example of the rapid and repeated evolution of latitudinal and altitudinal clines in cyanogenesis in white clover, Trifolium repens L. Recognizing that adaptive clines represent trade-offs in the selective benefits of traits at different ends of a geographical transect, these researchers focus on whether evidence for selection can be found at regional (coarse) and local (fine) scales. After adjusting for population genetic patterns generated by demographic processes, Kooyers and Olsen provide evidence that the cyanogenesis cline is adaptive across a transect from Louisiana to Wisconsin, USA. Within local populations, divergent selection on coupling dominant and recessive alleles that underlie cyanogenesis is predicted to drive populations to gametic phase disequilibrium (LD), a pattern that has been found in several other studies reviewed by Kooyers and Olsen. The absence of LD within any sampled populations in this study leads the authors to suggest that selective patterns within these clines may be more complex than previously proposed, perhaps even following theoretical predictions of a geographic mosaic.     

Evolution of additive and nonadditive genetic variance in development time along a cline in Drosophila serrata

Abstract. Latitudinal clines provide natural systems that may allow the effect of natural selection on the genetic variance to be determined. Ten clinal populations of Drosophila serrata collected from the eastern coast of Australia were used to examine clinal patterns in the trait mean and genetic variance of the life-history trait egg-to-adult development time. Development time significantly lengthened from tropical areas to temperate areas. The additive genetic variance for development time in each population was not associated with latitude but was associated with the population mean development time. Additive genetic variance tended to be larger in populations with more extreme development times and appeared to be consistent with allele frequency change. In contrast, the nonadditive genetic variance was not associated with the population mean but was associated with latitude. Levels of nonadditive genetic variance were greatest in the region of the cline where the gradient in the change in mean was greatest, consistent with Barton's (1999) conjecture that the generation of linkage disequilibrium may become an important component of the genetic variance in systems with a spatially varying optimum.     

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation—or clines—have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by‐product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.     

The use of AFLP to find an informative SNP: genetic differences across a migratory divide in willow warblers

We used the amplified fragment length polymorphism (AFLP) method to obtain genetic markers distinguishing two subspecies of willow warblers Phylloscopus trochilus that have different migratory behaviours but are not differentiated in mitochondrial DNA or at several microsatellite loci. With the inverse-polymerase chain reaction (PCR) approach we converted a dominant AFLP-marker to a codominant single nucleotide polymorphism (SNP). Across Scandinavia we typed 621 birds at the SNP locus AFLP-WW1 and we found a sigmoid change in allele frequencies centred around 62 degrees latitude. North of the latitudinal cline was a west-east cline. Both clines are narrower than one would expect from dispersal distances in willow warblers, which suggests that these are maintained by selection. The latitudinal cline at the locus AFLP-WW1 is paralleled by changes in several other traits, all of which might be maintained by a single selective force. The most plausible selection factor that we have identified is selection against hybrids because of inferior migratory behaviour. The selective force maintaining the east-west cline is less obvious. We discuss alternatives to the selection scenario, involving colonization history and asymmetric gene flow.     

Genetic Drift in a Cline

A model is developed of genetic drift in a cline maintained by spatially varying natural selection and local dispersal of individuals. The model is analyzed by an approximation scheme which is valid for weak selection and small migration rates. The results, which are based on numerical iterations of the approximate equations, are that the cline is less steep than predicted on the basis of the deterministic theory but that for weak selection the correlation between random fluctuations in neighboring colonies is approximately the same as in models of migration and drift in the absence of selection.     

Comparing clines on molecular and phenotypic traits in hybrid zones: a window on tension zone models

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.     

Parapatric diversification after post-glacial range expansion in the gall fly Urophora cardui (Tephritidae)

Aim Primary and secondary genetic clines in post-glacial colonized regions have different implications for biogeographic distributions and the origin of species. Primary clines arise in situ after colonization as adaptive responses to environmental gradients, while secondary clines are caused by contact between vicariant lineages. Here we analyse primary versus secondary origin of a genetic cline in the tephritid fly Urophora cardui in Jutland, Denmark, in a post-glacial landscape. Location Western Palaearctic. Methods Phylogeographic and demographic analyses of U. cardui based on mitochondrial DNA (mtDNA) genealogies, hierarchical genetic variance tests based on allozymes and distribution analysis of a rare allele from the Jutland cline. Results There was no phylogeographic divergence between the Jutland population of U. cardui north of the cline and neighbouring western European regional populations, which all shared the common western European mtDNA haplotype H1. At nuclear loci, by contrast, the North Jutland population was diverged above the mean level of divergence among regional populations and had no loss of genetic variation. A rare allozyme allele that was frequent in the cline area (up to 45%) and was missing north of the cline also occurred at low frequency (0–14%) elsewhere in the sampling range. Shallow phylogeographic divergence was observed between Russian and western European populations and between English and continental populations. Main conclusions The genetic variation patterns support primary cline evolution and parapatric divergence in Jutland following a demographic expansion of a western European ancestral source population of U. cardui, and suggest cryptic refugia and/or selection in other European population assemblages. The patterns of intra-specific regional divergence are discussed with respect to the interpretation of cryptic refugia in Europe after the most recent ice age.     

Variation in mitochondrial DNA in a cline of allele frequencies and shell phenotype in the dog-whelk Nucella lapillus (L.)

Genetic constitution in the intertidal gastropod Nucella lapillus (L.) influences shell shape, growth rate and physiology. Clinal variation in these traits along a 5 km stretch of coastline in south Devon can be related to environmental variation in temperature and desiccation stress. We have examined mtDNA variation along this shore to investigate whether the cline represents primary or secondary contact. Two distinct mtDNA haplotypes were found which exhibit coincident step clines with karyotypic, allozymic and phenotypic variation and covary with the environmental pressures of temperature and desiccation. These results are interpreted in the context of the wider scale distribution of genetic and phenotypic variation in N. lapillus. It is suggested that the shore studied may represent one of a number of regions of secondary contact within a mosaic hybrid zone in N. lapillus , where coadapted phenotypic variation correlates with habitat and the position of the clines represents an environmental transition.     

Genetic Drift in Clines Which Are Maintained by Migration and Natural Selection

Genetic drift will cause a migration-selection cline to wobble about its expected position. A rough linear approximation is developed, valid when local populations are large. This is used to calculate effects of genetic drift on clines in a stepping-stone model with abrupt and with gradual changes of selection coefficients at a single haploid locus. Among the quantities calculated are measures of slope, standardized variation of gene frequencies around their expected values, and correlation among neighboring populations with respect to deviations from the expected gene frequencies. These quantities appear to be primarily functions of Ns and Nm for a given pattern of selection. Computer simulation gives rough confirmation of these results. Standardized variances of gene frequencies and correlation of neighbors differ along the cline in the case of smooth changes in selection. In no case is pathological behavior of gene frequency deviations found near the boundaries of selective regions. Local behavior of gene frequences of nearby colonies is approximately predicted by a simple adaptation of the stepping-stone theory of Kimura and Weiss. Approximate measures of the lateral variation of the midpoint of a cline and the probability of non-monotonicity are also calculated and discussed.     

Different cell size and cell number contribution in two newly established and one ancient body size cline of Drosophila subobscura

Latitudinal genetic clines in body size occur in many ectotherms including Drosophila species. In the wing of D. melanogaster, these clines are generally based on latitudinal variation in cell number. In contrast, differences in wing area that evolve by thermal selection in the laboratory are in general based on cell size. To investigate possible reasons for the different cellular bases of these two types of evolutionary response, we compared the newly established North and South American wing size clines of Drosophila subobscura. The new clines are based on latitudinal variation in cell area in North America and cell number in South America. The ancestral European cline is also based on latitudinal variation in cell number. The difference in the cellular basis of wing size variation in the American clines, which are roughly the same age, together with the similar cellular basis of the new South American cline and the ancient European one, suggest that the antiquity of a cline does not explain its cellular basis. Furthermore, the results indicate that wing size as a whole, rather than its cellular basis, is under selection. The different cellular bases of different size clines are most likely explained either entirely by chance or by different patterns of genetic variance--or its expression--in founding populations.     

Little and large: body size and genetic clines in a New Zealand gecko (Woodworthia maculata) along a coastal transect

Clinal variation can result from primary differentiation or secondary contact and determining which of these two processes is responsible for the existence of a cline is not a trivial problem. Samples from a coastal transect of New Zealand geckos (Woodworthia maculata) identified for the first time a body size cline 7-10 km wide. The larger geckos are almost twice the mass of the small adult geckos. Clines in allele and haplotype frequency were found at two of the four genetic loci examined. Estimated width of the morphological cline was concordant with neither the narrower mtDNA cline (3-7 m) nor the wider nuclear cline (RAG-2; 34-42 km), and cline centers were not coincident. Although the body size cline is narrow compared to the entire range of the species, it is 2-3 orders of magnitude greater than estimates of dispersal distance per generation for these geckos. No evidence of assortative mating, nor of hybrid disadvantage was identified, thus there is little evidence to infer that endogenous selection is maintaining a hybrid zone. We cannot distinguish secondary contact from primary origin of this body size cline but conclude that secondary contact is likely due to the occurrence of mtDNA haplotypes from three distinct clades within the coastal transect and the presence of two frequency clines within this region.     

Geographic distributions of Idh -1 alleles in a cricket are linked to differential enzyme kinetic performance across thermal environments

Background  

Geographic clines within species are often interpreted as evidence of adaptation to varying environmental conditions. However, clines can also result from genetic drift, and these competing hypotheses must therefore be tested empirically. The striped ground cricket, Allonemobius socius, is widely-distributed in the eastern United States, and clines have been documented in both life-history traits and genetic alleles. One clinally-distributed locus, isocitrate dehydrogenase (Idh -1), has been shown previously to exhibit significant correlations between allele frequencies and environmental conditions (temperature and rainfall). Further, an empirical study revealed a significant genotype-by-environmental interaction (GxE) between Idh -1 genotype and temperature which affected fitness. Here, we use enzyme kinetics to further explore GxE between Idh -1 genotype and temperature, and test the predictions of kinetic activity expected under drift or selection.     

Inference of selection and stochastic effects in the house mouse hybrid zone

We explored the transition of 13 X-linked markers across two separate portions of the house mouse hybrid zone, asking whether such a comparison can distinguish the effects of selection from random factors. A heuristic search in the likelihood landscape revealed more complex likelihood profiles for data sampled in two-dimensional (2D) space relative to data sampled along a linear transect. Randomized resampling of localities analyzed for individual loci showed that deletion of sites away from the zone center can decrease cline width estimates whereas deletion of sites close to the center can significantly increase the width estimates. Deleting localities for all loci resulted in wider clines if the number of samples from the center was limited. The results suggest that, given the great variation in width estimates resulting from inclusion/exclusion of sampling sites, the geographic sampling design is important in hybrid zone studies and that our inferences should take into account measures of uncertainty such as support intervals. The comparison of the two transects indicates cline widths are narrower for loci in the central part of the X chromosome, suggesting selection is stronger in this region and genetic incompatibilities may have at least partly common architecture in the house mouse hybrid zone.     

Adaptation to a steep environmental gradient and an associated barrier to gene exchange in Littorina saxatilis

Steep environmental gradients offer important opportunities to study the interaction between natural selection and gene flow. Allele frequency clines are expected to form at loci under selection, but unlinked neutral alleles may pass easily across these clines unless a generalized barrier evolves. Here we consider the distribution of forms of the intertidal gastropod Littorina saxatilis, analyzing shell shape and amplified fragment length polymorphism (AFLP) loci on two rocky shores in Britain. On the basis of previous work, the AFLP loci were divided into differentiated and undifferentiated groups. On both shores, we have shown a sharp cline in allele frequencies between the two morphs for differentiated AFLP loci. This is coincident with a habitat transition on the shore where the two habitats (cliff and boulder field) are immediately contiguous. The allele frequency clines coincide with a cline in shell morphology. In the middle of the cline, linkage disequilibrium for the differentiated loci rises in accordance with expectation. The clines are extremely narrow relative to dispersal, probably as a result of both strong selection and habitat choice. An increase in F(ST) for undifferentiated AFLPs between morphs, relative to within-morph comparisons, is consistent with there being a general barrier to gene flow across the contact zone. These features are consistent either with an episode of allopatric divergence followed by secondary contact or with primary, nonallopatric divergence. Further data will be needed to distinguish between these alternatives.     

Genetic architecture of a body colour cline in Drosophila americana

Phenotypic variation within a species is often structured geographically in clines. In Drosophila americana, a longitudinal cline for body colour exists within North America that appears to be due to local adaptation. The tan and ebony genes have been hypothesized to contribute to this cline, with alleles of both genes that lighten body colour found in D. americana. These alleles are similar in sequence and function to the allele fixed in D. americana's more lightly pigmented sister species, Drosophila novamexicana. Here, we examine the frequency and geographic distribution of these D. novamexicana‐like alleles in D. americana. Among alleles from over 100 strains of D. americana isolated from 21 geographic locations, we failed to identify additional alleles of tan or ebony with as much sequence similarity to D. novamexicana as the D. novamexicana‐like alleles previously described. However, using genetic analysis of 51 D. americana strains derived from 20 geographic locations, we identified one new allele of ebony and one new allele of tan segregating in D. americana that are functionally equivalent to the D. novamexicana allele. An additional 5 alleles of tan also showed marginal evidence of functional similarity. Given the rarity of these alleles, however, we conclude that they are unlikely to be driving the pigmentation cline. Indeed, phenotypic distributions of the 51 backcross populations analysed indicate a more complex genetic architecture, with diversity in the number and effects of loci altering pigmentation observed both within and among populations of D. americana. This genetic heterogeneity poses a challenge to association studies and genomic scans for clinal variation, but might be common in natural populations.