Skip to main content
  • Research article
  • Open access
  • Published:

Schizophrenia and vitamin D related genes could have been subject to latitude-driven adaptation

Abstract

Background

Many natural phenomena are directly or indirectly related to latitude. Living at different latitudes, indeed, has its consequences with being exposed to different climates, diets, light/dark cycles, etc. In humans, one of the best known examples of genetic traits following a latitudinal gradient is skin pigmentation. Nevertheless, also several diseases show latitudinal clinals such as hypertension, cancer, dismetabolic conditions, schizophrenia, Parkinson's disease and many more.

Results

We investigated, for the first time on a wide genomic scale, the latitude-driven adaptation phenomena. In particular, we selected a set of genes showing signs of latitude-dependent population differentiation. The biological characterization of these genes showed enrichment for neural-related processes. In light of this, we investigated whether genes associated to neuropsychiatric diseases were enriched by Latitude-Related Genes (LRGs). We found a strong enrichment of LRGs in the set of genes associated to schizophrenia. In an attempt to try to explain this possible link between latitude and schizophrenia, we investigated their associations with vitamin D. We found in a set of vitamin D related genes a significant enrichment of both LRGs and of genes involved in schizophrenia.

Conclusions

Our results suggest a latitude-driven adaptation for both schizophrenia and vitamin D related genes. In addition we confirm, at a molecular level, the link between schizophrenia and vitamin D. Finally, we discuss a model in which schizophrenia is, at least partly, a maladaptive by-product of latitude dependent adaptive changes in vitamin D metabolism.

Background

During recent decades the study of human evolution has been of increasing interest, due also to the large amount of data now available, concerning populations geographically and widely distributed. Meanwhile, new applications of evolutionary biology to medical problems are being discovered at an accelerating rate [1].

The genetic diversity is one of the most important instruments available to understand our evolutionary history. Differences among individuals belonging to the same population are generally smaller than those of individuals belonging to different populations [2]. This is due to both demographic history and selection that shaped the genome to adapt the genome of different populations to the experienced environment. All loci in a population share the same demographic history, thus they are expected to show similar patterns of variation. But it is a matter of fact that substantial differences in the level of among-population differences across loci exist [3]. It is therefore reasonable to assume that loci showing more among-population differentiation than the rest might identify regions of the genome that have been subjected to selection [4, 5].

In 1966 Cavalli-Sforza suggested using measures of population divergence to detect natural selection [6] and Lewontin and Krakauer proposed using for this purpose, Wright's F-statistics [7]. In particular, fixation index (FST) is one of the most widely used descriptive measures of population and evolutionary genetics [8]. FST is directly related to the variance in allele frequency among populations. A small FST value means that the allele frequencies among populations are similar, whereas a large FST value implies that the allele frequencies are different. If natural selection favours one allele among others at a particular locus and only for some particular populations, the corresponding FST for that locus will be larger than in any other locus where genetic drift alone is at work.

Among the environmental factors that strongly influenced our evolutionary history, geographical latitude deserves particular attention. Latitude, indeed, severely affects many natural phenomena such as climate, flora and fauna, light-dark cycle, and all of them, in turn, have an impact on many aspects of our life. For sake of brevity hereafter we refer to all these phenomena simply as "latitude".

Genetic traits following a latitudinal gradient have been observed for several polymorphisms in humans as well as in natural populations of model organisms like Drosophila and Arabidopsis thaliana [9–12]. The best known example of this kind of spatial variation in Homo Sapiens is skin pigmentation. The clinal gradation of skin colouration is correlated with UV radiation levels and represents a compromising solution to the conflicting physiological requirements of photoprotection and vitamin D UV-dependent synthesis [13]. The latter is very important since vitamin D is involved in many health outcomes (e.g. cardiovascular diseases, rickets, pelvic deformities, infections, etc.).

A possible influence of latitude on the circadian phenotype has also been suggested [14, 15]. Circadian rhythm is a ubiquitous feature of living systems. Daylight hours vary with latitude and seasons therefore adaptability of circadian clocks is of fundamental importance for the adaptation of organisms to the alternating light/dark cycles.

Another example of spatial variation is human body size and shape. These phenotypes show a correlation with climate (that in turn has a strong relationship with latitude) suggesting also for humans the adaptation to the classical ecological rules that individuals living in colder regions are bulkier and have shorter limb lengths [16, 17].

Nevertheless, it is also possible that an allelic variant increasing the fitness of individuals at particular latitudes will no longer be advantageous or even increase the risk for some pathologies at different latitudes. In particular, several diseases show latitudinal clinals. A well known example regards sodium homeostasis. It have been postulated that in hot climate regions, genetic variants inducing enhanced sodium retention were positively selected. This adaptive process would allow a proper vascular tone and salt storage in conditions of excessive sweating. These same variants, when carried by individuals migrated in colder climates (i.e. African American), would increase the risk for sodium retention-related hypertension [18]. Supporting this hypothesis, several studies reported a strong correlation between latitude and the frequencies of hypertension susceptibility variants [18, 19].

In addition, there is growing evidence supporting the idea that these diseases are frequently due to a negative by-product of adaptive changes during human evolution [20]. This is the case when natural selection favours a vital phenotype at the price of predisposing to some other pathology that, for example, do not directly affect the reproduction or are characterized by a late onset-age. Indeed, contrasting forces often affect the outcome of natural selection. For instance, depigmentation is crucial for vitamin D synthesis at higher latitudes but it also exposes a higher risk for skin cancer. As most individuals do not develop cancer until they are past their reproductive age, from an evolutionary point of view, skin cancer represents a less powerful selective force than vitamin D availability in serum [21].

Many other common diseases like different types of cancer, dismetabolic conditions, schizophrenia, Parkinson's disease, etc. have an incidence following a latitudinal gradient [22–25]. However, the relative importance of variation of environmental exposures or genetic predisposition is not yet fully defined. In some cases, a genetic adaptation has been suggested [20, 22], even if the direct target of this process is often unclear.

In this paper, we investigated the latitude-driven adaptation phenomena, for the first time, on a wide genomic scale. In particular, we selected a set of SNPs and genes showing signs of latitude-dependent population differentiation, and by a biological characterization of the genes, we found enrichment for neural-related processes. In light of this result, we investigated whether genes associated to pathological phenotypes, namely psychiatric and neurological diseases, were enriched for Latitude-Related Genes (LRGs). Remarkably, we found a strong enrichment of LRGs in the set of genes associated with schizophrenia. In an attempt to try to explain this possible link between latitude and schizophrenia, we investigated their association with vitamin D, which had been previously associated, separately, to both of them. Our findings suggest a molecular link among latitude, schizophrenia and vitamin D.

Results

A set of SNPs showing high levels of latitude-dependent population differentiation was selected by using a two-step approach. Our starting point consists of genotype data concerning about 660,000 SNP loci of 938 unrelated individuals from 51 populations of the Human Genome Diversity Panel [26].

The first step was the estimation of the population differentiation level of each SNP. After the exclusion for minor allele frequency and for SNPs falling in intergenic regions, we obtained a set of 224,501 SNPs. For all of them we calculated the fixation index (FST) according to the Weir and Cockerham estimator [27], and to select SNPs with high levels of population differentiation, we extracted those at the top of the empirical distribution of FST values. Because of the differences existing between the distribution of FST values for autosomic and X-linked SNPs [28], these two sets were handled separately. In particular, we selected 22,132 autosomic and 459 X-linked SNPs falling in the top 10% of their own distributions (namely with FST value greater than 0.153 and 0.262, respectively).

In the second step, we computed for each SNP the absolute value of correlation between the frequency of the ancestral allele in the 51 populations and the absolute value of geographical latitude of the population location. We again handled the distributions of correlation values of autosomic and X-linked SNPs separately because of their differences. Moreover, since population sizes were very different, we computed the correlation by taking into account the number of individuals in each population. We selected the autosomic and X-linked SNPs with absolute value of correlation greater than 0.567 and 0.575 respectively, corresponding to the highest 10% of their respective distributions. We finally obtained two sets of 2193 autosomic and 46 X-linked SNPs corresponding to 1307 and 29 unique genes. Hereafter, we denote by LRGs (Latitude-Related Genes) this set of genes whose SNPs showed both high FST (mean = 0.230) and high latitude correlation (mean = 0.616) values (Additional file 1).

High FST values can be produced both by selection and by demography. The effect of population histories is potentially the major confounding factor in the interpretation of genetic differences among populations. To roughly estimate the extent of phenomena that could be ascribed to natural selection in our LRG set, we used an approach similar to that proposed by Barreiro et al. [29]. We found in LRG an enrichment of non-synonymous SNP of about 30% (p = 0.016, Fisher's exact test; Additional file 2). This finding seems to indicate that, at least in part, LRGs set is enriched for genes under selective pressure.

As a further test for the procedure, we chose a well-known latitude dependent phenomenon such as skin pigmentation and we checked whether LRGs were statistically overrepresented in the set of genes associated to this phenotype. As expected a significant enrichment was found (7 out of 24 genes in common, Fisher's exact test p-value 0.0003).

To characterize the set of LRGs we used two methods, tissue localization and functional characterization. Firstly, we explored the tissue localization of proteins encoded by these genes. To achieve this, we used the Database for Annotation, Visualization and Integrated Discovery (DAVID) focusing on the "Uniprot Tissue" list, a curated list of localization based on literature mining. We found that genes expressed in brain and brain-related tissues were significantly enriched in LRGs, accounting for more than a half of them (Table 1).

Table 1 Enrichment in tissue for LRGs computed using the DAVID's Uniprot Tissue category

LRGs were also functionally characterized by looking for overrepresentation of Gene Ontology (GO) annotation terms [30]. Because of the high redundancy of GO, we used the Model-based Gene-Set Analysis (MGSA) method included in Ontologizer 2.0 [31, 32]. This promising and novel approach analyses all categories together by embedding them in a Bayesian network. Differently from other methods, it provides for each term a marginal posterior probability that reflects a measure of certainty in its involvement in the process. Following the authors' recommendation, we repeated the analysis 20 times in order to see whether the reported marginal probabilities of the top terms fluctuated. We found that two terms out of the three showing a posterior probability above 0.5 consistently among the 20 runs were related to neural processes (Table 2).

Table 2 Enrichment for GO terms by LRGs

All these results suggested a further investigation about a possible relationship between LRGs and genes involved in neuropsychiatric diseases related with latitude. Indeed, several neurological diseases were previously described to have a latitude-shaped incidence and/or prevalence [23, 25, 33]. To perform this task we compared the list of LRGs with publicly available collections of genes involved in schizophrenia, multiple sclerosis, Parkinson's and Alzheimer's disease [34–37]. While there is a weak or non significant overlap with genes related to multiple sclerosis, Parkinson's and Alzheimer's disease, we found a significant enrichment of 85 LRGs in genes related to schizophrenia (Fisher's exact test, Bonferroni adjusted p-value 1.6 × 10-5; Table 3, Figure 1 and Additional file 3). We found a similar enrichment by comparing the list of LRGs with other lists of schizophrenia related genes (Additional file 2).

Table 3 Enrichment of neuropsychiatric diseases lists by LRGs
Figure 1
figure 1

Overlaps among latitude, vitamin D and schizophrenia related genes. In each circle is reported the number of genes present in the list. For the two-way intersections is reported the size and the significance of the overlap (Fisher's exact test).

We then investigated for possible latitude-dependent biological mechanisms linking latitude to neural development. An important factor hypothesized to be both latitude dependent and neural development-related is vitamin D. We checked for an enrichment of LRGs in Vitamin D Related (VDR) genes by manually creating a list of 943 genes that broadly comprised the most important processes in which it is involved, since a comprehensive list is not yet present in the literature (Additional file 4). To achieve this, we merged: a list of 6 genes implied in the metabolism of vitamin D obtained from Reactome, a list of 26 genes in the pathway of the control of the expression by vitamin D receptor from Biocarta and a list of 911 genes from a large-scale identification of 1,25(OH)2D3 target genes by Wang and colleagues [38].

We computed the overlap between LRGs and VDR genes and found a significant overlap of 97 genes (p-value 3.5 × 10-8, Fisher's exact test; Figure 1 and Additional file 3). This result suggests that VDR genes show signature of latitude-dependent population differentiation. Also the overlap of 70 genes between VDR and schizophrenia related genes was significant (p-value 1.4 × 10-6, Fisher's exact test; Figure 1 and Additional file 3), confirming the role of vitamin D in schizophrenia pathogenesis.

Finally, we found 9 genes (SMARCA2, MITF, DLGAP1, MAGI1, IL4R, NTRK3, RUNX1, PPP3CA, and INPP4B) in common among those ones related to latitude, vitamin D and schizophrenia (Figure 1). We checked whether or not in these 9 genes belonged SNPs, selected by our procedure, that were previously studied in relationship to either schizophrenia or vitamin D related phenotypes. One SNP resulted from the analysis, rs3793490 (FST = 0.202, correlation = 0.626), an intronic SNP of the SMARCA2 gene. In Figure 2 is reported its alleles geographic distribution (A) and the values of cross-Population Extended Haplotype Homozygosity (XP-EHH) (B) of the genomic region. This test detects alleles that have risen to high frequency rapidly, enough that long-range association with nearby polymorphisms (the long-range haplotype) have not been eroded by recombination. The analysis showed a strong sign of recent selective pressure.

Figure 2
figure 2

Alleles geographic distribution (A) and cross-population extended haplotype homozygosity (B) for the SNP rs3793490. Images were readapted from the "HGDP Selection Browser". The red line marks the SNP's position inside the SMARCA2 gene region.

Discussion

Many natural phenomena are directly or indirectly related to latitude. Living at different latitudes has consequences in being generally exposed to different climates, diets, light/dark cycles, etc. Therefore, it is reasonable to presume that exposure of individuals to different latitudes could have shaped genetic background as a result of the adaptation process. Indeed, relationships between allelic frequencies of specific genes and latitude have been identified in plants [39–41] and animals [42, 43].

Previous studies in humans, at specific loci, have found evidence of correlation between allelic frequencies and latitude of sampled populations. For example allelic frequencies of a genetic polymorphism in the prion gene PRNP showed a clear correlation with the latitude within Europe [44] and one in the ACP1 worldwide [45].

To explore possible genetic adaptations to latitude, within this work we defined a set of latitude-related genes (LRGs) following a two step approach. Firstly, we identified SNPs with a high level of population differentiation (FST) with the aim to enrich for variants under selective pressure. From these we then extracted those SNPs showing high values of correlation of allelic frequencies with the geographical latitude (Additional file 1). To the best of our knowledge this is the first search at a wide genomic level for loci showing latitude-dependent populations differentiation.

Both functional characterization and expression localization of LRGs resulted in a strong enrichment of neural-related processes (Tables 1 and 2). The relationship between neural development and latitude is partially known. In particular, there exists evidence of a latitude correlation of some physiological neural-related phenotypes. It was reported for humans that cranial capacity is different between populations, probably, as results of adaptation for brain cooling and that the craniofacial diversity results from the tissues of neural crest origin [17]. It is worth stressing that although population differences are observed in these phenotypes, no relationship exists with mental functioning. Also, for several pathological neural phenotypes there has been described previously a relationship with latitude. For multiple sclerosis and schizophrenia there was a latitudinal variation in incidence and prevalence [25, 33]. An association between mortality related to Parkinsonism and birthplace geographical latitude was also found [46].

When we compared the list of LRGs with those containing genes associated to neurological and psychiatric diseases, we found a vanishing enrichment of LRGs in genes related to multiple sclerosis, Alzheimer's and Parkinson's disease (Table 3). This is not surprising. Concerning Alzheimer's and Parkinson's diseases, it should be noted that the relationships with latitude reported in the literature are not largely confirmed as in the case of the other two diseases. Multiple sclerosis has a strong immunity component in the etiology. Since genes involved in immunity related processes usually exhibit low levels of FST [4, 28], they are likely to be excluded by our procedure, in this way weakening the enrichment.

In contrast, we found a strong significant enrichment of LRGs in genes associated with schizophrenia (Table 3, Figure 1 and Additional file 3). The correlation of schizophrenia prevalence with latitude is both large and robust [47]. In previous surveys it was noted that there is a strong tendency for schizophrenia prevalence to increase with increasing latitude [25, 48]. The present result also agrees with that obtained in one of our previous works, based on different data and different statistical approaches. In this previous study, we explored for signs of natural selection in genes associated to complex diseases, and the strongest level of population differentiation was observed in genes associated to psychiatric diseases [28]. In addition, testing for recent selective sweeps in human populations, Crespi and colleagues found significant evidence for adaptive evolution of several genes underlying schizophrenia [20]. Our results may suggest that a latitude-related adaptation occurred for some schizophrenia associated genes, but to which extent this phenomenon is related to the latitude-shaped prevalence of schizophrenia is of course still an open question.

Searching for a mechanism that could connect latitude and schizophrenia at a molecular level, we reasoned that a well-known phenomenon separately linked to both of them is vitamin D. Indeed, vitamin D is essential for normal growth, calcium absorption and skeletal development. The cutaneous synthesis of vitamin D is a function of skin pigmentation and of the solar zenith angle which, in turn, depends on latitude, season, and time of day [49]. With the same dietary intake, the most important determinant for vitamin D levels is considered to be where individuals live, because of the dependence on geographical location of the availability of UV radiation for vitamin D synthesis. Ancestral populations who migrated out of Africa, moving from south to north were exposed to less incident sunlight. It is commonly accepted that, in these conditions, depigmentation was favoured [50]. This adaptation process was compulsory, since pelvic deformities due to vitamin D deficiency could prevent normal childbirth, but it is reasonable to suggest that also other molecular mechanisms apart from depigmentation could have been subject to selective pressure in order to increase the levels of vitamin D in the body.

To test this link at a genomic level, we created a hand-curated list of genes that, at different levels, are related to vitamin D (see Results; Additional file 4). We found among vitamin D related genes a significant enrichment for genes latitude-related (Figure 1 and Additional file 3). This result may suggest that vitamin D related genes could have been subject, at least in part, to a latitude-dependent adaptation, supporting the hypothesis that selection also acted on other mechanisms different from depigmentation.

On the other hand, the relationship between vitamin D and schizophrenia is well established. Vitamin D receptors were found in most tissues other than those classically involved in the vitamin D action (bone, gut, kidney, etc.). In particular, receptors for vitamin D are widely distributed in the nervous system and vitamin D has been recently implicated in brain function [51]. There is growing evidence that low vitamin D levels adversely impact on brain development [52]. In mice, it has been suggested that changes in brain development induced by prenatal vitamin D deficiency lead to specific functional alterations in hippocampal synaptic plasticity [53]. Also in humans, maternal vitamin D insufficiency has been associated with childhood rickets and longer term problems including schizophrenia [54]. In addition, epidemiological data suggest that vitamin D deficiency may be associated with increased risks of mental health disorders such as schizophrenia [55]. Despite the consistence of this relationship, to our knowledge studies connecting schizophrenia and vitamin D at a molecular level are not yet available.

We found a significant overlap of 70 genes between schizophrenia and vitamin D related genes (Figure 1 and Additional file 3). Our analysis provides the first hint, at a genomic level, of the existence of a relationship between them. The largest part of vitamin D related genes in our list is made by genes differentially expressed in epithelial cells after treatment with the biological active form of vitamin D [56]. It is therefore possible that the same genes could be regulated by vitamin D levels also in neurons during brain development that, in turn, have been associated to schizophrenia [47].

According to our hypothesis that vitamin D could be the link between latitude and schizophrenia, we focused on the 9 LRGs present in both the lists of vitamin D and schizophrenia related genes (Figure 1). Among these genes, that were related to schizophrenic phenotypes, several of them were also previously described in association with bone development. However, no evidence of selective adaptation was present in the literature until now.

For example, one of these genes is the neurotrophic tyrosine kinase receptor type 3 (NTRK3) which encodes a member of the NTRK family. These neurotrophins (NTs) receptors are best known for their role in the differentiation and survival of various types of neurons [57]. Gene expression of NTRK3 has been reported to be reduced in patients with schizophrenia [58, 59]. Furthermore, it has been suggested that the NTRK3 gene influences hippocampal function and may modify the risk of schizophrenia [60]. Nevertheless, NTs and their receptors are also produced by a growing list of non-neuronal cells [61] including osteoblastic cell lines [62]. NTs receptors were observed in the bone forming area during fracture healing. NTRK3 was observed in osteoblast like cells and hypertrophic chondrocytes [63].

Another example is PPP3CA (protein phosphatase 3 catalytic subunit alpha isoform). This gene, also known as calcineurin A alpha, acts as a calcium-sensor and regulator of calcium homeostasis. For this reason, it shapes calcium and cyclic AMP dependent processes like synaptic activity, receptor desensitization, cell survival and neuroplasticity. Expression of PPP3CA, previously described by microarrays analysis from cortical [64] and rat hippocampus tissues [65], was found down regulated in schizophrenic anterior temporal lobe [66]. On the other hand, PPP3CA is also expressed in osteoclasts, playing a role in the regulation of bone resorption and its deletion results in osteoporosis [67, 68].

At least in one case, we also found independent evidence of an adaptive process, reasonably suggesting a molecular mechanism. This is the case of the SMARCA2 (SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily a, member 2) gene. The protein encoded by SMARCA2 (also known as BRM) is part of the large ATP-dependent chromatin remodelling complex SNF/SWI, which is required for transcriptional activation of genes normally repressed by chromatin. Mammalian SWI/SNF complex actually consists of a small series of compositionally distinct assemblies distinguished by the presence of alternative subunits. The complexes contain either one of two closely related alternative ATPases: Brahma (BRM i.e. SMARCA2) or Brahma-related gene 1 (BRG1). The combinatorial assembly of these complexes could account for the specificity of their functions in different tissues and development phases. SWI/SNF components and DNA replication-related factors form, in turn, a human multiprotein complex (WINAC) that directly interacts with vitamin D receptor [69]. WINAC and vitamin D receptor are targeted to vitamin D responsive promoters in the absence of ligand to both positively and negatively regulated genes. WINAC may rearrange the nucleosome array around the positive and negative vitamin D responsive elements (VDREs), thereby facilitating the coregulatory complexes access for further transcription control. Subsequent binding of coregulators requires ligand binding. Several studies have revealed that one family of the SWI/SNF complexes based on the BRG1 and BRM ATPases has particular critical dosage-dependent roles in the development of the nervous system [70, 71].

One of the five SNPs with the highest levels of both population differentiation and correlation with latitude that we identified at this locus was rs3793490. In a previous study this SNP was associated with SMARCA2 expression levels in the human brain [72]. In particular, the T variant, whose frequency increases in a latitude dependent way, is associated with low SMARCA2 expression levels. Much evidence suggests that low levels of SMARCA2 may play a role in the pathophysiology of schizophrenia. SMARCA2 knockout mice showed impaired social interaction and prepulse inhibition. In the mouse brain, psychotogenic drugs lowered SMARCA2 expression while antipsychotic drugs increased it [72]. On the other hand, in MC3T3-E1 cell line, deficiency of SMARCA2 results in an accelerated rate of mineralization with higher levels of expression of osteogenic markers [73]. In addition, the most prominent phenotype of SMARCA2 null mice is a larger (about 14% more than normal) body size with a disproportionately increased bone and muscle mass [74]. Putting together all this data seems to suggest that SMARCA2 deficiency is associated with different effects on bones and neurons. It is possible to speculate that a selective pressure could have favoured a haplotype containing the T allele in an environment with low vitamin D availability, for its positive effect on the bone phenotype. In turn, the low SMARCA2 expression levels linked to this variant could be responsible for the increased risk of schizophrenia. The XP-EHH analysis that we performed (Figure 2) confirms the presence of a strong recent positive selection in the SMARCA2 locus. Interestingly, the homozygosity geographical pattern detected by XP-EHH recognises the effect that we found with correlation analysis. These results agree with a previously described hypothesis of schizophrenia being, at least partly, a maladaptive by-product of adaptive changes during human evolution [20].

In Figure 1, besides the described overlaps among the lists used, large non-intersecting areas exist. This can be partially due to technical limitations of our work. Anyway, these areas of non-intersection were, at least in part, expected. As stated, latitude acts directly or indirectly, on a wide variety of phenomena, therefore it is expected that its effects are not limited just to vitamin D and neural development. For similar reasons, we expected that also the link between latitude and schizophrenia cannot be completely explained by vitamin D. For example, in the intersection between LRGs and genes related with schizophrenia (Additional file 3) there is a gene involved in the circadian rhythms, TIMELESS.

Circadian rhythm consists of light and dark phases which coincide with the phases of the solar day and that is, obviously, correlated with the different photoperiods existing at different latitudes. TIMELESS is required for normal progression of S-phase and is involved in the circadian rhythm autoregulatory loop. Associations between circadian gene polymorphisms, including TIMELESS, and some mental disorders have been found, including schizophrenia [75]. In addition, expression of TIMELESS was investigated in the pitcher-plant mosquito, Wyeomyia smithii, and was found to vary with latitude of origin. This suggests that other mechanisms should be taken into account [76].

It is worth stressing that this work represents just an initial exploration of this complex problem and can suffer from some limitations. The first aspect to take into account is the arbitrary choice of the thresholds used as inclusion criteria in the list of LRGs. However, it should be underlined that our aim was just to obtain a list enriched for genes showing both high levels of population differentiation and correlation with latitude. Another important aspect is that the largest part of our analysis is based on hand-curated lists. Nevertheless, it should be noted that these lists are widely used [77, 78]. The only exception is the list of vitamin D related genes that we had to build ex novo since no others were present in literature. In addition, a cause-effect relationship can never be conclusively established based only on a reciprocal enrichment between sets of genes. Moreover, alternative explanations of the observed enrichments can be imagined. For example, we cannot exclude the presence of an unknown common factor, different from vitamin D, linking together schizophrenia and latitude. A further limit of our approach concerns the gene level analysis. In fact, focusing on genes rather than on SNPs, does not allow to conclusively assess whether variants related to latitude coincide with those involved in the studied phenotypes. In other words, even if a gene is associated to both latitude adaptation and schizophrenia, we cannot exclude that this is due to functionally independent variants. For this reason our conclusions should be corroborated by further and more detailed studies. Finally, we are aware that population differentiation is influenced by demographic history and thus it cannot be straightforwardly interpreted as a sign of natural selection. Some approaches have been proposed to evaluate the impact of natural selection on population differentiation. Under the assumption of neutrality, any set of SNPs, even if classified according to their physical location and functional impact, should show the same degree of population differentiation. According to Barreiro and colleagues, any deviation from this expectation should be attributable to selection [29]. In particular, the authors observed that variants leading to amino-acid changes (non-synonymous mutations) are overrepresented among SNPs showing high level of FST. They interpreted this excess as the result from the action of natural selection. We used a similar approach to analyse our data. We found an excess of non-synonymous polymorphisms in LRGs. This result seems to suggest the appreciable presence of genes under selective pressure in LRGs.

During the preparation of the manuscript, an international team of researchers presented the first detailed analysis of the draft sequence of the Neanderthals' genome [79]. In particular, they showed the presence of interbreeding between Neanderthals and Homo Sapiens occurred after the Out-of-Africa migration. Affecting mainly the non-African populations, this genetic flow could mimic a latitudinal effect on the Africa-Europe axis and therefore potentially influences our conclusions. First of all, LRGs are selected according to a worldwide latitude correlation, which is only partially due to the Africa-Europe axis. This is the case, for example, of the rs3793490 SNP (Figure 2) where it is also clear a latitudinal effect out of this axis. In addition, Europeans and Asians share only 1% to 4% of their nuclear DNA with Neanderthals suggesting a limited impact, if any, on our conclusions. Finally, we observed that approximately 10% of LRGs are present in the list of genes that Green et al. showed to have evolved recently in our lineage after we split from Neanderthals, and thus, at least in these cases, we can confidently exclude the influence of interbreeding.

Conclusions

We found that frequencies of genetic variants involved in brain development, schizophrenia and vitamin D-related processes significantly vary with latitude. We also found a molecular relationship of schizophrenia with both latitude and vitamin D. At least in one case, our results suggest that schizophrenia could be a maladaptive by-product of latitude dependent adaptive changes in vitamin D metabolism. The great importance of latitude in human evolutionary history suggests that also other pathologies should be explored.

Methods

Data

The whole analysis is based on genotypes data of 660,918 single-nucleotide polymorphisms (SNPs) in samples from the Human Genome Diversity Panel (HGDP-CEPH), which represents 938 unrelated individuals from 51 populations from sub-Saharan Africa, North Africa, Europe, the Middle East, South/Central Asia, East Asia, Oceania, and the Americas. We removed the two small heterogeneous Southern Bantu populations as in [80]. Genotypic data was retrieved from http://hagsc.org/hgdp/[81] while geographical information was obtained from http://www.cephb.fr/en/hgdp/.

For all the SNPs, we computed the allelic frequencies within each population by using the R package "genetics" version 1.3.4. Additional SNP information about physical positions and SNP-gene mapping were obtained from dbSNP build 129 http://www.ncbi.nlm.nih.gov/projects/SNP. Data from the HGDP-CEPH project and dbSNP were merged in a local MySQL database.

After excluding SNPs with minor allele frequencies less than 5% in all of the populations, we obtained a set of 655,810 SNPs. Since we were interested in performing this study at a gene level, we retained only intragenic SNPs obtaining a final set of 224,501 SNPs. In particular, we considered all SNPs within 2 kb of a gene, according to dbSNP classification. We also excluded the Y-linked SNPs because of scarce numbers.

Genes underlying skin pigmentation were obtained from a study by Myles and colleagues [82]. Genes related to neuropsychiatric diseases were obtained from a set of publicly available comprehensive, uniform and regularly updated database of genes considered involved in schizophrenia, multiple sclerosis, Parkinson's and Alzheimer's disease [34–37]. These databases collect genes related to these phenotypes resulting from different approaches (association, genome-wide, candidate, etc.). We used the January 30th 2010 update of these databases, containing 892, 197, 511 and 622 genes, respectively. The list of vitamin D related genes was manually created by merging three different lists of genes related in different ways to vitamin D (Additional file 4). The first list is the "Vitamin D (calciferol) metabolism" pathway by Reactome (REACT_13523.2). The second one was extracted from the Biocarta's pathway "Control of the expression by vitamin D receptor" (h_vdrPathway) as present on the Cancer Genome Anatomy Project http://cgap.nci.nih.gov/. Finally, candidate transcriptional target genes of vitamin D were obtained from a set of genes differentially expressed in SCC25 cells treated with 1,25(OH)2D3 by a genome-wide microarray analysis [56].

A common problem in using lists of genes is that different types of gene identifiers are used. The effect is that usually only part of the genes is recognized. Therefore, each tool that we used lost a variable number of genes during the analyses. Nevertheless, to obtain the best performances from each tool we decided to provide them with the whole lists. In particular, DAVID recognized 1207 identifiers and MGSA recognized 1101 ids. In addition, overlaps between lists of genes were computed by using MatchMiner web tool, which also takes into account gene aliases and/or different types of gene identifiers [83]. This tool recognized 885 schizophrenia, 178 multiple sclerosis, 490 Parkinson's disease, 618 Alzheimer's disease and 1254 latitude related genes.

Statistical analysis

For all SNPs we calculated the fixation index (FST) according to the Weir and Cockerham estimator [27] using a previously developed Perl script [28]. We then computed the absolute value of the correlation between the frequency in each population of the ancestral allele and the absolute value of the latitude of the population weighted by its sample size. This was to avoid the overweighting of allelic frequencies inside small populations and was implemented simply by using the weighted mean, standard deviation and covariance.

Statistical significance of overlaps was estimated by using Fisher's exact test, considering 21463 genes (number of distinct gene symbols present in the "refGene" track of UCSC Genome Browser) as background population. All statistical analyses were carried out with R ver. 2.10.1 [84]. The whole study was conducted considering a p-value of 0.001 as statistical significance threshold.

To check for a potential bias toward larger genes (since they can contain more genotyped SNPs), we applied the non -parametric Mann-Whitney test to compare the difference among the lengths of LRGs with respect to the remaining genes in the "refGene" track of UCSC Genome Browser. There is no statistical evidence for a difference (p = 0.07), with median length of LRGs slightly smaller than median of the others.

Biological characterization

The analysis of gene expression localization was performed by using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 and the "UP_TISSUE" category [85, 86]. The "Uniprot tissue" (UP_TISSUE) list is based on literature mining and reports for each gene in which tissues it has been found to be expressed, by using a curated vocabulary.

Overrepresentation of GO terms was assessed by using Ontologizer 2.0 and the Model-based Gene Set Analysis (MGSA) method. The idea behind this approach is to estimate the marginal posterior probability of a term being enriched by using a Bayesian Network. The greater probability for a GO term being near to 1, the higher is the certainty of its involvement in the process. All parameters were left to "auto", allowing the system to automatically estimate the optimal a priori probabilities and the false positive and false negative rates. Probabilities are estimated by MGSA using 107 Markov-Chain Monte Carlo (MCMC) steps. Since MCMC is not guaranteed to converge in any a priori defined number of steps, we repeated the analysis 20 times and retained only terms having a marginal posterior probability greater than 0.5 in every run, as recommended by authors.

References

  1. Nesse RM, Bergstrom CT, Ellison PT, Flier JS, Gluckman P, Govindaraju DR, Niethammer D, Omenn GS, Perlman RL, Schwartz MD, Thomas MG, Stearns SC, Valle D: Evolution in health and medicine Sackler colloquium: Making evolutionary biology a basic science for medicine. Proc Natl Acad Sci USA. 2010, 107 (Suppl 1): 1800-1807. 10.1073/pnas.0906224106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Rosenberg NA, Pritchard JK, Weber JL, Cann HM, Kidd KK, Zhivotovsky LA, Feldman MW: Genetic structure of human populations. Science. 2002, 298: 2381-2385. 10.1126/science.1078311.

    Article  CAS  PubMed  Google Scholar 

  3. Weir BS, Cardon LR, Anderson AD, Nielsen DM, Hill WG: Measures of human population structure show heterogeneity among genomic regions. Genome Res. 2005, 15: 1468-1476. 10.1101/gr.4398405.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Akey JM, Zhang G, Zhang K, Jin L, Shriver MD: Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 2002, 12: 1805-1814. 10.1101/gr.631202.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Xue Y, Zhang X, Huang N, Daly A, Gillson CJ, Macarthur DG, Yngvadottir B, Nica AC, Woodwark C, Chen Y, Conrad DF, Ayub Q, Mehdi SQ, Li P, Tyler-Smith C: Population differentiation as an indicator of recent positive selection in humans: an empirical evaluation. Genetics. 2009, 183: 1065-1077. 10.1534/genetics.109.107722.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Cavalli-Sforza LL: Population structure and human evolution. Proc R Soc Lond B Biol Sci. 1966, 164: 362-379. 10.1098/rspb.1966.0038.

    Article  CAS  PubMed  Google Scholar 

  7. Lewontin RC, Krakauer J: Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics. 1973, 74: 175-195.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. Holsinger KE, Weir BS: Genetics in geographically structured populations: defining, estimating and interpreting F(ST). Nat Rev Genet. 2009, 10: 639-650. 10.1038/nrg2611.

    Article  CAS  PubMed  Google Scholar 

  9. Sezgin E, Duvernell DD, Matzkin LM, Duan Y, Zhu C, Verrelli BC, Eanes WF: Single-locus latitudinal clines and their relationship to temperate adaptation in metabolic genes and derived alleles in Drosophila melanogaster. Genetics. 2004, 168: 923-931. 10.1534/genetics.104.027649.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Verrelli BC, Eanes WF: Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster. Genetics. 2001, 157: 1649-1663.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. Balasubramanian S, Sureshkumar S, Agrawal M, Michael TP, Wessinger C, Maloof JN, Clark R, Warthmann N, Chory J, Weigel D: The PHYTOCHROME C photoreceptor gene mediates natural variation in flowering and growth responses of Arabidopsis thaliana. Nat Genet. 2006, 38: 711-715. 10.1038/ng1818.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Caicedo AL, Stinchcombe JR, Olsen KM, Schmitt J, Purugganan MD: Epistatic interaction between Arabidopsis FRI and FLC flowering time genes generates a latitudinal cline in a life history trait. Proc Natl Acad Sci USA. 2004, 101: 15670-15675. 10.1073/pnas.0406232101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Jablonski NG, Chaplin G: The evolution of human skin coloration. J Hum Evol. 2000, 39: 57-106. 10.1006/jhev.2000.0403.

    Article  CAS  PubMed  Google Scholar 

  14. Pereira DS, Tufik S, Louzada FM, Benedito-Silva AA, Lopez AR, Lemos NA, Korczak AL, D'Almeida V, Pedrazzoli M: Association of the length polymorphism in the human Per3 gene with the delayed sleep-phase syndrome: does latitude have an influence upon it?. Sleep. 2005, 28: 29-32.

    PubMed  Google Scholar 

  15. Cruciani F, Trombetta B, Labuda D, Modiano D, Torroni A, Costa R, Scozzari R: Genetic diversity patterns at the human clock gene period 2 are suggestive of population-specific positive selection. Eur J Hum Genet. 2008, 16: 1526-1534. 10.1038/ejhg.2008.105.

    Article  CAS  PubMed  Google Scholar 

  16. Roberts DF: Body weight, race and climate. Am J Phys Anthropol. 1953, 11: 533-558. 10.1002/ajpa.1330110404.

    Article  CAS  PubMed  Google Scholar 

  17. Irmak MK, Korkmaz A, Erogul O: Selective brain cooling seems to be a mechanism leading to human craniofacial diversity observed in different geographical regions. Med Hypotheses. 2004, 63: 974-979. 10.1016/j.mehy.2004.05.003.

    Article  CAS  PubMed  Google Scholar 

  18. Young JH, Chang YC, Kim JD, Chretien J, Klag MJ, Levine MA, Ruff CB, Wang N, Chakravarti A: Differential susceptibility to hypertension is due to selection during the out-of-Africa expansion. PLoS Genet. 2005, 1: e82-10.1371/journal.pgen.0010082.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Thompson EE, Kuttab-Boulos H, Witonsky D, Yang L, Roe BA, Di Rienzo A: CYP3A variation and the evolution of salt-sensitivity variants. Am J Hum Genet. 2004, 75: 1059-1069. 10.1086/426406.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Crespi B, Summers K, Dorus S: Adaptive evolution of genes underlying schizophrenia. Proc Biol Sci. 2007, 274: 2801-2810. 10.1098/rspb.2007.0876.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Parra EJ: Human pigmentation variation: evolution, genetic basis, and implications for public health. Am J Phys Anthropol. 2007, 85-105. 10.1002/ajpa.20727. Suppl 45

  22. Hancock AM, Witonsky DB, Gordon AS, Eshel G, Pritchard JK, Coop G, Di Rienzo A: Adaptations to climate in candidate genes for common metabolic disorders. PLoS Genet. 2008, 4: e32-10.1371/journal.pgen.0040032.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Eisen A, Calne D: Amyotrophic lateral sclerosis, Parkinson's disease and Alzheimer's disease: phylogenetic disorders of the human neocortex sharing many characteristics. Can J Neurol Sci. 1992, 19: 117-123.

    CAS  PubMed  Google Scholar 

  24. Karami S, Boffetta P, Stewart P, Rothman N, Hunting KL, Dosemeci M, Berndt SI, Brennan P, Chow W, Moore LE: Occupational sunlight exposure and risk of renal cell carcinoma. Cancer. 2010

    Google Scholar 

  25. Kinney DK, Teixeira P, Hsu D, Napoleon SC, Crowley DJ, Miller A, Hyman W, Huang E: Relation of schizophrenia prevalence to latitude, climate, fish consumption, infant mortality, and skin color: a role for prenatal vitamin d deficiency and infections?. Schizophr Bull. 2009, 35: 582-595. 10.1093/schbul/sbp023.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, Myers RM: Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008, 319: 1100-1104. 10.1126/science.1153717.

    Article  CAS  PubMed  Google Scholar 

  27. Weir BS, Cockerham CC: Estimating F-Statistics for the Analysis of Population Structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.

    Article  Google Scholar 

  28. Amato R, Pinelli M, Monticelli A, Marino D, Miele G, Cocozza S: Genome-wide scan for signatures of human population differentiation and their relationship with natural selection, functional pathways and diseases. PLoS ONE. 2009, 4: e7927-10.1371/journal.pone.0007927.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Barreiro LB, Laval G, Quach H, Patin E, Quintana-Murci L: Natural selection has driven population differentiation in modern humans. Nat Genet. 2008, 40: 340-345. 10.1038/ng.78.

    Article  CAS  PubMed  Google Scholar 

  30. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Bauer S, Grossmann S, Vingron M, Robinson PN: Ontologizer 2.0--a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics. 2008, 24: 1650-1651. 10.1093/bioinformatics/btn250.

    Article  CAS  PubMed  Google Scholar 

  32. Bauer S, Gagneur J, Robinson PN: GOing Bayesian: model-based gene set analysis of genome-scale data. Nucleic Acids Res. 2010

    Google Scholar 

  33. Taylor BV, Lucas RM, Dear K, Kilpatrick TJ, Pender MP, van der Mei IA, Chapman C, Coulthard A, Dwyer T, McMichael AJ, Valery PC, Williams D, Ponsonby A: Latitudinal variation in incidence and type of first central nervous system demyelinating events. Mult Scler. 2010

    Google Scholar 

  34. Allen NC, Bagade S, McQueen MB, Ioannidis JPA, Kavvoura FK, Khoury MJ, Tanzi RE, Bertram L: Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet. 2008, 40: 827-834. 10.1038/ng.171.

    Article  CAS  PubMed  Google Scholar 

  35. Lill C, McQueen MB, Roehr J, Bagade S, Schjeide B, Zipp F, Bertram L: The MSGene Database. 2010

    Google Scholar 

  36. Lill C, Bagade S, McQueen MB, Roehr J, Kavvoura FK, Schjeide B, Allen NC, Tanzi RE, Khoury MJ, Ioannidis JPA, Bertram L: The PDGene Database. 2010

    Google Scholar 

  37. Bertram L, McQueen MB, Mullin K, Blacker D, Tanzi RE: Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet. 2007, 39: 17-23. 10.1038/ng1934.

    Article  CAS  PubMed  Google Scholar 

  38. Wang T, Tavera-Mendoza LE, Laperriere D, Libby E, MacLeod NB, Nagai Y, Bourdeau V, Konstorum A, Lallemant B, Zhang R, Mader S, White JH: Large-scale in silico and microarray-based identification of direct 1,25-dihydroxyvitamin D3 target genes. Mol Endocrinol. 2005, 19: 2685-2695. 10.1210/me.2005-0106.

    Article  CAS  PubMed  Google Scholar 

  39. Cober ER, Morrison MJ: Regulation of seed yield and agronomic characters by photoperiod sensitivity and growth habit genes in soybean. Theor Appl Genet. 2010, 120: 1005-1012. 10.1007/s00122-009-1228-6.

    Article  CAS  PubMed  Google Scholar 

  40. Zhang D, Zhang H, Wang M, Sun J, Qi Y, Wang F, Wei X, Han L, Wang X, Li Z: Genetic structure and differentiation of Oryza sativa L. in China revealed by microsatellites. Theor Appl Genet. 2009, 119: 1105-1117. 10.1007/s00122-009-1112-4.

    Article  CAS  PubMed  Google Scholar 

  41. Guo Z, Song Y, Zhou R, Ren Z, Jia J: Discovery, evaluation and distribution of haplotypes of the wheat Ppd-D1 gene. New Phytol. 2009

    Google Scholar 

  42. Van 't Land J, Van Putten WF, Villarroel H, Kamping A, Van Delden W: Latitudinal variation for two enzyme loci and an inversion polymorphism in Drosophila melanogaster from Central and South America. Evolution. 2000, 54: 201-209.

    PubMed  Google Scholar 

  43. Johnsen A, Fidler AE, Kuhn S, Carter KL, Hoffmann A, Barr IR, Biard C, Charmantier A, Eens M, Korsten P, Siitari H, Tomiuk J, Kempenaers B: Avian Clock gene polymorphism: evidence for a latitudinal cline in allele frequencies. Mol Ecol. 2007, 16: 4867-4880. 10.1111/j.1365-294X.2007.03552.x.

    Article  CAS  PubMed  Google Scholar 

  44. Lucotte G, Mercier G: The population distribution of the Met allele at the PRNP129 polymorphism (a high risk factor for Creutzfeldt-Jakob disease) in various regions of France and in West Europe. Infect Genet Evol. 2005, 5: 141-144. 10.1016/j.meegid.2004.08.002.

    Article  CAS  PubMed  Google Scholar 

  45. Greene LS, Bottini N, Borgiani P, Gloria-Bottini F: Acid phosphatase locus 1 (ACP1): Possible relationship of allelic variation to body size and human population adaptation to thermal stress-A theoretical perspective. Am J Hum Biol. 2000, 12: 688-701. 10.1002/1520-6300(200009/10)12:5<688::AID-AJHB14>3.0.CO;2-C.

    Article  PubMed  Google Scholar 

  46. Treves TA, de Pedro-Cuesta J: Parkinsonism mortality in the US, 1. Time and space distribution. Acta Neurol Scand. 1991, 84: 389-397. 10.1111/j.1600-0404.1991.tb04975.x.

    Article  CAS  PubMed  Google Scholar 

  47. McGrath J, Saha S, Chant D, Welham J: Schizophrenia: a concise overview of incidence, prevalence, and mortality. Epidemiol Rev. 2008, 30: 67-76. 10.1093/epirev/mxn001.

    Article  PubMed  Google Scholar 

  48. Torrey EF: Prevalence studies in schizophrenia. Br J Psychiatry. 1987, 150: 598-608. 10.1192/bjp.150.5.598.

    Article  CAS  PubMed  Google Scholar 

  49. Lu Z, Chen TC, Zhang A, Persons KS, Kohn N, Berkowitz R, Martinello S, Holick MF: An evaluation of the vitamin D3 content in fish: Is the vitamin D content adequate to satisfy the dietary requirement for vitamin D?. J Steroid Biochem Mol Biol. 2007, 103: 642-644. 10.1016/j.jsbmb.2006.12.010.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Jablonski NG, Chaplin G: The evolution of human skin coloration. J Hum Evol. 2000, 39: 57-106. 10.1006/jhev.2000.0403.

    Article  CAS  PubMed  Google Scholar 

  51. Kuningas M, Mooijaart SP, Jolles J, Slagboom PE, Westendorp RGJ, van Heemst D: VDR gene variants associate with cognitive function and depressive symptoms in old age. Neurobiol Aging. 2009, 30: 466-473. 10.1016/j.neurobiolaging.2007.07.001.

    Article  CAS  PubMed  Google Scholar 

  52. Kesby JP, Cui X, Ko P, McGrath JJ, Burne THJ, Eyles DW: Developmental vitamin D deficiency alters dopamine turnover in neonatal rat forebrain. Neurosci Lett. 2009, 461: 155-158. 10.1016/j.neulet.2009.05.070.

    Article  CAS  PubMed  Google Scholar 

  53. Grecksch G, Rüthrich H, Höllt V, Becker A: Transient prenatal vitamin D deficiency is associated with changes of synaptic plasticity in the dentate gyrus in adult rats. Psychoneuroendocrinology. 2009, 34 (Suppl 1): S258-264. 10.1016/j.psyneuen.2009.07.004.

    Article  CAS  PubMed  Google Scholar 

  54. Holmes VA, Barnes MS, Alexander HD, McFaul P, Wallace JMW: Vitamin D deficiency and insufficiency in pregnant women: a longitudinal study. Br J Nutr. 2009, 102: 876-881. 10.1017/S0007114509297236.

    Article  CAS  PubMed  Google Scholar 

  55. Lucas RM, Repacholi MH, McMichael AJ: Is the current public health message on UV exposure correct?. Bull World Health Organ. 2006, 84: 485-491. 10.2471/BLT.05.026559.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Wang T, Tavera-Mendoza LE, Laperriere D, Libby E, MacLeod NB, Nagai Y, Bourdeau V, Konstorum A, Lallemant B, Zhang R, Mader S, White JH: Large-scale in silico and microarray-based identification of direct 1,25-dihydroxyvitamin D3 target genes. Mol Endocrinol. 2005, 19: 2685-2695. 10.1210/me.2005-0106.

    Article  CAS  PubMed  Google Scholar 

  57. Levi-Montalcini R, Hamburger V: Selective growth stimulating effects of mouse sarcoma on the sensory and sympathetic nervous system of the chick embryo. J Exp Zool. 1951, 116: 321-361. 10.1002/jez.1401160206.

    Article  CAS  PubMed  Google Scholar 

  58. Weickert CS, Ligons DL, Romanczyk T, Ungaro G, Hyde TM, Herman MM, Weinberger DR, Kleinman JE: Reductions in neurotrophin receptor mRNAs in the prefrontal cortex of patients with schizophrenia. Mol Psychiatry. 2005, 10: 637-650. 10.1038/sj.mp.4001678.

    Article  CAS  PubMed  Google Scholar 

  59. Schramm M, Falkai P, Feldmann N, Knable MB, Bayer TA: Reduced tyrosine kinase receptor C mRNA levels in the frontal cortex of patients with schizophrenia. Neurosci Lett. 1998, 257: 65-68. 10.1016/S0304-3940(98)00807-6.

    Article  CAS  PubMed  Google Scholar 

  60. Otnaess MK, Djurovic S, Rimol LM, Kulle B, Kähler AK, Jönsson EG, Agartz I, Sundet K, Hall H, Timm S, Hansen T, Callicott JH, Melle I, Werge T, Andreassen OA: Evidence for a possible association of neurotrophin receptor (NTRK-3) gene polymorphisms with hippocampal function and schizophrenia. Neurobiol Dis. 2009, 34: 518-524. 10.1016/j.nbd.2009.03.011.

    Article  CAS  PubMed  Google Scholar 

  61. Shibayama E, Koizumi H: Cellular localization of the Trk neurotrophin receptor family in human non-neuronal tissues. Am J Pathol. 1996, 148: 1807-1818.

    PubMed Central  CAS  PubMed  Google Scholar 

  62. Nakanishi T, Takahashi K, Aoki C, Nishikawa K, Hattori T, Taniguchi S: Expression of nerve growth factor family neurotrophins in a mouse osteoblastic cell line. Biochem Biophys Res Commun. 1994, 198: 891-897. 10.1006/bbrc.1994.1127.

    Article  CAS  PubMed  Google Scholar 

  63. Asaumi K, Nakanishi T, Asahara H, Inoue H, Takigawa M: Expression of neurotrophins and their receptors (TRK) during fracture healing. Bone. 2000, 26: 625-633. 10.1016/S8756-3282(00)00281-7.

    Article  CAS  PubMed  Google Scholar 

  64. Hakak Y, Walker JR, Li C, Wong WH, Davis KL, Buxbaum JD, Haroutunian V, Fienberg AA: Genome-wide expression analysis reveals dysregulation of myelination-related genes in chronic schizophrenia. Proc Natl Acad Sci USA. 2001, 98: 4746-4751. 10.1073/pnas.081071198.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Eastwood SL, Salih T, Harrison PJ: Differential expression of calcineurin A subunit mRNA isoforms during rat hippocampal and cerebellar development. Eur J Neurosci. 2005, 22: 3017-3024. 10.1111/j.1460-9568.2005.04518.x.

    Article  PubMed  Google Scholar 

  66. Martins-de-Souza D, Gattaz WF, Schmitt A, Rewerts C, Marangoni S, Novello JC, Maccarrone G, Turck CW, Dias-Neto E: Alterations in oligodendrocyte proteins, calcium homeostasis and new potential markers in schizophrenia anterior temporal lobe are revealed by shotgun proteome analysis. J Neural Transm. 2009, 116: 275-289. 10.1007/s00702-008-0156-y.

    Article  CAS  PubMed  Google Scholar 

  67. Sun L, Zhu L, Zaidi N, Yang G, Moonga BS, Abe E, Iqbal J, Epstein S, Blair HC, Huang CL, Zaidi M: Cellular and molecular consequences of calcineurin A alpha gene deletion. Ann N Y Acad Sci. 2007, 1116: 216-226. 10.1196/annals.1402.050.

    Article  CAS  PubMed  Google Scholar 

  68. Awumey EM, Moonga BS, Sodam BR, Koval AP, Adebanjo OA, Kumegawa M, Zaidi M, Epstein S: Molecular and functional evidence for calcineurin-A alpha and beta isoforms in the osteoclast: novel insights into cyclosporin A action on bone resorption. Biochem Biophys Res Commun. 1999, 254: 248-252. 10.1006/bbrc.1998.9785.

    Article  CAS  PubMed  Google Scholar 

  69. Kitagawa H, Fujiki R, Yoshimura K, Mezaki Y, Uematsu Y, Matsui D, Ogawa S, Unno K, Okubo M, Tokita A, Nakagawa T, Ito T, Ishimi Y, Nagasawa H, Matsumoto T, Yanagisawa J, Kato S: The chromatin-remodeling complex WINAC targets a nuclear receptor to promoters and is impaired in Williams syndrome. Cell. 2003, 113: 905-917. 10.1016/S0092-8674(03)00436-7.

    Article  CAS  PubMed  Google Scholar 

  70. Khavari PA, Peterson CL, Tamkun JW, Mendel DB, Crabtree GR: BRG1 contains a conserved domain of the SWI2/SNF2 family necessary for normal mitotic growth and transcription. Nature. 1993, 366: 170-174. 10.1038/366170a0.

    Article  CAS  PubMed  Google Scholar 

  71. Lessard J, Wu JI, Ranish JA, Wan M, Winslow MM, Staahl BT, Wu H, Aebersold R, Graef IA, Crabtree GR: An essential switch in subunit composition of a chromatin remodeling complex during neural development. Neuron. 2007, 55: 201-215. 10.1016/j.neuron.2007.06.019.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Koga M, Ishiguro H, Yazaki S, Horiuchi Y, Arai M, Niizato K, Iritani S, Itokawa M, Inada T, Iwata N, Ozaki N, Ujike H, Kunugi H, Sasaki T, Takahashi M, Watanabe Y, Someya T, Kakita A, Takahashi H, Nawa H, Muchardt C, Yaniv M, Arinami T: Involvement of SMARCA2/BRM in the SWI/SNF chromatin-remodeling complex in schizophrenia. Hum Mol Genet. 2009, 18: 2483-2494. 10.1093/hmg/ddp166.

    Article  CAS  PubMed  Google Scholar 

  73. Flowers S, Nagl NG, Beck GR, Moran E: Antagonistic roles for BRM and BRG1 SWI/SNF complexes in differentiation. J Biol Chem. 2009, 284: 10067-10075. 10.1074/jbc.M808782200.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Reyes JC, Barra J, Muchardt C, Camus A, Babinet C, Yaniv M: Altered control of cellular proliferation in the absence of mammalian brahma (SNF2alpha). EMBO J. 1998, 17: 6979-6991. 10.1093/emboj/17.23.6979.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. Lamont EW, Legault-Coutu D, Cermakian N, Boivin DB: The role of circadian clock genes in mental disorders. Dialogues Clin Neurosci. 2007, 9: 333-342.

    PubMed  Google Scholar 

  76. Mathias D, Jacky L, Bradshaw WE, Holzapfel CM: Geographic and developmental variation in expression of the circadian rhythm gene, timeless, in the pitcher-plant mosquito, Wyeomyia smithii. J Insect Physiol. 2005, 51: 661-667. 10.1016/j.jinsphys.2005.03.011.

    Article  CAS  PubMed  Google Scholar 

  77. Palsdottir A, Helgason A, Palsson S, Bjornsson HT, Bragason BT, Gretarsdottir S, Thorsteinsdottir U, Olafsson E, Stefansson K: A Drastic Reduction in the Life Span of Cystatin C L68Q Carriers Due to Life-Style Changes during the Last Two Centuries. PLoS Genet. 2008, 4: 10.1371/journal.pgen.1000099.

    Google Scholar 

  78. Huffaker SJ, Chen J, Nicodemus KK, Sambataro F, Yang F, Mattay V, Lipska BK, Hyde TM, Song J, Rujescu D, Giegling I, Mayilyan K, Proust MJ, Soghoyan A, Caforio G, Callicott JH, Bertolino A, Meyer-Lindenberg A, Chang J, Ji Y, Egan MF, Goldberg TE, Kleinman JE, Lu B, Weinberger DR: A primate-specific, brain isoform of KCNH2 affects cortical physiology, cognition, neuronal repolarization and risk of schizophrenia. Nat Med. 2009, 15: 509-518. 10.1038/nm.1962.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  79. Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz MH, Hansen NF, Durand EY, Malaspinas A, Jensen JD, Marques-Bonet T, Alkan C, Prüfer K, Meyer M, Burbano HA, Good JM, Schultz R, Aximu-Petri A, Butthof A, Höber B, Höffner B, Siegemund M, Weihmann A, Nusbaum C, Lander ES, Russ C, Novod N, Affourtit J, Egholm M, Verna C, Rudan P, Brajkovic D, Kucan Z, Gusic I, Doronichev VB, Golovanova LV, Lalueza-Fox C, de la Rasilla M, Fortea J, Rosas A, Schmitz RW, Johnson PLF, Eichler EE, Falush D, Birney E, Mullikin JC, Slatkin M, Nielsen R, Kelso J, Lachmann M, Reich D, Pääbo S: A draft sequence of the Neandertal genome. Science. 2010, 328: 710-722. 10.1126/science.1188021.

    Article  CAS  PubMed  Google Scholar 

  80. Balloux F, Handley LL, Jombart T, Liu H, Manica A: Climate shaped the worldwide distribution of human mitochondrial DNA sequence variation. Proc Biol Sci. 2009, 276: 3447-3455. 10.1098/rspb.2009.0752.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  81. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, Myers RM: Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008, 319: 1100-1104. 10.1126/science.1153717.

    Article  CAS  PubMed  Google Scholar 

  82. Myles S, Somel M, Tang K, Kelso J, Stoneking M: Identifying genes underlying skin pigmentation differences among human populations. Hum Genet. 2007, 120: 613-621. 10.1007/s00439-006-0256-4.

    Article  CAS  PubMed  Google Scholar 

  83. Bussey KJ, Kane D, Sunshine M, Narasimhan S, Nishizuka S, Reinhold WC, Zeeberg B, Ajay W, Weinstein JN: MatchMiner: a tool for batch navigation among gene and gene product identifiers. Genome Biol. 2003, 4: R27-10.1186/gb-2003-4-4-r27.

    Article  PubMed Central  PubMed  Google Scholar 

  84. R Development Core Team: R: A Language and Environment for Statistical Computing. 2009, Vienna, Austria

    Google Scholar 

  85. Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57. 10.1038/nprot.2008.211.

    Article  CAS  Google Scholar 

  86. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Michelle Kutzner for her help in the revision of the manuscript. RA is a recipient of a fellowship by Doctorate of Computational Biology and Bioinformatics, University ''Federico II'', Naples.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Roberto Amato.

Additional information

Authors' contributions

RA, MP, GM and SC participated in the design of the study and performed the statistical analysis. AM, GM and SC conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2010_1566_MOESM1_ESM.XLS

Additional file 1: LRGs. Excel spreadsheet containing the list of Latitude Related Genes with FST and correlation value for each SNP. (XLS 312 KB)

12862_2010_1566_MOESM2_ESM.PDF

Additional file 2: supp_data. Additional information about the enrichment of SNP classes and the enrichment of LRGs in various lists of schizophrenia related genes. (PDF 319 KB)

Additional file 3: overlaps. Excel spreadsheet containing the list of genes in each two-ways intersection. (XLS 47 KB)

Additional file 4: VDRGs. Excel spreadsheet containing the list of vitamin D related genes. (XLS 74 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Amato, R., Pinelli, M., Monticelli, A. et al. Schizophrenia and vitamin D related genes could have been subject to latitude-driven adaptation. BMC Evol Biol 10, 351 (2010). https://doi.org/10.1186/1471-2148-10-351

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-10-351

Keywords