Latitudinal clines: an evolutionary view on biological rhythms

Properties of the circadian and annual timing systems are expected to vary systematically with latitude on the basis of different annual light and temperature patterns at higher latitudes, creating specific selection pressures. We review literature with respect to latitudinal clines in circadian phenotypes as well as in polymorphisms of circadian clock genes and their possible association with annual timing. The use of latitudinal (and altitudinal) clines in identifying selective forces acting on biological rhythms is discussed, and we evaluate how these studies can reveal novel molecular and physiological components of these rhythms.     

Chromosomal inversion polymorphisms and adaptation

Chromosomal inversion polymorphisms continue to be identified from an increasing number of populations of insects, plants, bacteria and humans. In the fruit fly Drosophila, chromosomal polymorphisms were used in classic studies of natural selection. Recent molecular genetic studies suggest that inversion polymorphisms are dynamical systems. These studies also indicate patterns of disequilibrium and variation that are consistent with co-adapted gene complexes. Although these complexes have yet to be identified, recent studies have identified traits, such as body size, that are linked to inversion polymorphisms. Selection acting on these polymorphisms is strong because latitudinal clines in inversion frequency become re-established rapidly after a new continent is colonized. A combined molecular and phenotypic approach is helping to identify the role of inversion polymorphisms in adaptive divergence, but the genes responsible for associations between traits and inversion polymorphisms remain to be identified.     

Estimates of gene flow from rare alleles in natural populations of medfly Ceratitis capitata (Diptera: Tephritidae)

Gene flow based on the spatial distribution of rare alleles at 25 gene loci was estimated in 15 populations of Ceratitis capitata (Wiedemann) from different parts of the world. Estimates of Nm, the number of migrants exchanged per generation among populations in different regions of the world, appeared to be quite similar, ranging from 3.36 in tropical Africa to 2.94 in the New World and 2.72 in Mediterranean basin populations. This suggests that gene flow among neighbouring populations of medfly is quite extensive. The genetic differentiation in American, Mediterranean and African populations was related to major climatic differences between North and South. These differences arise mainly from five loci that showed gene frequency patterns suggestive of latitudinal clines in allele frequencies. The clinal variation was such that tropical-subtropical populations were more heterozygous than temperate populations. It was concluded that gene flow, counteracting the forces of natural selection and genetic drift, determines the extent to which geographical populations of C. capitata are differentiated.     

Patterns of male sterility in a grasshopper hybrid zone imply accumulation of hybrid incompatibilities without selection

It is now widely accepted that post-zygotic reproductive isolation is the result of negative epistatic interactions between derived alleles fixed independently at different loci in diverging populations (the Dobzhansky-Muller model). What is less clear is the nature of the loci involved and whether the derived alleles increase in frequency through genetic drift, or as a result of natural or sexual selection. If incompatible alleles are fixed by selection, transient polymorphisms will be rare and clines for these alleles will be steep where divergent populations meet. If they evolve by drift, populations are expected to harbour substantial genetic variation in compatibility and alleles will introgress across hybrid zones once they recombine onto a genetic background with which they are compatible. Here we show that variation in male sterility in a naturally occurring Chorthippus parallelus grasshopper hybrid zone conforms to the neutral expectations. Asymmetrical clines for male sterility have long tails of introgression and populations distant from the zone centre show significant genetic variation for compatibility. Our data contrast with recent observations on 'speciation genes' that have diverged as a result of strong natural selection.     

Clines in clock genes: fine-tuning circadian rhythms to the environment

The dissection of the circadian clock into its molecular components represents the most striking and well-studied example of a gene regulatory network underlying a complex behavioural trait. By contrast, the evolutionary analysis of the clock has developed more slowly. Here we review studies that have surveyed intraspecific clock gene variation over large geographical areas and have discovered latitudinal clines in gene frequencies. Such spatial patterns traditionally suggest that natural selection shapes genetic variation, but it is equally possible that population history, or a mixture of demography and selection, could contribute to the clines. We discuss how population genetics, together with functional assays, can illuminate these possible cases of natural selection in Drosophila clock genes.     

The Genetic Content of Chromosomal Inversions across a Wide Latitudinal Gradient

There is increasing evidence regarding the role of chromosomal inversions in relevant biological processes such as local adaptation and speciation. A classic example of the adaptive role of chromosomal polymorphisms is given by the clines of inversion frequencies in Drosophila subobscura, repeatable across continents. Nevertheless, not much is known about the molecular variation associated with these polymorphisms. We characterized the genetic content of ca. 600 individuals from nine European populations following a latitudinal gradient by analysing 19 microsatellite loci from two autosomes (J and U) and the sex chromosome (A), taking into account their chromosomal inversions. Our results clearly demonstrate the molecular genetic uniformity within a given chromosomal inversion across a large latitudinal gradient, particularly from Groningen (Netherlands) in the north to Málaga (Spain) in the south, experiencing highly diverse environmental conditions. This low genetic differentiation within the same gene arrangement across the nine European populations is consistent with the local adaptation hypothesis for th evolutionof chromosomal polymorphisms. We also show the effective role of chromosomal inversions in maintaining different genetic pools within these inverted genomic regions even in the presence of high gene flow. Inversions represent thus an important barrier to gene flux and can help maintain specific allelic combinations with positive effects on fitness. Consistent patterns of microsatellite allele-inversion linkage disequilibrium particularly in loci within inversions were also observed. Finally, we identified areas within inversions presenting clinal variation that might be under selection.     

Interactions of polymorphisms in different clock genes associated with circadian phenotypes in humans

Several studies have shown that mutations and polymorphisms in clock genes are associated with abnormal circadian parameters in humans and also with more subtle non-pathological phenotypes like chronotypes. However, there have been conflicting results, and none of these studies analyzed the combined effects of more than one clock gene. Up to date, association studies in humans have focused on the analysis of only one clock gene per study. Since these genes encode proteins that physically interact with each other, combinations of polymorphisms in different clock genes could have a synergistic or an inhibitory effect upon circadian phenotypes. In the present study, we analyzed the combined effects of four polymorphisms in four clock genes (Per2, Per3, Clock and Bmal1) in people with extreme diurnal preferences (morning or evening). We found that a specific combination of polymorphisms in these genes is more frequent in people who have a morning preference for activity and there is a different combination in individuals with an evening preference for activity. Taken together, these results show that it is possible to detect clock gene interactions associated with human circadian phenotypes and bring an innovative idea of building a clock gene variation map that may be applied to human circadian biology.     

Antigen and haplotype frequencies at three human leucocyte antigen loci (HLA-A, -B, -C) in the Pawaia of Papua New Guinea

The genetic profile of the Pawaia, a seminomadic, linguistic isolate from the highlands fringe of Papua New Guinea, is described in terms of antigen and haplotype frequencies at three class I human leucocyte antigen loci (HLA-A, -B, and -C). The Pawaia, like other Papua New Guinea populations, exhibit restricted polymorphisms at all three loci studied, both in the number of alleles segregating and in the level of average heterozygosity. An extremely high frequency (52.9%) of HLA-B27, the antigen implicated in the pathogenesis of seronegative arthropathies, was found. A novel HLA-C locus specificity, CNG, resulting probably from a gene duplication event, was also observed in significant numbers. Although the gene frequency comparisons suggest their strong affinities with the highlanders, the Pawaia haplotypes reveal significant admixture from other neighbouring groups as well. The usefulness of HLA haplotypes in tracing the movements of human populations in the New Guinea area is discussed.     

Classic clover cline clues

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

Single-locus latitudinal clines and their relationship to temperate adaptation in metabolic genes and derived alleles in Drosophila melanogaster

We report a study in Drosophila melanogaster of latitudinal clines for 23 SNPs embedded in 13 genes (Pgi, Gapdh1, UGPase, Pglym78, Pglym87, Eno, Men, Gdh, Sod, Pgk, Mdh1, TreS, Treh) representing various metabolic enzymes. Our samples are from 10 populations spanning latitude from southern Florida to northern Vermont. Three new clines with latitude were detected. These are the amino acid polymorphisms in the NAD-dependent glutamate dehydrogenase (Gdh) and trehalase (Treh) genes, and a silent site polymorphism in the UDP-glucose pyrophosphorylase gene (UGPase). The result, when combined with the overall incidence and pattern of reports for six other genes (Adh, Gpdh, Pgm, G6pd, 6Pgd, Hex-C), presents a picture of latitudinal clines in metabolic genes prevalent around the branch point of competing pathways. For six of the seven amino acid polymorphisms showing significant latitudinal clines in North America, the derived allele is the one increasing with latitude, suggesting temperate adaptation. This is consistent with a model of an Afrotropical ancestral species adapting to temperate climates through selection favoring new mutations.     

Re-establishment of clinal variation in flowering time among introduced populations of purple loosestrife (Lythrum salicaria, Lythraceae)

Range expansion during biological invasion requires that invaders adapt to geographical variation in climate, which should yield latitudinal clines in reproductive phenology. We investigated geographic variation in life history among 25 introduced populations of Lythrum salicaria, a widespread European invader of North American wetlands. We detected a strong latitudinal cline in initiation of flowering and size at flowering, which paralleled that reported among native populations. Plants from higher latitudes flowered earlier and at a smaller size than those from lower latitudes, even when raised in a uniform glasshouse. Early flowering was associated with greatly reduced reproductive output, but this was not associated with latitudinal variation in abundance, and probably did not result from a genetic correlation between time to and size at flowering. As introduction to North America c. 200 years ago, L. salicaria has re-established latitudinal clines in life history, probably as an evolutionary response to climatic selection.     

Patterns of variation within and between Greek populations of Ceratitis capitata suggest extensive gene flow and latitudinal clines

Four natural Greek populations of Mediterranean fruit fly, Ceratitis capitata (Wiedemann), was studied for genetic variability at 25 enzyme loci. The comparison of polymorphism within and between populations shows two loci with high between-population heterozygosity (HT) and very high fixation index (F(ST)) values, suggesting the presence of balancing selection. The gradual decline of common allele frequency of the polymorphic loci tested indicated that latitudinal clines are present in Greece. Indirect estimates of gene flow based both on Wright's method (Nm*) and on the Slatkin's method (Nm*), which depends on the frequencies of rare alleles found in only one population, revealed a substantial amount of gene flow (Nm = 3.493 and Nm* = 3.197). These estimates of gene flow may well explain why the "introduced" Greek populations of C. capitata, in spite of their low genetic variability, display the same polymorphic loci. Gene flow in combination with natural selection and genetic drift may have played an important role to genetic differentiation in this species in Greece.     

Avian Clock gene polymorphism: evidence for a latitudinal cline in allele frequencies

In comparison with most animal behaviours, circadian rhythms have a well-characterized molecular genetic basis. Detailed studies of circadian clock genes in 'model' organisms provide a foundation for interpreting the functional and evolutionary significance of polymorphic circadian clock genes found within free-living animal populations. Here, we describe allelic variation in a region of the avian Clock orthologue which encodes a functionally significant polyglutamine repeat (ClkpolyQcds), within free-living populations of two passerine birds, the migratory bluethroat (Luscinia svecica) and the predominantly nonmigratory blue tit (Cyanistes caeruleus). Multiple ClkpolyQcds alleles were found within populations of both species (bluethroat: 12 populations, 7 alleles; blue tit: 14 populations, 9 alleles). Some populations of both species were differentiated at the ClkpolyQcds locus as measured by F(ST) and R(ST) values. Among the blue tit, but not bluethroat populations, we found evidence of latitudinal clines in (i) mean ClkpolyQcds repeat length, and (ii) the proportions of three ClkpolyQcds genotype groupings. Parallel analyses of microsatellite allele frequencies, which are considered to reflect selectively neutral processes, indicate that interpopulation allele frequency variation at the ClkpolyQcds and microsatellite loci does not reflect the same underlying demographic processes. The possibility that the observed interpopulation ClkpolyQcds allele frequency variation is, at least in part, maintained by selection for microevolutionary adaptation to photoperiodic parameters correlated with latitude warrants further study.     

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long‐standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion‐size and inversion‐shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.     

Latitudinal clines in<Emphasis Type="Italic">Drosophila melanogaster</Emphasis>: Body size,allozyme frequencies,inversion frequencies,and the insulin-signalling pathway

Many latitudinal clines exist inDrosophila melanogaster: in adult body size, in allele frequency at allozyme loci, and in frequencies of common cosmopolitan inversions. The question is raised whether these latitudinal clines are causally related. This review aims to connect data from two very different fields of study, evolutionary biology and cell biology, in explaining such natural genetic variation inD. melanogaster body size and development time. It is argued that adult body size clines, inversion frequency clines, and clines in allele frequency at loci involved in glycolysis and glycogen storage are part of the same adaptive strategy. Selection pressure is expected to differ at opposite ends of the clines. At high latitudes, selection onD. melanogaster would favour high larval growth rate at low temperatures, and resource storage in adults to survive winter. At low latitudes selection would favour lower larval critical size to survive crowding, and increased male activity leading to high male reproductive success. Studies of the insulin-signalling pathway inD. melanogaster point to the involvement of this pathway in metabolism and adult body size. The genes involved in the insulin-signalling pathway are associated with common cosmopolitan inversions that show latitudinal clines. Each chromosome region connected with a large common cosmopolitan inversion possesses a gene of the insulin transmembrane complex, a gene of the intermediate pathway and a gene of the TOR branch. The hypothesis is presented that temperateD. melanogaster populations have a higher frequency of a ’thrifty’ genotype corresponding to high insulin level or high signal level, while tropical populations possess a more ’spendthrift’ genotype corresponding to low insulin or low signal level.     

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.     

Genetic adaptation of the human circadian clock to day-length latitudinal variations and relevance for affective disorders

Background

The temporal coordination of biological processes into daily cycles is a common feature of most living organisms. In humans, disruption of circadian rhythms is commonly observed in psychiatric diseases, including schizophrenia, bipolar disorder, depression and autism. Light therapy is the most effective treatment for seasonal affective disorder and circadian-related treatments sustain antidepressant response in bipolar disorder patients. Day/night cycles represent a major circadian synchronizing signal and vary widely with latitude.

Results

We apply a geographically explicit model to show that out-of-Africa migration, which led humans to occupy a wide latitudinal area, affected the evolutionary history of circadian regulatory genes. The SNPs we identify using this model display consistent signals of natural selection using tests based on population genetic differentiation and haplotype homozygosity. Signals of natural selection driven by annual photoperiod variation are detected for schizophrenia, bipolar disorder, and restless leg syndrome risk variants, in line with the circadian component of these conditions.

Conclusions

Our results suggest that human populations adapted to life at different latitudes by tuning their circadian clock systems. This process also involves risk variants for neuropsychiatric conditions, suggesting possible genetic modulators for chronotherapies and candidates for interaction analysis with photoperiod-related environmental variables, such as season of birth, country of residence, shift-work or lifestyle habits.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0499-7) contains supplementary material, which is available to authorized users.     

Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster

Determining the genetic basis of environmental adaptation is a central problem of evolutionary biology. This issue has been fruitfully addressed by examining genetic differentiation between populations that are recently separated and/or experience high rates of gene flow. A good example of this approach is the decades-long investigation of selection acting along latitudinal clines in Drosophila melanogaster. Here we use next-generation genome sequencing to reexamine the well-studied Australian D. melanogaster cline. We find evidence for extensive differentiation between temperate and tropical populations, with regulatory regions and unannotated regions showing particularly high levels of differentiation. Although the physical genomic scale of geographic differentiation is small--on the order of gene sized--we observed several larger highly differentiated regions. The region spanned by the cosmopolitan inversion polymorphism In(3R)P shows higher levels of differentiation, consistent with the major difference in allele frequencies of Standard and In(3R)P karyotypes in temperate vs. tropical Australian populations. Our analysis reveals evidence for spatially varying selection on a number of key biological processes, suggesting fundamental biological differences between flies from these two geographic regions.