Allele Frequency Clines Research Articles

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species.

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow-leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid-capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single-nucleotide polymorphisms (SNPs) for FST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.

Environmental Correlation Analysis for Genes Associated with Protection against Malaria.

Genome-wide searches for loci involved in human resistance to malaria are currently being conducted on a large scale in Africa using case-control studies. Here, we explore the utility of an alternative approach—“environmental correlation analysis, ECA,” which tests for clines in allele frequencies across a gradient of an environmental selection pressure—to identify genes that have historically protected against death from malaria. We collected genotype data from 12,425 newborns on 57 candidate malaria resistance loci and 9,756 single nucleotide polymorphisms (SNPs) selected at random from across the genome, and examined their allele frequencies for geographic correlations with long-term malaria prevalence data based on 84,042 individuals living under different historical selection pressures from malaria in coastal Kenya. None of the 57 candidate SNPs showed significant (P < 0.05) correlations in allele frequency with local malaria transmission intensity after adjusting for population structure and multiple testing. In contrast, two of the random SNPs that had highly significant correlations (P < 0.01) were in genes previously linked to malaria resistance, namely, CDH13, encoding cadherin 13, and HS3ST3B1, encoding heparan sulfate 3-O-sulfotransferase 3B1. Both proteins play a role in glycoprotein-mediated cell-cell adhesion which has been widely implicated in cerebral malaria, the most life-threatening form of this disease. Other top genes, including CTNND2 which encodes δ-catenin, a molecular partner to cadherin, were significantly enriched in cadherin-mediated pathways affecting inflammation of the brain vascular endothelium. These results demonstrate the utility of ECA in the discovery of novel genes and pathways affecting infectious disease.

Inference of biogeographical ancestry across central regions of Eurasia.

The inference of biogeographical ancestry (BGA) can provide useful information for forensic investigators when there are no suspects to be compared with DNA collected at the crime scene or when no DNA database matches exist. Although public databases are increasing in size and population scope, there is a lack of information regarding genetic variation in Eurasian populations, especially in central regions such as the Middle East. Inhabitants of these regions show a high degree of genetic admixture, characterized by an allele frequency cline running from NW Europe to East Asia. Although a proper differentiation has been established between the cline extremes of western Europe and South Asia, populations geographically located in between, i.e, Middle East and Mediterranean populations, require more detailed study in order to characterize their genetic background as well as to further understand their demographic histories. To initiate these studies, three ancestry informative SNP (AI-SNP) multiplex panels: the SNPforID 34-plex, Eurasiaplex and a novel 33-plex assay were used to describe the ancestry patterns of a total of 24 populations ranging across the longitudinal axis from NW Europe to East Asia. Different ancestry inference approaches, including STRUCTURE, PCA, DAPC and Snipper Bayes analysis, were applied to determine relationships among populations. The structure results show differentiation between continental groups and a NW to SE allele frequency cline running across Eurasian populations. This study adds useful population data that could be used as reference genotypes for future ancestry investigations in forensic cases. The 33-plex assay also includes pigmentation predictive SNPs, but this study primarily focused on Eurasian population differentiation using 33-plex and its combination with the other two AI-SNP sets.

Genetic footprints reveal geographic patterns of expansion in Fennoscandian red foxes.

Population expansions of boreal species are among the most substantial ecological consequences of climate change, potentially transforming both structure and processes of northern ecosystems. Despite their importance, little is known about expansion dynamics of boreal species. Red foxes (Vulpes vulpes) are forecasted to become a keystone species in northern Europe, a process stemming from population expansions that began in the 19th century. To identify the relative roles of geographic and demographic factors and the sources of northern European red fox population expansion, we genotyped 21 microsatellite loci in modern and historical (1835-1941) Fennoscandian red foxes. Using Bayesian clustering and Bayesian inference of migration rates, we identified high connectivity and asymmetric migration rates across the region, consistent with source-sink dynamics, whereby more recently colonized sampling regions received immigrants from multiple sources. There were no clear clines in allele frequency or genetic diversity as would be expected from a unidirectional range expansion from south to north. Instead, migration inferences, demographic models and comparison to historical red fox genotypes suggested that the population expansion of the red fox is a consequence of dispersal from multiple sources, as well as in situ demographic growth. Together, these findings provide a rare glimpse into the anatomy of a boreal range expansion and enable informed predictions about future changes in boreal communities.

The effective founder effect in a spatially expanding population.

The gradual loss of diversity and the establishment of clines in allele frequencies associated with range expansions are patterns observed in many species, including humans. These patterns can result from a series of founder events occurring as populations colonize previously unoccupied areas. We develop a model of an expanding population and, using a branching process approximation, show that spatial gradients reflect different amounts of genetic drift experienced by different subpopulations. We then use this model to measure the net average strength of the founder effect, and we demonstrate that the predictions from the branching process model fit simulation results well. We further show that estimates of the effective founder size are robust to potential confounding factors such as migration between subpopulations. We apply our method to data from Arabidopsis thaliana. We find that the average founder effect is approximately three times larger in the Americas than in Europe, possibly indicating that a more recent, rapid expansion occurred.

Individual performance in relation to cytonuclear discordance in a northern contact zone between long-toed salamander (Ambystoma macrodactylum) lineages.

Cytonuclear discordance in contact zones between related lineages is common, with mitochondrial clines often being displaced from clines in nuclear allele frequency. Proposed explanations for such a pattern include adaptive introgression of mtDNA or a neutral wake of mtDNA being left behind following hybrid zone movement. However, studies investigating these hypotheses are rare. Our previous survey of genetic variation in the long-toed salamander (Ambystoma macrodactylum) highlighted a potential case of cytonuclear discordance between two lineages in western Canada. Here, we use additional markers and samples to clarify the extent of this discordance. We simultaneously assess the feeding performance of individuals in a common environment to test for an association between mitotype and individual performance. The genetic results confirm a general pattern of cytonuclear discordance in the focal region. However, we also observed more limited introgression of a diagnostic nuclear marker. Intriguingly, although there were differences in individual performance associated with the transition between mitotypes, these differences were not fully explained by mitotype. Instead, the lowest performance was observed in individuals demonstrating the greatest mismatch between mtDNA and all nuclear markers, suggesting the potential for cytonuclear incompatibilities to be acting. These results highlight the complexity of understanding the causes and consequences of mtDNA introgression and cytonuclear discordance in contact zones.

Spatial Genetic Structure Patterns of Phenotype-Limited and Boundary-Limited Expanding Populations: A Simulation Study

Range expansions may create a unique spatial genetic pattern characterized by alternate genetically homogeneous domains and allele frequency clines. Previous attempts to model range expansions have mainly focused on the loss of genetic diversity during expansions. Using individual-based models, we examined spatial genetic patterns under two expansion scenarios, boundary-limited range expansions (BLRE) and phenotype-limited range expansions (PhLRE). Our simulation revealed that the genetic diversity within populations lost quickly during the range expansion, while the genetic difference accumulated between populations. Consequently, accompanying the expansions, the overall diversity featured a slow decrease. Specifically, during BLREs, high speed of boundary motion facilitated the maintenance of total genetic diversity and sharpened genetic clines. Very slight constraints on boundary motion of BLREs drastically narrowed the homogeneous domains and increased the allele frequency fluctuations from those levels exhibited by PhLREs. Even stronger constraints, however, surprisingly brought the width of homogeneous domains and the allele frequency fluctuations back to the normal levels of PhLREs. Furthermore, high migration rates maintained a higher total genetic diversity than low ones did during PhLREs. Whereas, the total genetic diversities during BLREs showed a contrary pattern: higher when migration was low than those when migration was high. Besides, the increase of migration rates helped maintain a greater number of homogeneous domains during PhLREs, but their effects on the number of homogeneous domains during BLREs were not monotonous. Previous studies have showed that the homogenous domains can merge to form a few broad domains as the expansion went on, leading to fewer homogeneous domains. Our simulations, meanwhile, revealed that the range expansions could also rebuild homogeneous domains from the clines during the range expansion. It is possible that that the number of homogeneous domains was determined by the interaction of merging and newly emerging homogeneous domains.

Morphological Clines and Weak Drift along an Urbanization Gradient in the Butterfly, Pieris rapae

Urban areas are increasing globally, providing opportunities for biodiversity researchers to study the process in which species become established in novel, highly disturbed habitats. This ecological process can be understood through analyses of morphological and genetic variation, which can shed light on patterns of neutral and adaptive evolution. Previous studies have shown that urban populations often diverge genetically from non-urban source populations. This could occur due to neutral genetic drift, but an alternative is that selection could lead to allele frequency changes in urban populations. The development of genome scan methods provides an opportunity to investigate these outcomes from samples of genetic variation taken along an urbanization gradient. Here we examine morphological variation in wing size and diversity at neutral amplified fragment length polymorphisms in the butterfly Pieris rapae L. (Lepidoptera, Pieridae) sampled from the center to the periphery of Marseille. We utilize established and novel environmental correlation approaches to scan genetic variation for evidence of selection. We find significant morphological differences in urban populations, as well as weak genetic structure and decreased genetic diversity in urban versus non-urban sites. However, environmental correlation tests provide little support for selection in our dataset. Our comparison of different methods and allele frequency clines suggests that loci identified as significant are false positives. Although there is some indication that selection may be acting on wing size in urban butterflies, genetic analyses suggest P. rapae are undergoing neutral drift.

Genomic consequences of multiple speciation processes in a stick insect

Diverse geographical modes and mechanisms of speciation are known, and individual speciation genes have now been identified. Despite this progress, genome-wide outcomes of different evolutionary processes during speciation are less understood. Here, we integrate ecological and spatial information, mating trials, transplantation data and analysis of 86 130 single nucleotide polymorphisms (SNPs) in eight populations (28 pairwise comparisons) of Timema cristinae stick insects to test the effects of different factors on genomic divergence in a system undergoing ecological speciation. We find patterns consistent with effects of numerous factors, including geographical distance, gene flow, divergence in host plant use and climate, and selection against maladaptive hybridization (i.e. reinforcement). For example, the number of highly differentiated ‘outlier loci’, allele-frequency clines and the overall distribution of genomic differentiation were recognizably affected by these factors. Although host use has strong effects on phenotypic divergence and reproductive isolation, its effects on genomic divergence were subtler and other factors had pronounced effects. The results demonstrate how genomic data can provide new insights into speciation and how genomic divergence can be complex, yet predictable. Future work could adopt experimental, mapping and functional approaches to directly test which genetic regions are affected by selection and determine their physical location in the genome.

Sequence variation patterns along a latitudinal cline in Scots pine (Pinus sylvestris): signs of clinal adaptation?

Adaptive clinal variation is abundant in nature, and its genetic basis is of great interest. The polygenic nature of this variation poses a challenge for identifying the causative loci underlying these adaptations. Here, we have examined the patterns of sequence variation in ten candidate genes for timing of bud set and cold tolerance, traits that form strong latitudinal clines in Scots pine (Pinus sylvestris). A set of reference genes was used for comparison and to infer a simplified demographic background model with approximate Bayesian computation methods. Against the resulting bottleneck model, the neutrality tests show little signs for local adaptation, but species-wide selection is suggested in some of the studied genes. In line with the theoretical expectations, we see little evidence for adaptive differentiation with F ST methods. Instead, allele frequency clines were found in three genes (dhn1, ftl2, and prr1). Our results add to the discussion on which molecular signals best characterize a gene subject to clinally varying selection. This will be especially relevant when these kinds of observations can be examined in parallel with association study results.

Evidence for strong assortative mating, limited gene flow, and strong differentiation across the blue‐footed/Peruvian booby hybrid zone in northern Peru

Hybrid zones represent natural laboratories in which the processes of divergence and genetic isolation can be examined. The generation and maintenance of a hybrid zone requires mispairing and successful reproduction between organisms that differ in one or more heritable traits. Understanding the dynamics of hybridization between two species requires an understanding of the extent to which they have diverged genetically, the frequency of misparing and hybrid production, and the extent of introgression. Three hundred and twenty one blue‐footed Sula nebouxii and Peruvian S. variegata boobies from the eastern tropical Pacific Ocean were analyzed using 19 putatively neutral genetic markers to evaluate interspecific differentiation, to classify morphological hybrids using Bayesian assignments, and to characterize hybridization using cline theory and Bayesian assignments. The species were well differentiated at mitochondrial and nuclear microsatellites, the hybrid zone was bimodal (contained a high frequency of each parental species but a low frequency of hybrids), and morphologically intermediate individuals were most likely F1 hybrids resulting from mating between female Peruvian boobies and male blue‐footed boobies. Clines in allele frequency could be constrained to share a common geographic centre but could not be constrained to share a common width. Peruvian and blue‐footed boobies hybridize infrequently, potentially due to strong premating reproductive isolation; however, backcrossing appears to facilitate introgression from blue‐footed to Peruvian boobies in this hybrid system.

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.

Population genetic structure of Penaeus monodon, in relation to monsoon current patterns in Southwest, East and Andaman coastal waters of India

The black tiger shrimp ( Penaeus monodon), a commercially important penaeid species, is widely distributed across the Indo-Pacific region. Genetic diversity in P. monodon collected from eight geographical regions in Southwest, East and Andaman coastal waters of India ( N = 418) was investigated using 10 polymorphic microsatellite loci. Average observed heterozygosity at sampled loci were high, ranging from 0.643 (Coromandel Coast) to 0.753 (South Andaman). Pairwise F ST (ranged from 0.005 to 0.078) and R ST (ranged from 0.005 to 0.171) estimates revealed surprisingly strong and statistically significant genetic structure among tiger shrimp populations. A synthetic map generated by multidimensional scaling shows an apparent cline in allele frequencies paralleling the roughly circular flow of surface currents in the Bay of Bengal. Significant heterozygote deficiencies were noted in most population samples at most loci. Andaman Island sites showed the highest diversity. Recognition of high genetic diversity and distinct population structuring of P. monodon in Indian seas has important implications for future domestication of this species in India, for two reasons: identification of the best wild founding stocks for aquaculture and, subsequently, the potential impacts of release of domesticates to the wild, either accidentally or deliberately (i.e. for stock enhancement).

Polymorphism in the couch potato gene clines in eastern Australia but is not associated with ovarian dormancy in Drosophila melanogaster

Natural selection can generate parallel latitudinal clines in traits and gene frequencies across continents, but these have rarely been linked. An amino acid (isoleucine to lysine, or I462K) polymorphism of the couch potato (cpo) gene in Drosophila melanogaster is thought to control female reproductive diapause cline in North America (Schmidt et al. 2008, Proc Natl Acad Sci USA, 105, 16207-16211). Here, we show that under standard diapause-inducing conditions (12 °C and short photoperiod) (Saunders et al. 1989, Proc Natl Acad Sci USA, 86, 3748-3752), egg maturation in Australian flies is delayed, but not arrested at previtellogenic stages. At 12 °C, the phenotypic distribution in egg development was bimodal at stages 8 and 14 and showed a strong nonlinear pattern on the east coast of Australia, with incidence of egg maturation delay (ovarian dormancy) increasing both toward tropical and temperate climates. Furthermore, we found no evidence for an association between the cpo I462K polymorphism and ovarian dormancy at either 12 or 10 °C (when egg maturation was often delayed at stage 7). Owing to strong linkage disequilibrium, the latitudinal cline in cpo allele frequencies was no longer evident once variation in the In(3R)P inversion polymorphism was taken into account. Our results suggest that the standard diapause-inducing conditions (12 °C and short photoperiod) were not sufficient to cause the typical previtellogenic developmental arrest in Australian flies and that the cpo I462K polymorphism does not explain the observed delay in egg development. In conclusion, ovarian dormancy does not show a simple latitudinal cline, and the lack of cpo-dormancy association suggests a different genetic basis to reproductive dormancy in North America and Australia.

Genetic Differentiation, Clinal Variation and Phenotypic Associations With Growth Cessation Across the Populus tremula Photoperiodic Pathway

Perennial plants monitor seasonal changes through changes in environmental conditions such as the quantity and quality of light. To ensure a correct initiation of critical developmental processes, such as the initiation and cessation of growth, plants have adapted to a spatially variable light regime and genes in the photoperiodic pathway have been implicated as likely sources for these adaptations. Here we examine genetic variation in genes from the photoperiodic pathway in Populus tremula (Salicaceae) for signatures diversifying selection in response to varying light regimes across a latitudinal gradient. We fail to identify any loci with unusually high levels of genetic differentiation among populations despite identifying four SNPs that show significant allele frequency clines with latitude. We do, however, observe large covariance in allelic effects across populations for growth cessation, a highly adaptive trait in P. tremula. High covariance in allelic effects is a signature compatible with diversifying selection along an environmental gradient. We also observe significantly higher heterogeneity in genetic differentiation among SNPs from the photoperiod genes than among SNPs from randomly chosen genes. This suggests that spatially variable selection could be affecting genes from the photoperiod pathway even if selection is not strong enough to cause individual loci to be identified as outliers. SNPs from three genes in the photoperiod pathway (PHYB2, LHY1, and LHY2) show significant associations with natural variation in growth cessation. Collectively these SNPs explain 10-15% of the phenotypic variation in growth cessation. Covariances in allelic effects across populations help explain an additional 5-7% of the phenotypic variation in growth cessation.

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 )–i nL. 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 2 = 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.

Connecting genetic variation to phenotypic clines

From early allozyme work to recent genome-wide scans, many studies have reported associations between molecular markers and latitude. These geographic patterns are tantalizing because they hint at the possibility of identifying specific mutations responsible for climatic adaptation. Unfortunately, few studies have done so because these exciting first glances often prove extremely challenging to follow up. Many difficulties can hinder connecting genetic and phenotypic variation in this context, and without such links, distinguishing the action of spatially varying selection from the other evolutionary processes capable of generating these patterns can be quite thorny. Nevertheless, two papers in this issue report excellent progress in overcoming these obstacles and provide persuasive evidence supporting the involvement of specific natural variants in clinal adaptation of Drosophila melanogaster populations (Fig. 1). In the first paper, Paaby et al. (2010) describe replicated allele frequency clines for a coding polymorphism in the Insulin-like Receptor (InR) gene on two continents, findings that strongly point to selection acting at this locus and that likely reflect life history adaptation. McKechnie et al. (2010) report compelling functional evidence that cis-regulatory variation in the Dca (drosophila cold acclimation) gene contributes to an adaptive cline in wing size. Notably, these papers employ largely alternative and complementary approaches, and together they exemplify how diverse strategies may be interwoven to draw convincing connections between genotype, phenotype, and evolutionary process.

Identification of a candidate adaptive polymorphism forDrosophilalife history by parallel independent clines on two continents

Life history traits are critical components of fitness and frequently reflect adaptive responses to environmental pressures. However, few genes that contribute to natural life history variation have been identified. Insulin signalling mediates the determination of life history traits in many organisms, and single gene manipulation in Drosophila melanogaster suggests that individual genes in the pathway have the potential to produce major effects on these quantitative traits. We evaluated allelic variation at two insulin signalling genes, the Insulin-like Receptor (InR) and its substrate, chico, in natural populations of D. melanogaster. We found different patterns of variation: InR shows evidence of positive selection and clines in allele frequency across latitude; chico exhibits neutral patterns of evolution. The clinal patterns at InR are replicated between North America and Australia, showing striking similarity in the distribution of specific alleles and the rate at which allele frequencies change across latitude. Moreover, we identified a polymorphism at InR that appears to be functionally significant and consistent with hypothetical patterns of selection across geography. This polymorphism provides new characterization of genic regions of functionality within InR, and is likely a component in a suite of genes and traits that respond adaptively to climatic variation.

The genetic structure of populations of an invading pest fruit fly, Bactrocera tryoni, at the species climatic range limit

Previous population genetic studies of the Queensland fruit fly, Bactrocera tryoni Froggatt (Diptera: Tephritidae), in its central range have shown barely detectable genetic differentiation across distances of almost 3000 km (F(st)=0.003). In this study, we investigated the genetic structuring of southern border populations of B. tryoni, in the region extending from the central population to the recently colonized southern range limit. The expectation was that marginal populations would be small, fragmented population sinks, with local adaptation limited by gene flow or drift. Unexpectedly, we found that the population at the southern extreme of the range was a source population, rather than a sink, for the surrounding region. This was shown by assignment testing of recent outbreaks in an adjoining quarantine area and by indirect migration estimates. Furthermore, populations in the region had formed a latitudinal cline in microsatellite allele frequencies, spanning the region between the central population and the southern range limit. The cline has formed within 250 generations of the initial invasion and appears stable between years. We show that there is restricted gene flow in the region and that effective population sizes are of the order of 10(2)-10(3). Although the cline may result from natural selection, neutral evolutionary processes may also explain our findings.

A Decreasing Gradient of 374F Allele Frequencies in the Skin Pigmentation Gene SLC45A2, from the North of West Europe to North Africa

The 374F mutation in the SLC45A2 gene, encoding the membrane-associated transporter protein (MATP) that plays an important role in melanin synthesis, has been suggested to be associated with skin color in Caucasians. In this study, the detailed distribution of the 374F allele has been investigated in 2063 unrelated subjects from 18 European and three North African populations. The highest allele frequency is observed in Denmark (0.980), and the lowest frequencies are observed in Tunisia (0.610) and in Morocco (0.691). A significant latitudinal cline in 374F allele frequencies was observed, ranging from the north of West Europe to North Africa (r = 0.869). The results confirm that the distribution of the 374F allele may reflect the ultraviolet radiation level and can be associated with skin color variation in these regions.