Inversions fail to account for allozyme clines

Allozyme and inversion data from natural populations of Drosophila melanogaster from the eastern United States were analyzed to determine whether the clines at allozyme loci are due to nonrandom associations with common cosmopolitan inversions. All inversions show strong clines. Clines were large and significant for half of the eight allozyme loci. An analysis of the contribution of inversions to clines of allozyme genes revealed three outcomes: the inversion cline (1) enhanced the allozyme cline, but was only partly responsible, (2) reduced the allozyme cline, and (3) had no effect. The allozyme clines were mainly determined by the pattern of allele frequencies within the chromosomal arrangements. Consequently, it was concluded that allozyme clines would exist in the absence of inversion clines.     

Genetic Differentiation between Geographically Distant Populations of DROSOPHILA MELANOGASTER

We have studied allozyme variation at 26 gene loci in nine populations of Drosophila melanogaster originating on five different continents. The distant populations show significant genetic differentiation. However, only half of the loci studied have contributed to this differentiation; the other half show identical patterns in all populations. The genetic differentiation in North American, European and African populations is correlated with the major climatic differences between north and south. These differences arise mainly from seven loci that show gene-frequency patterns suggestive of latitudinal clines in allele frequencies. The clinal variation is such that subtropical populations are more heterozygous than temperate populations. These results are discussed in relation to the selectionist and neutralist hypotheses of genetic variation in natural populations.     

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.     

Cline coupling and uncoupling in a stickleback hybrid zone

Strong ecological selection on a genetic locus can maintain allele frequency differences between populations in different environments, even in the face of hybridization. When alleles at divergent loci come into tight linkage disequilibrium, selection acts on them as a unit and can significantly reduce gene flow. For populations interbreeding across a hybrid zone, linkage disequilibria between loci can force clines to share the same slopes and centers. However, strong ecological selection on a locus can also pull its cline away from the others, reducing linkage disequilibrium and weakening the barrier to gene flow. We looked for this “cline uncoupling” effect in a hybrid zone between stream resident and anadromous sticklebacks at two genes known to be under divergent natural selection (Eda and ATP1a1) and five morphological traits that repeatedly evolve in freshwater stickleback. These clines were all steep and located together at the top of the estuary, such that we found no evidence for cline uncoupling. However, we did not observe the stepped shape normally associated with steep concordant clines. It thus remains possible that these clines cluster together because their individual selection regimes are identical, but this would be very surprising given their diverse roles in osmoregulation, body armor, and swimming performance.     

Evidence for a link between local and seasonal cycles in gene frequencies and latitudinal gene clines in a cyclic parthenogen

In an earlier study (Rhombergh et al., Can. J. Genet. Cytol. 27: 224-232, 1985) of natural populations of the cyclic parthenogenetic Rose aphids, Macrosiphum rosae, 6 out of 31 loci were found to be polymorphic and one locus (Esterase-4) showed cyclic seasonal changes in gene and genotypic frequencies. Assuming that the Est-4 polymorphism was balanced and due to some climatic factor that varies seasonally, and realizing that most environmental factors that vary seasonally also vary latitudinally, we predicted existence of a latitudinal gene cline at this locus. In the present study we surveyed four polymorphic loci (chosen to be used as markers) in six geographic populations spanning over 1200 km between the United States and Canada and found all four loci to have latitudinal clines. We think that the gene clines are due to a latitudinal cline in the degree of advancement of local populations through the seasonal cycle, and have called such a pattern a 'seasonal phase cline'. The results are discussed in relation to the temporal instability of local patterns and persistence of genetic variability on the large scale in aphids. It is argued that population structure of aphids makes retention of selectively neutral or weakly selected polymorphisms difficult.     

Genetic population structure of the large blue butterfly Maculinea alcon in Denmark

We have investigated the genetic population structure within and the genetic differentiation between local populations of the large blue butterfly Maculinea alcon throughout the Jutland peninsula. Samples were collected as eggs on foodplants (Gentiana pneumonanthe), and reared to 4th instar caterpillars in the laboratory. A significant excess of homozygotes was found for all the investigated allozyme loci in most of the populations. A North-South cline was observed for the allele frequencies at some of the loci and for several linkage groups. Because some of the allele frequency clines were parallel to clines in adult morphological variation, we interpret our results as evidence for the co-existence of at least two gene pools within the Danish Maculinea alcon populations. Multilocus electrophoretic data revealed highly positive but variable FST values, which under this scenario would reflect varying frequencies of the Maculinea gene pools across the local populations. The significantly positive FIS values indicate that these gene pools are at least partly reproductively isolated (Wahlund effect). The co-occurrence of several Maculinea alcon gene pools on many local sites in Jutland is of great importance for conservation of the fragmented Maculinea populations. Our results show that there is probably more Maculinea biodiversity to conserve than was previously thought, and suggests that extant populations are more fragmented and vulnerable than counts of flying adults or eggs on foodplants indicate.     

B chromosome polymorphism in maize landraces: adaptive vs. demographic hypothesis of clinal variation

Cytogenetic analysis of maize landraces from northwestern Argentina has revealed an altitudinal cline in the mean number of B chromosomes (B's) per plant, with cultivars growing at higher altitudes exhibiting a higher number of B's. Altitudinal and longitudinal clines are frequently interpreted as evidence of selection, however, they can also be produced by the interplay between drift and spatially restricted gene flow or by admixture between previously isolated populations that have come into secondary contact. Here, we test the adaptive significance of the observed altitudinal gradient by comparing the levels of differentiation in the mean number of B's to those obtained from 18 selectively neutral loci [simple sequence repeats (SSRs)] among seven populations of the cline. The adequacy of alternative genetic-differentiation measures was determined, and associations between cytogenetic, genetic, and altitudinal distances were assessed by means of matrix- correspondence tests. No evidence for association between pairwise F(ST) and altitudinal distance or B-chromosome differentiation was found. The contrasting pattern of altitudinal divergence between the mean number of B's per plant and the genetic differentiation at SSR loci indicates that demographic processes cannot account for the observed levels of divergence in the mean number of B's.     

Biochemical genetics of Fundulus heteroclitus (L.). I. Temporal and spatial variation in gene frequencies of Ldh-B,Mdh-A,Gpi-B,and Pgm-A

Natural populations of Fundulus heteroclitus show extensive spatial variation in gene frequencies at four unlinked polymorphic loci. Large clinal changes in gene frequencies were found for Ldh-B, Mdh-A, and Gpi-B, whereas the spatial variation for the Pgm-B locus was small. Since the geographical area over which these clines are found is characterized by a steep thermal gradient, the clines in gene frequency are correlated with a directional change in mean water temperature. Maximum gene diversity of these four loci was correlated with annual fluctuations in water temperature. Temporal stability of the allelic frequencies was established for a 2–4 year period.This work was supported by NSF Grants GB37548 and DEB 76-19877 and by a grant from the National Geographic Society. A. R. P. is an NIH trainee supported by a training grant (No. HD00139) to the Department of Biology.This is contribution No. 960 of the Department of Biology, Johns Hopkins University.     

Steep,coincident, and concordant clines in mitochondrial and nuclear‐encoded genes in a hybrid zone between subspecies of Atlantic killifish,Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus‐specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non‐neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non‐neutral nuclear clines, SNPs in two nuclear‐encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin‐binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non‐unimodal, suggesting selection against advanced‐generation hybrids, possibly due to effects on processes involved in oxidative metabolism.     

Gene Flow and Selection in a Two-Locus System

A model of gene flow and selection in two linked loci is analyzed. The problems considered are the effects of linkage on the clines in frequencies at the two loci and the role of gene flow in producing linkage disequilibrium between the loci. Also, the possible significance of linkage as a mechanism for permitting a population of "track" spatial changes in the environment is considered. The results are that when the recombination fraction between the loci is of the same order of magnitude as the selection coefficients or smaller, then linkage is important in determining the gene frequencies and a substantial amount of linkage disequilibrium is present in the cline. Depending on the spatial pattern of selection on the two loci, linkage can either decrease or increase a population's response to local selection.     

Inferences of selection and migration in the Danish house mouse hybrid zone

We analysed the patterns of allele frequency change for ten diagnostic autosomal allozyme loci in the hybrid zone between the house mouse subspecies Mus musculus domesticus and M. m. musculus in central Jutland. After determining the general orientation of the clines of allele frequencies, we analysed the cline shapes along the direction of maximum gradient. Eight of the ten clines are best described by steep central steps with coincident positions and an average width of 8.9 km (support limits 7.6–12.4) flanked by tails of introgression, indicating the existence of a barrier to gene flow and only weak selection on the loci studied. We derived estimates of migration from linkage disequilibrium in the centre of the zone, and by applying isolation by distance methods to microsatellite data from some of these populations. These give concordant estimates of σ =  0.5–0.8 km generation     . The barrier to gene flow is of the order of 20 km (support limits 14–28), and could be explained by selection of a few per cent at 43–120 underdominant loci that reduces the mean fitness in the central populations to 0.45. Some of the clines appear symmetrical, whereas others are strongly asymmetrical, and two loci appear to have escaped the central barrier to gene flow, reflecting the differential action of selection on different parts of the genome. Asymmetry is always in the direction of more introgression into musculus , indicating either a general progression of domesticus into the musculus territory, possibly mediated by differential behaviour, or past movement of the hybrid zone in the opposite direction, impeded by potential geographical barriers to migration in domesticus territory.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 593–616.     

The deviation from linkage equilibrium with multiple loci varying in a stepping-stone cline

The balance between the creation of associations between alleles at different loci by immigration and the convergence to linkage equilibrium due to the recombination process is studied in a theoretical model. The geographical structure of the model is a stepping-stone chain of populations linking two genetically constant source populations. The model assumes an arbitrary number of autosomal loci and considers genetic variation (two alleles at each locus) that is not subject to natural selection. The gene frequencies at each locus will then show a linear cline through the stepping-stone chain of populations. The deviation from linkage equilibrium through the stepping-stone cline is characterized by an equation for linear measures that provide the linkage disequilibrium measures for a given set of loci in terms of the gene frequencies and the linkage disequilibria in the source populations and in terms of the linkage disequilibrium measures through the cline for lower numbers of loci. Numerical examples of this iterative solution are given, and it is shown that the build-up of the higher order Bennett-disequilibria through the cline is considerably more pronounced than the build-up of two-locus disequilibria.     

Evidence for a link between local and seasonal cycles in gene frequencies and latitudinal gene clines in a cyclic parthenogen

In an earlier study (Rhombergh et al., Can. J. Genet. Cytol. 27: 224–232, 1985) of natural populations of the cyclic parthenogenetic Rose aphids, Macrosiphum rosae, 6 out of 31 loci were found to be polymorphic and one locus (Esterase-4) showed cyclic seasonal changes in gene and genotypic frequencies. Assuming that the Est-4 polymorphism was balanced and due to some climatic factor that varies seasonally, and realizing that most environmental factors that vary seasonally also vary latitudinally, we predicted existence of a latitudinal gene cline at this locus. In the present study we surveyed four polymorphic loci (chosen to be used as markers) in six geographic populations spanning over 1200 km between the United States and Canada and found all four loci to have latitudinal clines. We think that the gene clines are due to a latitudinal cline in the degree of advancement of local populations through the seasonal cycle, and have called such a pattern a seasonal phase cline. The results are discussed in relation to the temporal instability of local patterns and persistence of genetic variability on the large scale in aphids. It is argued that population structure of aphids makes retention of selectively neutral or weakly selected polymorphisms difficult.     

Nonrandom larval dispersal can steepen marine clines

Sharp and stable clinal variation is enigmatic when found in species with high gene flow. Classical population genetic models treat gene flow as a random homogenizing force countering local adaptation across habitat discontinuities. Under this view, dispersal over large spatial scales will lower the effectiveness of adaptation by natural selection at finer spatial scales. Thus, random gene flow will create a shallow phenotypic cline across an ecotone in response to a steep selection gradient. In sedentary marine species that disperse primarily as larvae, nonrandom dispersal patterns are expected due to coastal hydrodynamics. Surprisingly sharp phenotypic and genotypic clines have been documented in marine species with high gene flow. We are interested in the extent to which nonrandom dispersal could accentuate such clines. We model a linear species range in which populations have stable and uniform densities along a selection gradient; in contrast to random dispersal, convergent advection of larvae can amplify phenotypic differentiation if coupled with a semipermeable dispersal barrier in the convergence zone. The migration load caused by directional dispersal pushes the phenotypic mean away from the local trait optimum in downstream populations, that is, near the convergence zone. A dispersal barrier is possible as a result of colliding currents if the water and larvae are mostly displaced offshore, away from suitable settlement habitat. Disjunctions in a quantitative trait were enlarged in the convergence zone by faster current flows or a more complete dispersal barrier. With advection of larvae per generation one-third as far as the average dispersal distance by diffusion, convergence on a dispersal barrier with 40% permeability generated a trait disjunction across the convergence zone of two phenotypic standard deviations. Without directional dispersal, similar clines also developed across a habitat gap, where population density was low, or across dispersal barriers with less than 1% permeability. These findings suggest that the types of hydrographic phenomena often associated with marine transition zones can strongly affect the balance between gene flow and selection and generate surprisingly steep clines given the large-scale gene flow expected from larvae.     

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.     

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.     

Y-chromosome polymorphisms and the origins of the European gene pool

Gradients of allele frequencies have long been considered the main genetic characteristic of the European population, but mitochondrial DNA diversity seems to be distributed differently. One Alu insertion (YAP), five tetranucleotide (DYS19, DYS389B, DYS390, DYS391 and DYS393) and one trinucleotide (DYS392) microsatellite loci of the Y chromosome were analysed for geographical patterns in 59 European populations. Spatial correlograms showed clines for most markers, which paralleled the gradients previously observed for two RFLP polymorphisms. Effective separation times between populations were estimated from genetic distances at microsatellite loci. Even after correcting for the possible effects of continuous local gene flow, the most distant Indo-European-speaking populations seem to have separated no more than 7000 years ago. The clinal patterns and the estimated, recent separation times between populations jointly suggest that Y-chromosome diversity in Europe largely reflects a directional demic expansion, which is unlikely to have occurred before the Neolithic period.     

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.     

Inference from Clines Stabilized by Frequency-Dependent Selection

Frequency-dependent selection against rare forms can maintain clines. For weak selection, s, in simple linear models of frequency-dependence, single locus clines are stabilized with a maximum slope of between square root of s/square root of 8 sigma and square root of s/square root of 12 delta, where sigma is the dispersal distance. These clines are similar to those maintained by heterozygote disadvantage. Using computer simulations, the weak-selection analytical results are extended to higher selection pressures with up to three unlinked genes. Graphs are used to display the effect of selection, migration, dominance, and number of loci on cline widths, speeds of cline movements, two-way gametic correlations ("linkage disequilibria"), and heterozygote deficits. The effects of changing the order of reproduction, migration, and selection, are also briefly explored. Epistasis can also maintain tension zones. We show that epistatic selection is similar in its effects to frequency-dependent selection, except that the disequilibria produced in the zone will be higher for a given level of selection. If selection consists of a mixture of frequency-dependence and epistasis, as is likely in nature, the error made in estimating selection is usually less than twofold. From the graphs, selection and migration can be estimated using knowledge of the dominance and number of genes, of gene frequencies and of gametic correlations from a hybrid zone.     

Estimates of Selection and Gene Flow from Measures of Cline Width and Linkage Disequilibrium in Heliconius Hybrid Zones

Hybrid zones can yield estimates of natural selection and gene flow. The width of a cline in gene frequency is approximately proportional to gene flow (sigma) divided by the square root of per-locus selection (square root of s). Gene flow also causes gametic correlations (linkage disequilibria) between genes that differ across hybrid zones. Correlations are stronger when the hybrid zone is narrow, and rise to a maximum roughly equal to s. Thus cline width and gametic correlations combine to give estimates of gene flow and selection. These indirect measures of sigma and s are especially useful because they can be made from collections, and require no field experiments. The method was applied to hybrid zones between color pattern races in a pair of Peruvian Heliconius butterfly species. The species are Müllerian mimics of one another, and both show the same changes in warning color pattern across their respective hybrid zones. The expectations of cline width and gametic correlation were generated using simulations of clines stabilized by strong frequency-dependent selection. In the hybrid zone in Heliconius erato, clines at three major color pattern loci were between 8.5 and 10.2 km wide, and the pairwise gametic correlations peaked at R approximately 0.35. These measures suggest that s approximately 0.23 per locus, and that sigma approximately 2.6 km. In erato, the shapes of the clines agreed with that expected on the basis of dominance. Heliconius melpomene has a nearly coincident hybrid zone. In this species, cline widths at four major color pattern loci varied between 11.7 and 13.4 km. Pairwise gametic correlations peaked near R approximately 1.00 for tightly linked genes, and at R approximately 0.40 for unlinked genes, giving s approximately 0.25 per locus and sigma approximately 3.7 km. In melpomene, cline shapes did not perfectly fit theoretical shapes based on dominance; this deviation might be explained by long-distance migration and/or strong epistasis. Compared with erato, sample sizes in melpomene are lower and the genetics of its color patterns are less well understood. In spite of these problems, selection and gene flow are clearly of the same order of magnitude in the two species. The relatively high per locus selection coefficients agree with "major gene" theories for the evolution of Müllerian mimicry, but the genetic architecture of the color patterns does not. These results show that the genetics and evolution of mimicry are still only sketchily understood.