LOCAL ADAPTATION MAINTAINS CLINAL VARIATION IN MELANIN‐BASED COLORATION OF EUROPEAN BARN OWLS (TYTO ALBA)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.     

Genetic drift widens the expected cline but narrows the expected cline width

Random genetic drift shifts clines in space, alters their width, and distorts their shape. Such random fluctuations complicate inferences from cline width and position. Notably, the effect of genetic drift on the expected shape of the cline is opposite to the naive (but quite common) misinterpretation of classic results on the expected cline. While random drift on average broadens the overall cline in expected allele frequency, it narrows the width of any particular cline. The opposing effects arise because locally, drift drives alleles to fixation--but fluctuations in position widen the expected cline. The effect of genetic drift can be predicted from standardized variance in allele frequencies, averaged across the habitat: . A cline maintained by spatially varying selection (step change) is expected to be narrower by a factor of √1- relative to the cline in the absence of drift. The expected cline is broader by the inverse of this factor. In a tension zone maintained by underdominance, the expected cline width is narrower by about 1- relative to the width in the absence of drift. Individual clines can differ substantially from the expectation, and we give quantitative predictions for the variance in cline position and width. The predictions apply to clines in almost one-dimensional circumstances such as hybrid zones in rivers, deep valleys, or along a coast line and give a guide to what patterns to expect in two dimensions.     

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.     

Population structure in the barn swallow,Hirundo rustica: a comparison between neutral DNA markers and quantitative traits

Local adaptation to variable environments can generate clinal variation in morphological traits. Alternatively, similar patterns of clinal variation may be generated simply as a result of genetic drift/migration balance. Teasing apart these different processes is a continuing focus in evolutionary ecology. We compare genetic differentiation at molecular loci and quantitative traits to analyse the effect of these different processes in a morphological latitudinal cline of the barn swallow, Hirundo rustica, breeding across Europe. The results obtained show no structuring at neutral microsatellite loci, which contrasts with positive structuring at five quantitative morphometric traits. This supports the hypothesis that the observed morphometric cline in barn swallows is the result of selection acting in a spatially heterogeneous environment. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 306–314.     

Divergent tail and throat ornamentation in the barn swallow across the Japanese islands

The geographic variations in male ornamentation provide insights into how different populations reach a different mean trait value under opposing forces of natural and sexual selection. Although the latitudinal cline of the elongated tail streamer, a sexually selected trait in the European subspecies of the barn swallow Hirundo rustica rustica, is a classic example, it has recently been shown that other subspecies of swallows have different targets of sexual selection. Here, we studied the latitudinal cline of ornamentation in the Asian subspecies, H. r. gutturalis, in which not the tail length but the white tail spot and red throat patch are important sexually selected traits. After controlling for covariates, the size of the white tail spot increased with latitude, while the size of the red throat patch decreased with latitude. On the other hand, we could not find any clear pattern regarding the elongated tail streamer, measured as fork depth. The divergent ornamentation across populations could be explained by latitudinal clines of sexually selected advantages of each ornament.     

An analysis of selection on a colour polymorphism in the northern leopard frog

In this study, we investigated the role of selection in the maintenance of a dorsal colour polymorphism in natural populations of the northern leopard frog, Rana pipiens. We determined genetic structure both spatially and temporally from a suite of putatively neutral molecular markers and tested whether or not the colour locus exhibited patterns of genetic variation that differed from those of the neutral loci. Spatial genetic structure at the colour locus was indistinguishable from structure at neutral loci [95% confidence intervals of F(ST) (neutral) = (0.07, 0.35), F(ST) (colour locus) = 0.114]. In the temporal analysis, we found that the variance among populations in the change in allele frequency at the colour locus (equal to 0.004) lies within the 95% confidence intervals for the variance among populations in changes in allele frequencies at neutral loci. In light of our inability to show evidence for the selective maintenance of the colour polymorphism, we used computer simulations to infer the power of our spatial and temporal techniques to detect selection. The computer simulations showed that although the strength of selection (s) would need to be relatively strong to have been detected by the temporal approach (s = 0.1-0.4, depending on the model tested), the spatial analysis would have detected all but weak selection (s = 0.01-0.04, depending on the model tested). This study illustrates the importance of using a locus comparison approach to detect evidence for selective maintenance before conducting studies to measure the selective mechanisms maintaining a polymorphism.     

Geographic variation in sexual dimorphism in the barn owl Tyto alba: a role for direct selection or genetic correlation?

Geographic variation in sexually selected traits is commonly attributed to geographic variation in the net benefit accrued from bearing such traits. Although natural and sexual selection are potentially important in shaping geographic variation, genetic constraints may also play a role. Although a genetic correlation between two traits may itself be the outcome of natural or sexual selection, it may indirectly reinforce the establishment and maintenance of cline variation with respect to one particular trait when across the cline different values of other traits are selected. Using the barn owl Tyto alba, a species in which the plumage of females is more reddish‐brown and more marked with black spots than that of males, I report results that are consistent with the hypothesis that both direct selection and genetic constraints may help establish and maintain cline variation in sexual dichromatism. In this species, inter‐individual variation in plumage coloration and spottiness has a genetic basis, and these traits are not sensitive to the environment. Data, based on the measurement of skin specimens, is consistent with the hypothesis that the stronger European cline variation in male spottiness than in female spottiness depends on the combined effects of (1) the similar cline variation in male and female plumage coloration and (2) the more intense phenotypic correlation between plumage coloration and spottiness in males (darker birds are more heavily spotted in the two sexes, but especially males) which is a general feature among the globally distributed barn owls. In northern Europe, male and female T. a. guttata are reddish‐brown and heavily spotted, and in southern Europe male and female T. a. alba are white, but only females display many spots. Here, I discuss the relative importance of direct selection, genetic correlation and the post‐ice age invasion of Europe by T. alba, in generating sex‐specific cline variation in plumage spottiness and non‐sex‐specific cline variation in plumage coloration.     

Spatial variation in the temporal change of male and female melanic ornamentation in the barn owl

Because the magnitude of selection can vary between sexes and in space and time, sexually antagonistic selection is difficult to demonstrate. In a Swiss population of barn owls (Tyto alba), a heritable eumelanic colour trait (size of black spots on ventral feathers) was positively selected with respect to yearling survival only in females. It remains unclear whether the absence of negative selection in males is typical in this species. To tackle this issue indirectly, we measured the size of black spots in 1733 skin specimens collected by museums from 1816 to 2001 in seven European countries and in the Middle-East. The temporal change in spot size was sex- and country-specific. In males, spots became smaller particularly in three countries (Middle-East, Italy and Switzerland). In females, the size of spots increased significantly in two countries (UK and Spain) and decreased in two others (Germany and Switzerland). Because migration and phenotypic plasticity cannot explain these results, selection is the most likely cause. The weaker temporal change in spot size in females than males may be because of the combined effect of strong genetic correlation between the sexes and stronger negative selection in males than positive selection in females. We thus suggest that in the barn owl, spot size (or genetically correlated traits) is sexually antagonistically selected and that its pattern of selection may account for the maintenance of its variation and sexual dimorphism.     

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.     

Mate‐choice copying: A fitness‐enhancing behavior that evolves by indirect selection

A spatially explicit, individual‐based simulation model is used to study the spread of an allele for mate‐choice copying (MCC) through horizontal cultural transmission when female innate preferences do or do not coevolve with a male viability‐increasing trait. Evolution of MCC is unlikely when innate female preferences coevolve with the trait, as copier females cannot express a higher preference than noncopier females for high‐fitness males. However, if a genetic polymorphism for innate preference persists in the population, MCC can evolve by indirect selection through hitchhiking: the copying allele hitchhikes on the male trait. MCC can be an adaptive behavior—that is, a behavior that increases a population's average fitness relative to populations without MCC—even though the copying allele itself may be neutral or mildly deleterious.     

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.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

Can Diversifying Selection Be Distinguished from History in Geographic Clines? A Population Genomic Study of Killifish (Fundulus heteroclitus)

A common geographical pattern of genetic variation is the one-dimensional cline. Clines may be maintained by diversifying selection across a geographical gradient but can also reflect historical processes such as allopatry followed by secondary contact. To identify loci that may be undergoing diversifying selection, we examined the distribution of geographical variation patterns across the range of the killifish (Fundulus heteroclitus) in 310 loci, including microsatellites, allozymes, and single nucleotide polymorphisms. We employed two approaches to detect loci under strong diversifying selection. First, we developed an automated method to identify clinal variation on a per-locus basis and examined the distribution of clines to detect those that exhibited signifcantly steeper slopes. Second, we employed a classic -outlier method as a complementary approach. We also assessed performance of these techniques using simulations. Overall, latitudinal clines were detected in nearly half of all loci genotyped (i.e., all eight microsatellite loci, 12 of 16 allozyme loci and 44% of the 285 SNPs). With the exception of few outlier loci (notably mtDNA and malate dehydrogenase), the positions and slopes of Fundulus clines were statistically indistinguishable. The high frequency of latitudinal clines across the genome indicates that secondary contact plays a central role in the historical demography of this species. Our simulation results indicate that accurately detecting diversifying selection using genome scans is extremely difficult in species with a strong signal of secondary contact; neutral evolution under this history produces clines as steep as those expected under selection. Based on these results, we propose that demographic history can explain all clinal patterns observed in F. heteroclitus without invoking natural selection to either establish or maintain the pattern we observe today.     

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.     

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.     

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.     

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.     

Temporal variation in glucocorticoid levels during the resting phase is associated in opposite way with maternal and paternal melanic coloration

Sex‐dependent selection can help maintain sexual dimorphism. When the magnitude of selection exerted on a heritable sex trait differs between the sexes, it may prevent each sex to reach its phenotypic optimum. As a consequence, the benefit of expressing a sex trait to a given value may differ between males and females favouring sex‐specific adaptations associated with different values of a sex trait. The level of metabolites regulated by genes that are under sex‐dependent selection may therefore covary with the degree of ornamentation differently in the two sexes. We investigated this prediction in the barn owl, a species in which females display on average larger black spots on the plumage than males, a heritable ornament. This melanin‐based colour trait is strongly selected in females and weakly counter‐selected in males indicating sex‐dependent selection. In nestling barn owls, we found that daily variation in baseline corticosterone levels, a key hormone that mediates life history trade‐offs, covaries with spot diameter displayed by their biological parents. When their mother displayed larger spots, nestlings had lower corticosterone levels in the morning and higher levels in the evening, whereas the opposite pattern was found with the size of paternal spots. Our study suggests a link between daily regulation of glucocorticoids and sex‐dependent selection exerted on sexually dimorphic melanin‐based ornaments.     

Isolation by distance,resistance and/or clusters? Lessons learned from a forest‐dwelling carnivore inhabiting a heterogeneous landscape

Landscape genetics provides a valuable framework to understand how landscape features influence gene flow and to disentangle the factors that lead to discrete and/or clinal population structure. Here, we attempt to differentiate between these processes in a forest‐dwelling small carnivore [European pine marten (Martes martes)]. Specifically, we used complementary analytical approaches to quantify the spatially explicit genetic structure and diversity and analyse patterns of gene flow for 140 individuals genotyped at 15 microsatellite loci. We first used spatially explicit and nonspatial Bayesian clustering algorithms to partition the sample into discrete clusters and evaluate hypotheses of ‘isolation by barriers’ (IBB). We further characterized the relationships between genetic distance and geographical (‘isolation by distance’, IBD) and ecological distances (‘isolation by resistance’, IBR) obtained from optimized landscape models. Using a reciprocal causal modelling approach, we competed the IBD, IBR and IBB hypotheses with each other to unravel factors driving population genetic structure. Additionally, we further assessed spatially explicit indices of genetic diversity using sGD across potentially overlapping genetic neighbourhoods that matched the inferred population structure. Our results revealed a complex spatial genetic cline that appears to be driven jointly by IBD and partial barriers to gene flow (IBB) associated with poor habitat and interspecific competition. Habitat loss and fragmentation, in synergy with past overharvesting and possible interspecific competition with sympatric stone marten (Martes foina), are likely the main factors responsible for the spatial genetic structure we observed. These results emphasize the need for a more thorough evaluation of discrete and clinal hypotheses governing gene flow in landscape genetic studies, and the potential influence of different limiting factors affecting genetic structure at different spatial scales.     

Clines in polygenic traits

This article outlines theoretical models of clines in additive polygenic traits, which are maintained by stabilizing selection towards a spatially varying optimum. Clines in the trait mean can be accurately predicted, given knowledge of the genetic variance. However, predicting the variance is difficult, because it depends on genetic details. Changes in genetic variance arise from changes in allele frequency, and in linkage disequilibria. Allele frequency changes dominate when selection is weak relative to recombination, and when there are a moderate number of loci. With a continuum of alleles, gene flow inflates the genetic variance in the same way as a source of mutations of small effect. The variance can be approximated by assuming a Gaussian distribution of allelic effects; with a sufficiently steep cline, this is accurate even when mutation and selection alone are better described by the 'House of Cards' approximation. With just two alleles at each locus, the phenotype changes in a similar way: the mean remains close to the optimum, while the variance changes more slowly, and over a wider region. However, there may be substantial cryptic divergence at the underlying loci. With strong selection and many loci, linkage disequilibria are the main cause of changes in genetic variance. Even for strong selection, the infinitesimal model can be closely approximated by assuming a Gaussian distribution of breeding values. Linkage disequilibria can generate a substantial increase in genetic variance, which is concentrated at sharp gradients in trait means.