Alexander V. Igoshin
Konstantin V. Gunbin2,3,4
Nikolay S. Yudin- 1Sector of the Genetics of Industrial Microorganisms, The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch, The Russian Academy of Sciences, Novosibirsk, Russia
- 2Center of Brain Neurobiology and Neurogenetics, The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch, The Russian Academy of Sciences, Novosibirsk, Russia
- 3V. Zelman Institute for Medicine and Psychology Novosibirsk State University, Novosibirsk, Russia
- 4Center for Mitochondrial Functional Genomics, Institute of Living Systems, Immanuel Kant Baltic Federal University, Kaliningrad, Russia
- 5Laboratory of Livestock Molecular Genetics and Breeding, The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch, The Russian Academy of Sciences, Novosibirsk, Russia
- 6Laboratory of Human Molecular Genetics, The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch, The Russian Academy of Sciences, Novosibirsk, Russia
Dispersal of Homo sapiens across the globe during the last 200,000 years was accompanied by adaptation to local climatic conditions, with severe winter temperatures being probably one of the most significant selective forces. The TRPM8 gene codes for a cold-sensing ion channel, and adaptation to low temperatures is the major determinant of its molecular evolution. Here, our aim was to search for signatures of cold climate adaptation in TRPM8 gene using a combined data set of 19 populations of East Asian ancestry from the 1000 Genomes Project and Human Genome Diversity Project. As a result, out of a total of 60 markers under study, none showed significant association with the average winter temperatures at the locations of the studied populations considering the multiple testing thresholds. This might suggest that the principal mode of TRPM8 evolution may be different from widespread models, where adaptive alleles are additive, dominant or recessive, at least in populations with the predominant East Asian component. For example, evolution by means of selectively preferable epistatic interactions among amino acids may have taken place. Despite the lack of strong signals of association, however, a very promising single nucleotide polymorphism (SNP) was found. The SNP rs7577262 is considered the best candidate based on its allelic correlations with winter temperatures, signatures of selective sweep and physiological evidences. The second top SNP, rs17862920, may participate in adaptation as well. Additionally, to assist in interpreting the nominal associations, the other markers reached, we performed SNP prioritization based on functional evidences found in literature and on evolutionary conservativeness.
Introduction
Recent paleoanthropological evidences show the presence of anatomically modern humans in Africa as early as 300 kya (Hublin et al., 2017), with the earliest known “Out of Africa” migration event dating back to 200 kya (Hershkovitz et al., 2018). Dispersal of Homo sapiens across the globe during the last 200,000 years was accompanied by adaptation to local environments. Spatial variations in selective pressures have ultimately led to observable geographic distribution of many physiological and anatomical traits in present-day humans. For example, the low level of UV radiation at higher latitudes is now considered to be the major cause of evolution of depigmented skin (Jablonski and Chaplin, 2010).
Since the mid-2000s, there has been significant progress in genotyping technologies, followed by publicly available databases of human genetic variation. This circumstance helped the population geneticists to discover signatures of human local adaptation from genome-wide genotyping data. Human microevolution driven by the action of low temperatures has long been attracting attention of the scientific community. Now, a number of studies have been dedicated to this issue both at the genome level (e.g., Hancock et al., 2011a; Cardona et al., 2014; Valverde et al., 2015) and at the level of selected regions or genes (Hancock et al., 2008; Ohashi et al., 2011; Hancock et al., 2011b; Sazzini et al., 2014; Quagliarello et al., 2017).
Probably the best-known gene in terms of its possible role in adaptation to cold climate is TRPM8 located on human chromosome 2. This gene codes for ion channel functioning as a thermal sensor, detecting temperatures in the range from 15 to 30°C (Fernández et al., 2011). There are evidences supporting its physiological role in thermoregulation, and in fact, TRPM8 is the only well-established cold receptor in mammals (Bautista et al., 2007; Colburn et al., 2007; Dhaka et al., 2007). Besides these, there are data on associations of its single nucleotide polymorphisms (SNPs) with sensitivity to cold (Kozyreva et al., 2011), the respiratory system response to cooling (Kozyreva et al., 2014), blood lipids, and anthropometric parameters in humans (Potapova et al., 2014). The TRPM8 gene was suggested to underlie genetic adaptation to cold in ground squirrel and hamster (Matos-Cruz et al., 2017), sheep (Fariello et al., 2014; Liu et al., 2016), and humans (Cardona et al., 2014; Key et al., 2018). According to modern views, adaptation to low temperatures is the major determinant of TRPM8 molecular evolution (Myers et al., 2009; Majhi et al., 2015).
To our knowledge, a study by Key and colleagues (2018) is the only one to use environmental data to search for signatures of cold climate adaptation in the TRPM8 gene. The authors used latitudes and annual average temperatures at the locations of the populations of the Old World as predictors for SNP allele frequency distributions. They found evidences that SNP rs10166942 had undergone climate-mediated selection, which raised its derived allele frequency from south to north.
In our opinion, focusing on closely related populations is preferable to using large population sets for the following reasons. First, the ability to survive in a severely cold climate is supposed to be highly polygenic, as many biological processes like vasoconstriction, nonshivering thermogenesis, regulation of adipocyte differentiation, and thermoception are expected to be involved. It is known that the adaptation to cold can be associated with quite different genetic bases (Yudin et al., 2017). Because different branches of Homo sapiens are likely to have had distinct genetic background before and during the process of climate-driven selection, it is possible that in phylogenetically distant groups, adaptation may have recruited different genes. Second, even in the case when selection is acting on the same gene, variants involved in adaptation may differ in different branches. Our supposition is supported by the example of variants associated with lactase persistence. Thus, within European populations, the activity of the lactase enzyme in adulthood is connected with the C/T-13910 variant in the enhancer region of the LCT gene, whereas in sub-Saharan Africa, this trait is mainly correlated with the presence of the G/C-14010 mutation (Tishkoff et al., 2007). Therefore, it would be sensible to search for microevolution in clusters of related populations.
For the above reasons, the aim of this study was to search for signatures of adaptation to low temperatures in the TRPM8 gene under various null hypotheses of population structure and dynamics using a combined data set of 19 populations of the East Asian ancestry from Human Genome Diversity Project (HGDP) and 1000 Genomes (1000G) Project with the assistance of environmental correlation analysis techniques. Locations of chosen populations are characterized by a large range of average winter temperatures (−37–+27°C), implying substantial differences in selection pressures.
Materials and Methods
Genotypic and Environmental Data
In this study, we used genotypic data on 656 individuals from 19 HGDP (Cann et al., 2002) and 1000G Project (1000 Genomes Project Consortium et al., 2010) populations (Supplementary Table 1) having predominantly the East Asian genetic component. Data on SNPs belonging to TRPM8 gene were obtained from NCBI dbSNP (https://www.ncbi.nlm.nih.gov/snp/), resulting in 60 polymorphic markers (minor allele frequency >0.01) being at the intersection of HGDP and 1000G sets. The missing genotypes in HGDP data were imputed using fastPHASE v.1.4.8 software (Scheet and Stephens, 2006) with default parameters. Genotypic information from HGDP and VCF formats was combined by using a self-made Python 3 script, so that inconsistency of DNA strands between databases (if that was the case) was resolved using 1000G VCF as a reference. Besides TRPM8 SNPs, we used 5 Mb regions upstream and downstream this gene (1,309 markers at r2 < 0.7) to infer the phylogenetic tree used in PGLS, to estimate the covariance matrix used in Bayenv2-BLM and Bayenv2-SRC, to correct for background levels of the population structure in LFMM and to make population inferences in BayScEnv more precise. Also, this set of SNPs was used to construct a null distribution for empirical p-value calculation. Information on latitudes and longitudes for 1000G populations was taken from Key et al. (2018). Latitudes and longitudes for HGDP populations were taken from Cann et al. (2002), and average winter temperature values were obtained from ClimateCharts.net database (https://climatecharts.net) using the corresponding coordinates. We believe the average winter temperature to be a more pertinent predictor for distribution of cold-adaptive alleles than the annual average temperature, as regions with a continental climate may have cold winters and hot summers.
Construction of a Phylogenetic Tree and Statistical Analysis
As our sample consisted of phylogenetically close populations, we first performed conventional Spearman’s rank correlation test not accounting for the sample structure (Spearman, 1904) at the population (i.e., using allele frequencies) and individual’s (i.e., using allele dosages) levels.
PGLS analysis (R package ‘ape,’ Paradis et al., 2004) was carried out at the individual’s level using the simplest Brownian motion model. We selected this type of analysis because it is an opposite alternative (phyletic evolution) to the conventional Spearman’s rank correlation test. The phylogenetic tree used in this test was reconstructed with IQ-TREE v. 1.5.5 subprogram ModelFinder (Kalyaanamoorthy et al., 2017) based on the best nucleotide substitution model.
Our primary aim was to test the association between the climatic factor and allele frequencies. For this purpose, we chose two independent approaches: the Bayesian linear model from Bayenv2 (further referred to as Bayenv2-BLM) software (Günther and Coop, 2013) and LFMM (Frichot et al., 2013), each building a regression model relating allele frequencies to environmental values. To minimize the problem of false-positive associations between allele frequencies and environment because of the population structure, the above methods take into account allele frequency correlations across populations while performing the analysis by various ways.
In addition, we used Spearman’s rank correlation test from Bayenv2 (further referred to as Bayenv2-SRC) that uses allele frequencies standardized to have no covariance. It is less powerful than Bayenv2-BLM but more robust to outliers and can detect monotonic relationships.
BayScEnv test (de Villemereuil and Gaggiotti, 2015) was used as an alternative to Bayenv2-BLM and LFMM. This method assumes that all populations are independent and exchange genes through the limited migrant pool; it includes a locus-specific effect unrelated to the environmental variable, taking into consideration locus-specific deviations from a neutral model. BayScEnv software was also used to calculate Fst distances for each locus averaged over populations (Supplementary Table 2).
In addition to the correlation techniques, we tried the XP-CLR test (Chen et al., 2010) as a complementary approach. This test is designed for detecting selective sweeps on the basis of joint modeling of the multilocus allele frequency differentiation between two populations. The method does not require information on the ancestral/derived status at each SNP (Chen et al., 2010; Vatsiou et al., 2016).
For more details on phylogenetic tree construction and statistical analysis, please see the Supplementary Material.
SNP Prioritization
We prioritized SNPs based on three types of evidences found in literature, “association with trait relevant to survival in a cold climate,” “evidences for cold-mediated selection,” and “association with any other phenotype or risk.” These categories were given weights 3, 2, and 1, respectively, and a score for each SNP was summarized (Supplementary Table 3). The key assumption behind this prioritization is that because of the pleiotropic nature of TRPM8 gene, allelic substitutions having any functional manifestation may potentially have more chances of affecting survival in cold climate conditions than those not having any known effects.
Additionally, we obtained the PhyloP100way vertebrate conservation score for each SNP (Supplementary Table 4) from UCSC Genome Browser (Casper et al., 2018). Currently, it is commonly thought that the genetic drift plays a minor role in the evolution of conservative sites, and relatively rare allele replacements occurring therein are mostly driven by positive selection (Andolfatto, 2005; Cai et al., 2009; Halligan et al., 2011; Bazykin and Kondrashov, 2012). Therefore, significant allelic correlation with environmental gradient supported by a high conservation score promises to be the true sign of local adaptation.
Results
Contrary to what we had expected, only three SNPs out of a total of 60 markers under study showed nominally significant association with the average winter temperatures at the locations of the studied populations by any two types of analysis (the results of tests carried out using default/recommended parameters are shown in Table 1; for full results, see Supplementary Tables 5 and 6). When considering the multiple testing threshold, however, none of them is significant (adjusted p values not shown).
Table 1 Genic and upstream TRPM8 variants showing nominally significant (in bold) association in at least two correlation tests with default/recommended parameters (K = 2 for LFMM and pi = 0.1/p = 0.5 for BayScEnv).
SNP rs11682848 was previously reported as associated with the prognosis of colorectal cancer (Walther, 2010). Interestingly, such connection of climate-associated loci with cancer has already been noticed by other researchers (Hancock et al., 2011a). Furthermore, it has been recently shown by combination of 247 genome-wide association studies that cold selected genes are enriched with cancer-associated genes (Voskarides, 2018). Ironically, SNP rs11682848 has the lowest conservation score among 60 markers under study. This means that either rs11682848 is being a false-positive finding or being linked to some functional variant.
SNP rs17862920 has evidences of associations with migraine susceptibility (Freilinger et al., 2012; Meng et al., 2018). The rs17862920-С allele predisposing to migraine is more prevalent in northern latitudes. Also, rs17862920 has been shown to be associated with sensing cold pain in Finnish and Norwegian individuals, with C allele carriers being more susceptible (Kaunisto et al., 2013). Migraine has been reported to be related to increased pain perception of nonnoxious cold temperatures (Burstein et al., 2000). Unlike rs11682848, PhyloP100way conservation score for rs17862920 is more promising and has a rank of 20/60 while still being negative. It is possible that allele substitution in rs17862920 has a functional effect. Thus, rs17862920 was predicted to regulate TRPM8 transcription by TFsearch and GoldenPath in F-SNP bioinformatics tool (Ghosh et al., 2013).
As for rs6723922, this SNP is a genetic risk factor for severe cutaneous adverse drug reactions (Park et al., 2018). One could hypothesize that there is a certain mechanism underlying both the altered cold sensation and the increased cutaneous susceptibility to chemicals. It could be no surprise given that the TRPM8 channel is activated by a variety of chemical ligands (Beccari et al., 2017). Like rs11682848, SNP rs6723922 has a low conservation score (the rank of 55/60), implying conclusions for this marker similar to those for rs11682848.
The results of the XP-CLR test are more encouraging (Supplementary Figure S1). It appears that there is a pronounced trend for several pairs of populations to show the signature of a selective sweep 10 ± 6 Kb upstream from the TRPM8 gene. The direction of selection in this region is seen when reversing tested and reference populations in pairs (e.g., compare “JPT vs. KHV” and “KHV vs. JPT”). The strongest XP-CLR peaks within this putative sweep are mainly located near rs10929317 and rs7577262 SNP loci. The former was removed from the analysis because of high LD with rs17862920: r2 = 0.966/D’ = 0.995 in East Asian populations (LDlink tool; Machiela and Chanock, 2015) and is therefore expected to be as significant as rs17862920. The latter was used in the control set of 1,309 markers. Surprisingly, this SNP demonstrates significant association with our climatic variable in almost all of the correlation tests (Table 1). Also, rs7577262 has been reported to be associated with susceptibility to migraine (Anttila et al., 2013) and blood pressure response to the cold pressor test (He et al., 2013).
Discussion
A variety of facts have led us to think of TRPM8 gene as being under intense positive selection. We expected that the large amount of SNPs in TRPM8 would demonstrate strong signals because of being under selection immediately or being linked to some causal variants. However, this is not the case. Furthermore, SNPs detected do not pass the corrected threshold, considering multiple testing. Among possible explanations are the following hypotheses:
1. The predominant mode of TRPM8 evolution may be different from the widespread models exploited by our tests, where adaptation is assumed to be mediated by selecting alleles with additive or at least recessive/dominant trait coding. However, it is suggested that epistatic interactions may play a role in the evolution of thermoTRP channels (Saito and Tominaga, 2017). Therefore, conformational epistasis-based evolution, where some epistatic interactions among amino acids are preferred, might have resulted in the inability of approaches we used to detect strong signatures of selection in the TRPM8 gene.
2. The level of allele frequency variation (see Supplementary Table 7 for minor allele frequency distributions) in the populations studied is not sufficient to robustly discriminate the loci under selection. Thus, the averaged Fst distances for populations used by Key et al. (2018) are much higher than those for our data set (Supplementary Table 2).
3. Only 5 out of 19 locations of populations from our data set have the average winter temperature below −10°C. It is possible that the underrepresentation of northern populations in this study might lead to insufficient signal strength. Further accumulation of open access data on genetic variation in the north would help in detecting loci under selection.
Despite the lack of strong signals of association, however, a very promising candidate SNP was found. SNP rs7577262 is 7.1 kb upstream of the transcription start site for TRPM8 mRNA, implying its possible involvement in transcriptional regulation. In addition to correlations and signatures of sweep, physiological data contribute equally to the evidences in favor of selection acting on rs7577262. The rs7577262-G allele is associated with a higher blood pressure response to the cold pressor test (He et al., 2013). It is known that the blood pressure response to the cold pressor test primarily stems from alpha-adrenergically mediated peripheral vasoconstriction (Leppäluoto and Hassi, 1991; Larra et al., 2015), which is, in turn, one of the basic mechanisms of cold adaptation (Daanen and Lichtenbelt, 2016). Given that this allele is more prevalent in northern latitudes (Figure 1), its adaptive role may be assumed.
Figure 1 Geographic distribution of allele frequencies for rs7577262 polymorphism in populations of East Asian ancestry. Average winter temperatures at the locations of the populations studied are shown. 1000 Genomes population: CHB, Han Chinese in Beijing, China; JPT, Japanese in Tokyo, Japan; CDX, Chinese Dai in Xishuangbanna, China; CHS, Han Chinese South, China; KHV, Kinh in Ho Chi Minh City, Vietnam.
Another SNP, rs17862920, is linked with rs7577262 (r2 = 0.59 in East Asians). Probably, this is the reason for the correlation the former demonstrates. Both SNPs are risk loci for migraine. At the same time, it has been mentioned above that rs17862920 is associated with sensing cold pain in Finnish and Norwegian individuals, with C allele carriers (more prevalent in northern latitudes) being more susceptible. It can be assumed that both loci are independently involved in adaptation to low temperatures. In that case, however, the adaptive role of rs17862920-C allele is hard to explain. The possible mechanism of differential survival might be avoidance of potentially lethal hypothermia by those harboring С allele.
As for rs6723922 and rs11682848 loci, none of them shows any sign of selective sweep in the XPCLR test. Probably, those are false positives, or at least, linked to an unobserved variant under selection. It is also worth noting that SNP evidence scores for these loci are quite low.
In addition to the search for signatures of selection, we would like to note some details on the BayScEnv test not published anywhere (as far as we know).
Changing model parameters drastically affects the output in BayScEnv (see Supplementary Table 6). For example, significant results (q value <0.05) were obtained when using model parameters pi = 0.5/p = 0.1 (SNP rs11682848 being significant) or pi = 0.9/p = 0.1 (18/60 SNPs being significant). At the same time, empirical p values are more stable. Thus, we suggest using them in hypothesis-driven studies (with default model parameters) of local adaptation and choosing the significance threshold based on expert’s opinions rather than relying on FDR outputs.
Counterintuitively, a reduction in the number of tests in BayScEnv does not lead to a greater number of statistically significant FDR outputs (the posterior error probability and the q value). Furthermore, in our case, given parameters pi = 0.9/ p = 0.1, 18 out of a total of 60 SNPs reach significance level when analyzing 1,369 markers, whereas none is significant when using 60 SNPs. This discrepancy might be explained by less precise constructing a null model of population structure.
Conclusions
Several lines of evidence point to possible involvement of rs7577262 in cold adaptation. This SNP is considered the best candidate based on its allelic correlations with winter temperatures, signatures of selective sweep and physiological evidences. The second top SNP, rs17862920, may participate in adaptation as well. As for rs6723922 and rs11682848 loci, these appear to be false positives or at least linked to some unobserved selected variant.
Author Contributions
NY and MV conceived the project. KG supervised the project. AI and KG processed and analyzed the data. AI, KG, and NY drafted the manuscript.
Funding
This study was supported by budget from project No. 0324-2019-0041 of the Federal Research Center «Institute of Cytology and Genetics» SB RAS (ICG SB RAS).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The Common Use Center «Bioinformatics» (ICG SB RAS) is gratefully acknowledged for providing computer facilities.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00759/full#supplementary-material
References
B1000 Genomes Project Consortium, Abecasis, G. R., Altshuler, D., Auton, A., Brooks, L. D., Durbin, R. M., et al. (2010). A map of human genome variation from population-scale sequencing. Nature 467, 1061–1073. doi: 10.1038/nature09534
Andolfatto, P. (2005). Adaptive evolution of non-coding DNA in Drosophila. Nature 437, 1149–1152. doi: 10.1038/nature04107
Anttila, V., Winsvold, B. S., Gormley, P., Kurth, T., Bettella, F., McMahon, G., et al. (2013). Genome-wide meta-analysis identifies new susceptibility loci for migraine. Nat. Genet. 45, 912–917. doi: 10.1038/ng.2676
Bautista, D. M., Siemens, J., Glazer, J. M., Tsuruda, P. R., Basbaum, A. I., Stucky, C. L., et al. (2007). The menthol receptor TRPM8 is the principal detector of environmental cold. Nature 448, 204–208. doi: 10.1038/nature05910
Bazykin, G. A., Kondrashov, A. S. (2012). Major role of positive selection in the evolution of conservative segments of Drosophila proteins. Proc. Biol. Sci. 279, 3409–3417. doi: 10.1098/rspb.2012.0776
Beccari, A. R., Gemei, M., Lo Monte, M., Menegatti, N., Fanton, M., Pedretti, A., et al. (2017). Novel selective, potent naphthyl TRPM8 antagonists identified through a combined ligand- and structure-based virtual screening approach. Sci. Rep. 7, 10999. doi: 10.1038/s41598-017-11194-0
Burstein, R., Yarnitsky, D., Goor-Aryeh, I., Ransil, B. J., Bajwa, Z. H. (2000). An association between migraine and cutaneous allodynia. Ann. Neurol. 47, 614–624. doi: 10.1002/1531-8249(200005)47:5<614::AID-ANA9>3.0.CO;2-N
Cai, J. J., Macpherson, J. M., Sella, G., Petrov, D. A. (2009). Pervasive hitchhiking at coding and regulatory sites in humans. PLoS Genet. 5, e1000336. doi: 10.1371/journal.pgen.1000336
Cann, H. M., de Toma, C., Cazes, L., Legrand, M. F., Morel, V., Piouffre, L., et al. (2002). A human genome diversity cell line panel. Science 296, 261–262. doi: 10.1126/science.296.5566.261b
Cardona, A., Pagani, L., Antao, T., Lawson, D. J., Eichstaedt, C. A., Yngvadottir, B., et al. (2014). Genome-wide analysis of cold adaptation in indigenous Siberian populations. PLoS One 9, e98076. doi: 10.1371/journal.pone.0098076
Casper, J., Zweig, A. S., Villarreal, C., Tyner, C., Speir, M. L., Rosenbloom, K. R., et al. (2018). The UCSC genome browser database: 2018 update. Nucleic Acids Res. 46, D762–D769. doi: 10.1093/nar/gkx1020
Chasman, D. I., Schürks, M., Anttila, V., de Vries, B., Schminke, U., Launer, L. J., et al. (2011). Genome-wide association study reveals three susceptibility loci for common migraine in the general population. Nat. Genet. 43, 695–698. doi: 10.1038/ng.856
Chen, H., Patterson, N., Reich, D. (2010). Population differentiation as a test for selective sweeps. Genome Res. 20, 393–402. doi: 10.1101/gr.100545.109
Colburn, R. W., Lubin, M. L., Stone, D. J., Jr., Wang, Y., Lawrence, D., D’Andrea, M. R., et al. (2007). Attenuated cold sensitivity in TRPM8 null mice. Neuron 54, 379–386. doi: 10.1016/j.neuron.2007.04.017
Daanen, H. A., Lichtenbelt, W. D. (2016). Human whole body cold adaptation. Temperature 3, 104–118. doi: 10.1080/23328940.2015.1135688
de Villemereuil, P., Gaggiotti, O. E. (2015). A new FST-based method to uncover local adaptation using environmental variables. Methods Ecol. Evol. 6, 1248–1258. doi: 10.1111/2041-210X.12418
Dhaka, A., Murray, A. N., Mathur, J., Earley, T. J., Petrus, M. J., Patapoutian, A. (2007). TRPM8 is required for cold sensation in mice. Neuron 54, 371–378. doi: 10.1016/j.neuron.2007.02.024
Fariello, M. I., Servin, B., Tosser-Klopp, G., Rupp, R., Moreno, C., International Sheep Genomics Consortium, et al. (2014). Selection signatures in worldwide sheep populations. PLoS One 9, e103813. doi: 10.1371/journal.pone.0103813
Fernández, J. A., Skryma, R., Bidaux, G., Magleby, K. L., Scholfield, C. N., McGeown, J. G., et al. (2011). Voltage- and cold-dependent gating of single TRPM8 ion channels. J. Gen. Physiol. 137, 173–195. doi: 10.1085/jgp.201010498
Freilinger, T., Anttila, V., de Vries, B., Malik, R., Kallela, M., Terwindt, G. M., et al. (2012). Genome-wide association analysis identifies susceptibility loci for migraine without aura. Nat. Genet. 44, 777–782. doi: 10.1038/ng.2307
Frichot, E., Schoville, S. D., Bouchard, G., François, O. (2013). Testing for associations between loci and environmental gradients using latent factor mixed models. Mol. Biol. Evol. 30, 1687–1699. doi: 10.1093/molbev/mst063
Ghosh, J., Pradhan, S., Mittal, B. (2013). Genome-wide-associated variants in migraine susceptibility: a replication study from North India. Headache 53, 1583–1594. doi: 10.1111/head.12240
Günther, T., Coop, G. (2013). Robust identification of local adaptation from allele frequencies. Genetics 195, 205–220. doi: 10.1534/genetics.113.152462
Halligan, D. L., Oliver, F., Guthrie, J., Stemshorn, K. C., Harr, B., Keightley, P. D. (2011). Positive and negative selection in murine ultra-conserved noncoding elements. Mol. Biol. Evol. 28, 2651–2660. doi: 10.1093/molbev/msr093
Hancock, A. M., Clark, V. J., Qian, Y., Di Rienzo, A. (2011b). Population genetic analysis of the uncoupling proteins supports a role for UCP3 in human cold resistance. Mol. Biol. Evol. 28, 601–614. doi: 10.1093/molbev/msq228
Hancock, A. M., Witonsky, D. B., Alkorta-Aranburu, G., Beall, C. M., Gebremedhin, A., Sukernik, R., et al. (2011a). Adaptations to climate-mediated selective pressures in humans. PLoS Genet. 7, e1001375. doi: 10.1371/journal.pgen.1001375
Hancock, A. M., Witonsky, D. B., Gordon, A. S., Eshel, G., Pritchard, J. K., Coop, G., et al. (2008). Adaptations to climate in candidate genes for common metabolic disorders. PLoS Genet. 4, e32. doi: 10.1371/journal.pgen.0040032
He, J., Kelly, T. N., Zhao, Q., Li, H., Huang, J., Wang, L., et al. (2013). Genome-wide association study identifies 8 novel loci associated with blood pressure responses to interventions in Han Chinese. Circ. Cardiovasc. Genet. 6, 598–607. doi: 10.1161/CIRCGENETICS.113.000307
Hershkovitz, I., Weber, G. W., Quam, R., Duval, M., Grün, R., Kinsley, L., et al. (2018). The earliest modern humans outside Africa. Science 359, 456–459. doi: 10.1126/science.aap8369
Hublin, J. J., Ben-Ncer, A., Bailey, S. E., Freidline, S. E., Neubauer, S., Skinner, M. M., et al. (2017). New fossils from Jebel Irhoud, Morocco and the pan-African origin of homo sapiens. Nature 546, 289–292. doi: 10.1038/nature22336
Jablonski, N. G., Chaplin, G. (2010). Human skin pigmentation as an adaptation to UV radiation. Proc. Natl. Acad. Sci. U. S. A. 107 Suppl 2, 8962–8968. doi: 10.1073/pnas.0914628107
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A., Jermiin, L. S. (2017). Model finder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285
Kaunisto, M., Holmström, E., Anttila, V., Kallela, M., Hämäläinen, E., Stubhaug, A., et al (2013). TRPM8 variants that protect from migraine make individuals less susceptible to cold pain. In: Presented at the 63rd Annual Meeting of The American Society of Human Genetics. http://www.ashg.org/2013meeting/abstracts/fulltext/f130121596.htm.
Key, F. M., Abdul-Aziz, M. A., Mundry, R., Peter, B. M., Sekar, A., D’Amato, M., et al. (2018). Human local adaptation of the TRPM8 cold receptor along a latitudinal cline. PLoS Genet. 14, e1007298. doi: 10.1371/journal.pgen.1007298
Kozyreva, T. V., Tkachenko, E., Potapova, T. A., Voevoda, M. I. (2014). Respiratory system response to cooling in subjects with single nucleotide polymorphism rs11562975 in gene of thermosensitive TRPM8 ion channel. Fiziol Cheloveka 40, 94–98. doi: 10.1134/S0362119714020108
Kozyreva, T. V., Tkachenko, E., Potapova, T. A., Romashchenko, A. G., Voevoda, M. I. (2011). Relationship of single-nucleotide polymorphism rs11562975 in thermo-sensitive ion channel TRPM8 gene with human sensitivity to cold and menthol. Fiziol Cheloveka 37, 71–76. doi: 10.1134/S0362119711020101
Larra, M. F., Schilling, T. M., Röhrig, P., Schächinger, H. (2015). Enhanced stress response by a bilateral feet compared to a unilateral hand Cold Pressor Test. Stress 18, 589–596. doi: 10.3109/10253890.2015.1053452
Leppäluoto, J., Hassi, J. (1991). Human physiological adaptations to the arctic climate. Arctic 44, 139–145. doi: 10.14430/arctic1530
Liu, Z., Ji, Z., Wang, G., Chao, T., Hou, L., Wang, J. (2016). Genome-wide analysis reveals signatures of selection for important traits in domestic sheep from different ecoregions. BMC Genomics 17, 863. doi: 10.1186/s12864-016-3212-2
Machiela, M. J., Chanock, S. J. (2015). LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 31, 3555–3557. doi: 10.1093/bioinformatics/btv402
Majhi, R. K., Saha, S., Kumar, A., Ghosh, A., Swain, N., Goswami, L., et al. (2015). Expression of temperature-sensitive ion channel TRPM8 in sperm cells correlates with vertebrate evolution. PeerJ. 3, e1310. doi: 10.7717/peerj.1310
Matos-Cruz, V., Schneider, E. R., Mastrotto, M., Merriman, D. K., Bagriantsev, S. N., Gracheva, E. O. (2017). Molecular prerequisites for diminished cold sensitivity in ground squirrels and hamsters. Cell Rep. 21, 3329–3337. doi: 10.1016/j.celrep.2017.11.083
Meng, W., Adams, M. J., Hebert, H. L., Deary, I. J., McIntosh, A. M., Smith, B. H. (2018). A genome-wide association study finds genetic associations with broadly-defined headache in UK Biobank (N = 223,773). E Bio. Med. 28, 180–186. doi: 10.1016/j.ebiom.2018.01.023
Myers, B. R., Sigal, Y. M., Julius, D. (2009). Evolution of thermal response properties in a cold-activated TRP channel. PLoS One 4, e5741. doi: 10.1371/journal.pone.0005741
Ohashi, J., Naka, I., Tsuchiya, N. (2011). The impact of natural selection on an ABCC11 SNP determining earwax type. Mol. Biol. Evol. 28, 849–857. doi: 10.1093/molbev/msq264
Paradis, E., Claude, J., Strimmer, K. (2004). APE: analyses of phylogenetics and evolution in r language. Bioinformatics 20, 289–290. doi: 10.1093/bioinformatics/btg412
Park, H. W., Kim, S. H., Chang, Y. S., Kim, S. H., Jee, Y. K., Lee, A. Y., et al. (2018). The Fas signaling pathway is a common genetic risk factor for severe cutaneous drug adverse reactions across diverse drugs. Allergy Asthma Immunol. Res. 10, 555–561. doi: 10.4168/aair.2018.10.5.555
Potapova, T. A., Babenko, V. N., Kobzev, V. F., Romashchenko, A. G., Maksimov, V. N., Voevoda, M. I. (2014). Associations of cold receptor TRPM8 gene single nucleotide polymorphism with blood lipids and anthropometric parameters in Russian population. Bull. Exp. Biol. Med. 157, 757–761. doi: 10.1007/s10517-014-2660-4
Quagliarello, A., De Fanti, S., Giuliani, C., Abondio, P., Serventi, P., Sarno, S., et al. (2017). Multiple selective events at the PRDM16 functional pathway shaped adaptation of western European populations to different climate conditions. J. Anthropol. Sci. 95, 235–247. doi: 10.4436/JASS.95011
Saito, S., Tominaga, M. (2017). Functional diversity and evolutionary dynamics of thermoTRP channels. Cell Calcium 57, 214–221. doi: 10.1016/j.ceca.2014.12.001
Sazzini, M., Schiavo, G., De Fanti, S., Martelli, P. L., Casadio, R., Luiselli, D. (2014). Searching for signatures of cold adaptations in modern and archaic humans: hints from the brown adipose tissue genes. Heredity 113, 259–267. doi: 10.1038/hdy.2014.24
Scheet, P., Stephens, M. (2006). A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am. J. Hum. Genet. 78, 629–644. doi: 10.1086/502802
Spearman, C. (1904). The proof and measurement of association between two things. Am. J. Psychol. 15, 72–101. doi: 10.2307/1412159
Tishkoff, S. A., Reed, F. A., Ranciaro, A., Voight, B. F., Babbitt, C. C., Silverman, J. S., et al. (2007). Convergent adaptation of human lactase persistence in Africa and Europe. Nat. Genet. 39, 31–40. doi: 10.1038/ng1946
Valverde, G., Zhou, H., Lippold, S., de Filippo, C., Tang, K., LópezHerráez, D., et al. (2015). A novel candidate region for genetic adaptation to high altitude in Andean populations. PLoS One 10, e0125444. doi: 10.1371/journal.pone.0125444
Vatsiou, A. I., Bazin, E., Gaggiotti, O. E. (2016). Detection of selective sweeps in structured populations: a comparison of recent methods. Mol. Ecol. 25, 89–103. doi: 10.1111/mec.13360
Voskarides, K. (2018). Combination of 247 genome-wide association studies reveals high cancer risk as a result of evolutionary adaptation. Mol. Biol. Evol. 35, 473–485. doi: 10.1093/molbev/msx305
Walther, A. (2010). Germline determinants of outcome and risk in colorectal cancer. Ph.D. thesis. London: University College London.
Keywords: TRPM8, environmental correlation analysis, SNP, cold adaptation, East Asian ancestry
Citation: Igoshin AV, Gunbin KV, Yudin NS and Voevoda MI (2019) Searching for Signatures of Cold Climate Adaptation in TRPM8 Gene in Populations of East Asian Ancestry. Front. Genet. 10:759. doi: 10.3389/fgene.2019.00759
Received: 30 October 2018; Accepted: 17 July 2019;
Published: 23 August 2019.
Edited by:
Ancha Baranova, George Mason University, United StatesReviewed by:
Margarida Matos, University of Lisbon, PortugalToni Gossmann, University of Sheffield, United Kingdom
Copyright © 2019 Igoshin, Gunbin, Yudin and Voevoda. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexander V. Igoshin, aWdvc2hpbkBiaW9uZXQubnNjLnJ1